Search Results - EGA European Genome-Phenome Archive

Patient-Derived Breast Cancer Organoid Study

This study uses normal and tumor tissues from breast cancer patients to derive 3D cell culture models called organoids. Organoid technology is an innovative platform that can be used to study a patient's tumor in a dish ex-vivo. As part of this study we have comprehensively characterized these models using next-generation sequencing for: 1) pathogenic SNVs in a targeted gene panel 2) copy number alterations 3) transcriptome analysis using RNA sequencing and 4) single-cell RNA sequencing

Study phs002722

MutWP5__CRUK_Mutographs_of_CancerBreastCancer_Mastectomy__Exome_

Sequencing of LCM-derived microbiopsies from 30 women who underwent mastectomies due to a breast cancer diagnosis. Goal to assess the mutational burden, spectrum, and clonal dynamics within the tissue. Exome data will be used as a driver and clonality screen, highly clonal or driver-containing samples will subsequently be sent for whole-genome sequencing. Results from this portion of the study will be compared to women who had cosmetic breast reduction surgeries and those with germline BRCA 1/2 mutations.

Study EGAS00001003319

Investigating_the_genetics_of_immunity_against_Salmonella_in_humans

This study entails whole genome sequencing of an interleukin (IL)-12 b-1 receptor-deficient individual who presented with a chronic systemic Salmonella Enteritidis infection that did not resolve with standard IFNg and antibiotic treatment. Whole genome sequencing of the patient’s parents are also included.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 EGAS00001001664

Genomic Signatures Reveal DNA Damage Response Deficiency in Brain Metastases of Colorectal Cancer

Brain metastasis (BM) of colorectal cancer (CRC) is rare but lethal, and lacks effective therapies or a good understanding of its genomic landscapes. We conduct an analysis of whole-exome sequencing (WES) and whole-genome sequencing (WGS) data on 19 trios of patient-matched BMs, primary CRC tumors, and adjacent normal tissue to (1) document genomic signatures during the evolution of CRC BM; (2) identify clinically actionable targets for the BM treatment; and (3) identify potential drivers in CRC BM.

Study EGAS00001003659

Interplay of somatic alterations and immune infiltration modulates response to PD-1 blockade in advanced clear cell renal cell carcinoma -- CA209-010

Immune checkpoint inhibitors targeting the PD-1 pathway have transformed the management of many advanced malignancies, including clear cell renal cell carcinoma (ccRCC), but the drivers and resistors of PD-1 response remain incompletely elucidated. Here, we analyzed 592 tumors collected from advanced ccRCC patients enrolled in prospective clinical trials of treatment with PD-1 blockade (or mTOR inhibition as control arm) by whole-exome and RNA-sequencing, integrated with immunofluorescence analysis, to define the somatic alteration landscape of late-stage ccRCC and to uncover the immunogenomic determinants of therapeutic response. While conventional genomic markers (tumor mutation burden, neoantigen load) and degree of CD8+ T cell infiltration were not associated with clinical response, we discovered numerous chromosomal alterations in advanced ccRCC associated with response or resistance to PD-1 blockade. These advanced tumors were highly CD8+ T cell infiltrated, with only 22% and 5% with an immune desert and immune excluded phenotype, respectively. Our analysis revealed that CD8+ infiltrated tumors are depleted of favorable PBRM1 mutations and are enriched for unfavorable chromosomal losses of 9p21.3 when compared to non-infiltrated tumors. These data demonstrate how the interplay of somatic alterations and immunophenotypes impacts therapeutic efficacy.

Dataset EGAD00001006026

Interplay of somatic alterations and immune infiltration modulates response to PD-1 blockade in advanced clear cell renal cell carcinoma -- CA209-009

Immune checkpoint inhibitors targeting the PD-1 pathway have transformed the management of many advanced malignancies, including clear cell renal cell carcinoma (ccRCC), but the drivers and resistors of PD-1 response remain incompletely elucidated. Here, we analyzed 592 tumors collected from advanced ccRCC patients enrolled in prospective clinical trials of treatment with PD-1 blockade (or mTOR inhibition as control arm) by whole-exome and RNA-sequencing, integrated with immunofluorescence analysis, to define the somatic alteration landscape of late-stage ccRCC and to uncover the immunogenomic determinants of therapeutic response. While conventional genomic markers (tumor mutation burden, neoantigen load) and degree of CD8+ T cell infiltration were not associated with clinical response, we discovered numerous chromosomal alterations in advanced ccRCC associated with response or resistance to PD-1 blockade. These advanced tumors were highly CD8+ T cell infiltrated, with only 22% and 5% with an immune desert and immune excluded phenotype, respectively. Our analysis revealed that CD8+ infiltrated tumors are depleted of favorable PBRM1 mutations and are enriched for unfavorable chromosomal losses of 9p21.3 when compared to non-infiltrated tumors. These data demonstrate how the interplay of somatic alterations and immunophenotypes impacts therapeutic efficacy.

Dataset EGAD00001006027

Interplay of somatic alterations and immune infiltration modulates response to PD-1 blockade in advanced clear cell renal cell carcinoma -- CA209-025

Immune checkpoint inhibitors targeting the PD-1 pathway have transformed the management of many advanced malignancies, including clear cell renal cell carcinoma (ccRCC), but the drivers and resistors of PD-1 response remain incompletely elucidated. Here, we analyzed 592 tumors collected from advanced ccRCC patients enrolled in prospective clinical trials of treatment with PD-1 blockade (or mTOR inhibition as control arm) by whole-exome and RNA-sequencing, integrated with immunofluorescence analysis, to define the somatic alteration landscape of late-stage ccRCC and to uncover the immunogenomic determinants of therapeutic response. While conventional genomic markers (tumor mutation burden, neoantigen load) and degree of CD8+ T cell infiltration were not associated with clinical response, we discovered numerous chromosomal alterations in advanced ccRCC associated with response or resistance to PD-1 blockade. These advanced tumors were highly CD8+ T cell infiltrated, with only 22% and 5% with an immune desert and immune excluded phenotype, respectively. Our analysis revealed that CD8+ infiltrated tumors are depleted of favorable PBRM1 mutations and are enriched for unfavorable chromosomal losses of 9p21.3 when compared to non-infiltrated tumors. These data demonstrate how the interplay of somatic alterations and immunophenotypes impacts therapeutic efficacy.

Dataset EGAD00001006029

TRIGR collaboration between University of Helsinkin and Children's Mercy Research Institute

DAC for collaborations between by the University of Helsinki and Children's Mercy Hospital, such as the Trial to Reduce IDDM in the Genetically at Risk (TRIGR), used to study the evolving epigenomics of immune cells in type 1 diabetes at single nuclei resolution

Dac EGAC50000000543

RRBS - MBD4-deficient AML

Reduced Representation Bisulfite Sequencing for WEHI-AML-1 and WEHI-AML-2. RRBS libraries were made with the NuGEN Ovation RRBS Methyl-Seq System. Bisulfite conversion was performed with the Qiagen Epitect kit. Sequencing was performed on an Illumina HiSeq2500.

Dataset EGAD00001003567

H3Africa TrypanoGEN Main WGS

This dataset is for TrypanoGEN Phase 1: Variant discovery, and includes 233 samples sequenced to approximately 10X coverage. Samples are from Guinea, Cote D’Ivoire, DRC and Uganda using Illumina HiSeq 2500.

Dataset EGAD00001005076

Martin_Ravenscroft-AUS1

Whole genome sequencing data on D19-0702 (AUS1), presented in Martin et al. 2020 (AUS1). WGS (Illumina HiSeq) was performed at Kinghorn Centre for Clinical Genetics, Garvan Institute of Medical Research. Data was analyzed using the Seave bioinformatic analysis pipeline (https://www.seave.bio).

Dataset EGAD00001006276

Fecal whole metagenomic shotgun sequencing data.

Human fecal WMS data from patients treated with combined anti-CTLA-4 and anti-PD-1 immunotherapy for advanced melanoma.

Dataset EGAD00001006734

HV31 - MGI standard short-read sequencing

MGI standard short-read sequencing data for individual HV31 generated using DNA from peripheral blood mononuclear cells, to a sequencing depth of ~57×.

Dataset EGAD00001007044

HV31 - MGI CoolMPS short-read sequencing

MGI CoolMPS short-read sequencing data for individual HV31 generated using DNA from peripheral blood mononuclear cells, to a sequencing depth of ~57×.

Dataset EGAD00001007048

Treatment arm information and sample sheet

Table of treatment arm information for the 418 RNAseq evaluable population.

Dataset EGAD00001009502

Genomic and transcriptomic determinants of host susceptibility, protection, and viral mutation in experimental SARS CoV-2 infection: RNA (2025-08-11)

Single-cell profiling of healthy adult volunteers that were inoculated with SARS-CoV-2. . This dataset contains all the data available for this study on 2025-08-11.

Dataset EGAD00001015679

A study of the immune system in patients with peripheral inflammatory neuropathy (CIDP): RNA adult (2025-10-02)

Leveraging single-cell sequencing technologies to shed light on the immune aetiology of chronic inflammatory demyelinating polyradiculoneuropathy (CIDP) . This dataset contains all the data available for this study on 2025-10-02.

Dataset EGAD00001015725

Transcriptomic profiling of the Enteric Nervous System in Hirschsprung Disease (2025-07-31)

This study involves atlasing the development of the postnatal gut nervous system in order to elucidate the pathogenic mechanisms of Hirschsprung disease. . This dataset contains all the data available for this study on 2025-07-31.

Dataset EGAD00001015667

MASS_Pilot: Muscle and Ageing Science Study - Collaboration with Wellcome Sanger Institute to characterise skeletal muscle ageing using Human Cell Atlas approaches (2025-07-31)

Transcriptomic analysis of skeletal muscle, to understand the mechanisms of muscle ageing. . This dataset contains all the data available for this study on 2025-07-31.

Dataset EGAD00001015670

Additional RNAseq data of metastatic breast cancer samples

Paired-end RNA sequencing data from 33 metastatic breast cancer samples. Sequencing was performed on the Illumina NovSeq 6000 using NEBNext Single Cell and Low Input RNA Library Prep Kit for Illumina.

Dataset EGAD00001015678

Ischemia Reperfusion Responses in Human Lung Transplants at the Single Cell Resolution

Ischemia reperfusion is an unavoidable step of organ transplantation. Development of therapeutics for lung injury during transplantation has proved challenging; understanding lung injury from human data at the single cell resolution is required to accelerate the development of therapeutics. Donor lung biopsies from six human lung transplant cases were collected at the end of cold preservation and 2-hour reperfusion and underwent single cell RNA sequencing. Donor and recipient origin of cells from the reperfusion timepoint were deconvolved. Gene expression profiles were (1) compared between each donor cell type between timepoints and (2) compared between donor and recipient cells. Inflammatory responses from donor lung macrophages were found after reperfusion with upregulation of multiple cytokines and chemokines, especially IL-1 and IL-1. Significant inflammatory responses were found in alveolar epithelial cells (featured by CXCL8) and lung endothelial cells (featured by IL-6 upregulation). Different inflammatory responses were noted between donor and recipient monocytes and CD8+ T cells. The inflammatory signals and differences between donor and recipient cells observed provide insight into the cellular and molecular mechanisms of ischemia reperfusion induced lung injury. Further investigations may lead to the development of novel targeted therapeutics.

Study EGAS50000000490

Single-cell RNA sequencing of pediatric Hodgkin Lymphoma to study the inhibition of T cell subtypes

Here, we applied single-cell RNA-sequencing (scRNA-seq) on isolated HRS cells and the immune cells from the same cHL tumors. This allowed us to identify genes of cell surface proteins that are consistently overexpressed in HRS cells and can potentially be used as targets for antibody-drug conjugates or CAR T cells. Finally, we identify potential interactions by which HRS cells inhibit T cells, among which the Galectin-1/CD69 and HLA-DRA/LAG3 interactions. However, high levels of inter-patient heterogeneity of the interaction strength were observed. In conclusion, this study identifies new potential therapeutic targets for cHL and highlights the importance of studying heterogeneity when identifying therapy targets.

Study EGAS50000000432

Lebanon_HighCov_seq

We propose to investigate the genomics of a Lebanese sample. we plan (1) genotyping (2.5M) of 120 Lebanese individuals, choosing this platform to allow comparison with the 1000 Genomes and African Genome Variation project populations, (2) 8x sequencing of 100 individuals to provide unbiased population-level sequence information including allele frequencies, and (3) 30x sequencing of four individuals to allow inferences such as from PSMC/MSMC to be obtained. The current form is for the high-coverage sequencing. 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 EGAS00001002085

Lebanon_LowCov_seq

we propose to investigate the genomics of a Lebanese sample. we plan (1) genotyping (2.5M) of 120 Lebanese individuals, choosing this platform to allow comparison with the 1000 Genomes and African Genome Variation project populations, (2) 8x sequencing of 100 individuals to provide unbiased population-level sequence information including allele frequencies, and (3) 30x sequencing of four individuals to allow inferences such as from PSMC/MSMC to be obtained. The current form is for the low-coverage sequencing. 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 EGAS00001002084

ChIP-seq for 10 samples

This dataset contains ChIP-seq data generated from two patient-derived xenograft (PDX) models. For each model, sequencing data are provided from input control (n=1), IgG control (n=2), and FOXF1 immunoprecipitation (n=2). All the 10 samples are in bam file type.

Dataset EGAD50000001786

Sequencing data for oesophageal and related samples - Alex Frankell et al (RNA)

Data supporting: "The landscape of selection in 551 Esophageal Adenocarcinomas defines genomic biomarkers for the clinic." Frankell et al. RNAseq (BAM files) 116 tumours

Dataset EGAD00001004423

Mapping genetic variants underlying gene regulation in healthy intestinal cell types to identify novel IBD drug targets

Biopsies from the terminal ileum and rectum of healthy individuals are digested on ice to single cells and processed for single-cell RNA-sequencing (10X Genomics and Illumina) 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/.

Dataset EGAD00001006331

DETECT-A NGS Data Batch 1

This batch is a subset of the full DETECT-A dataset, containing 43963 fastq files generated from 3024 subjects. All sequencing was conducted using Illumina HiSeq 4000 and Illumina MiSeq platforms. Note that the division into batches follows no specific criteria and that the sequencing data for each subject has multiple files which may span multiple datasets. Thus, for a comprehensive analysis, it is recommended to request access to all datasets that comprise this study.

Dataset EGAD00001008415

iPSC and iNeuron RNAseq, chip-seq and single cell CRISPR activation

Dataset includes 1) RNAseq for Bob_Ngn2 cell differentiation from iPSC to iNeuron stage (time point 0, 24, 48 and 96 hours). 2) Bob_Ngn2 cell line chip-seq (H3K4me3, H3K4me1, H3K27me3, H3K9me3, H3K27ac, H3K36me3) for both iPSC and iNeuron stage with 3 replicates at each time point. 3) Single cell CRISPR activation experiment with 96 endogenous genes.

Dataset EGAD00001010050

Sequencing data for oesophageal / related samples - Kazachenka et al (DNA)

Data supporting: "The transcriptional landscape of endogenous retroelements delineates esophageal adenocarcinoma subtypes" WGS for 452 samples

Dataset EGAD00001011095

DETECT-A NGS Data Batch 6

This batch is a subset of the full DETECT-A dataset, containing 9 fastq files generated from 8 subjects. All sequencing was conducted using Illumina HiSeq 4000 and Illumina MiSeq platforms. Note that the division into batches follows no specific criteria and that the sequencing data for each subject has multiple files which may span multiple datasets. Thus, for a comprehensive analysis, it is recommended to request access to all datasets that comprise this study.

Dataset EGAD00001015458

DETECT-A NGS Data Batch 4

This batch is a subset of the full DETECT-A dataset, containing 25210 fastq files generated from 3040 subjects. All sequencing was conducted using Illumina HiSeq 4000 and Illumina MiSeq platforms. Note that the division into batches follows no specific criteria and that the sequencing data for each subject has multiple files which may span multiple datasets. Thus, for a comprehensive analysis, it is recommended to request access to all datasets that comprise this study.

Dataset EGAD00001015360

DETECT-A NGS Data Batch 5

This batch is a subset of the full DETECT-A dataset, containing 12740 fastq files generated from 3211 subjects. All sequencing was conducted using Illumina HiSeq 4000 and Illumina MiSeq platforms. Note that the division into batches follows no specific criteria and that the sequencing data for each subject has multiple files which may span multiple datasets. Thus, for a comprehensive analysis, it is recommended to request access to all datasets that comprise this study.

Dataset EGAD00001015426

WXS dataset of Oncogenic and immunological targets for matched therapy of pediatric blood cancer patients: Dutch iTHER study experience

Dataset EGAD00001015638

RNA-seq dataset of Oncogenic and immunological targets for matched therapy of pediatric blood cancer patients: Dutch iTHER study experience

Dataset EGAD00001015639

Resuscitation Outcomes Consortium (ROC) Cardiac Epidemiologic Registry (Cardiac Epistry) Versions 1 and 2 (ROC-Cardiac Epistry 1 and 2-BioLINCC)

Data Access NOTE: Please refer to the “Authorized Access” section below for information about how access to the data from this accession differs from many other dbGaP accessions. Available Data: The ROC Cardiac Epistry versions 1 and 2 include all cardiac arrest cases entered into the ROC database from December 1, 2005 to April 1, 2011. ROC Cardiac Epistry 3 includes cardiac arrest cases captured between 2011 and 2015, and introduced significant changes in how and which data were collected.For ROC traumatic injury Epistry data, please see: ROC-Trauma Epistry Objectives: To build a prospective population-based registry of participants with out-of-hospital cardiac arrest responded to by Emergency Medical Services (EMS). Specific aims: to establish whether the results of Resuscitation Outcomes Consortium (ROC) trials can be generalized to the larger population of participants that experience cardiac arrest;to more fully establish the burden of cardiac arrest; and to examine the relationships between variation in EMS structure and process, regional and periodic factors, and participant outcomes.Background: Cardiac arrest is a common, serious, debilitating and costly public health problem. Although there has been a steady decline in morbidity and mortality from most cardiovascular diseases, high mortality rates for out-of-hospital cardiac arrests continue to pose a challenge for healthcare providers and a significant public health burden. The Resuscitation Outcomes Consortium (ROC) was established in 2004 to conduct clinical research in the areas of cardiopulmonary arrest and life-threatening traumatic injury with the overall goal of improving resuscitation outcomes. Participant and care characteristics can predict favorable outcomes in cardiac arrests, but there is still a wide variation in outcomes that is not well understood. EMS factors such as service level, number of responding providers, use of procedures or drugs in the field, training, quality assurance/feedback, and response time intervals also vary significantly by region. Variations in geographic, socioeconomic and periodic factors may also be associated with differences in outcomes.Prior to ROC Cardiac Epistry, there were no North American population-based registries for out-of-hospital cardiac arrests. Therefore there was a need for standardized data collection of out-of-hospital cardiac arrests in diverse geographic locations in order to identify the independent effects of prognostic or treatment factors accounting for variations in survival. Participants: The registry included 109,326 cardiac arrest events from 264 EMS agencies transporting to 287 acute care hospitals from the following regional centers: Birmingham, Alabama; Dallas, Texas; Iowa City, Iowa; Milwaukee, Wisconsin; Pittsburgh, Pennsylvania; Portland, Oregon; San Diego, California; Seattle/King County, Washington; Ottawa, Ontario; Toronto, Ontario; and Vancouver, British Columbia.Design: ROC Epistry collected standardized data regarding episode-specific factors, participant demographics, clinical information, pre-hospital interventions and disposition, hospital information, and participant outcome for all out-of-hospital cardiac arrests in the ROC regions. Each ROC site had to ensure capture of all eligible cases within the EMS service areas. Out-of-hospital data were extracted from existing databases whenever possible and augmented with targeted review of EMS reports. Hospital data were abstracted directly from the hospital file in most cases. Alternative methods included linkage to death registries and obituaries if the death occurred within 30 days. Sites submitted data using a web-based interface or batch uploads. Participants were not contacted directly.

Study phs003803

FEGA and European GDI: working together to improve human health

The Federated EGA is a global resource for discovery of and access to sensitive human omics and associated data consented for secondary use, through a network of human data repositories to accelerate biomedical research and improve human health. The Federated EGA network was launched in September 2022 with five inaugural nodes, and since 2023 seven operational nodes can share data across national borders in adherence to European and national laws. A few weeks after FEGA's official launch, in November 2022, the European Genomic Data Infrastructure (GDI) project was kicked-off . This European Commission co-funded project, coordinated by ELIXIR, is aimed to deliver federated, sustainable and secure data infrastructure to access genomic and related phenotypic and clinical data across Europe. This project supports the aim of the 1+MG initiative (25 EU countries, Norway and UK) to enable personalised medicine and health through a shared framework and infrastructure for securely accessing and integrating high quality genomic data and other health data across borders.1+MG will be an integral component of the European Health Data Space (EHDS) for secondary use (Healthdata@EU) as an authorised participant. How are Federated EGA and GDI similar? Given their shared visions, the Federated EGA and European GDI networks have a lot in common. To begin, they share the same overall goal of establishing networks of “nodes” that host sensitive human data within a jurisdiction and connecting these nodes in a global network to support data discovery and promote human genome data access for research and healthcare. FEGA and GDI are both initially focused on enabling access to human genomic data, while FEGA ambition is to expand the scope to clinical research data and other omics data in the future. FEGA and GDI are both built on open and interoperable software solutions, a subset of which are based on the LocalEGA components. FEGA and GDI implementation solutions are based on international community standards, for example those developed by the Global Alliance for Genomics and Health, which contributes to making them interoperable. Both FEGA and GDI allow for the nodes to make use of any solution that is fully compatible. As illustrated in Figure 1, a final key point is the significant overlap of institutions involved in envisioning and operating FEGA and GDI nodes, with the clear wish to keep them interoperable in the future. What distinguishes Federated EGA from GDI? Despite being largely similar, there are some differences between the FEGA and GDI networks, which we aim to clarify in this post. The first difference is the governance model that the nodes will operate to comply with national laws and GDPR. In the FEGA network, nodes have taken inspiration from the EGA data access model, where the infrastructure is a data processor of the hosted data and data controllership remains with the originating Data Access Committees (i.e. Data controllers) for each dataset. On the other hand, in the GDI network the controllership of datasets will be transferred to a 1+MG European Digital Infrastructure Consortium (EDIC) legal entity created by the Member States (MS), who will make data access decisions, with data holders having veto powers for their datasets. Importantly, FEGA nodes have the flexibility to choose another model to fit with their data protection framework, including becoming data controllers for their datasets. The second difference is the inclusion criteria for data. While FEGA nodes are designed to accept almost any type of omics data in need of control access (e.g. genomics, transcriptomics, genotyping, single cells sequencing, patient-tracked metagenomics), GDI nodes are initially, but not only, focused on accepting whole genome and exome sequencing data and affiliated data from sources such as (i) the Genome of Europe use case of the 1+MG initiative which specifically aims to fulfil the mission of building “a European network of national genomic reference cohorts of at least 500,000 citizens”) (ii) data collection of other types of genomic data identified by countries through the 1+MG dashboard and (iii) genomic data coming from data holders which would need to fulfil EHDS requirements. The third difference is the maturity of the software stacks. The FEGA network provides a set of software for data and metadata submission, storage, permissions management, and file distribution (the LocalEGA package). GDI is building a complete set of open source reference for the five functionalities covering the full data life cycle (a few more compared to the LocalEGA) including federated processing (analytics, AI/ML), which is still under active development. Notably, one of the GDI functionalities - storage and interfaces - can be satisfied by using the LocalEGA storage solution, highlighting the ability of FEGA and GDI to be interoperable. In both FEGA and GDI, nodes are allowed to use alternative solutions as long as they are fully interoperable with the networks. Figure 1 provides a simple overview of these described commonalities and differences. Figure1: schematic overview of the main commonalities and difference among a FEGA and a GDI node Can the same institute run a FEGA and a GDI node at the same time? The answer is Yes! We believe this is entirely possible and encouraged. As long as 1+MG requirements are continuously fulfilled in the implementation. The same trained personnel could operate the infrastructure and leverage the same national funding to run a FEGA and a GDI node. Several European nodes, especially GDI vanguard nodes like Norway, Sweden, Finland and Spain are following this model. The idea is that a lot of the work can be reused, given the overlapping scope and the interoperable technology. Thus, the hosted datasets can be discovered via the two catalogues, perhaps accessible under the same or different governance models. From a FEGA node perspective, the GDI datasets could simply be a subset of the datasets hosted in the FEGA node, which have their specific governance as any other. What's next? So, the amazing teams of people building the Federated EGA network and the European GDI network are working together to create infrastructures able to provide secure access to human genomic and associated data around Europe. Because we all know how much human genomics can improve healthcare and precision medicine. And we all want to collaborate to make it happen. This article has been reviewed collectively by members of the EGA and the GDI coordination team.

Blog fega-and-gdi

Resident memory CD8 T cells persist for years in human small intestine

Resident memory CD8 T cells (Trm) have been shown to provide effective protective responses in the small intestine (SI) in mice. A better understanding of the generation and persistence of SI CD8 Trm cells in humans may have implications for intestinal immune-mediated diseases and vaccine development. Analyzing normal and transplanted human SI we demonstrated that the majority of SI CD8 T cells were bona fide CD8 Trm cells that survived for over 1 year in the graft. Intraepithelial and lamina propria CD8 Trm cells showed a high clonal overlap and a repertoire dominated by expanded clones, conserved both spatially in the intestine and over time. Functionally, lamina propria CD8 Trm cells were potent cytokine-producers, exhibiting a polyfunctional (IFN-γ+ IL-2+ TNF-α+) profile, and efficiently expressed cytotoxic mediators after stimulation. These results suggest that SI CD8 Trm cells could be relevant targets for future oral vaccines and therapeutic strategies for gut disorders.

Study EGAS00001003676

Comprehensive Genomic Characterization of Translocation Renal Cell Carcinoma

Xp11 translocation renal cell carcinoma (tRCC) is a rare, female-predominant cancer driven by a fusion between the TFE3 gene on chromosome Xp11.2 and a partner gene on either chrX or an autosome. It remains unknown what types of rearrangements underlie TFE3 fusions, whether fusions can arise from both the active (chrXa) and inactive X (chrXi) chromosomes, and whether TFE3 fusions from chrXi translocations account for the female predominance of tRCC. To address these questions, we performed haplotype-specific analyses of chrX rearrangements in tRCC whole genomes. We show that TFE3 fusions universally arise as reciprocal translocations and that oncogenic TFE3 fusions can arise from chrXi:autosomal translocations. Female-specific chrXi:autosomal translocations result in a 2:1 female-to-male ratio of TFE3 fusions involving autosomal partner genes and account for the female predominance of tRCC. Our results highlight how X chromosome genetics uniquely constrains somatic chrX alterations and underlies cancer sex differences.This version of the dbGaP study provides a comprehensive set of whole genome data and tables with clinical information collected for the 16 patients (15 pateints as tRCC and 1 patient as ccRCC) in the study cohort. The data also contains 12 RNA-seq samples used in the study and analyses. We anticipate that this data will considerably advance our molecular and genetic understanding of oncogenic translocations in TFE3-tRCC involving inactive chromosome X and its role in cancer sex bias.

Study phs003008

Field studies of Cryptosporidiosis and Enteropathogens in Bangladesh

Cryptosporidiosis causes severe diarrhea in infants in the developing world. There is no vaccine to prevent it, and little in the way of treatment. This study on Bangladeshi urban slum children aims to support the design of a vaccine, both by determining how the immune system protects from infection and by identifying the genotypes of the parasite that should be included in a vaccine, as well as aid in development of therapies by identifying human genes that control infection. The primary objective of the study was to determine the incidence and contribution to disease of the different species and genotypes of cryptosporidia. Secondary objectives were designed to determine acquired immune response to cryptosporidiosis and identify human genes that influence susceptibility to cryptosporidiosis. This is an observational study. Children were recruited from Mirpur Dhaka slum (Cohort 1) and rural Mirzapur (Cohort 2) and followed for cryptosporidium infection longitudinally from birth through age 4 years in Cohort 1 and birth through age 2 years in Cohort 2. Biweekly household visits for diarrheal surveillance were made in addition to anthropometric measurements of mother and child, blood samples collected two times each year from the child, a work-up of diarrheal stools and non-diarrheal surveillance stools for cryptosporidium and other enteropathogens, and blood and breast milk samples from the mother. A detailed description of the study design and procedures can be obtained from publication: Kevin L Steiner, et al., 2018, PMID: 29897482.

Study phs001665

Human genomic and phenotypic synthetic data for the study of rare diseases

The purpose of this dataset is to facilitate development of technical implementations for rare disease data integration, analysis, discovery, and federated access. This synthetic dataset includes clinical and genomic data from 6 rare disease cases. It consists of 18 whole genomes (6 index cases with their parents) which have genetic background based on public human data sequenced in the context of the Illumina Platinum initiative (Eberle, MA et al. (2017)) and made available by the HapMap project (https://www.genome.gov/10001688/international-hapmap-project). In each of the cases, real causative variants correlating with the phenotypic data provided were spiked-in. The cases included in this synthetic dataset correspond to the following type of disorders: CASE 1- Congenital myasthenic syndrome (Autosomal Dominant -de novo variant) CASE 2- Macular dystrophy (Autosomal Dominant) CASE 3- Muscular dystrophy (Autosomal Recessive-compound heterozygous variants) CASE 4- Mitochondrial disorder (Autosomal Recessive-consanguineous case - homozygous variant) CASE 5- Breast cancer (Autosomal Dominant) CASE 6- Similar as case 1 for patient matchmaking tests: Congenital myasthenic syndrome (Autosomal Dominant-de novo variant) For each case you will be able to download the following data: clinical information (phenopackets per individual and pedigree per family), raw genomic data (FASTQ and BAMs) and processed genomic data (vcfs). When using the data, the following should be acknowledged: the RD-Connect GPAP (https://platform.rd-connect.eu/), EC H2020 project EJP-RD (grant # 825575), EC H2020 project B1MG (grant # 951724) and Generalitat de Catalunya VEIS project (grant # 001-P-001647).

Study EGAS00001005702

BCR Signaling in human BM PC

The functions of bone marrow plasma cells (BMPC) beyond antibody production are not fully elucidated and distinct subsets of BMPC suggest potential different functions. Phenotypic differences were identified for human BMPC depending on CD19 expression. Since CD19 is a co-stimulatory molecule of the B-cell-receptor (BCR), and IgA+ and IgM+ BMPC express the BCR on their surface, we here studied whether CD19 expression affects cellular responses, such as BCR signaling and the expression of checkpoint molecules. We analyzed 132 BM samples from individuals undergoing routine total hip arthroplasty. We found that both CD19+ and CD19- BMPC expressed BCR signaling molecules. Notably, the BCR-associated kinase spleen tyrosine kinase (SYK) including pSYK was higher expressed in CD19+BMPC compared to CD19- BMPC. BCR stimulation also resulted in increased kinase phosphorylation downstream of the BCR while expression of CD19 remained stable afterwards. Interestingly, the BCR response was restricted to IgA+ BMPC independently of CD19 expression. With regard to the expression of checkpoint molecules, CD19- BMPC expressed higher levels of co-inhibitory molecule programmed cell death protein-1 (PD-1) than CD19+ BMPC. IgA+ BMPC characteristically upregulated PD-1 upon BCR stimulation in contrast to other PC subsets and inhibition of the kinase SYK abrogated PD-1 upregulation. In contrast, expression of PD-1 ligand, B and T lymphocyte attenuator (BTLA) and CD28 did not change upon BCR activation of IgA+ BMPC. Here, we identify a distinct functiotype of IgA+ BMPC that is independent of the phenotypic heterogeneity of the subsets according to their CD19 expression. The data suggest that IgA+ BMPC underlie different regulatory principles and/or exert distinct regulatory functions.

Study EGAS00001004948

Human breast transcriptome analysis

Obesity is a risk factor for the development of post-menopausal breast cancer and is associated with a worse breast cancer prognosis. Here we carried out transcriptome profiling on nontumorous breast tissue to define the molecular changes that occur related to excess body fat and adipose inflammation. A second goal was to determine if blood biomarkers reflect molecular changes in the breast. Increases in a variety of pathways linked to tumor development and progression were found. Moreover, increased expression of several drug targets including aromatase, the rate-limiting enzyme for estrogen biosynthesis, and PD-1 was observed. Various blood biomarkers that are associated with either elevated risk of breast cancer or worse prognosis were altered and correlated with molecular changes in the breast. Collectively, this study helps to explain both the link between obesity and breast cancer and the utility of blood biomarkers for determining risk and prognosis.

Study EGAS00001004665

Common origin and somatic mutation patterns of composite lymphomas and leukemias

The dataset contains six composite lymphomas and leukemias cases (four composite B-cell lymphomas (case 1-4), and two combinations of B-cell and T-cell lymphoma (case 5+6)). Two tumor and one non-tumor (NTC) sample were sequenced from each case after library prep with NEBNext Ultra DNA Library Prep Kit, Exome capturing was done with NimbleGen SeqCap EZ Choice kit, and sequenced on a Illumina HiSeq 2000 HO and given as fastq files. WGA (Qiagen REPLI-g kit) was used for case 2 and 3. Case 1-4 were microdisected, case 5+6 were flow-sorted.

Dataset EGAD50000001486

RNASeq from PC12 cell lines (derived from PPGL) cultivated under 21% (normoxia) and 1% (hypoxia) oxygen conditions.

RNASeq fastq files from PC12 cells (derived from rat PPGL) were cultivated in normoxia conditions (37˚C, 5% CO2 and 21% O2 balanced with N2) or in hypoxia conditions in a hypoxic incubator (37˚C, 5% CO2 and 1% O2 balanced with N2), for short time (12h, 24h, and 48h) or prolonged time (36 days). Samples at each timepoint were cultured in triplicate (total 24 paired end fastq files). Sequencing was performed with the Illumina NovaSeq6000.

Dataset EGAD50000001200

Multi-region RNA-Seq data of 10 neuroblastoma cases

To investigate intratumour heterogeneity and to better understand tumour evolution in neuroblastoma, we have performed a multi-region RNA sequencing on a total of 51 spatially separated tumor samples from 9 primary neuroblastomas (2 low-risk, 1 intermediate-risk and 6 high-risk) and 1 relapsed neuroblastoma. We also assessed the impact of chemotherapy on the clonal expansion by sequencing tumour regions from one medium risk and one high-risk tumour for which we had matched samples obtained at diagnosis and after chemotherapy.

Dataset EGAD00001008133

eMERGE Network Phase III Clinical Sequencing: eMERGEseq Panel

The Electronic Medical Records and Genomics (eMERGE) Network is a National Institutes of Health (NIH)-organized and funded consortium of U.S. medical research institutions. The primary goal of the eMERGE Network is to develop, disseminate, and apply approaches to research that combine biorepositories with electronic medical record (EMR) systems for genomic discovery and genomic medicine implementation research. eMERGE was announced in September 2007 and began its third phase in September 2015. eMERGE III consists of nine study sites, two central sequencing and genotyping facilities, and a coordinating center. eMERGE Phase III aims to: 1) sequence and assess the phenotypic implication of rare variants in a custom designed eMERGEseq panel consisting of 109 genes (including 56 ACMG actionable finding list genes and the top 6 genes from each site relevant to their specific aims), as well as approximately 1400 SNPs; 2) assess the phenotypic implications of these variants by developing, validating and implementing new phenotype algorithms, 3) integrate genetic variants into EMRs to inform clinical care; and 4) create community resources. Included in this study are: ~24,000 eMERGE participants from 10 eMERGE III study sites. Corresponding demographics, body mass index measurements. Top PheWAS codes generated from a collated list of ICD codes from all study sites. Study sites and participants include: Cincinnati Children's Hospital Medical Center (CCHMC): Cincinnati Children's Hospital Medical Center (CCHMC) is a not-for-profit hospital and research center pioneering breakthrough treatments, providing outstanding family-centered patient care and training healthcare professionals for the future, and dedicated to improving health and welfare of children and to the shared purpose of discovery and practical application of new genomic information to the ordinary care of children. We bring a comprehensive electronic health record (EPIC), a deidentified i2b2 data warehouse of 680K patient records, a biobank with >261,000 consents that allow return of results to >84,000 patients and guardians who have provided DNA samples, and hundreds of faculty and senior staff who make genomics or informatics an active focus of their research. CCHMC will help the eMERGE III Steering Committee identify genes for the eMERGE III targeted sequencing panel, provide 3,000 DNA samples from CCHMC patients to be sequenced, review targeted gene panels from clinical care at CCHMC for somatic mosaicism and reinterpretation, and further develop and disseminate a software workflow suite for sequence analysis. We will also extend our work generating phenotype algorithms using heuristic and machine learning methods to many new childhood diseases. We will develop tools to evaluate adolescent return of results preferences, examine the ethical and legal obligations and potential to reanalyze results, and develop clinical decision support for phenotyping, test ordering, and returning sequencing results. Children's Hospital of Philadelphia (CHOP): The Center for Applied Genomics (CAG) is a specialized Center of Emphasis at the Children's Hospital of Philadelphia (CHOP), and one of the world's largest genetics research programs, with to state-of-the-art high-throughput sequencing and genotyping technology. Our primary goal is to translate basic research findings to medical innovations. We aim to develop new and better ways to diagnose and treat children affected by rare and complex medical disorders, including asthma, autism, epilepsy, pediatric cancer, learning disabilities, and a range of rare diseases. Ultimately, our objective is to generate new diagnostic tests and to guide physicians to the most appropriate therapies. Participants were recruited from the CAG biorepository (n>450,000), specifically from >100,000 CHOP pediatric patients and family members, which is enriched for rare-diseases (n>12,000). Center for Applied Genomics, The Children's Hospital of Philadelphia We gratefully thank all the children and their families who enrolled in this study, and all individuals who donated blood samples for research purposes. Genotyping for this project was performed at the Center for Applied Genomics and supported by an Institutional Development Award from The Children's Hospital of Philadelphia. Sequencing was supported by the National Institutes of Health through an award from the National Human Genome Research Institute's Electronic Medical Records and Genomics (eMERGE) program (U01HG008684). Columbia University: The goal of the Columbia eMERGE III project is to develop methods for integrating genomic data in EHRs and to study the impact of such genomic informatics interventions on the health of a diverse, underserved urban adult English- and Spanish-speaking patient population in Northern Manhattan served by Columbia University Medical Center/New York-Presbyterian Hospital system. The study group is 2500 patients recruited from diverse clinics and community outreach centers of self-reported White (~61%), Asian (~11%), African-American (~11%), American Indian/Alaska Native (<1%) racial and Hispanic (~33%) ethnic backgrounds. There are two subgroups in the study cohort - a retrospective group (N=1052) that includes patients from oncology and nephrology clinics, and a prospective one (N=1448) that includes healthy individuals as well as participants with diverse medical conditions. Confirmed pathogenic variants in 70 selected genes will be returned to participants and their healthcare providers through the EHR integration. Participants are able to choose the results they receive and will have the freedom to meet with a genetic counselor and a geneticist to review results. The impact of genetic testing on clinical care is determined by periodic monitoring of EHRs. Geisinger: Samples and phenotype data in this study were provided by the Geisinger MyCode® Community Health Initiative. Participants are recruited across the Geisinger System via online consents or in-person consents at a hospital or clinic visit. Enrollment is ongoing with over 100,000 individuals currently consented. Partners Healthcare (Harvard University): The Partners HealthCare Biobank is a large research program designed to help researchers understand how people's health is affected by their genes, lifestyle, and environment. This large research data and sample repository provides access to high-quality, consented blood samples to help foster research, advance our understanding of the causes of common diseases, and advance the practice of medicine. For the Partners research community (Massachusetts General Hospital and Brigham and Women's Hospital), the Biobank provides: Banked samples (plasma, serum, and DNA) collected from consented patients Blood samples that were discarded after clinical testing in the Crimson Cores maintained in the Brigham and Women's Hospital and Massachusetts General Hospital Pathology Departments Sample handling and preparation services Link to the biobank data to the Partners Research Patient Data Registry (RPDR) a research instance of our electronic clinical chart Data access through our research portal. To date, over 70,000 Partners patients have given their consent to enroll, give a blood sample, receive research results and agreed to be re-contacted for additional research studies. The Biobank has enabled Partners investigators to compete for nationally recognized grants in personalized medicine such as a clinical electronic Medical Records and Genomics network (eMERGE) site and the national All of US program. The Biobank currently supports over 120 Partners investigators and over 130 million dollars in NIH research. Kaiser Permanente Washington/ (KPWA) / University of Washington (UW): KPWA participants were enrolled in the eMERGE Network through the Northwest Institute of Genetic Medicine (NWIGM) biorepository, and provided the appropriate consent to receive clinically relevant genetic results (N=2,500.) NWIGM is based at the University of Washington and co-managed by the University of Washington and KPWA. The purpose of the NWIGM biorepository is to build infrastructure and resources to carry out a broad range of future genetic research. KPWA members enrolled in the biorepository are asked to provide informed consent to providing a DNA sample for storage in the NWIGM biorepository. The consent is purposefully broad to serve the dual purpose of reducing the burden on researchers who wish to use this biorepository and the IRB committees who will be responsible for reviewing these requests in the future. Participants were eligible if aged 50 - 65 years old at the time of their enrollment into the NWIGM repository, living, enrolled in KPWA's integrated group practice, and had completed an online Health Risk Appraisal. The selection algorithm was based on several data sources from the EHR at KPWA. 1) Demographics - participants with self-reported race as Asian ancestry were prioritized and selected to enrich for non-European ancestry. The KPWA eMERGE cohort includes N=1,245 members of Asian ancestry. 2) Participants were also selected for a history of colorectal cancer (N=1,255), in order to allow us to enrich germline pathogenic variants. Mayo Clinic: The Return of Actionable Variants Empirical (RAVE) Study was approved by the Mayo Clinic IRB. We recruited 2537 participants from Mayo Clinic biobanks in Rochester, MN, who had hypercholesterolemia or colon polyps, thereby enriching for Familial hypercholesterolemia (FH) and monogenic causes of colorectal cancer (CRC). Additional eligibility criteria were: 1) residents of Southeast MN who were alive and aged 18-70 years; 2) LDL-C level >155 or >120 mg/dl while on lipid-lowering therapy; 3) no known cause of secondary hyperlipidemia; and 4) no cognitive impairment or dementia that would compromise their ability to give written informed consent. Based on these criteria, we identified 5270 eligible patients and obtained informed consent from 3030 participants. Recruitment was conducted in waves and utilized mailed recruitment packets consisting of a study brochure, a written informed consent form, a baseline psychosocial questionnaire, and a return postage-paid envelope. DNA of 2537 participants was sent for CLIA-certified targeted sequencing of 109 genes including genes associated with FH and CRC. Targeted sequencing and genotyping was performed in a Central Laboratory Improvement Amendment (CLIA)-certified laboratory. Northwestern University: Samples and data used in this study were obtained from patients from Northwestern Medicine, an integrated healthcare system, formed through a partnership of Northwestern Memorial HealthCare and Northwestern University Feinberg School of Medicine. Participants include a retrospective cohort from the Northwestern Pharmacogenomics Study, funded through the eMERGE II project, NHGRI (3U01HG006388-02S1) and a prospective cohort from the Genetic Testing and Your Health Study, funded through the eMERGE III project, NHGRI (U01HG008673). Patients were eligible to participate if they were18 years or older and see a physician at Northwestern Medicine. Patients consented to genetic testing and to allow their results to be placed in their electronic medical record. Vanderbilt University Medical Center: Vanderbilt University Medical Center (VUMC) participants were enrolled in the eMERGE Network through the Vanderbilt Genome-Electronic Records (VGER) project. Patients were provided the appropriate consent to receive clinically relevant genetic results (N=2,700). Participants were eligible if aged 21 or over, had a healthcare provider at VUMC, and visited the provider at least 3 times in the past 3 years. Meharry Medical College: Inclusion of ethnic groups in genomic research is critical to identify possible reasons for health disparities. African-Americans are being enrolled in various outpatient clinics of Nashville General Hospital at Meharry, an inner city hospital primary serving a poorer patient group. A total of 500 African Americans with four cancer types demonstrating health disparities in this population - prostate, colon, breast, lung are identified and approached by clinical research coordinators. The purpose of the study is to determine if any genetic information can be identified from these patients who have or are at high risk of one of these disparate cancers. All participants provide written informed consent and HIPAA authorization to provide blood samples for broad research use and permission to access data in their hospital electronic medical record for research now and in the future. An extensive demographic profile is obtained and entered into a REDCap database. Blood samples are obtained for a panel of alleles from extracted DNA at Baylor. In addition, de-identified coded samples are processed and stored in a central biorepository for further DNA, RNA and proteomic analyses. The survey and phlebotomy are performed at the time of the initial contact and agreement to participate. Nearly all patients approached willingly agree to participate for potential benefit to themselves, family members, or humankind. Little concern is voiced of providing samples for genetic analysis. Study investigators will share results with the participants and providers if testing does not indicate high risk. Results indicating increased risk or actionable alleles for the patient and/or family will be returned by a genetic counselor. Monitoring of the patients' health in this cohort will continue to be followed in the EMR to identify any future associations that might explain health disparities in African Americans. Proposals will be reviewed from investigators to study the genetic or proteomic samples as well as the clinical and demographic information in the repository. Please note that this version of the dataset has a handful of mismatches between genotyped and provided sex. Data with the following IDs should be removed prior to analysis: 420252874213744142412243424569384245694642672223

Study phs001616

RNA-Seq Brest Patient-derived Tumor Organoids

This dataset includes 14 RNASeq samples from 8 patients. For some patients, organoids at different passages have been profiled. They can be identified by their alias, composed of the following structure: <subjectId>_O<organoid_number>_<passage>_RNA. For example, ICSBCS002_O1_1_RNA indicates the earlier passage of organoid 1 for subject ICSBCS002, and ICSBCS002_O1_2_RNA indicates the later passage. Additionally, WCM2137_O1_RNA and WCM2137_O2_RNA indicate two organoids (O1 and O2) derived from different tissues. More information about the passages and tissues can be found in the sample information table.

Dataset EGAD50000000960

Molecular Characterization of Prostate Cancer Specimens by Bulk and Single Cell Analysis

Version 1: We performed single-cell RNA-sequencing (scRNA-seq) of fresh biopsies from metastatic castration resistant prostate cancer. Additionally, where available, we performed matched bulk whole exome sequencing (WXS; tumor and normal paired) and whole transcriptome sequencing for the same individuals. We also performed bulk whole exome (WXS) and RNA-sequencing (RNA-Seq) from FFPE pre-treatment biopsies from patients with clinically localized high-risk prostate cancer.Version 2: We added data from a randomized phase 2 trial of radium-223 with or without pembrolizumab in patients with Metastatic Castration-Resistant Prostate Cancer (mCRPC) to bone. After enrollment, study participants underwent bone biopsy and were randomized 2:1 to receive radium-223 at 55 kBq/kg every 4 weeks + pembrolizumab at 200 mg every 3 weeks (Arm A, R223+P) or radium-223 at 55 kBq/kg every 4 weeks alone (Arm B, R223). Stratification factors were alkaline phosphatase ≥220 vs. Data generated includes bulk whole exome sequencing (WXS), ultra-low pass whole genome sequencing (WGS) at 0.1X coverage, and targeted-capture sequencing. The ultra-low pass sequencing was used as part of the selection criteria for the target capture sequencing. Note: Two subjects in this study are linked to phs001577 and phs000915 respectively.

Study phs001988

Genomic evolution and transcriptional changes in the evolution of prostate cancer into neuroendocrine and ductal carcinoma types (RNAseq)

Prostate cancer is typically of acinar adenocarcinoma type but can occasionally present as neuro-endocrine and/or ductal type carcinoma. These are associated with clinically aggressive disease and the former often arises on a background of androgen deprivation therapy, although it can al-so arise de novo. Two prostate cancer cases were sequenced by exome capture from archival tis-sue. Case 1 was de novo small cell neuroendocrine carcinoma and ductal adenocarcinoma with 3 longitudinal samples over 5 years. Case 2 was a single time point after development of treatment related neuroendocrine prostate carcinoma. Case 1 showed whole genome doubling in all sam-ples and focal amplification of AR in all samples except the first time point. Phylogenetic analysis revealed a common ancestry for the ductal and small cell carcinoma. Case 2 showed 13q loss (involving RB1) in both adenocarcinoma and small cell carcinoma regions, and 3p gain, 4p loss, 17p loss (involving TP53) in the latter. By using highly curated samples, we demonstrate for the first time that small cell neuroendocrine and ductal prostatic carcinoma can have a common an-cestry and the process of evolution over time. We highlight whole genome doubling in a patient with prostate cancer relapse, reinforcing its poor prognostic nature.

Study EGAS00001007428

15388 results for "1 mil coins fc 26 Visit Buyfc26coins.com for latest FC 26 coins news..buoc"

in 147.92 milliseconds.