Search Results - EGA European Genome-Phenome Archive

Systematic comparative analysis of single-nucleotide variants detection methods from single-cell RNA sequencing data

Systematic interrogation of single nucleotide variations (SNVs) is one of the most promising approaches to delineate the cellular heterogeneity and phylogenetic relationships at the single cell level. While SNV detection from abundant single cell RNA sequencing (scRNA-seq) data is applicable and cost-effective in identifying expressed variants, inferring sub-clones, and deciphering genotype-phenotype linkages, there is a lack of computational methods specifically developed for SNV calling in scRNA-seq. Although variant callers for bulk RNA-seq have been sporadically used in scRNA-seq, the performances of different tools have not been assessed. Here, we performed a systematic comparison of seven tools including SAMtools, the GATK pipeline, CTAT, FreeBayes, MuTect2, Strelka2 and VarScan2, using both simulation and scRNA-seq datasets, and identified multiple elements influencing their performance. Our study provided the first benchmarking to evaluate the performances of different SNV detection tools for scRNA-seq data.

Study EGAS00001003883

CRISPR transduction of iPS cells

Genome engineering using CRISPR/Cas9 technology enables simple, efficient and precise genomic modifications in human cells. Conventional immortalized cell lines can be easily edited or screened using genome-wide libraries with lentiviral transduction. However, cell types derived from the differentiation of induced Pluripotent Stem Cells (iPSC), which often represent more relevant, patient-derived models for human pathology, are much more difficult to engineer as CRISPR/Cas9 delivery to these differentiated cells can be inefficient and toxic. Here, we present an efficient, lentiviral transduction protocol for delivery of CRISPR/Cas9 to macrophages derived from human iPSC with efficiencies close to 100%. We demonstrate CRISPR/Cas9 knockouts for three non-essential proof-of-concept genes - HPRT1, PPIB and CDK4. We then scale the protocol and validate for a genome-wide pooled CRISPR/Cas9 loss-of-function screen. This methodology enables, for the first time, systematic exploration of macrophage involvement in immune responses, chronic inflammation, neurodegenerative diseases and cancer progression, using efficient genome editing techniques.

Study EGAS00001005102

The genomic landscape of early stage ovarian high grade serous carcinoma

To determine whether early-stage ovarian high grade serous carcinoma (HGSC) represents a distinct genomic entity, we collected samples from 43 patients with stage IA-IIA HGSC to identify potential differences in short genomic variants and copy number aberrations, and compared them to a cohort of 52 late-stage (stage IIIC-IV) cases. We found no significant differences in somatic mutations or focal copy number alterations between early-stage and late-stage cohorts. There was, however, a significant difference in both ploidy and copy number signature exposure between early and late-stage samples, with higher ploidy and signature 4 exposure in late-stage cases. Unsupervised hierarchical clustering revealed three clusters, which were prognostic. Together, our data suggest that early and late-stage HGSC share fundamental genomic features, but that late- stage disease appears distinct from early-stage, with evidence of whole genome duplication that may provide evolutionary benefit.

Study EGAS00001005567

Integrative molecular analysis of pediatric Anaplastic large cell lymphoma reveals subtypes with distinct immune suppression signatures.

Anaplastic large cell lymphoma (ALCL) is a peripheral T-cell lymphoma accounting for 10–15% of all childhood lymphomas. While more than 90% of the ALCL cases contain ALK-rearrangement, these tumors possess significant inter-tumor molecular heterogeneity that contributes to distinct morphologic differences and clinical impact. To gain insight into the molecular heterogeneity within ALK+ ALCL, we performed whole-exome sequencing, RNA-sequencing, and methylome analysis of 42 primary pediatric ALK+ ALCL patients. Our data showed that ALK+ALCLs was subclassified into two subtypes based on ALK gene expression, methylation profiles, and somatic mutation patterns. ALK-low samples had more highly methylated gene regions while enriched with immune infiltration. ALK-high samples were enriched with somatic copy number alteration and high expression of MYC and PD-L1. These data indicate that ALK expression, somatic mutations, and aneuploidy status are negatively associated with immune infiltration, stipulating that ALK expression status might be associated with resistance to immune checkpoint inhibitors.

Study EGAS00001004189

Exome sequencing of familial high-grade serous ovarian carcinoma reveals heterogeneity for rare candidate susceptibility genes

High-grade serous ovarian carcinoma (HGSOC) has a significant hereditary component, approximately half of which cannot be explained by known genes. To discover genes, we analyse germline exome sequencing data from 516 BRCA1/2-negative women with HGSOC, focusing on genes enriched with rare, protein-coding loss-of-function (LoF) variants. Overall, there is a significant enrichment of rare protein-coding LoF variants in the cases (p < 0.0001, chi-squared test). Only thirty-four (6.6%) have a pathogenic variant in a known or proposed predisposition gene. Few genes have LoF mutations in more than four individuals and the majority are detected in one individual only. Forty-three highly-ranked genes are identified with three or more LoF variants that are enriched by three-fold or more compared to GnomAD. These genes represent diverse functional pathways with relatively few involved in DNA repair, suggesting that much of the remaining heritability is explained by previously under-explored genes and pathways.

Study EGAS00001004235

A microRNA-signature that correlates with cognition and is a target against cognitive decline

While some individuals age without pathological memory impairments, others develop age-associated cognitive diseases. Since changes in cognitive function develop slowly over time in these patients, they are often diagnosed at an advanced stage of molecular pathology, a time point when causative treatments fail. Thus, there is great need for the identification of inexpensive and minimal invasive approaches that could be used for screening with the aim to to identify individuals at risk for cognitive decline that can then undergo further diagnostics and eventually stratified therapies. In this study we use an integrative approach combining the analysis of human data and mechanistic studies in model systems to identify a circulating 3-microRNA signature that reflect key processes linked to neural homeostasis and informs about cognitive status. We furthermore provide evidence that expression changes of this signature represent multiple mechanisms de-regulated in the aging and diseased brain and are a suitable target for RNA therapeutics.

Study EGAS00001005627

Whole genome sequencing data for 10 hepatocellular carcinomas (HCC) and matched non-tumor liver tissues + optical mapping data for 4 HCC and 3 matched non-tumor liver tissues.

Oncogene activation leads to replication stress and promotes genomic instability. Here we combine optical mapping and whole genome sequencing to explore in-depth the nature of structural variants (SVs) induced by replication stress in cyclin-activated hepatocellular carcinomas (CCN-HCC). In addition to classical tandem duplications, CCN-HCC display frequent intra- and inter-chromosomal templated insertion cycles (TIC) likely resulting from template switching events. Template switching preferentially involves active topologically associated domains that are close in the 3D genome organization. Template sizes depend on the type of cyclin activation and are coordinated within each TIC. Replication stress induces continuous accumulation of SVs during CCN-HCC progression, fostering the acquisition of new driver alterations and large-scale copy-number changes at TIC borders. Together, this analysis sheds light on the mechanisms, dynamics and consequences of SV accumulation in tumors with oncogene-induced replication stress.

Study EGAS00001005629

Single-cell RNA-seq of bronchoalveolar lavage (BAL) fluid of late stage severe COVID-19 patients

COVID-19-induced ‘acute respiratory distress syndrome’ (ARDS) is associated with prolonged respiratory failure and high mortality, but the mechanistic basis of lung injury remains incompletely understood. Here, we analyzed pulmonary immune responses and lung pathology in two cohorts of patients with COVID-19 ARDS using functional single cell genomics, immunohistology and electron microscopy. We describe an accumulation of CD163-expressing monocyte-derived macrophages, that acquired a profibrotic transcriptional phenotype during COVID-19 ARDS. Gene set enrichment and computational data integration revealed a significant similarity between COVID-19-associated macrophages and profibrotic macrophage populations identified in idiopathic pulmonary fibrosis. COVID-19 ARDS was associated with clinical, radiographic, histopathological, and ultrastructural hallmarks of pulmonary fibrosis. Exposure of human monocytes to SARS-CoV-2, but not Influenza A virus or viral RNA analogs, was sufficient to induce a similar profibrotic phenotype in vitro. In conclusion, we demonstrate that SARS-CoV-2 triggers profibrotic macrophage responses and pronounced fibroproliferative ARDS.

Study EGAS00001005634

Genetic and immune landscape evolution defines subtypes of MMR deficient colorectal cancer

Mismatch repair deficient colorectal cancers have high mutation loads and many respond to immune checkpoint-inhibitors. We investigated how genetic and immune landscapes co-evolve in these tumors. All cases had high truncal mutation loads. Driver aberrations showed a clear hierarchy despite pervasive intratumor heterogeneity: Those in WNT/βCatenin, mitogen-activated protein kinase and TGFβ receptor family genes were almost always truncal. Immune evasion drivers were predominantly subclonal and showed parallel evolution. Pan-tumor evolution, subclonal evolution, and evolutionary stasis of genetic immune evasion drivers defined three MMRd CRC subtypes with distinct T-cell infiltrates. These immune evasion drivers have been implicated in checkpoint-inhibitor resistance. Clonality and subtype assessments are hence critical for predictive immunotherapy biomarker development. Cancer cell PD-L1 expression was conditional on loss of the intestinal homeobox transcription factor CDX2. This explains infrequent PD-L1 expression by cancer cells and likely contributes to the high recurrence risk of MMRd CRCs with impaired CDX2 expression.

Study EGAS00001005769

Circulating tumor cells for comprehensive and multiregional non-invasive genetic characterization of multiple myeloma (arrays set)

Multiple myeloma (MM) patients undergo repetitive bone marrow (BM) aspirates for genetic characterization. Circulating tumor cells (CTCs) are detectable in peripheral blood of virtually all MM cases and are prognostic, but their applicability for non-invasive screening has been poorly investigated. Here, we used next-generation flow cytometry to isolate matched CTCs and BM tumor cells from 42 patients and compare their genetic profile. In 8 cases, tumor cells from extramedullary plasmacytomas were also sorted and whole-exome sequencing performed in the three spatially-distributed clones. Noteworthy, ≥84% mutations present in BM and extramedullary clones were detectable in CTCs. Concordance between BM tumor cells and CTCs was high for chromosome arm-level copy-number alterations (≥94%) though not for translocations (43%). All high-risk genetic abnormalities except one t(4;14) were detected in CTCs whenever present in BM tumor cells. Altogether, these results support CTCs for non-invasive risk-stratification of MM patients based on their numbers and genetic profile.

Study EGAS00001004314

A molecular signature for IL-10-producing Th1 cells in protozoan parasitic diseases

Control of the intracellular parasites that cause malaria and visceral leishmaniasis (VL) is dependent on the generation of pro-inflammatory, IFNγ+ Tbet+ CD4+ T (Th1) cells by infected hosts. Immunoregulatory IL-10 produced by Th1 cells serves to mitigate subsequent inflammation and related disease pathology. However, these IL-10-producing Th1 (Tr1) cells can also not only promote parasite persistence, but may impair immunity to re-infection and potential vaccine efficacy. Here, we identify molecular and phenotypic signatures that distinguish Th1 cells from Tr1 cells in both experimental VL, caused by infection of C57BL/6 mice with the human parasite Leishmania donovani, and in Plasmodium falciparum-infected humans participating in controlled human malaria infection (CHMI) studies. Such characterisations allow for the better understanding of Tr1 cell development and function in the context of these parasitic diseases, and for the identification of targets for immune modulation to improve anti-parasitic immunity.

Study EGAS00001004454

Combined clinical and gene expression score identifies high-risk individuals among follicular lymphoma patients on immunotherapy

Follicular lymphoma (FL) is usually an indolent malignant B-cell lymphoma, but at least 30% of the patients experience transformation to Diffuse Large B-cell Lymphoma (DLBCL) with inferior survival. Early identification of patients with poor prognosis is essential for tailoring treatment. Here, we describe a network-based integrative analysis of 168 FL-biopsies from 127 patients with long clinical follow-up and transformation status. A distinct shift in expression towards a DLBCL expression phenotype is observed in BTK associated genes in cases with transformation. This shift is present years prior to transformation. An integrated score combining the BTK score and a clinical index is proposed, termed the BTK-FLIPI score. It identifies patients (16%) with very poor prognosis, with median time to progression 21.7 months (95% CI: 14.3-42.6) versus 70.9 months (95% CI: 41.5-95.0) for the remainder. The results suggest BTK as a potential target for treatment of high-risk FL-patients.

Study EGAS00001002566

The evolution of hematopoietic cells under cancer therapy

Chemotherapies may increase mutagenesis of healthy cells and change the selective pressures in tissues thus influencing their evolution. However, their contributions to the mutation burden and clonal expansions of healthy somatic tissues are not clear. Here, exploiting the mutational footprint of some chemotherapies, we explore their influence on the evolution of hematopoietic cells. Cells of Acute Myeloid Leukemia (AML) secondary to treatment with platinum-based drugs showed the mutational footprint of these drugs, indicating that non-malignant blood cells received chemotherapy mutations. No trace of 5-fluorouracil (5FU) mutational signature was found in AML secondary to exposure to 5FU, suggesting that cells establishing the AML were quiescent during treatment. Using the platinum-based mutational signature as a barcode, we determined that secondary AMLs began their clonal expansion after the start of the cytotoxic treatment. Its absence in clonal hematopoiesis cases is consistent with a clonal expansion that predates the exposure to the cytotoxic agent.

Study EGAS00001005234

Myoepithelial progenitors as founder cells of hyperplastic human breast lesions upon PIK3CA transformation

The myoepithelial (MEP) lineage of human breast comprises bipotent and multipotent progenitors in ducts and terminal duct lobular units (TDLUs). We here assess whether this heterogeneity impacts on oncogenic PIK3CA transformation. Single cell RNA sequencing (scRNA-seq) and multicolor imaging reveal that terminal ducts represent the most enriched source of cells with ductal MEP markers including a-smooth muscle actin (a-SMA), keratin K14, K17 and CD200. Furthermore, we find neighboring CD200high and CD200low progenitors within terminal ducts. When sorted and kept in ground state conditions, their CD200low and CD200high phenotypes are preserved. Upon differentiation, progenitors remain multipotent and bipotent, respectively. Immortalized progenitors are transduced with shp53 and mutant PIK3CA. Upon transplantation, CD200low MEP progenitors distinguish from CD200high by the formation of multilayered structures with a hyperplastic inner layer of luminal epithelial cells. We suggest a model with spatially distributed MEP progenitors as founder cells of biphasic breast lesions with implications for early detection and prevention strategies.

Study EGAS00001005933

Gain of function DNMT3A mutations cause microcephalic dwarfism and hypermethylation of Polycomb-regulated regions

DNA methylation and Polycomb are key factors in the establishment of vertebrate cellular identity and fate. Here we report de novo missense mutations in DNMT3A, encoding the DNA methyltransferase DNMT3A, that cause microcephalic dwarfism, a hypocellular disorder of extreme global growth failure. Substitutions in the PWWP domain abrogate binding to the histone modifications H3K36me2/3, and alter DNA methylation in patient cells. Polycomb-associated DNA methylation canyons/valleys, hypomethylated domains encompassing developmental genes, become methylated with concomitant depletion of H3K27me3 and H3K4me3 bivalent marks. Such de novo DNA methylation occurs during differentiation of Dnmt3aW326R pluripotent cells in vitro, and is also evident in Dnmt3aW326R/+ dwarf mice. We therefore propose that the interaction of the DNMT3A PWWP domain with H3K36me2/3 normally limits DNA methylation of polycomb-marked regions. Our findings implicate the interplay between DNA methylation and polycomb at key developmental regulators as a determinant of organism size in mammals.

Study EGAS00001003232

Fibroblast heterogeneity drives metastatic spread in breast cancer through distinct mechanisms

The major cause of death among breast cancer patients is metastatic dissemination, suggesting that understanding the mechanisms driving this process remains an unmet medical need. Here, by concomitantly using 5 markers, we identify 4 Carcinoma-Associated Fibroblasts (CAF) subsets in primary tumors and corresponding metastatic lymph nodes (LN). Two CAF subsets, called CAF-S1 and CAF-S4, correlate with the proportion of cancer cells in invaded LN and exert distinct functions in tumor cell spread. CAF-S1 promote cancer cell migration and initiates an epithelial to mesenchymal transition through a CXCL12 / TGFb positive loop. For their part, the highly contractile CAF-S4 foster cancer cell motility and invasion in 3-dimensional matrix via Notch pathway activation. Finally, CAF-S1 and CAF-S4 accumulation in LN is a new prognostic factor at diagnostic, independent of breast cancer subtypes and LN status, strengthening the validity of our findings in human breast cancer pathophysiology.

Study EGAS00001003238

Homologous recombination DNA repair deficiency and PARP inhibition activity in primary triple negative breast cancer

Triple negative breast cancer (TNBC) encompasses molecularly different subgroups, with a subgroup harboring evidence of defective homologous recombination (HR) DNA repair. Within a phase 2 window clinical trial, RIO trial (EudraCT 2014-003319-12), we investigate the activity of PARP inhibitors in 43 patients with untreated TNBC. The primary endpoint, decreased Ki67, occurs in 12% of TNBC. Multiple different approaches identify HR deficiency in the majority of TNBC, 69% with the mutational-signature-based HRDetect assay. Cancers with mutational signatures of HR deficiency have a functional defect in HR. Following rucaparib treatment there is no association of Ki67 change with HR deficiency. In contrast, early circulating tumor DNA dynamics identify activity of rucaparib in cancers with mutation-signatures of HR deficiency, with rucaparib inducing expression of interferon response genes in HR deficient cancers. The majority of TNBCs have a defect in DNA repair, identifiable by mutational signature analysis, that may be targetable with PARP inhibitors.

Study EGAS00001004405

Genome-wide cell-free DNA mutational integration enables ultra-sensitive cancer monitoring

In many areas of oncology, we lack sensitive tools to track low burden disease. While cell-free DNA (cfDNA) shows promise in detecting cancer mutations, we found that the combination of low tumor fraction (TF) and limited number of DNA fragments restricts low disease burden monitoring through the prevailing deep targeted sequencing paradigm. We reasoned that breadth may supplant depth of sequencing to overcome the barrier of cfDNA abundance. Whole genome sequencing (WGS) of cfDNA allowed ultra-sensitive detection, capitalizing on the cumulative signal of thousands of somatic mutations observed in solid malignancies, with TF detection sensitivity as low as 10-5. The WGS approach enabled dynamic tumor burden tracking and post-operative residual disease detection, associated with adverse outcome. Thus, we present an orthogonal framework for cfDNA cancer monitoring via genome wide mutational integration, enabling ultra-sensitive detection, overcoming the limitation of cfDNA abundance, and empowering treatment optimization in low-disease burden oncology care

Study EGAS00001004406

SARS‐CoV‐2 receptor ACE2 and TMPRSS2 are primarily expressed in bronchial transient secretory cells

The SARS‐CoV‐2 pandemic affecting the human respiratory system severely challenges public health and urgently demands for increasing our understanding of COVID‐19 pathogenesis, especially host factors facilitating virus infection and replication. SARS‐CoV‐2 was reported to enter cells via binding to ACE2, followed by its priming by TMPRSS2. Here, we investigate ACE2 and TMPRSS2 expression levels and their distribution across cell types in lung tissue (twelve donors, 39,778 cells) and in cells derived from subsegmental bronchial branches (four donors, 17,521 cells) by single nuclei and single cell RNA sequencing, respectively. While TMPRSS2 is strongly expressed in both tissues, in the subsegmental bronchial branches ACE2 is predominantly expressed in a transient secretory cell type. Interestingly, these transiently differentiating cells show an enrichment for pathways related to RHO GTPase function and viral processes suggesting increased vulnerability for SARS‐CoV‐2 infection. Our data provide a rich resource for future investigations of COVID‐19 infection and pathogenesis.

Study EGAS00001004419

Somatic_mutation_profiling_of_intestinal_crypts_from_IBD_Full_STDY

Following on from our pilot study we will investigate further the role that somatic mutations in intestinal stem cells play in IBD. We will use laser capture microdissection to ascertain epithelial cells from an individual crypts/vilus. We will extract DNA from these and use a new low-input material library prep methods to make accurate libraries for sequencing. We will perform 15X whole-genome sequencing of crypts biopsied from patients. Biopsies have been taken from inflammed and noninflammed patches to compare mutation rates and mutational signatures. We will also test if crypts from inflammed patches of a single patient are clonal (either within a given inflammed patch or across inflammed patches). We will quantify the rate of cancer driver mutations in inflammed and non-inflammed crypts. 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 EGAS00001003249

Custom long non-coding RNA capture enhances detection sensitivity in different human sample types.

Long non-coding RNAs (lncRNAs) are a heterogeneous group of transcripts that lack protein coding potential and display regulatory functions in various cellular processes. As a result of their cell- and cancer-specific expression patterns, lncRNAs have emerged as potential diagnostic and therapeutic targets. The accurate characterization of lncRNAs in bulk transcriptome data remains challenging due to their low abundance compared to protein coding genes. To tackle this issue, we describe a unique short-read custom lncRNA capture sequencing approach that relies on a comprehensive set of 565,878 capture probes for 49,372 human lncRNA genes. This custom lncRNA capture approach was evaluated on various sample types ranging from artificial high-quality RNA mixtures to more challenging formalin-fixed paraffin-embedded tissue and biofluid material. The custom enrichment approach allows the detection of a more diverse repertoire of lncRNAs, with better reproducibility and higher coverage compared to classic total RNA-sequencing.

Study EGAS00001005418

Detection of cell-free DNA fragmentation and copy number alterations in cerebrospinal fluid from glioma patients

Glioma is difficult to detect or characterise using current liquid biopsy approaches. Detection of cell-free tumor DNA (cftDNA) in cerebrospinal fluid (CSF) has been proposed as an alternative to detection in plasma. We used shallow whole-genome sequencing (sWGS, at a coverage of < 0.4x) of cell-free DNA from the CSF of 13 patients with primary glioma to determine somatic copy number alterations and DNA fragmentation patterns. This allowed us to determine the presence of cftDNA in CSF without any prior knowledge of point mutations present in the tumor. We also showed that the fragmentation pattern of cell-free DNA in CSF is different from that in plasma. This low-cost screening method provides information on the tumor genome and can be used to target those patients with high levels of cftDNA for further larger-scale sequencing, such as by whole-exome and whole-genome sequencing.

Study EGAS00001003255

CCL22 chemokine mutations drive natural killer cell lymphoproliferative disease by biasing GPCR signaling

Chronic lymphoproliferative disorder of natural killer cells (CLPD-NK) is characterized by clonal expansion of NK cells of incompletely understood genomic basis. Here, we report somatic mutations in the chemokine gene CCL22 as the hallmark of a distinct subset of CLPD-NK. CCL22 mutations were enriched at highly conserved residues, were mutually exclusive of STAT3 mutations, and were associated with gene expression similar to normal CD16dim/CD56bright NK cells. Mechanistically, the mutations resulted in ligand-biased chemokine receptor signaling, with decreased internalization of the G-protein-coupled receptor for CCL22, CCR4, via impaired β-arrestin recruitment. This resulted in increased cell chemotaxis in vitro, bidirectional crosstalk with the hematopoietic microenvironment, and enhanced NK cell proliferation in vivo in transgenic human IL-15 mice. Somatic CCL22 mutations illustrate a novel mechanism of tumor formation in which gain-of-function chemokine mutations promote tumorigenesis by biased G-protein-coupled receptor signaling and dysregulation of microenvironmental crosstalk.

Study EGAS00001006009

ENCORAFENIB COMBINED WITH BINIMETINIB FOR BRAFV600E-MUTATED RELAPSED/REFRACTORY MULTIPLE MYELOMA: THE PHASE II GMMG-BIRMA TRIAL (Hipo_K08K)

Activating BRAF mutations are rare events in multiple myeloma but have been shown to be promising therapeutic targets in small case studies. This multicenter phase II trial assessed the efficacy and safety of the BRAF/MEK inhibitors, encorafenib and binimetinib, in patients with relapsed/refractory multiple myeloma (RRMM) carrying an activating BRAFV600E mutation. Twelve RRMM patients with a median of five prior lines of therapy were enrolled. The study met its primary endpoint with 10 patients achieving at least a partial response according to IMWG criteria (overall response rate 83.3%). Responses occurred rapidly within the first cycle and deepened over time. The median progression-free survival (PFS) was 5.6 months and overall survival was 55% at 24 months. Genomic profiling revealed RAS mutations and structural variants involving the BRAF locus to be drivers of resistance to BRAF/MEK inhibition in RRMM. This trial demonstrates the value of translationally driven clinical trials in selected patient populations.

Study EGAS00001005973

22929 results for "ega"

in 30.88 milliseconds.