Search Results - EGA European Genome-Phenome Archive

Epithelioid haemangioendothelioma (EHE) case series from the Stafford Fox Rare Cancer Program

Nine EHE cases with comprehensive histologic and molecular profiling from the Walter and Eliza Hall Institute of Medical Research Stafford Fox Rare Cancer Program were collated via nation-wide referral to the Australian Rare Cancer Portal. Molecular analysis of RNA, DNA or both was possible in seven of nine cases. Survival was highly variable and unrelated to molecular profile. DNA sequences for two cases are deposited here.

Study EGAS00001007474

Genome-wide association study of esophageal squamous cell cancer identifies shared and distinct risk variants in African and Chinese populations

Study EGAS00001007477

Whole Genome Sequence and RNASeq Samples for Lung Cancer

Short Read Whole Genome Sequence and RNASeq Samples for Lung Cancer

Study EGAS00001007832

Discrepancies in Tumour Mutation Burden (TMB) reporting from sequential Endobronchial ultrasound trans bronchial needle aspiration (EBUS TBNA) samples within single lymph node stations

Introduction: Tumour Mutation Burden (TMB) is a potential biomarker for immune therapies. Identifying factors that may affect TMB prediction is key to its utility as biomarker for treatment options. Here we investigated parameters that might affect TMB using cytology smears obtained from endobronchial ultrasound transbronchial needle aspiration (EBUS TBNA)-sampled malignant lymph nodes, which are often the only diagnostic samples collected in patients with advanced lung cancer. Methods: Individual Diff-Quik cytology smears were prepared for each needle pass from the same enlarged malignant mediastinal lymph nodes. DNA extracted from smears of each pass underwent sequencing using a large gene panel sequencing (TSO500, Illumina). TMB was estimated for each individual needle pass using the TSO500 Local App v. 2.0 (Illumina). Results: Twenty patients had two or more Diff-Quik smears (total 45 smears) which passed sequencing quality control. Average TMB for all smears was 8.7 ± 5.0 mutations per megabase (Mb). Sixteen of the 20 patients had paired samples with minimal differences in TMB score (average difference 1.3 ± 0.85). Paired samples from 13 patients had concordant TMB in terms of scores being below or above a threshold of 10 mutations/Mb. Samples from seven patients had discrepant TMB, with four cases in particular exhibiting an average difference of 11.3 ± 2.7 mutations/Mb. Factors affecting TMB calling involved tumour content of the samples, the amount of DNA used in sequencing as well as bone fide intra tumour heterogeneity between paired samples. Conclusions: TMB assessment is feasible from EBUS-TBNA derived Diff-Quik slides from a single needle pass. Repeated samples of one or more lymph nodes have minimal variation in TMB in most cases (80%), suggesting accuracy of the test overall. To enhance accuracy of the result, care should be taken in relation to the tumour content of the cytology smear. Further, clinicians should be aware that even slight differences in lymph node site selection could give rise to variable genomics and discrepant TMB due to tumour heterogeneity. Therefore, it would be reasonable for multi-lymph node site analysis of TMB to be performed using EBUS TBNA aspirates from different easily accessible lymph nodes. Further investigations need to determine whether these samples should be combined as one sample or tested separately.

Study EGAS00001007484

Clinical and genetic factors associated with tumor response to neoadjuvant (chemo)radiotherapy, survival and recurrence risk in rectal cancer

Rectal cancer poses challenges in preoperative treatment response, with up to 30% achieving a complete response (CR). Personalized treatment relies on accurate identification of responders at diagnosis. This study aimed to unravel CR determinants, overall survival (OS), and time to recurrence (TTR) using clinical and targeted sequencing data. Analyzing 402 patients undergoing preoperative treatment, tumor stage, size, and treatment emerged as robust response predictors. CR rates were higher in smaller, early-stage, and intensively treated tumors. Targeted sequencing analyzed 216 cases, while 120 patients provided hotspot mutation data. KRAS mutation dramatically reduced CR odds by over 50% (odds ratio [OR]=0.3 in the targeted sequencing and OR=0.4 hotspot cohorts, respectively). In contrast, SMAD4 and SYNE1 mutations were associated with higher CR rates (OR=6.0 and 6.8, respectively). Favorable OS was linked to younger age, CR, and low baseline carcinoembryonic antigen levels. Notably, CR and an APC mutation increased TTR, while a BRAF mutation negatively affected TTR. Beyond tumor burden, SMAD4 and SYNE1 mutations significantly influenced CR. KRAS mutations independently correlated with radiotherapy resistance, and BRAF mutations heightened recurrence risk. Intriguingly, nonresponding tumors with initially small sizes carried a higher risk of recurrence. These findings offer insights into rectal cancer treatment response, guiding personalized therapeutic strategies. By uncovering factors impacting CR, OS, and TTR, this study underscores the importance of tailored approaches for rectal cancer patients. These findings, based on extensive analysis and mutation data, pave the way for personalized interventions, optimizing outcomes in the challenges of rectal cancer preoperative treatment.

Study EGAS00001007501

RNASeq Neoantigen Immunogenicity Landscape

Immune checkpoint blockade (ICB) therapy is a cornerstone of oncologic treatment for patients with advanced stage non-small cell lung cancer (NSCLC) and other malignancies. Neoantigen immunoediting drives ICB efficacy, yet the fundamental physiochemical characteristics of neoantigens and how neoantigen immunogenicity shapes treatment response remains poorly understood. To help address these questions, a prospective clinical trial of NSCLC patients treated with nivolumab was conducted. We assessed genomic alterations in tumors from 58 patients and performed large-scale neoepitope immunogenicity analyses, before and during treatment (CheckMate153, CA209-153). Tumors were analyzed by whole-exome and transcriptome sequencing. In responding patients, loss of mutation and neoantigen burden early during therapy associated with clinical benefit. We evaluated the immunogenicity of 1,453 candidate neoantigens and identified 502 neopeptides that bound to MHC I and 196 neopeptides that were immunogenically recognized by T cells in the setting of nivolumab treatment. These T cell reactive neoantigens were differentially present in clonal populations that underwent distinctive evolutionary trajectories across responders and nonresponders. Mapping these neoantigens to tumor clonal dynamics and clinical response revealed strong selection against immunogenic neoantigen harboring clones compared to non-immunogenic clones. Using this large collection of neoantigens, we identified position specific amino acid features related to immunogenicity, which we used to develop and validate an immunogenicity score. Changes in the genomic and neoantigen immunogenicity landscapes were accompanied by temporal changes in the tumor microenvironment. Nivolumab-induced microenvironmental evolution in NSCLC shared some similarities with that in melanoma, yet critical differences in immunologic programs were apparent from comparative network analysis between tumor types. This study provides unprecedented molecular portraits of the neoantigen landscapes underlying nivolumab’s mechanism of action.

Study EGAS00001007509

RNASeq files for GenomePaint paper

RNASeq files for GenomePaint paper titled "Exploration of coding and non-coding variants in cancer using GenomePaint."

Dataset EGAD00001006680

Machine learning guided signal enrichment for ultrasensitive plasma tumor burden monitoring

In solid tumor oncology, circulating tumor DNA (ctDNA) is poised to transform care through accurate assessment of minimal residual disease (MRD) and therapeutic response monitoring. To overcome the sparsity of ctDNA fragments in low tumor fraction (TF) settings and increase MRD sensitivity, we previously leveraged genome-wide mutational integration through plasma whole genome sequencing (WGS). We now introduce MRD-EDGE, a composite machine learning-guided WGS ctDNA single nucleotide variant (SNV) and copy number variant (CNV) detection platform designed to increase signal enrichment. MRD-EDGESNV uses deep learning and a ctDNA-specific feature space to increase SNV signal-to-noise enrichment in WGS by 300X compared to previous WGS error suppression. MRD-EDGECNV also reduces the degree of aneuploidy needed for ultrasensitive CNV detection through WGS from 1 Gb to 200 Mb, thereby expanding its applicability to a wider range of solid tumors. We harness the improved performance to identify MRD following surgery in multiple cancer types, track changes in tumor burden in response to neoadjuvant immunotherapy in non-small cell lung cancer (NSCLC) and demonstrate ctDNA shedding in precancerous colorectal adenomas. Finally, the radical signal-to-noise enrichment in MRD-EDGESNV enables de novo mutation calling in melanoma and small cell lung cancer (SCLC) without matched tumor, yielding clinically informative TF monitoring for patients on immune checkpoint inhibition (ICI).

Study EGAS00001007545

Immune and clinicopathological features predict HER2-positive breast cancer prognosis in the neoadjuvant NeoALTTO and CALGB 40601 trials

The identification of prognostic drivers in patients receiving neoadjuvant therapy is crucial for treatment optimization in HER2-positive breast cancer, with the immune microenvironment being a key factor. Here, we investigated the complexity of B and T cell receptor (BCR and TCR) repertoires, reflecting antigen-specific immune response, in the context of two phase III trials, namely NeoALTTO and CALGB 40601, evaluating neoadjuvant paclitaxel with trastuzumab and/or lapatinib in women affected by HER2-positive breast cancer. BCR/TCR features were heterogeneous, and BCR measures describing clonal expansion demonstrated independent prognostic value. A model developed in NeoALTTO and validated in CALGB 40601 predicted event-free survival (EFS) by integrating hormone receptor and clinical nodal status, treatment response (breast pCR)breast pathological complete response (pCR) as well as tumor-infiltrating lymphocytes and BCR evenness, measuring the equality of BCR clonal distribution. The HER2-EveNT score derived from the model allows the identification of patients with 5-year EFS > 90%, and, in those not achieving pCR, of a subgroup of immune-enriched tumors with an excellent outcome despite residual disease.

Study EGAS00001007563

Idiosyncratic and generic single nuclei and spatial transcriptional patterns in papillary and anaplastic thyroid cancers

Thyroid cancer is the most common endocrine malignancy. This dataset encompasses two types of thyroid cancer : anaplastic which is the most de-differentiated and aggressive one, and papillary which is the most common one. We profiled 14 patients, including 10 papillary and 4 anaplastic thyroid carcinomas, using both single nuclei RNA sequencing and spatial transcriptomics to link single cell resolution RNA sequencing with tissue morphology and better understand inter and intratumoral thyroid cancer heterogeneity.

Study EGAS00001007574

DERMATLAS__Leiomyoma_WES

Virtually no other tumour type is associated with so many different forms as skin cancer. Histologically, tumours of the skin may arise from epithelium, including epidermis, hair follicle, sebaceous or sweat gland, melanocytes, dermal-associated mesenchymal structures or tissue resident immune cells, making for a diversity of clinical presentations. Importantly, many skin tumour types have an extremely poor prognosis. Many skin tumour subtypes have never undergone molecular profiling, or if they have, targeted sequencing has been used and the number of cases analyzed has been so limited that deriving firm conclusions about the profile of driver genes, DNA mutational signatures and germline alleles has not been possible For 70 key skin tumour subtypes defined by the World Health Organization, that have not undergone extensive genetic analysis previoulsy, we propose to perform whole exome and transcriptome sequencing, from a range of body sites, to build a genomic atlas of dermatological tumours, including detailed maps of SNVs, copy number alterations, genome-wide methylation and expression profiles. The cases in this cohort are leiomyomas.

Study EGAS00001007629

germline variants in children with hematological cancer

Trio parent-child whole exome sequencing has been used to identify cancer predisposing syndrome (CPS) in our pediatric patients with hematological cancer, and it was identified in 3.8% of the total number of patients(131).

Study EGAS00001006907

DERMATLAS__Leiomyoma_RNAseq

Virtually no other tumour type is associated with so many different forms as skin cancer. Histologically, tumours of the skin may arise from epithelium, including epidermis, hair follicle, sebaceous or sweat gland, melanocytes, dermal-associated mesenchymal structures or tissue resident immune cells, making for a diversity of clinical presentations. Importantly, many skin tumour types have an extremely poor prognosis. Many skin tumour subtypes have never undergone molecular profiling, or if they have, targeted sequencing has been used and the number of cases analyzed has been so limited that deriving firm conclusions about the profile of driver genes, DNA mutational signatures and germline alleles has not been possible For 70 key skin tumour subtypes defined by the World Health Organization, that have not undergone extensive genetic analysis previoulsy, we propose to perform whole exome and transcriptome sequencing, from a range of body sites, to build a genomic atlas of dermatological tumours, including detailed maps of SNVs, copy number alterations, genome-wide methylation and expression profiles. The cases in this cohort are leiomyomas.

Study EGAS00001007630

DERMATLAS__Leiomyosarcoma_RNAseq

Virtually no other tumour type is associated with so many different forms as skin cancer. Histologically, tumours of the skin may arise from epithelium, including epidermis, hair follicle, sebaceous or sweat gland, melanocytes, dermal-associated mesenchymal structures or tissue resident immune cells, making for a diversity of clinical presentations. Importantly, many skin tumour types have an extremely poor prognosis. Many skin tumour subtypes have never undergone molecular profiling, or if they have, targeted sequencing has been used and the number of cases analyzed has been so limited that deriving firm conclusions about the profile of driver genes, DNA mutational signatures and germline alleles has not been possible For 70 key skin tumour subtypes defined by the World Health Organization, that have not undergone extensive genetic analysis previoulsy, we propose to perform whole exome and transcriptome sequencing, from a range of body sites, to build a genomic atlas of dermatological tumours, including detailed maps of SNVs, copy number alterations, genome-wide methylation and expression profiles. The cases in this cohort are leiomyosarcomas.

Study EGAS00001007631

Buccal Sample Methylation for Breast Cancer Detection

This study encompasses buccal sample DNA methylation (IlluminaMethylationEPIC array) from breast cancer cases and healthy age-matched controls.

Study EGAS00001007658

Unraveling mutagenic processes influencing the tumor mutational patterns of individuals with Constitutional Mismatch Repair Deficiency

Constitutional mismatch repair deficiency (CMMRD), caused by bi-allelic germline variants in one of the mismatch repair (MMR) genes, is a childhood cancer predisposition syndrome that often results in the development of multiple tumors early in life. A better understanding of mutational processes driving subsequent CMMRD tumors could advance optimal treatment. Therefore, we performed a genomic characterization of 41 tumors from 17 individuals with CMMRD. Mutational patterns in these tumors were found to be influenced by multiple factors, including the affected MMR genes, tumor types and somatic polymerase proofreading mutations. Temozolomide treatment left prominent mutational signatures in two second primary hematologic malignancies. Furthermore, a novel indel signature was found in 54% of the tumors, characterized by one base pair cytosine insertions in cytosine homopolymers. In conclusion, the analysis of sequential tumors in individuals with CMMRD reveals diverse mutational patterns influenced by the underlying mutated MMR gene, tumor type and treatment history.

Study EGAS00001007660

DNA-methylation variability in normal mucosa of patients with adenomatous polyps: a marker of field cancerization

A study linking DNA methylation modifications in the normal mucosa of the right colon in individuals with adenomatous polyp formation

Study EGAS00001007666

Evaluation of EBUS-TBNA Aspirates from Advanced NSCLC for Comprehensive Sequencing Platforms Including Whole Genome Sequencing

Endobronchial ultrasound-guided transbronchial needle aspiration (EBUS-TBNA) is often the only source of tumor tissue from patients with advanced, inoperable lung cancer. EBUS-TBNA aspirates are used for the diagnosis, staging, and genomic testing to inform therapy options. Here we extracted DNA and RNA from 220 EBUS-TBNA aspirates to evaluate their suitability for whole genome (WGS), whole exome (WES), and comprehensive panel sequencing. For a subset of 40 cases, the same nucleic acid extraction was sequenced using WGS, WES, and the TruSight Oncology 500 assay. Genomic features were compared between sequencing platforms and compared with those reported by clinical testing. A total of 204 aspirates (92.7%) had sufficient DNA (100 ng) for comprehensive panel sequencing, and 109 aspirates (49.5%) had sufficient material for WGS. Comprehensive sequencing platforms detected all seven clinically reported tier 1 actionable mutations, an additional three (7%) tier 1 mutations, six (15%) tier 2–3 mutations, and biomarkers of potential immunotherapy benefit (tumor mutation burden and microsatellite instability). As expected, WGS was more suited for the detection and discovery of emerging novel biomarkers of treatment response. WGS could be performed in half of all EBUS-TBNA aspirates, which points to the enormous potential of EBUS-TBNA as source material for large, well-curated discovery-based studies for novel and more effective predictors of treatment response. Comprehensive panel sequencing is possible in the vast majority of fresh EBUS-TBNA aspirates and enhances the detection of actionable mutations over current clinical testing.

Study EGAS00001007708

Discovery of immunotherapy targets for pediatric solid and brain tumors by exon-level expression

Immunotherapy with chimeric antigen receptor T cells for pediatric solid and brain tumors is constrained by available targetable antigens. Cancer-specific exons present a promising reservoir of targets; however, these have not been explored and validated systematically in a pan-cancer fashion. To identify cancer specific exon targets, we analyze 1,532 RNA-seq datasets from 16 types of pediatric solid and brain tumors for comparison with normal tissues using a newly developed workflow. We find 2,933 exons in 157 genes encoding proteins of the surfaceome or matrisome with high cancer specificity either at the gene (n=148) or the alternatively spliced isoform (n=9) level. Expression of selected alternatively spliced targets, including the EDB domain of fibronectin 1, and gene targets, such as COL11A1, are validated in pediatric patient derived xenograft tumors. We generate T cells expressing chimeric antigen receptors specific for the EDB domain or COL11A1 and demonstrate that these have antitumor activity. The full target list, explorable via an interactive web portal (https://cseminer.stjude.org/), provides a rich resource for developing immunotherapy of pediatric solid and brain tumors using gene or AS targets with high expression specificity in cancer.

Study EGAS00001007766

To_profile_the_landscape_of_sebaceous_tumours___RNA

Sebaceous tumours are a rare cutaneous cancer with potential for aggressive behaviour. However, limited information is available on these cancers with few published cases. Here we wish to exome sequence these cancers to define the first genomic landscape for this malignancy. We will extract RNA from formalin-fixed, paraffin-embedded (FFPE) cores. Cores may be obtained from lesional and non-lesional tissues of primaries as well as matching metastases. The extracted RNA will be used for RNA sequencing.

Study EGAS00001007803

MethylationEPIC BeadChip samples of pre- and post-5ZA-treated head and neck cancer patients refractory to anto-PD-1 therapy

Background Immune checkpoint blockade (ICB) is an effective treatment in a subset of patients diagnosed with head and neck squamous cell carcinoma (HNSCC); however, most patients are refractory. Methods In a non-randomized, open-label phase 1b clinical trial, participants with recurrent and/or metastatic (R/M) HNSCC were treated with low-dose 5’azacytidine (5’aza) daily for either 5- or 10-days in combination with fixed doses of durvalumab and tremelimumab after progression on ICB therapy. The primary objective was to assess the biologically effective dose of 5’aza as determined by molecular changes in paired baseline and on-treatment tumor biopsies; the secondary objective was safety. Results Thirty-eight percent (3/8) of participants with evaluable paired tissue samples had a greater than 2-fold increase from baseline in interferon-gamma (IFN-ɣ) signature and CD274 (programmed cell death protein 1 ligand, PD-L1) expression within the tumor microenvironment (TME), which was associated with increased CD8+ T cell infiltration and decreased infiltration of CD4+ T regulatory cells. The mean neutrophil-to-lymphocyte ratio (NLR) decreased by greater than 50%, from 14.2 (SD 22.6) to 6.9 (SD 5.2). Median overall survival (OS) was 16.3 months (95% CI 1.9, NA) and 2-year OS rate was 24.7% (95% CI: 4.5%, 53.2%) with 58% (7/12) of treated participants demonstrated prolonged OS of greater than 12 months. Conclusion Our findings suggest that low-dose 5’aza can reprogram systemic host immune responses and the local TME to increase IFN-ɣ and PD-L1 expression. The increased expression of these established biomarkers correlated with prolonged OS upon re-challenge with ICB.

Study EGAS00001007998

Phosphoproteomics adds value to treatment recommendations in molecular tumor boards

Precision oncology approaches employing genomics-guided targeted therapies for individual patients have provided significant survival benefits in several cancer types. However, low response rates in most solid malignancies, many patients without actionable genomic lesions, and increasing evidence that non-genomic mechanisms may play an important role in tumors indicate that genomics alone is often insufficient to inform and guide the clinical care of patients. Here, we show for the first time that comprehensive (phospho)proteome profiling is feasible and informative in a real-world prospective precision oncology setting. We developed a novel tumor pathway activity (TUPAC) scoring methodology that provides a holistic perspective on dysregulated receptor tyrosine kinase (RTK) signaling in individual patients. Based on 919 tumor tissue profiles of patients enrolled in the national NCT/DKTK MASTER study and the INFORM registry trial, we illustrate how TUPAC scoring can be used to uncover individual tumor biology in rare cancers and demonstrate that TUPAC methodology identifies EGFR-driven sarcomas despite the absence of genetic EGFR alterations. Overall, our work demonstrates the utility of the additional phosphoproteome data layer to enhance therapeutic recommendations in molecular tumor boards.

Study EGAS00001007891

The University of Hong Kong Intestinal Metaplasia Organoids Study

Background: Gastric intestinal metaplasia (IM) is a pre-cancerous stage spanning a morphological spectrum that is poorly represented by human cell line models. Objective: We aim to establish and characterize human IM cell models to better understand IM progression along the cancer spectrum. Design: A large human gastric IM organoid (IMO) cohort (n=28), their clonal derivatives and normal gastric organoids (n=42) for comparison were established. Comprehensive multi-omics profiling and functional characterization were performed. Results: Single-cell transcriptomes revealed IMO cells spanning a spectrum from hybrid gastric/intestinal to advanced intestinal differentiation. Their lineage trajectories connected different cycling and quiescent stem and progenitors, highlighting differences in gastric to IM transition and the potential origin of IM from STMN1 cycling isthmus stem cells. Hybrid IMOs showed impaired differentiation potential, high lineage plasticity beyond gastric or intestinal fates, and reactivation of a fetal gene program. Cell populations in gastric IM and cancer tissues were highly similar to those derived from IMOs and exhibited a fetal signature. Genomically, IMOs showed elevated mutation burden, frequent chromosome 20 gain, and epigenetic de-regulation of many intestinal and gastric genes. Functionally, IMOs were FGF10 independent and showed downregulated FGFR2. Several IMOs exhibited a cell-matrix adhesion independent (CMi) subpopulation that displayed chromosome 20 gain but lacked key cancer driver mutations, potentially representing the earliest neoplastic precursor of IM-induced gastric cancer. Conclusions: Overall, our IM organoid biobank captured the heterogeneous nature of IM, revealing mechanistic insights on IM pathogenesis and progression, offering an ideal platform for studying early gastric neoplastic transformation and chemoprevention.

Study EGAS00001007899

Organoid drug profiling identifies methotrexate as a therapy for SCCOHT, a rare pediatric cancer

Small cell carcinoma of the ovary, hypercalcemic type (SCCOHT) is a rare and lethal tumor in adolescent and young adult patients. Currently, there is no standard-of-care treatment for these patients. Reliable models that represent this disease and can be used for translational research are scarce. To model SCCOHTs, we have established eight patient-derived tumoroid lines from tumor lesions of three patients with SCCOHT. The tumoroids recapitulate genomic and transcriptomic characteristics of the corresponding patient tumors and capture intrapatient tumor heterogeneity. Organoid drug profiling using a library of 153 clinical compounds identified methotrexate as an effective and selective drug against SCCOHTs with a clinically relevant IC50 of 35 nmol/L. RNA sequencing demonstrated that methotrexate induced TP53 pathway activation and apoptosis. These data underscore that organoid technology can support the design of therapeutic strategies for rare cancers.

Study EGAS00001007911

Integrated genomic analysis for HCC

Liver cancer is a common malignancy and is the third cause of cancer-related death in the world and hepatocellular carcinoma (HCC) is the major type of liver cancer, accounting for ~90% of cases. Only 20 ~ 30% of patients can benefit from combined immunotherapy and targeted therapy due to the high inter-tumor heterogeneity of HCC. This work aims to identify homogeneous and robust molecular subtypes in HCC based on a large, homogenous and well-annotated cohort. We performed whole genome/exome sequencing (WGS/WES) and/or RNA sequencing (RNA-seq) on 529 HCC from 461 patients collected mainly in France. Based on the great consistence between genomic and transcriptomic data, we identified 9 robust HCC subtypes mainly based on driver mutations. And we further characterized HCC subtypes using genomic, transcriptomic and clinicopathological features. Among 9 subtypes, 5 subtypes belonged to chromosome instable tumors, while 3 subtypes belonged to chromosome stable tumors. In addition, our subtypes were associated with prognosis and we also found different distributions of subtypes across various features, such as etiology and gender. Our study provided a robust molecular classification based on a large HCC dataset from Western country, it improves our understanding of the mechanisms of carcinogenesis and will facilitate development of genome-based precision medicine in HCC.

Study EGAS00001007957

7201 results for "canc*"

in 17.45 milliseconds.