Search Results - EGA European Genome-Phenome Archive

High Resolution Maps of the HeLa 3D Genome Using Hi-C

We use in situ Hi-C to probe the three-dimensional architecture of genomes, constructing haploid and diploid maps of nine cell types (including HeLa). The densest, in human lymphoblastoid cells, contains 4.9 billion contacts, achieving 1-kilobase resolution. We find that genomes are partitioned into contact domains (median length, 185kb), which are associated with distinct patterns of histone marks and segregate into six subcompartments. We identify ~10,000 loops. These loops frequently link promoters and enhancers, correlate with gene activation, and show conservation across cell types and species. Loop anchors typically occur at domain boundaries and bind CTCF. CTCF sites at loop anchors occur predominantly (>90%) in a convergent orientation, with the asymmetric motifs 'facing' one another. The inactive X-chromosome splits into two massive domains and contains large loops anchored at CTCF-binding repeats. The data for the other (non-HeLa) cell types has been uploaded to GEO.

Study phs001010

Genomic Landscape of Apical Periodontitis

We undertook the first genome-wide association (GWAS) study of apical periodontitis, enrolling 932 patients with deep carious lesions at UTHealth School of Dentistry and the University of Pittsburgh School of Dental Medicine. We used a case-control design to compare those with apical periodontitis to those without. Genotype data was generated for all participants and imputed using the TOPMed reference panel. We performed a single-variant association testing, a subphenotype analysis, and a genetically regulated expression association analysis, and generated a polygenic risk score and performed phenome-wide association analysis. We identified eight near-genome-wide significant loci for our primary outcome, eight genes where genetically regulated expression was associated with apical periodontitis status, and 33 phecodes that were associated with the generated polygenic risk score. We also replicated three loci previously reported in the literature. Genotype, phenotype, and results from the primary association analysis will be available via dbGaP.

Study phs003252

Pooled Mutant KRAS-Targeted Long Peptide Vaccine Combined with Nivolumab and Ipilimumab for Patients with Resected MMR-p Colorectal and Pancreatic Cancer

We developed INtraCEllular Reverse Transcription with Sorting and sequencing (INCERTS), an approach to combine molecular indexing of receptor repertoires within intact cells and fluorescence-activated cell sorting (FACS). We demonstrate that INCERTS enables efficient processing of millions of cells from pooled human peripheral blood mononuclear cell (PBMC) samples while retaining robust association between T cell receptor (TCR) sequences and cellular phenotypes. We used INCERTS to discover antigen-specific TCRs from patients with cancer immunized with a novel mutant KRAS peptide vaccine. After ex vivo stimulation, 28 uniquely barcoded samples were pooled prior to FACS into peptide-reactive and non-reactive CD4+ and CD8+ populations. Combining complementary patient-matched single-cell RNA sequencing (scRNA-seq) data enabled retrieval of full-length, paired TCR alpha and beta chain sequences for future validation of therapeutic utility. Data available through dbGaP include scRNA-seq with matched TCR-seq.

Study phs003425

Efficient prediction of a spatial transcriptomics profile better characterizes breast cancer tissue sections without costly experimentation

Spatial transcriptomics is an emerging technology requiring costly reagents and considerable skills, limiting the identification of transcriptional markers related to histology. Here, we show that predicted spatial gene-expressions in unmeasured regions and tissues can enhance biologists�� histological interpretations. We developed the Deep learning model for Spatial gene Clusters and Expression, DeepSpaCE and confirmed its performance using the spatial-transcriptome profiles and immunohistochemistry images of consecutive human breast cancer tissue sections. For example, the predicted expression patterns of SPARC, an invasion marker, highlighted a small tumor-invasion region that is difficult to identify using raw data of spatial transcriptome alone because of a lack of measurements. We further developed semi-supervised DeepSpaCE using unlabeled histology images and increased the imputation accuracy of consecutive sections, enhancing applicability for a small sample size. Our method enables users to derive hidden histological characters via spatial transcriptome and gene annotations, leading to accelerated biological discoveries without additional experiments.

Study JGAS000290

DNA methylation and transposable element landscapes define human regulatory T cells in tissues and identify their blood recirculating counterpart

Regulatory T cells (Treg cells) in tissues can perform immunoregulatory and tissue-regenerative functions. Despite these important tasks, human-tissue Treg cells lack understanding, specifically, the epigenetic basis of their differentiation as well as re-circulating tissue-Treg cells are incompletely understood. Here, we performed genome-wide DNA methylation analysis of human Treg cells from skin and blood and embedded this in a multi-omic epigenetic framework including chromatin accessibility and gene-expression. We thereby revealed novel epigenetic programs defining human skin Treg cell differentiation and delineated how these are present on the different epigenetic levels. We discovered that transposable elements are an integral constituent of the demethylation landscape in skin Treg cells. Based on TCR-sequencing and DNA demethylation homologies, our data unveil that blood CCR8+ Treg cells are recirculatory human tissue-Treg cells. Our data provide novel insights in human tissue-Treg biology, required to harness these cells for therapeutic avenues.

Study EGAS50000000738

T cell receptor repertoire in cell-free DNA as a proxy for tumor infiltrates in patients treated with pembrolizumab

Profiling tumour and systemic T-cell receptor (TCR) repertoires can identify features associated with clinical benefit of cancer patients on immune checkpoint blockade therapy. This study reports on diversity and specificity analysis of TCR sequences from 59 tumor, 306 peripheral blood (PBMC), and 73 cell-free DNA (cfDNA) samples from 81 patients with solid tumors treated with pembrolizumab. This analysis reveals that head and neck carcinomas have significantly lower PBMC diversity and shorter persistence of pembrolizumab-induced diversification compared to other cancers. Using GLIPHII (Grouping Lymphocyte Interactions with Paratope Hotspots), we find that tumor-derived clonotypes form non-microbial specificity signatures that are shared across patients and cancer types. Tumor-derived clonotypes are enriched in cfDNA compared to PBMCs, despite the diversity of cfDNA TCR repertories being two orders of magnitude lower than PBMCs. TCRs shared between tumors and cfDNAs persist in PBMCs beyond 50 weeks of treatment, suggesting cfDNA captures clonotypes contributing to anti-tumor immune surveillance.

Study EGAS50000001306

Single cell sequencing data of PBMC and CSF from a cohort of Multiple Sclerosis patients and other neurological disease controls

Neuroinflammation is often characterised by immune cell infiltrates in the cerebrospinal fluid (CSF). Here we apply single-cell RNA sequencing to explore the functional characteristics of these cells in patients with various inflammatory, infectious and non-inflammatory neurological disorders. We show that CSF is distinct from the peripheral blood in terms of both cellular composition and gene expression. We report that the cellular and transcriptional landscape of CSF is altered in neuroinflammation, but is strikingly similar across different neuroinflammatory disorders. We find clonal expansion of CSF B and T cells in all disorders but most pronounced in inflammatory diseases, and we functionally characterise the transcriptional features of these cells. Finally, we explore the genetic control of gene expression in CSF lymphocytes. Our results highlight the common features of immune cells in the CSF compartment across diverse neurological diseases and may help to identify new targets for drug development or repurposing in Multiple Sclerosis.

Study EGAS50000000739

Human Brain Small Extracellular Vesicles Contain Selectively-Packaged, Full-Length mRNA

Brain cells release and take up small extracellular vesicles (sEVs) containing bioactive nucleic acids. sEV exchange is hypothesized to contribute to stereotyped spread of neuropathological changes in the diseased brain. We assessed mRNA from sEVs of non-diseased (ND) and Alzheimer’s disease (AD) human postmortem brain, using short- and long-read sequencing. sEV transcriptomes were distinct from bulk tissue, showing enrichment for multiple genes including mRNAs encoding ribosomal proteins and L1Hs transposable elements. AD versus ND sEVs showed enrichment of inflammation-related and depletion of synaptic signaling mRNAs. sEV mRNA from cultured murine primary neurons, astrocytes, or microglia showed similarities with human brain sEVs and revealed cell-type specific packaging. Nearly 80% of human brain sEV transcripts sequenced using long-read sequencing were full-length. Motif analyses of sEV-enriched full-length isoforms revealed RNA-binding proteins that may be associated with sEV loading. Collectively, we show that mRNA in brain sEVs is selectively-packaged and altered by disease state.

Study EGAS50000000029

Uterine_Atlas_Endometriosis

Uterine Atlas Endometriosis Endometriosis is a chronic inflammatory disease driven by oestrogen, affecting around 10% of women of child-bearing again world-wide. It is described by the presence of endometrium-like tissue outside of the uterus (ectopic endometriosis tissue) and debilitating chronic pain. There is a lack of understanding of endometriosis pathogenesis, as well as cellular composition of both the eutopic endometrium and ectopic endometriosis tissue of endometriosis patients.There are therefore two main aims for this study. Using single cell genomics and spatial methods we aim to: I) characterise endometrial cellular census in women with and without endometriosis II) characterise cellular census of ectopic endometriosis tissue Furthermore, we aim to establish endometrial organoids and endometriosis organoids from the samples collected and characterise these using single cell genomics and spatial methods. 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 EGAS00001004725

Transcriptome analysis in very preterm infants with chronic lung disease after birth

Bronchopulmonary dysplasia (BPD) is the most common complication of prematurity and is characterized by impaired alveolar and vascular development. The aim of the study was to identify early BPD-relevant pathways by applying whole genome transcriptional profiling combined with protein measurements from umbilical arterial blood samples of preterm infants with and without BPD. Only newborn infants <32 weeks gestational age (GA) were prospectively included in this study. Blood was obtained from an indwelling umbilical artery catheter at birth and after 72 hours after birth and subjected to gene expression analysis using the Codelink Human I 10k Bioarray and the the Codelink Human Whole Genome Bioarray. The study reveals a monocyte-centered immune response at birth characterizing preterm infants later developing BPD. The findings provide potential markers for early risk stratification in the preterm infant with respect to the development of chronic pulmonary disease and may subsequently enable the development of new treatment strategies in this high risk patient cohort.

Study EGAS00001002586

Single cell RNA sequencing and Whole Genome Sequencing on different cells from the same sample for a triple negative patient derived xenograft and ovarian cancer cell lines.

Measuring gene expression of genomically defined tumour clones at single cell resolution would associate functional consequences to somatic alterations, as a prelude to elucidating pathways driving cell population growth, resistance and relapse. In the absence of scalable methods to simultaneously assay DNA and RNA from the same single cell, independent sampling of cell populations for parallel measurement of single cell DNA and single cell RNA must be computationally mapped for genome-transcriptome association. Here we present clonealign, a robust statistical framework to assign gene expression states to cancer clones using single-cell RNA-seq and DNA-seq independently sampled from an heterogeneous cancer cell population. We apply clonealign to triple-negative breast cancer patient derived xenografts and high-grade serous ovarian cancer cell lines and discover clone-specific dysregulated biological pathways not visible using either DNA-Seq or RNA-Seq alone.

Study EGAS00001003387

Reference epigenomes generated as part of the International Human Epigenomics Consortium (IHEC)

The Centre for Epigenome Mapping Technologies (CEMT) is an epigenome sequencing platform funded by the Canadian Institutes of Health Research (CIHR) and Genome BC as part of the Canadian Epigenetics, Environment and Health Research Consortium (CEEHRC). CEMT aims to produce 100 human reference epigenomes, each comprising the following datasets: whole genome; whole genome bisulfite and oxidative bisulfite; transcriptome; microRNA; and chromatin immunoprecipitation (ChIP)-seq for the core chromatin marks (H3K4me1, H3K4me3, H3K9me3, H3K27me3, H3K36me3, and H3K27ac).CEEHRC is a full member of the International Human Epigenomics Consortium (IHEC), which aims to coordinate the production of 1,000 reference maps of human epigenomes for key cellular states relevant to health and diseases. As such, all datasets submitted by the CEMT platform conform to IHEC assay standards and metadata guidelines.

Study EGAS00001000552

DONSON encodes a novel replication fork protection factor mutated in microcephalic dwarfism.

To ensure efficient genome duplication, cells have evolved a multitude of factors that promote unperturbed DNA replication, and protect, repair and restart damaged forks. Here we identify DONSON as a novel fork protection factor, and report biallelic DONSON mutations in 29 individuals with microcephalic dwarfism. We demonstrate that DONSON is a component of the replisome that stabilises forks during normal genome replication. Loss of DONSON leads to severe replication-associated DNA damage arising from nucleolytic cleavage of stalled replication forks. Furthermore, ATR-dependent signalling in response to replication stress is impaired in DONSON-deficient cells, resulting in decreased checkpoint activity, and potentiating chromosomal instability. Hypomorphic mutations substantially reduce DONSON protein levels and impair fork stability in patient cells, consistent with defective DNA replication underlying the disease phenotype. In summary, we identify mutations in DONSON as a common cause of microcephalic dwarfism, and establish DONSON as a critical replication fork protein required for mammalian DNA replication and genome stability.

Study EGAS00001002224

X chromosomal genetic variants are associated with childhood obesity

Current genetic association studies are usually focused on autosomal variants only, and the sex chromosomes are often neglected. In recent years, a number of statistical techniques and strategies have been widely described making much easier overcoming X-chromosome technical hurdles and including this region within genetic studies. Tenomodulin (TNMD) is a Xq22 chromosome anti angiogenic locus which has been recently linked to different obesity-related phenotypes. These results have not been replicated to date. Given these facts, we have conducted a genetic association analysis in Spanish children population including seven TNMD SNPs as potential candidate markers for obesity and metabolic dysfunctions. Additionally genotypes for another locus located in the X chromosome, the SLC6A14, have been included in the dataset. A total of 915 DNA samples from 258 normal weight, 177 overweight and 480 obese Spanish children (438 males, 477 females) were genotyped for seven TNMD SNPs and one SLC6A14 SNP. Associations with anthropometric measurements and glucose metabolism were investigated.

Study EGAS00001002738

Exome Sequencing of Gastric Cancer

GASTRIC cancer (GC) is a heterogeneous disease with multiple environmental etiologies and alternative pathways of carcinogenesis. Beyond mutations in TP53, alterations in other genes or pathways account only for small subsets of the disease. We performed exome sequencing of 22 GCs and identified mutated genes and pathway alterations previously unreported: genes involved in chromatin modification were commonly mutated. A downstream validation study confirmed frequent inactivating mutations and/or protein deficiency of the SWI/SNF chromatin remodeling complex gene ARID1A in 83% of GCs with microsatellite instability, 73% of GCs with Epstein-Barr virus (EBV) infection, and 11% of non-EBV microsatellite stable GCs. The mutation spectrum of ARID1A differs between molecular subtypes and the mutation prevalence is negatively associated with mutations in TP53. Clinically, ARID1A alterations were associated with better prognosis in a stage-independent manner. These results revealed the genomic landscape and highlighted the importance of chromatin remodeling in the molecular taxonomy of GCs.

Study EGAS00001000153

Molecular Subtype-specific Biomarkers Improves Colorectal Cancer Prognostication

Colorectal cancer (CRC) is characterized by major inter-tumor diversity that complicates the prediction of disease and treatment outcomes. Recent efforts help resolve this by sub-classification of CRC according to natural molecular subtypes, however, this strategy is not yet able to provide clinicians with improved tools for decision-making. We here present an extended framework for CRC stratification that specifically aims to improve patient prognostication. Using transcriptional profiles from 1,100 CRCs, including a novel set of >300 samples, we identify novel cancer cell and tumor archetypes and suggest the tumor microenvironment as a major prognostic determinant that can be influenced by the microbiome. Notably, our subtyping strategy allowed identification of novel archetype-specific prognostic biomarkers that provided information beyond and independent of UICC-TNM staging, MSI-status and consensus molecular subtyping. The results illustrate that our extended subtyping framework, combining subtyping and subtype-specific biomarkers, could contribute to improved patient prognostication and may form a strong basis for future studies.

Study EGAS00001002376

Orthotopic Patient-Derived Xenografts of Pediatric Solid Tumors

Pediatric solid tumors arise from endodermal, ectodermal, or mesodermal lineages. Although the overall survival of children with solid tumors is 75%-80%, that of children with recurrent disease is below 30%. To capture the complexity and diversity of pediatric solid tumors and establish new models of recurrent disease, we developed a protocol to produce orthotopic patient-derived xenografts (O-PDXs) at diagnosis, recurrence, and autopsy. Tumor specimens were received from 168 patients, and 64 O-PDXs were established for 12 types of pediatric solid tumors. The origins of the O-PDX tumors were reflected in their gene-expression profiles and epigenomes. Genomic profiling of the tumors, including detailed clonal analysis, was performed to determine whether the clonal population in the xenograft recapitulated the patient’s tumor. We identified several drug vulnerabilities and showed that the combination of a WEE1 inhibitor (AZD1775), irinotecan, and vincristine can lead to complete response in multiple rhabdomyosarcoma O-PDX tumors in vivo.

Study EGAS00001002528

A biobank of patient-derived pediatric brain tumor models

Brain tumors are the leading cause of cancer-related death in children. Genomic studies have provided insights into molecular subgroups and oncogenic drivers of pediatric brain tumors that may lead to novel therapeutic strategies. To evaluate new treatments, better preclinical models adequately reflecting the biological heterogeneity are needed. Through the Children’s Oncology Group ACNS02B3 study, we have generated and comprehensively characterized 30 patient-derived orthotopic xenograft models and 7 cell lines representing 14 molecular subgroups of pediatric brain tumors. Patient-derived orthotopic xenograft models were found to be representative of the human tumors they were derived from in terms of histology, immunohistochemistry, gene expression, DNA methylation, copy number, and mutational profiles. In vivo drug sensitivity of targeted therapeutics was associated with distinct molecular tumor subgroups and specific genetic alterations. These models and their molecular characterization provide an unprecedented resource for the cancer community to study key oncogenic drivers and to evaluate novel treatment strategies.

Study EGAS00001002536

Human Pancreatic Beta Cell lncRNAs Control Cell-Specific Regulatory Networks

Recent studies have uncovered thousands of long non-coding RNAs (lncRNAs) in human pancreatic β cells. β cell lncRNAs are often cell type specific and exhibit dynamic regulation during differentiation or upon changing glucose concentrations. Although these features hint at a role of lncRNAs in β cell gene regulation and diabetes, the function of β cell lncRNAs remains largely unknown. In this study, we investigated the function of β cell-specific lncRNAs and transcription factors using transcript knockdowns and co-expression network analysis. This revealed lncRNAs that function in concert with transcription factors to regulate β cell-specific transcriptional networks. We further demonstrate that the lncRNA PLUTO affects local 3D chromatin structure and transcription of PDX1, encoding a key β cell transcription factor, and that both PLUTO and PDX1 are downregulated in islets from donors with type 2 diabetes or impaired glucose tolerance. These results implicate lncRNAs in the regulation of β cell-specific transcription factor networks.

Study EGAS00001002865

Single_cell_analysis_T_cell_activation

We isolated naïve and memory T cells from healthy donors and cultured them in resting and stimulated conditions. We performed scRNA sequencing to assess T cell actiavtaion and perform sceQTL analysis. Major challenges for translating GWAS variants to function include identification of the disease relevant cell types in which the associated variants function. We have shown that immune-disease variants are functional in CD4+ T cell activation, and understanding how GWAS variants affect the activation efficiency is critical for development of new therapies. To perform sceQTL analysis of naïve and memory T cell activation, here we stimulated T cells and profiled single cell transcriptome at different time points following stimulation. We mapped sceQTLs affecting cell activation. 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 EGAS00001003479

Genetic_profiling_of_mucosal_melanoma

Primary mucosal melanomas (MMs) arise from melanocytes located in mucosal membranes lining the respiratory, gastrointestinal and urogenital tracts. MMs frequently present late and have a poor prognosis; the 5-year survival rate is only 14%. MM makes up only ~1.4% of all melanomas and it is this rarity that makes knowledge of the genetic changes that contribute to its pathogenesis limited to a small number of exome/genome studies and other targeted studies. Thus to investigate the somatic alterations and mutation spectra in MM genomes, we have extracted genomic DNA from formalin-fixed, paraffin-embedded (FFPE) human MMs, and subjected them to whole exome sequencing. Given the propensity of MM to metastasize, we will also be sequencing metastatic MM lesions; primary and metastatic lesions from the same individual represent an excellent opportunity to identify potential drivers of metastasis in MM. Finally we will sequence ‘normal’ DNA from the same individual, where possible, to exclude germline variations.

Study EGAS00001001115

Adult granulosa cell tumor WGS data cohort with corresponding reference germline WGS data

Adult granulosa cell tumors (AGCTs) harbor a somatic FOXL2 c.402C>G mutation in ~95% of cases and are mainly surgically removed due to limited systemic treatment effect. In this study, potentially targetable genomic alterations in AGCTs were investigated by whole genome sequencing on 46 tumor samples and matched normal DNA. Copy number variant (CNV) analysis confirmed gain of chromosome 12 and 14, and loss of 22. Pathogenic TP53 mutations were identified in three patients with highest tumor mutational burden and mitotic activity, defining a high-grade AGCT subgroup. Within-patient tumor comparisons showed 29-80% unique somatic mutations per sample, suggesting tumor heterogeneity. FOXL2-wildtype AGCTs had DICER1, TERT(C228T) and TP53 mutations and similar CNV profiles as FOXL2-mutant tumors. Our study confirms that absence of the FOXL2 c.402C>G mutation does not exclude AGCT diagnosis. The lack of overlapping variants in targetable cancer genes, indicates the need for personalized treatment for AGCT patients

Study EGAS00001004249

A Universal Gut Metagenomic-Derived Signature Predicts Cirrhosis

Dysregulation of the gut microbiome has been implicated in the progression of nonalcoholic fatty liver disease (NAFLD) to advanced fibrosis and cirrhosis. To determine the diagnostic capacity of this association, stool microbiomes were compared across 163 well-characterized participants encompassing non-NAFLD controls, NAFLD-cirrhosis patients and their first-degree relatives. Interrogation of shotgun metagenomic and untargeted metabolomic profiles using the Random Forest machine learning algorithm and differential abundance analysis identified discrete metagenomic and metabolomic signatures that were similarly effective in detecting cirrhosis (diagnostic accuracy 0.91, AUC). Combining the metagenomic signature with age and serum albumin levels accurately distinguished cirrhosis in etiologically and genetically distinct cohorts from geographically separated regions. Additional inclusion of serum aspartate aminotransferase levels, which are increased in cirrhosis patients, enabled discrimination of cirrhosis from earlier stages of fibrosis. These findings demonstrate that a core set of gut microbiome species may offer universal utility as a non-invasive diagnostic test for cirrhosis.

Study EGAS00001004600

Genome-wide quantification of rare somatic mutations in normal human tissues using massively parallel sequencing

We present Bottleneck Sequencing System (BotSeqS), a nextgeneration sequencing method that simultaneously quantifies rare somatic point mutations across the mitochondrial and nuclear genomes. BotSeqS combines molecular barcoding with a simple dilution step immediately prior to library amplification. We use BotSeqS to show age and tissue-dependent accumulations of rare mutations and demonstrate that somatic mutational burden in normal tissues can vary by several orders of magnitude, depending on biologic and environmental factors. We further show major differences between the mutational patterns of the mitochondrial and nuclear genomes in normal tissues. Lastly, the mutation spectra of normal tissues were different from each other, but similar to those of the cancers that arose in them. This technology can provide insights into the number and nature of genetic alterations in normal tissues and can be used to address a variety of fundamental questions about the genomes of diseased tissues.

Study EGAS00001001838

DNA Methylation loss coupled with mitotic cell division promotes immune evasion of tumours with high mutation load

Mitotic cell division increases tumour mutation burden and copy number load with a positive and inverse correlation, respectively, with the clinical benefit of immunotherapy. Markers of cell division correlate also with genomic demethylation involving methylation loss in late-replicating partial methylation domains. Here we find that immunomodulatory pathway genes are concentrated in these domains and transcriptionally repressed in demethylated tumours with CpG island promoter hypermethylation. Global methylation loss correlated with immune evasion signatures independently of mutation burden and aneuploidy. Methylome data of our cohort (n = 60) and a published cohort (n = 81) in lung cancer and a melanoma cohort (n = 40) consistently demonstrated that genomic methylation alterations counteract the contribution of high mutation burden and increase immunotherapeutic resistance. Higher predictive power was observed for methylation loss than mutation burden. We also found that genomic hypomethylation correlates with the immune escape signatures of aneuploid tumours. Hence, DNA methylation alterations implicate epigenetic modulation as a combination regimens for precision immunotherapy.

Study EGAS00001003731

16787 results for "sec+and+unspecified+malignant+neoplasm"

in 19.17 milliseconds.