Search Results - EGA European Genome-Phenome Archive

Center for Inherited Disease Research (CIDR)-National Institute on Aging (NIA) Whole Exome Analysis of Ehlers-Danlos Syndrome

The Ehlers-Danlos syndromes (EDS) are a family of hereditary disorders of connective tissue characterized by joint hypermobility and skin involvement with soft, translucent, stretchy skin, and atrophic scarring. To date, over 12 different types of EDS are known, and the genes causing almost all of them have been found. However, the gene(s) underlying the most common of Ehlers-Danlos syndrome, known as the hypermobile type, remain elusive. The hypermobile type of Ehlers-Danlos syndrome is associated with multiple co-morbidities, including dysautonomia, chronic musculoskeletal pain, gastrointestinal dysmotility, mast cell activation and multiple neurologic complications. Affected persons often suffer from chronic fatigue and may be chronically disabled. While the precise frequency of the condition is not known, it may be as common as between 1/500 and 1/1000 persons. There is a significant overlap with the phenotype of fibromyalgia, and multiple persons diagnosed with fibromyalgia may well have the hypermobile type of Ehlers-Danlos syndrome. The phenotype of hypermobile EDS segregates as a single dominant trait with complete penetrance, variable expressivity, and a markedly evolving phenotype over time. Identification of the genetic cause of this hereditary disorder has the potential to inform rational therapy for this disabling condition. It is hoped that the insight into the pathogenesis of the hypermobile type of EDS will lead to a similar breakthrough in therapy for this disorder. Toward this end, whole exome sequencing (WES) of DNA from a total of 167 patients will be performed at the Center for Inherited Disease Research (CIDR) to identify genes underlying the hypermobile type of Ehlers-Danlos syndrome and those with EDS of an uncharacterized type.

Study phs001779

Clonal Lineage Tracing of Primary Human Cortical Progenitors, Cell Type Profiles in the Brain Vasculature and Genotyping Data for 22q11DS and Control iPS Lines

We sought to characterize the cell types produced by individual human cortical progenitors and describe their lineage relationships. We labeled individual human cortical cells with a massively complex lentivirus library encoding a GFP reporter gene and a molecular barcode that could be captured by standard 3' scRNA-seq methodologies. We then cultured the labeled human cells for 6 weeks in either a 1) mouse cortical astrocyte co-culture system or 2) in the mouse cortex as a xenograft. In a separate experiment, we maintained human cortical tissue slices in organotypic slice culture for 12 days. GFP+ cells were then isolated by FACS and then captured on the 10x chromium scRNA-seq platform. An additional sample was labeled by lentivirus and briefly cultured for 3 days prior to being mixed with mouse 3T3 cells in order to perform a cell mixing experiment ("barnyard") to benchmark our methodologies. We found that human individual cortical progenitors gave rise to both excitatory and inhibitory neurons as well as glia.We also sought to characterize the adult brain vasculature in both vascular malformations and non-malformed tissue. We use the 10x chromium scRNA-seq platform across 10 individuals (5 control and 5 arteriovenous malformations). We find broad disregulation across multiple cell types including the immune cells and find a specific monocyte subtype to contribute to a rupture phenotype. We also find a novel population of fibromyocytes in both conditions.Finally, we confirmed the genotype status of iPS lines that were derived from either karyotypically normal control individuals or from DiGeorge Syndrome patients carrying the 22q11.2 deletion using the Illumina Infinium Global Screening Array v3.0.

Study phs002624

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

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

Study EGAS00001002975

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

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

Study EGAS00001002430

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

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

Study EGAS00001002275

We evaluate the potential for routine WGS using ONT by sequencing the well-characterised reference sample NA12878 and the genome of an individual with ataxia-pancytopenia syndrome accompanied by severe immune dysregulation.

Whole-genome sequencing (WGS) is becoming widely used in clinical medicine in diagnostic contexts and to inform treatment choice. While current sequencing technologies have been extremely successful, their reliance on short read lengths necessarily involves some limitations: accurately assaying certain genomic regions and classes of variation can be problematic. Recent advances in throughput and cost have made WGS using the Oxford Nanopore Technologies (ONT) MinION long-read single-molecule sequencer a potential solution to these challenges. Here we evaluate the potential for routine WGS using ONT by sequencing the well-characterised reference sample NA12878 and the genome of an individual with ataxia-pancytopenia syndrome accompanied by severe immune dysregulation, to 82× and 30× respectively. For NA12878, we evaluated single-nucleotide variant (SNV) calls based on data from multiple base-calling algorithms. We demonstrate that phasing metrics from a novel, reference panel-free, long-read-based method can improve variant-calling performance from otherwise modest levels, resulting in a false discovery rate of 7.1% and false negative rate of 8.5%; remaining errors are concentrated near homopolymers and in regions of reduced sequencing coverage. In the clinical sample, we are able to identify and directly phase two non-synonymous de novo variants in SAMD9L, (OMIM #159550) inferring that they both lie on a common paternal haplotype by overlap with parental genotypes at nearby common variants. This work demonstrates that methodological innovation can substantially reduce variant-calling error rates in ONT data, and that with on-going improvements in throughput, base-calling and dedicated long-read-based SNV-calling methodology, ONT offers promise as an option for clinical WGS.

Study EGAS00001003469

A reference map of potential determinants for the human serum metabolome

The collection of metabolites circulating in the human blood, termed the serum metabolome, contains a plethora of biomarkers and causative agents. Although the origin of specific compounds is known, we have a poor understanding of the key determinants of most metabolites. Here, we measured the levels of 1251 circulating metabolites in serum samples from a healthy human cohort of 491 individuals, and devised machine learning algorithms to predict their levels in held-out subjects based on a comprehensive profile consisting of host genetics, gut microbiome, clinical parameters, diet, lifestyle, anthropometric measurements and medication data. Notably, we obtained statistically significant predictions for over 76% of the profiled metabolites. Despite using the strict out-of-sample prediction metric, which is a lower bound for the explained variance, diet and microbiome each explained hundreds of metabolites, with over 50% of the variance explained in some metabolites. We further validated the robustness of the microbiome related associations by showing a high replication rate in two geographically independent cohorts that were not available to us when developing the algorithms. We also demonstrate that some of these interactions are causal, as some metabolites we predicted to be positively associated with bread increased in level following a randomized clinical trial of bread intervention. Microbiome-explained metabolites were enriched with unnamed metabolites, and we devised an algorithm that accurately predicts their biological pathway, finding that they mainly include food components, aromatic amino acids and secondary bile acid derivatives. Overall, our results unravel potential determinants of over 800 metabolites, paving the way towards mechanistic understanding of alterations in metabolites under different conditions and to designing interventions for manipulating circulating metabolite levels.

Study EGAS00001004512

Genome-Wide associations of Lung Health Study (LHS)

The 'Genome-Wide Associations Environmental Interactions in the Lung Health Study' at Johns Hopkins University aims to test for association between lung function decline as a primary outcome associated with chronic obstructive pulmonary disease (COPD) using banked DNA and phenotype data on 4,287 European Americans from the longitudinal, multicenter Lung Health Study (LHS). The broad goals of the LungGO/ESP-GO falls into two general categories: (i) discovery of all variants (i.e., common and rare) in all protein-coding regions of the human genome (i.e., the exome) conferring risk to complex pulmonary diseases including COPD, in a subset of the LHS cohort. The Johns Hopkins University LHS cohort offers a unique opportunity to elucidate genetic variants that cause COPD. The Lung Health Study I was a randomized multicenter clinical trial with 5887 participants carried out from October 1986 to April 1994, designed to test the effectiveness of smoking cessation and bronchodilator administration in smokers aged 35 to 60 with mild lung function impairment. Participants were randomly assigned to one of three groups: usual care, who received no intervention smoking intervention with the inhaled bronchodilator ipratroprium bromide smoking intervention with an inhaled placebo. The effect of intervention was evaluated by the rate of decline of forced expiratory volume in one second (FEV1). For the GWAS, only the subset of European American LHS participants for whom lung function data from three time points or more are available. Thus, the GWAS represents 73% of the 5,887 volunteers who participated in the LHS study. Importantly, LHS subjects included had similar demographics (including age, gender and BMI) and rates of lung function decline (mean annual change in FEV1% predicted: -0.96 %/yr vs. -0.99 %/yr, p=0.57) compared with those not included in the GWAS, reflecting little selection bias for our primary outcome. They were, however, more likely to have quit smoking after 5 years. This study is part of the Gene Environment Association Studies initiative (GENEVA, http://www.genevastudy.org) funded by the trans-NIH Genes, Environment, and Health Initiative (GEI). The overarching goal is to identify novel genetic factors that contribute to lung function through large-scale genome-wide association studies of smokers enrolled in a multicenter clinical trial. Genotyping was performed at the Johns Hopkins University Center for Inherited Disease Research (CIDR). Data cleaning and harmonization were done at the GEI-funded GENEVA Coordinating Center at the University of Washington.

Study phs000335

North American Mitochondrial Disease Consortium Patient Registry and Biorepository

Mitochondrial diseases are caused by dysfunction of the mitochondria, which are specialized compartments that are present in every cell of the body except red blood cells. Mitochondria generate more than 90% of the energy that the body needs to sustain life and support growth. When they fail, less and less energy is generated within the cell. This injures the cell and can cause its death. If this process is repeated throughout the body, whole organ systems begin to fail, and the life of the person in whom this is happening is severely compromised. Mitochondrial diseases primarily affect children, but adult onset is becoming more and more common. Mitochondrial diseases are probably the most diverse human disorders at every level: clinical, biochemical, and genetic. Some affect only the nervous system but most affect many body systems, including the brain, heart, liver, skeletal muscles, kidney, and the endocrine and respiratory systems. Although mitochondrial disorders vary in severity, they are usually progressive, and often crippling. They can cause paralysis, seizures, mental retardation, dementia, hearing loss, blindness, weakness and premature death. Because of the range of symptoms and the frequent involvement of multiple body systems, mitochondrial diseases can be a great challenge to diagnose. Even when accurately diagnosed, they pose an even more formidable challenge to treat, as there are very few therapies and most are only partially effective. About this Study The first objective of this study is to establish a clinical registry of patients with suspected or confirmed mitochondrial diseases. We are collecting medical and family history, diagnostic test results, and prospective medical information for these patients and, using agreed procedures developed by the leading research clinicians in the field, providing, for the first time, standardized diagnoses of these complex disorders for the patients. The clinical information we collect from the participants will be used to learn about the spectrum of mitochondrial disorders and their prevalence. We will also develop studies which allow us to better understand how these diseases progress, which we do not understand well enough. When we begin clinical trials for mitochondrial diseases, patients enrolled in the registry who are identified as potentially eligible will be offered enrollment. Patients will only be included in studies if they give their consent in advance. The second objective of this study is to establish a biorepository for specimens and DNA from patients with mitochondrial diseases, in order to make materials easily available to consortium researchers.

Study phs001538

Cystic Fibrosis Nasal Epithelium Gene Expression by RNAseq

Cystic fibrosis (CF) is a monogenic disorder leading to multi-organ system dysfunction, with the primary cause of morbidity and mortality being lung disease. Mutations in CFTR result in defective ion channel transport leading to hyperviscous mucus, airway infection, inflammation and lung function decline. While the Phe508del mutation is the most common CFTR mutation in North American people with CF, it is well known that there is a broad range of lung disease severity even amongst Phe508del homozygotes, including twins/siblings. Heritability analyses showing the relative contribution of non-CFTR modifier genes in such differences in phenotype, have led to our genomic and transcriptomic studies of gene modifiers in CF. Complementary to two international CF GWAS studies, this study was performed to a) determine differential expression of genes relative to a quantitative lung disease severity phenotype score (KNoRMA; as used in CF GWAS studies) and b) to assess expression quantitative trait loci. To better understand genetic modification of CF lung disease severity, we analyzed gene expression (transcriptomic) profiles by mRNA sequencing from nasal mucosal biopsy specimens of 134 people with CF. Each study participants was characterized for CF lung disease severity phenotype by averaging three years of lung function data (forced expiratory volume in 1 second; FEV1), normalized to a CF-specific reference population matching for age, sex, and height. Participants completed nasal lavage and nasal curettage mucosal sampling, with directed sampling to the inferior nasal turbinate, a well validated surrogate for cellular characteristics and CFTR function of the lower CF airways. Participants were further invited to complete a second set of these sampling procedures for validation analysis. Specimen and concurrent clinical data collection occurred across four North American study sites and nasal mucosal specimens were subjected to RNASeq analysis. These data were to be made available in the dbGaP. Our association analyses of gene expression relative to quantitative CF lung disease phenotype, taken together with previously identified GWAS signals, identified biological pathways correlated with CF lung disease pathogenesis. Our further association analyses of gene expression relative to single nucleotide polymorphism variation at significant CF GWAS loci further identified potential mechanisms of genetic variants associated with CF lung disease. This resource is intended to aid researchers in characterizing airway mucosal gene expression from the CF airways for the purpose of gene expression profiling in CF and in muco-obstructive airways disease.

Study phs002254

Women's Environmental Cancer and Radiation Epidemiology (WECARE) Study

Study population The WECARE Study is an international multicenter population-based case-control study of cases with asynchronous contralateral breast cancer (CBC) and individually matched controls with unilateral breast cancer (UBC). Recruitment and data collection for the WECARE Study were conducted in two phases, herein referred to as WECARE I (2001–2004) and WECARE II (2009–2012). The study design of the first phase (WECARE I) has been described in detail elsewhere [Bernstein JL, et al 2004, PMID: 15084244]; the second phase (WECARE II) employed a nearly identical approach [Langballe R, et al 2016, PMID: 27400983]. Participants in each phase were identified through eight population-based cancer registries, including six in the United States that contribute to the National Cancer Institute's Surveillance, Epidemiology, and End Results (SEER) program: Los Angeles County Cancer Surveillance Program; Cancer Surveillance System of the Fred Hutchinson Cancer Center (Seattle region, WA); State Health Registry of Iowa; Cancer Surveillance Program of Orange County/San Diego-Imperial Organization for Cancer Control (Orange County/San Diego, CA); the Greater Bay Area Cancer Registry (San Francisco Bay Area region and Santa Clara region, CA); and the Sacramento and Sierra Center Registry (Sacramento region, CA). Participants were additionally identified using the Ontario Cancer Registry (Canada) and the Danish Breast Cancer Cooperative Group Registry, supplemented by data from the Danish Cancer Registry. The study protocol was approved by the institutional review or regulatory ethics boards at each study site. Data collection Study participants were interviewed by telephone using a structured questionnaire aimed at evaluating known or suspected breast cancer risk factors, including personal demographics, medical history, menstrual and reproductive history, family history of cancer, use of hormones, smoking, and alcohol intake. Risk factor status was assessed during the period prior to first diagnosis, as well as between first diagnosis and reference date (i.e., the at-risk period for CBC). Detailed data on treatment and tumor characteristics were obtained directly from cancer registry records or by abstracting medical records, including pathology and surgical reports, radiation oncology clinic notes, and systemic adjuvant treatment data. Complete medical treatment history information was collected, and for all women who received radiotherapy, the radiation dose to the contralateral breast was reconstructed using radiotherapy records and radiation measurements. Biospecimens were collected from study participants for genotyping: blood collected from WECARE I women; and saliva collected from WECARE II women.

Study phs003945

Exome sequencing reveals pathogenic variants in known and novel candidate genes for severe sperm motility disorders

STUDY QUESTION: What are the causative genetic variants in patients with male infertility due to severe sperm motility disorders? SUMMARY ANSWER: We identified high confidence disease-causing variants in multiple genes previously associated with severe sperm motility disorders in 9 out of 21 patients (43%). We also identified variants in novel candidate genes in 8 additional patients (38%). WHAT IS KNOWN ALREADY: Severe sperm motility disorders are a form of male infertility characterised by vital, but immotile sperm often in combination with a spectrum of structural abnormalities of the sperm flagellum. Examples of this spectrum of disorders include dysplasia of the fibrous sheath / multiple morphological abnormalities of the sperm flagella (DFS-MMAF) and severe forms of asthenozoospermia associated to non-specific structural defects. These disorders frequently recur in families and may include chronic respiratory disease. As such, genetic causes are suspected and many genes have recently been reported. Currently, depending on the clinical sub-categorisation, up to 50% of causality in patients with severe sperm motility disorders can be explained by pathogenic variants in at least 21 known genes. STUDY DESIGN, SIZE, DURATION: We performed exome sequencing in patients with severe sperm motility disorders from two different clinics, aiming to identify the underlying genetic defects and establishing a diagnostic yield in both. PARTICIPANTS/MATERIALS, SETTING, METHOD: Two groups of infertile men, one from Argentina (n= 9) and one from Australia (n=12), with clinically defined sperm motility disorders were included. All patients in the Argentine cohort were diagnosed with DFS-MMAF, based on transmission electron microscopy. Exome sequencing was performed in all 21 patients and variants with an allele frequency of <1% in the gnomAD population were prioritised and interpreted. MAIN RESULTS AND ROLE OF CHANCE: In 9/21 patients (43%), we identified pathogenic variants in known sperm motility disorder genes: CFAP43 (3 patients); CFAP44 (2 patients), QRICH2 (2 patients), DNAH1 (1 patient) and DNAH6 (1 patient). The diagnostic rate did not differ markedly between the Argentinian and the Australian cohort (44% and 42%, respectively). Furthermore, we identified patients with variants in the novel candidate sperm motility genes CFAP58, DNAH12, DRC1, MDC1, PACRG, SSPL2C and TPTE2. One patient presented with two possibly damaging variants in four candidate genes and it remains unclear which variants were responsible for the severe sperm motility defect in this patient. LIMITATIONS, REASONS FOR CAUTION: In this study, we described patients with potential bi-allelic variants in known and novel candidate genes for severe sperm motility disorders. Due to unavailability of parental DNA, we have not assessed de novo or maternally inherited dominant variants and could not phase the mutations to establish in all cases that the mutations occur on both alleles. WIDER IMPLICATIONS OF THE FINDINGS: Our results confirm an important role for five known genes for sperm motility and we demonstrate that exome sequencing is an effective method to diagnose patients with severe sperm motility defects (9/21 diagnosed; 43%). Furthermore, our analysis revealed seven novel candidate genes for severe sperm motility disorders. Genome-wide sequencing of additional patient cohorts and re-analysis of exome data of currently unsolved cases may reveal additional variants in these novel candidate genes.

Study EGAS00001005018

GabonDiv_SNP_genotyping

Genome-wide SNP genotyping data for 1,235 western Africans by Illumina HumanOmniExpress-12 array, used in the EGAS00001002078 study

Dataset EGAD00010001209

cqmuGWAS2

DNA for 2482 individuals from Chongqing was extract from peripheral blood and genotyped by Illumina Omni Zhonghua-8 version 2 gene chips.

Dataset EGAD00010001526

Genetic and transcriptomic landscape of DLBCL

Affymetrix SNP6.0 data for 341 DLBCL patients

Dataset EGAD00010001980

BelCovid_2_genotype

Cohort: Raw genotype files for BelCovid2 cohort. Genotype Chip: Illumina Global Screening Array-24 v3.0 + Multi-Disease beadchip genomic build: b37

Dataset EGAD00010002179

GEN_COVID_genotype

Cohort: Raw genotype files for GEN_COVID cohort. Genotype Chip: Illumina Global Screening Array-24 v3.0 + Multi-Disease beadchip genomic build: b37

Dataset EGAD00010002177

RealSeqS ovarian amplicon counts

Dataset contains text files that describe the chromosome, position, and read count for an amplicon using the RealSeqS assay.

Dataset EGAD50000000013

IMvigor130 heatmap Figure 2A

Source data for Figure.2A heatmap showing the gene expression across samples.

Dataset EGAD50000000139

Time course acetalax bisacodyl treatment

RNA-Seq raw data of PDX treated with acetalaxor bisacodyl for 4h00 to 24h00.

Dataset EGAD50000000875

HIV Viremic Non-Progressors (VNPs) and HIV Progressors Data Access Committee

Data Access Committe for the project: Analyzing exome sequencing in HIV Viremic Non-Progressors (VNPs) and HIV Progressors

Dac EGAC50000000062

A Risk Score Incorporating Low Pass Whole Genome Sequencing of Cell Free DNA from Patients Receiving CD19 CAR T-Cell Therapy for Large B-Cell Lymphoma

Study EGAS00001006308

12024 results for "coin buying fc 26 Visit Buyfc26coins.com for latest FC 26 coins news..2Dtf"

in 55.73 milliseconds.