Search Results - EGA European Genome-Phenome Archive

Autozygosity_pilot___British_Pakistani_from_Birmingham_2

The offspring of first cousin marriages have ~6% of their genome autozygous, i.e. homozygous identical by descent, or even more if there was further consanguinity in their ancestry. In the UK there are large populations with very high first cousin marriage rates of 20-50%. Sequencing the exomes of a sample of these individuals has the potential both to support genetic health programmes in these populations, and to provide genetic research information about rare loss of function mutations. This pilot study based on existing British-Pakistani cohort samples from Birmingham will identify homozygous individuals for almost all variants down to an allele frequency around 1%, plus individuals carrying hundreds of new homozygous rare loss-of-function variants, and will support development of community relations and ethics for a wider study currently being designed. The data deposited in the EGA consists of low coverage whole exome sequencing on these samples.

Study EGAS00001000567

Autozygosity_pilot___Pakistani_from_Birmingham

The offspring of first cousin marriages have ~6% of their genome autozygous, i.e. homozygous identical by descent, or even more if there was further consanguinity in their ancestry. In the UK there are large populations with very high first cousin marriage rates of 50-80%. Sequencing the exomes of a sample of these individuals has the potential both to support genetic health programmes in these populations, and to provide genetic research information about rare loss of function mutations. This pilot study based on existing British-Pakistani cohort samples from Birmingham will identify homozygous individuals for almost all variants down to an allele frequency around 1%, plus individuals carrying hundreds of new homozygous rare loss-of-function variants, and will support development of community relations and ethics for a wider study currently being designed. The data deposited in the EGA consist of low coverage whole exome sequencing on these samples.

Study EGAS00001000511

Autozygosity_pilot___QMUL

The offspring of first cousin marriages have ~6% of their genome autozygous, i.e. homozygous identical by descent, or even more if there was further consanguinity in their ancestry. In the UK there are large populations with very high first cousin marriage rates of 20-50%. Sequencing the exomes of a sample of these individuals has the potential both to support genetic health programmes in these populations, and to provide genetic research information about rare loss of function mutations. This pilot study based on existing British-Pakistani cohort samples will identify homozygous individuals for almost all variants down to an allele frequency around 1%, plus individuals carrying hundreds of new homozygous rare loss-of-function variants, and will support development of community relations and ethics for a wider study currently being designed. The data deposited in the EGA consists of low coverage whole exome sequencing on these samples.

Study EGAS00001000717

Whole Genome Sequencing to track subclonal heterogeneity in 18 samples from 3 Chronic Lymphocytic Leukemia patients subjected to repeated cycles of therapy.

This study took 18 samples from three patients with Chronic Lymphocytic Leukemia. The 3 CLL patients received multiple different treatments sequentially for a period of up to 7 years. Peripheral blood samples were taken at 5 specific time points during disease progression together with 1 matched buccal swab per patient. All the selected patients had unmutated IgVH status reflecting poor prognosis. However, none of the patients had complex genomic aberrations associated with aggressive disease. Peripheral blood samples and buccal smears were obtained from patients, who had given informed consent for sequential sample collection and analysis. Tumour samples were obtained before treatment, after treatment and at relapse and subjected to Ficoll separation. Lymphocyte cell purity was assessed by full blood count enumeration. Buccal smears were obtained and processed using Isohelix.

Study EGAS00001000885

Recapitulation of genetic predisposition to medulloblastoma in human neuroepithelial stem cells

Human neural stem cell cultures provide progenitor cells that are potential cells of origin for brain cancers. However, the extent to which genetic predisposition to tumor formation can be faithfully captured in stem cell lines is uncertain. Here we evaluated neuroepithelial stem (NES) cells, representative of cerebellar progenitors. We transduced NES cells with MYCN, observing medulloblastoma upon orthotopic implantation in mice. Significantly, transcriptomes and patterns of DNA methylation from xenograft tumors were globally more representative of human medulloblastoma compared to a MYCN-driven genetically engineered mouse model. Orthotopic transplantation of NES cells generated from Gorlin syndrome, patients predisposed to medulloblastoma due to germline mutated PTCH1, also generated medulloblastoma. We engineered candidate cooperating mutations in Gorlin NES cells, with mutation of DDX3X or loss of GSE1 both accelerating tumorigenesis. These findings demonstrate that human NES cells provide a potent experimental resource for dissecting genetic causation in medulloblastoma.

Study EGAS00001003620

Targeted sequencing of 12 genes in patients with HLH

Hemophagocytic lymphohistiocytosis (HLH) is a rapid-onset, potentially fatal hyperinflammatory syndrome. A prompt molecular diagnosis is crucial for appropriate clinical management. Here, we validated and prospectively evaluated a targeted high-throughput sequencing approach for HLH diagnostics. We have carried out high-throughput sequencing of 12 genes linked to HLH in 58 HLH patients. We achieved a diagnosis in 22 out of 58 patients (38%). Genetically undiagnosed HLH patients had a later age at onset and manifested higher frequencies of known secondary HLH triggers. Rare, putatively pathogenic monoallelic variants were identified in 9 patients. Almost all cases with reduced natural killer (NK) cell function received a diagnosis, but the majority of the prospective cases remain genetically unexplained, highlighting genetic heterogeneity and environmental impact within HLH. A complete understanding of the genetic susceptibility to HLH thus requires further in-depth investigations, including genome sequencing and detailed immunological characterization.

Study EGAS00001001605

RNA-seq consisting of FASTQ paired-end reads from cancer samples

Efforts to precisely identify tumor human leukocyte antigen presented peptides (HLAp) capable of mediating T cell based tumor rejection still face important challenges. Recent reports suggest that non-canonical cancer HLAp could be immunogenic but their identification requires highly sensitive and accurate mass-spectrometry (MS)-based proteogenomics approaches. Here, we present a MS-based analytical pipeline that can precisely characterize the non-canonical HLAp repertoire, incorporating whole exome sequencing, bulk and single cell transcriptomics, ribosome profiling, and a combination of two MS/MS search tools. This approach results in the accurate identification of hundreds of shared and tumor-specific non-canonical HLAp. Albeit often at low levels and in distinct subpopulations of cells, numerous non-canonical HLAp are shared across tumors. This analytical platform holds great promise for the discovery of novel cancer antigens for cancer immunotherapy.

Study EGAS00001003724

Cancer and germline exomes consisting of FASTQ paired-end reads from melanoma and lung cancer samples

Efforts to precisely identify tumor human leukocyte antigen presented peptides (HLAp) capable of mediating T cell based tumor rejection still face important challenges. Recent reports suggest that non-canonical cancer HLAp could be immunogenic but their identification requires highly sensitive and accurate mass-spectrometry (MS)-based proteogenomics approaches. Here, we present a MS-based analytical pipeline that can precisely characterize the non-canonical HLAp repertoire, incorporating whole exome sequencing, bulk and single cell transcriptomics, ribosome profiling, and a combination of two MS/MS search tools. This approach results in the accurate identification of hundreds of shared and tumor-specific non-canonical HLAp. Albeit often at low levels and in distinct subpopulations of cells, numerous non-canonical HLAp are shared across tumors. This analytical platform holds great promise for the discovery of novel cancer antigens for cancer immunotherapy.

Study EGAS00001003723

Profiling_molecular_heterogeneity_in_human_primary_microglia

In the brain the cells that control inflammation are called a type of white blood cell called microglia. Microglia are located throughout the brain and spinal cord and account for 10–15% of all cells found within the brain. As the resident white blood cells, they are the main active immune defence in the central nervous system (CNS). Microglia are part of an important class of cells known as macrophages that have two main states: M1 and M2. M1 cells are pro- inflammatory, leading to more inflammation, while M2 are anti-inflammatory, and drive wound healing. In this study, we will collect primary microglia from surgical biospies of 100 individuals. 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 EGAS00001002494

Single cell RNAseq of PBMC from bladder cancer patients

Frozen PBMC samples containing at least 1 million cells were thawed for 1 minute at 37C and washed twice with RPMI complete media (10% FBS with glutamate and Pen/Strep). All of the samples had >80% viable cells. Sample processing for single-cell RNA-seq was done using Chromium Single Cell 3’ Library and Gel bead kit v2 (PN-120237) following manufacturer’s user guide (CG00052, 10x Genomics, Pleasanton, CA). The total cell density was used to impute the volume of single cell suspension needed in the reverse transcription (RT) master mix, aiming to achieve ~ 6,000 cells per sample. cDNAs and libraries were prepared following manufacturer’s user guide (10x Genomics). cDNA amplification and indexed libraries were prepared using 12 and 14 cycles of PCR, respectively. Libraries were profiled, quantified, and sequenced as 5’ single-cell gene expression libraries.

Study EGAS00001004008

VIKING Health Study - Shetland

The Viking Health Study – Shetland, together with ORCADES collectively termed VIKING, aims to identify genes influencing risk factors for common diseases such as heart disease, diabetes, glaucoma and stroke. The genetic research cohort study in Shetland was established to identify genetic variants, both common and rare, increasing risk of disease. This isolated population has a number of characteristics, the most important being the very large number of distant relatives, which are favourable for the identification of rare variants influencing disease risk. This study describes PCR-free paired end whole genome sequencing (TruSeq DNA PCR-Free, Illumina) run on a HiSeqX platform at 30X coverage on 500 individuals from VIKING. Unrelated individuals from the largest families were selected first, followed by those from smaller families, until eventually related individuals were selected to best represent the variation in the full cohort.

Study EGAS00001003872

Splicing patterns in SF3B1 mutated uveal melanoma generate shared immunogenic tumor-specific neo-epitopes

Disruption of splicing patterns due to mutations of genes coding splicing factors in tumors represents a potential source of tumor neo-antigens, which would be both public (shared between patients) and tumor-specific (not expressed in normal tissues). In this study, we show that mutations of the splicing factor SF3B1 in uveal melanoma (UM) generate such immunogenic neo-antigens. Memory CD8+ T cells specific for these neo-antigens are preferentially found in 20% of UM patients bearing SF3B1 mutated tumors. Single cell analyses of neo-epitope specific T cells from the blood identified large clonal T cell expansions, with distinct effector transcription patterns. Some of these expanded TCRs are also present in the corresponding tumors. CD8+ T cell clones specific for the neo-epitopes specifically recognize and kill SF3B1-mutated tumor cells, supporting the use of this new family of neo-antigens as therapeutic targets.

Study EGAS00001005226

Accurate mapping of mitochondrial DNA deletions and duplications using deep sequencing

Deletions and duplications in mitochondrial DNA (mtDNA) cause mitochondrial disease and accumulate in conditions such as cancer and age-related disorders, but validated high-throughput methodology that can readily detect and discriminate between these two types of events is lacking. Here we present MitoSAlt, a computational method for accurate identification, quantification and visualization of mtDNA deletions and duplications from whole genome, whole exome or transcriptome sequencing data. MitoSAlt was tested on simulated sequencing reads and human patient samples with single deletions and duplications to verify its accuracy. Application to mouse models of mtDNA maintenance disease further demonstrated the ability to detect deletions and duplications even at low levels of heteroplasmy. MitoSAlt paves the way for simple and reliable determination of mtDNA deletions and duplications across a wide range of relevant conditions and available sequencing datasets.

Study EGAS00001004380

Single cell RNAseq of PBMC from RCC patients

Frozen PBMC samples containing at least 1 million cells were thawed for 1 minute at 37C and washed twice with RPMI complete media (10% FBS with glutamate and Pen/Strep). All of the samples had &gt;80% viable cells. Sample processing for single-cell RNA-seq was done using Chromium Single Cell 3’ Library and Gel bead kit v2 (PN-120237) following manufacturer’s user guide (CG00052, 10x Genomics, Pleasanton, CA). The total cell density was used to impute the volume of single cell suspension needed in the reverse transcription (RT) master mix, aiming to achieve ~ 6,000 cells per sample. cDNAs and libraries were prepared following manufacturer’s user guide (10x Genomics). cDNA amplification and indexed libraries were prepared using 12 and 14 cycles of PCR, respectively. Libraries were profiled, quantified, and sequenced as 5’ single-cell gene expression libraries.

Study EGAS00001004451

Epigenetic profiling of colorectal cancer initiating cells (CC-ICs) to identify bivalently marked genes (H3K4me3 and H3K27me3 ChIP-seq), and investigation of changes in transcriptome following EZH2 inhibition using RNA-seq.

Epigenetic regulation of transcription plays a crucial role in lineage commitment of embryonic stem cells. Promoters of key lineage-specific differentiation genes are found in a repressed bivalent state, having both activating H3K4me3 and repressive H3K27me3 histone marks, making them poised for transcription upon loss of H3K27me3 in response to environmental cues. Whether the tumour-initiating, self-renewing, cancer-initiating cells (C-ICs) have similar epigenetic regulatory mechanism that prevent lineage commitment is unknown. In order to investigate bivalently marked and repressed promoters, we used a patient-derived CC-IC enriched model to identify the changes in transcriptome following inhibition of EZH2, the H3K27 methyltransferase. We also performed ChIP-seq for H3K27me3 and H3K4me3 at baseline in order to identify repressed and bivalently marked promoters.

Study EGAS00001003003

SudanMitoSeq: Sudanese mitochondrial sequencing

In various contexts, mitochondrial function or dysfunction can be linked to mitochondrial genome variations. The use of mitochondrial genetics thus promises personalized diagnostics and treatments. In order to devise specific precision medicine approaches based on mitochondrial genetic variation and to test them within clinical studies, control data is indispensable. Such control data needs to comprehensively cover genetic variation commonly observed and thus expected. In the context of this study we assessed whether current, world-wide mitochondrial data sufficiently represents the region of North and East Africa, that is, whether current reference data is ready for precision medicine in this region. Towards this, we sequenced mitochondrial genomes of 159 Sudanese individuals provided as part of this EGA study and analyzed them together with various other, publicly available data concerning mitochondrial variation. For details, please refer to the publication.

Study EGAS00001005669

Hereditary Cancer Diagnostics with I2HCP gene panel

This dataset contains targetted DNA sequencing data generated in the context of genetic diagnostics of hereditary cancer. Target enrichment was performed with the I2HCP, a custom hereditary cancer gene panel based in Agilent SureSelect and developed by the Germans Trias i Pujol Research Institute (IGTP, Badalona, Catalonia, Spain) and the Catalan Institute of Oncology (ICO, L'Hospitalet de Llobregat, Catalonia, Spain). After the enrichment, 108 samples were sequenced in Illumina HiSeq machines and 130 in MiSeq machines. Some of the samples were sequenced on both platforms. The dataset also includes MLPA validation data for all identified copy-number alterations, annotated with exon resolution and information about negative MLPA results. All this data was generated in the context of routine diagnostics and compiled together with MPLA validation data for a genetic diagnostics oriented benchmark of germline CNV calling tools.

Study EGAS00001004316

Spatial atlas of clonal copy number alterations in co-existing benign and malignant tissue

Defining the transition from benign to malignant tissue is fundamental to improve early diagnosis of cancer. Here, we use a systematic approach to study spatial genome integrity in situ and describe previously unidentified clonal relationships. We employed spatially resolved transcriptomics to infer spatial copy number variations in >120,000 regions across multiple organs, in benign and malignant tissues. We demonstrate that genome-wide copy number variation reveals distinct clonal patterns within tumours and in nearby benign tissue using an organ-wide approach focused on the prostate. Our results suggest a model for how genomic instability arises in histologically benign tissue that may represent early events in cancer evolution. We highlight the power of capturing the molecular and spatial continuums in a tissue context and challenge the rationale for treatment paradigms, including focal therapy.

Study EGAS00001006124

Multi-omic Profiling of Central Nervous System Leukemia Identifies mRNA Translation as a Therapeutic Target

Central nervous system (CNS) dissemination of B-precursor acute lymphoblastic leukemia (B-ALL) has poor prognosis and remains a therapeutic challenge. Here we performed targeted DNA sequencing, transcriptional and proteomic profiling of paired leukemia infiltrating cells in the bone marrow (BM) and CNS of xenografts. Genes governing mRNA translation were upregulated in CNS leukemia, and subclonal genetic profiling confirmed this in both BM-concordant and BM-discordant CNS mutational populations. CNS leukemia cells were exquisitely sensitive to the translation inhibitor omacetaxine mepesuccinate, which reduced xenograft leptomeningeal disease burden. Proteomics demonstrated greater abundance of secreted proteins in CNS infiltrating cells including complement component 3 (C3) and drug targeting of C3 influenced CNS disease in xenografts. CNS infiltrating cells also exhibited selection for stemness traits and metabolic reprogramming. Overall, our study identifies therapeutic targeting of mRNA translation as a novel therapeutic approach for B-ALL leptomeningeal disease.

Study EGAS00001005647

Gene copy number variation in pediatric mental illness in a general population

We assessed the relationship of gene copy number variation (CNV) in mental health/neurodevelopmental traits and diagnoses, physical health, and cognitive biomarkers in a community sample of 7,100 unrelated European, and East Asian children and youth (Spit for Science). Diagnoses of mental health disorders were found in 17.5% of participants and 27.6% scored in the highest 10% on either or both ADHD and OCD trait measures. Clinically relevant CNVs were present in 3.9% of participants and were associated with elevated scores on a continuous measure of ADHD (p=5.0×10-3), on a cognitive biomarker of mental health (response inhibition (p=1.0×10-2)), and on prevalence of mental disorders (p=1.9×10-6, odds ratio: 3.09). With a rise of mental illness, our data establishes a baseline for delineating genetic contributors in paediatric-onset conditions.

Study EGAS00001006659

The impact of mutational clonality in predicting the response to anti-PD-L1/PD-L1 in advanced urothelial cancer

The advent of immune checkpoint inhibitors (ICI) has radically changed the management and therapeutic treatment landscape for patients with cancer. It has been shown that the use of these inhibitors can result in long-term complete remissions even in cases with advanced stages of the disease. However, only a small fraction of patients respond to the treatment. A better understanding of which factors could drive this clinical benefit is fundamental to building more accurate predictive biomarkers and developing alternative therapeutic options for patients who are unlikely to benefit. In order to address this issue, we have generated and analyzed whole exome sequencing (WES) and RNA sequencing (RNASeq) data from a cohort of metastatic urothelial carcinoma patients (N = 27) treated with ICI such as anti-PD-(L)1 monoclonal antibodies.

Study EGAS00001007086

Epithelial, fibroblast, myeloid, T cell, primary prostate cancer

The collection and use of tissue for this study had Melbourne Health institutional review board approval and patients provided written informed consent (Melbourne Health Local Project Number: 2016.087). Following the prostatectomy of 13 patients, ranging from 52 to 78 years of age and from CAPRA-S risk score of 0 (attributed to benign tissue samples, harvested from a site far from a low grade, low volume cancer) to 7 (Supplementary file 2), a four millimeter tissue core was collected from the prostate tumour site, conditional to histopathological verification66,67. If not otherwise specified, all procedures were carried out at 4 °C. Tissue blocks were washed in Phosphate-buffered saline (PBS) solution for 2 minutes and minced for 2 minutes with a scalpel. Homogenised tissue was added to a solution (total volume of 7 ml) composed by of 1 mg/ml collagenase IV (Worthington Biochemical Corp, USA), 0.02 mg/ml DNase 1 (New England Biolabs, USA), 0.2 mg/ml dispase (Merck, USA). The tissue homogenised was serially digested at 37 °C at 180 rpm, through three steps of 5, 10 and 10 minutes of duration, with the final 3 minutes dedicated to sedimentation at 0 rpm. After each digestion step, the supernatant was aspirated and filtered through a 70 μm strainer into a pre-chilled tube, diluting the solution with 15 ml of 2% bovine serum PBS to quench the enzymatic reaction. The resulting cumulative solution was then centrifuged at 1500 rpm for five minutes, with the supernatant collected and the cell pellet resuspended into 1 ml 2% PBS-serum prior to labelling (Fig. S1).

Dataset EGAD00001004948

H3K27ac ChIP-seq in a selected group of AML patients

ChIP-seq data were generated for a number of selected patients to investigate changes in enhancer and promoter regions. ChIP was performed as described previously with slight modifications27. Briefly, cells were crosslinked with 1% formaldehyde for 10 minutes at room temperature and the reaction was quenched with glycine at a final concentration of 0.125 M. Chromatin was sheared using the Covaris S220 focused-ultrasonicator to an average size of 250–350 bp. A total of 2.5 µg of antibody against H3K27ac (Abcam, ab4729) was added to sonicated chromatin of 2 × 106 cells and incubated overnight at 4 °C. Protein A sepharose beads (GE healthcare) were added to the ChIP reactions and incubated for 2 h at 4 °C. Beads were washed and chromatin was eluted. After crosslink reversal, RNase A and proteinase K treatment, DNA was extracted with the Monarch PCR & DNA Cleanup kit (NEB). Sequencing libraries were prepared with the NEBNext Ultra II DNA Library Prep Kit for Illumina (NEB) according to the manufacturer’s instructions. The quality of dsDNA libraries was analyzed using the High Sensitivity D1000 ScreenTape Kit (Agilent) and concentrations were assessed with the Qubit dsDNA HS Kit (Thermo Fisher Scientific). Libraries were single-end sequenced on a HiSeq 4000 (Illumina). ChIP-seq reads were aligned to the human reference genome build hg19 with bowtie

Dataset EGAD00001007582

PDX_models_from_Latin_America_WES_

Acral melanoma (AM) is a subtype of cutaneous melanoma that occurs on glabrous skin of the hands and feet (palms, soles and nail apparatus) and which is relatively more common in non-European descent populations such as those of Latin America. Compared to other subtypes, it is poorly studied and experimental models that faithfully reproduce human AM features are urgently needed. We report the generation and preliminary characterisation of a collection of patient-derived xenografts of acral melanomas (AM-PDX) from Latin-American patients. PDX came from Brazilian and Mexican patients cared for at the Instituto Nacional de Câncer (INCA) and the Instituto Nacional de Cancerologia (INCan). The majority of the samples are acral lentiginous melanomas, ulcerated, with > 4mm Breslow depth and > 2 mitosis / mm2. Samples were mechanically fragmented and subcutaneously implanted in NSG mice and then collected for molecular analysis.

Study EGAS00001005663

PDX_models_from_Latin_America_RNAseq

Acral melanoma (AM) is a subtype of cutaneous melanoma that occurs on glabrous skin of the hands and feet (palms, soles and nail apparatus) and which is relatively more common in non-European descent populations such as those of Latin America. Compared to other subtypes, it is poorly studied and experimental models that faithfully reproduce human AM features are urgently needed. We report the generation and preliminary characterisation of a collection of patient-derived xenografts of acral melanomas (AM-PDX) from Latin-American patients. PDX came from Brazilian and Mexican patients cared for at the Instituto Nacional de Câncer (INCA) and the Instituto Nacional de Cancerologia (INCan). The majority of the samples are acral lentiginous melanomas, ulcerated, with > 4mm Breslow depth and > 2 mitosis / mm2. Samples were mechanically fragmented and subcutaneously implanted in NSG mice and then collected for molecular analysis.

Study EGAS00001008150

PDX_models_from_Latin_America_Xenofiltered_WES

Acral melanoma (AM) is a subtype of cutaneous melanoma that occurs on glabrous skin of the hands and feet (palms, soles and nail apparatus) and which is relatively more common in non-European descent populations such as those of Latin America. Compared to other subtypes, it is poorly studied and experimental models that faithfully reproduce human AM features are urgently needed. We report the generation and preliminary characterisation of a collection of patient-derived xenografts of acral melanomas (AM-PDX) from Latin-American patients. PDX came from Brazilian and Mexican patients cared for at the Instituto Nacional de Câncer (INCA) and the Instituto Nacional de Cancerologia (INCan). The majority of the samples are acral lentiginous melanomas, ulcerated, with > 4mm Breslow depth and > 2 mitosis / mm2. Samples were mechanically fragmented and subcutaneously implanted in NSG mice and then collected for molecular analysis.

Study EGAS00001008231

PDX_models_from_Latin_America_RNAseq_Xenofiltered

Acral melanoma (AM) is a subtype of cutaneous melanoma that occurs on glabrous skin of the hands and feet (palms, soles and nail apparatus) and which is relatively more common in non-European descent populations such as those of Latin America. Compared to other subtypes, it is poorly studied and experimental models that faithfully reproduce human AM features are urgently needed. We report the generation and preliminary characterisation of a collection of patient-derived xenografts of acral melanomas (AM-PDX) from Latin-American patients. PDX came from Brazilian and Mexican patients cared for at the Instituto Nacional de Câncer (INCA) and the Instituto Nacional de Cancerologia (INCan). The majority of the samples are acral lentiginous melanomas, ulcerated, with > 4mm Breslow depth and > 2 mitosis / mm2. Samples were mechanically fragmented and subcutaneously implanted in NSG mice and then collected for molecular analysis.

Study EGAS00001008232

Transcriptome Changes in ASD (NRXN1α+/-) iPSC-derived neurons

NRXN1 deletions are identified as one of major rare risk factors for autism spectrum disorder (ASD) and other neurodevelopmental disorders. ASD has 30% co-morbidity with epilepsy and the latter is associated with excessive neuronal firing. NRXN1 encodes hundreds of presynaptic neuro-adhesion proteins categorized as NRXN1α/β/γ. Previous studies on cultured cells show that the short NRXN1β primarily exerts excitation effect, whereas the long NRXN1α which is more commonly deleted in patients involves in both excitation and inhibition. However, patient-derived models are essential for understanding functional consequences of NRXN1α deletions in human neurons Quantitative RNA sequencing was performed by BGI using 50bp single-end strategy with BGI-SEQ500 sequencing device on iPSC-derived cortical neurons from 6 control lines; (2VC1, 3VCX1, 3VC2, 4C3, 4CX1, NCRM-1) of 4 healthy donors, and 4 NRXN1α+/- lines; (ND1C1, ND2C11, ND2CX1, ND4-1C2) from 3 ASD patients.

Study EGAS00001005536

NanoString nCounter® PanCancer IO 360™on anti-PD1/anti-PD1+CTLA4 in patients with metastatic melanoma

Patients with advanced melanoma were treated with standard-of-care single agent anti-PD-1 (nivolumab or pembrolizumab) or a combined anti-PD-1+anti-CTLA-4 (ipilimumab) therapy. Patients were retrospectively identified, based on formalin-fixed paraffin-embedded (FFPE) tissue availability. Total RNA was isolated from macro-dissected FFPE tissue sections using the AllPrep DNA/RNA FFPE Kit (Qiagen) or High Pure FFPET RNA Isolation Kit (Roche) according to the manufacturer’s instructions. RNA quantity was assessed on Qubit, and RNA integrity was assessed using the Tapestation system (Agilent). Total RNA samples (60ng/ul, total 200ng) were used as input for the NanoString PanCancer IO360 Panel, run on the nCounter MAX/FLEX prep station and scanner. Samples were hybridised for 20hrs, and each cartridge contained a panel standard.

Study EGAS00001006977

Metabolic profiling of patient-derived organoids reveals nucleotide synthesis as a metabolic vulnerability in malignant rhabdoid tumors

Malignant rhabdoid tumor (MRT) is one of the most aggressive childhood cancers for which no effective treatment options are available. Reprogramming of cellular metabolism is an important hallmark of cancer, with various metabolism-based drugs being approved cancer treatment. In this study, we use patient derived tumor organoids (tumoroids) to map the metabolic landscape of several pediatric cancers. Combining gene expression analyses and metabolite profiling using mass spectrometry, we find nucleotide biosynthesis to be a particular vulnerability of MRT. Treatment of MRT tumoroids with de novo nucleotide synthesis inhibitors lowers nucleotide levels in MRT tumoroids and induces apoptosis. Lastly, we demonstrate invivo efficacy of methotrexate in an MRT PDX mouse model. Our study revealsnucleotide biosynthesis as an MRT-specific metabolic vulnerability, which can ultimately lead to better treatment options for children suffering from this lethal pediatric malignancy.

Study EGAS00001007877

Identification of the underlying causal variant in a multi-generational family with autosomal dominant common variable immunodeficiency

Common variable immunodeficiency (CVID) is the most common form of primary immunodeficiency with an estimated incidence of 1:10,000. It has been apparent for many years that CVID has a genetic component, occurs frequently in families and can have both a recessive or dominant mode of inheritance. In recent years, 4 genes underlying CVID have been identified; however, mutations within in them are estimated to account for no more than 10% of all cases of CVID. We have identified a multi-generational family with autosomal dominant CVID. Genome-wide linkage analysis has mapped the locus underlying CVID in this family to an approximately 9.2 Mb interval on chromosome 3q27.3-q29, between the markers D3S3570 and D3S1265. This locus is distinct from any of the previously mapped susceptibility loci suggesting a novel genetic variant is responsible for disease in this family. The aim of this study is to use exome sequencing of affected (n = 4) and unaffected (n = 4) individuals, in tandem with the available genetic mapping data, to identify the causal variant underlying CVID in this family.

Dataset EGAD00001000363

ScRNA-seq of human kidney immune cells of patients with ANCA-associated glomerulonephritis, Lupus Nephritis against a healthy nephrectomy control

ANCA-associated glomerulonephritis (AGN) associates with a high risk of end-stage kidneydisease. The role of kidney immune cells in local inflammation remains unclear. This study investigates kidney immune cell diversity and function. Kidney tissue from AGN patients (n=5) and a lupus nephritis (LN) patient (n=1) were aquired during a biopsy procedure for a clinical indication. Needle-core biopsies were obtained for histopathological examination, and an additional pass was performed to retrieve kidney tissue for scRNA-seq. Healthy kidney tissue (n=1) was obtained from a kidney that was surgically removed do tue due to a (non-invasive) papillary urothelial carcinoma. Immediately after collection, kidney tissue was processed into a single-cell suspension and sorted using a 4-color flow cytometry panel to isolate living, CD45+immune cells. To aid in the multi-omic characterization, surface markers and T and B cell repertoires were sequenced in 2 samples (1 AGN patient and the nephrectomy control). These samples were incubated with an oligo-antibody TotalSeq-C cocktail containing 130 unique cell surface antigens.

Dataset EGAD50000000229

Psoriasis Patient PBMC scRNA-seq data

This dataset contains raw and processed data of four scRNA-seq PBMC samples from four Psoriasis patients (Pso3 = male, Pso4 = female, Pso7 = female, Pso8 = male). Raw sample libraries were prepared with 10x 5'-kit v2 multiplexing the individual samples using TrueSeq barcoded antibodies to stain the cells. Resulting sequencing data consists of four FASTQs (2 x antibody multiplex to demultiplex the samples + 2 x transcriptome data containing the typical cell x UMI barcode data for count matrix generation). The raw sequencing data was processed with cellranger 6.0.1 with the GRCh38 3.0.0 reference and cell to sample assignment was subsequently done with couhto 1.1.0 (https://github.com/dmalzl/counhto) employing the same algorithm as cellranger multi (couldn't be used here though as multi does not handle 5'-chemistry data). Processed data contains all main outputs of cellranger count in raw text format or in h5 formats including the used feature reference for (i.e. antibody barcode to sample assignment) as well as the computed cell to sample assignment computed by counhto. For further details including all used code in the analysis of the data please consult https://doi.org/10.5281/zenodo.13846873

Dataset EGAD50000001102

Single-nuclei histone modification profiling of the adult human central nervous system unveils epigenetic memory of developmental programs

We profile the adult human CNS from distinct regions, sex and ages, for chromatin accessibility at the single-nuclei level. Tissue was collected from 20 different donors (10 male, 10 female) within the ages of 34 to 74 years old. Additionally we performed Single-cell nanoCUT&Tag H3K27ac and H3K27me3 profiling for 4 of the donors and multiome single cell assay (RNA-seq+ATAC-seq) for 2 of the donors. Each donor donated fresh frozen white matter from the following three tissue regions: primary motor cortex (Brodmann area 4, BA4), arbor vitae cerebelli (CB) and fasciculi cuneatus and gracilis from and cervical spinal cord (CSC). Tissue was processed semi-randomly, ensuring that each batch of experiments had a representation of both sexes and all three tissue regions. Completed libraries were again randomly multiplexed during sequencing to minimize batch effects. We also performed micro-C,2 replicates, of iPS human derived cell cultures of human oligodendrocyte precursor cells (hOPC).

Dataset EGAD50000000410

CAGE analysis for endometrial carcinoma

For endometrial cancer patients, lymphadenectomy is recommended to exclude rarely metastasized cancer cells. This procedure is performed even in patients with low risk of recurrence despite the risk of complications such as lymphedema. A method to accurately identify cases with no lymph node metastases (LN-) before lymphadenectomy is therefore highly required. We approached this clinical problem by examining primary lesions of endometrial cancers with CAGE (Cap Analysis Gene Expression), which quantifies promoter-level expression across the genome. Fourteen profiles delineated distinct transcriptional networks between LN?+?and LN- cases, within those classified as having the low or intermediate risk of recurrence.

Study JGAS000124

PD-associated regulatory variants in human dopaminergic neurons reveals modulators of SCARB2 and BAG3 expression

SNPs have been associated with the risk of Parkinson’s disease but not in regulatory regions. While those potentially disrupt TFBS altering gene expression and contributing to cis-regulatory variation. To study changes in A and B compartment interactions and in topologically associated domains (TADs), we sequenced: - 3D chromatin contact using LowC, 6 samples. 3 replicates for each condition: smNPC (controls) cells and differentiated neurons after 30 days. Resulting paired-end are provided as BAM files. In order to confirm the predicted impact of the PD-associated allele of rs144814361 on BAG3 promoter, we proceeded with genome editing of the TH-REP1 cell line by using prime editing to insert the “T” allele at the position chr10:119651405 in the BAG3 promoter. - Chromatin accessibility using ATAC-seq, 9 samples of derived cell line TH-REP1. 3 replicates for each condition: 1) iPSC Wild Type 2) SNP-BAG3 variant in iPSC 3) SNP-BAG3 variant in smNPC. BAG3 variant: heterozygous for the SNP variant: rs144814361 (chr10 119651405) Resulting sequencing are provided as paired-end FASTQ files.

Dataset EGAD50000002258

A WTCCC2 genome-wide association study for psychosis endophenotype (PE) in individuals from UK, Germany, Holland, Spain and Australia.

A WTCCC2 project genome-wide association study for psychosis endophenotype (PE) in 5831 individuals from UK, Germany, Holland, Spain and Australia, genotyped on the Affymetrix 6.0 array. Information about the family structure of the samples can be found in the .info file which accompanies the data. It should be noted that a number of samples in this study have no phenotype information including sex, case status and family membership. Details of the WTCCC2 analysis can be found in Bramon et al. [Biol Psychiatry. 2014 Mar 1;75(5):386-97].

Study EGAS00001000817

Benchmarking of ProSolo, a new probabilistic single nucleotide variant caller for single cell DNA sequencing data. This study provides the whole exome sequencing dataset used in this assessment.

Benchmarking of ProSolo, a new probabilistic single nucleotide variant caller for single cell DNA sequencing data that provides control over the false discovery rate of different single cell events at genomic sites (e.g. alternative allele presence or allele dropout). Here, we provide the newly generated whole exome sequencing data used in assessing ProSolo's performance. It contains whole exome sequencing of: a patient with a constitutional mismatch repair defect (MSH6-) and their parents and siblings, a granulocyte bulk sample and 5 single granulocytes sequenced after whole genome amplification.

Study EGAS00001004123

Clinical utility of combined low-pass whole genome and targeted sequencing in liquid biopsies for diagnosis and monitoring of pediatric solid tumors

We designed a liquid biopsy (LB) platform employing Low-Pass Whole Genome Sequencing (LP-WGS) and targeted sequencing of cell-free (cf) DNA isolated from plasma for detecting genome-wide copy number alterations and targeted gene fusions in pediatric solid tumors. A total of 143 plasma samples were analyzed from 19 controls (patients with non-oncologic diagnoses) and 73 patients, including 44 with bone or soft tissue sarcomas, 12 with renal tumors, 10 with germ cell tumors, 5 with hepatic tumors, and 2 with thyroid tumors.

Study EGAS00001006913

Oesophageal adenocarcinoma WGS from LUD2015-005 study (NCT02735239, EudraCT 2015-005298-19)

This dataset contains whole genome sequencing (WGS) generated from the inoperable cohorts of the LUD2015-005 study (NCT02735239, EudraCT 2015-005298-19). WGS data were generated with the aim to assess previously reported (e.g. tumour mutational burden) and novel predictive genomic markers for immunochemotherapy treatment in this setting. Using these data, mutation and copy number profiles were generated for the LUD2015-005 study, which were assessed for correlation with patient outcomes from this trial.

Dataset EGAD00001009400

Whole genome sequencing

Whole genome sequencing of high-grade serous ovarian cancer (HGSC) tumours and matched normals from 15 patients with homologous recombination deficiencies. The dataset includes fastq files from 56 HGSC tumours (1 primary, 1 relapse, 54 end-stage) and 15 matched normals. Sequence libraries were generated from tumour and matched normal genomic DNA using the KAPA HyperPrep PCR-free library preparation kit (Roche), or the Illumina TruSeq DNA Nano kit according to manufacturer’s instructions. Sequencing was carried out by the Kinghorn Centre for Clinical Genomics Sequencing Laboratory (Sydney, Australia) on the HiSeq X Ten System (Illumina) or by the Australian Genome Research Facility (Melbourne, Australia) on an Illumina NovaSeq to a minimum base coverage of 30-fold for normal DNA and 60-fold for tumour DNA samples.

Dataset EGAD00001009746

Single-cell genome-wide concurrent haplotyping and copy-number profiling through genotyping-by-sequencing

scGBS is a single-cell sequencing-based methodology to haplotype and copy-number profile single cells. Genomic size and complexity is reduced through restriction enzyme digestion and DNA is genotyped through sequencing of the restriction fragments. scGBS data serves as the input for haplarithmisis, an algorithm we previously developed for SNP array-based single-cell haplotyping (Zamani Esteki et al., 2015). We established technical parameters and developed an analysis pipeline enabling accurate concurrent haplotyping and copy-number profiling of single cells with the use of a HapMap cell line pedigree (7 single cells). A clinical validation of the methodology with a total of 14 single blastomeres and 3 trophectoderm samples biopsies from human preimplantation embryos for 6 PGT-M families were processed with scGBS and were previously haploptyped via SNP array.

Dataset EGAD00001007794

Single cell transcriptional evolution of myeloid leukaemia of Down syndrome – scRNA

Abstract: Children with Down syndrome have a 150-fold increased risk of developing myeloid leukaemia (ML-DS). Unusually for a childhood leukaemia, ML-DS arises from a preleukaemic state, termed transient abnormal myelopoiesis (TAM), and via a conserved sequence of mutations. Here, we examined the relationship between the genetic and transcriptional evolution of ML-DS from a rich collection of primary patient samples through single cell mRNA sequencing, complemented by phylogenetic analyses in progressive disease. We distilled the transcriptional consequence of each genetic step in the evolution of ML-DS, showing that TAM-defining GATA1 mutations account for most of the ML-DS transcriptome, including those of progressive disease. This transcriptional backbone may thus represent a common vulnerability in ML-DS. We extracted the transcriptional difference between TAM and ML-DS which, unexpectedly, reflected features shared across childhood leukaemia. Our approach delineates the transcriptional evolution of ML-DS and provides an analytical blueprint for distilling consequences of mutations directly from patient samples.

Dataset EGAD00001015452

POETIC clinical Trial Ki67 Bookend Breast Cancer Targeted Exome and RNAseq Project

Using RNAseq, we compare the 15% of poorest responders (PRs, n=177) as measured by proportional Ki67 changes after 2 weeks of neoadjuvant aromatase inhibitors to good responders (GRs, n=190) selected from the top 50% responders in the POETIC trial and matched for baseline Ki67 categories. In the POETIC trial, 4,480 postmenopausal women with primary ER+ BC were randomised 2:1 to receive either treatment with a non-steroidal AI (letrozole or anastrozole) for 2 weeks before and 2 weeks after surgery or to no perisurgical treatment. Only AI-treated patients with HER2- tumors, paired baseline and surgery Ki67 available, and baseline Ki67 immunohistochemistry (IHC) >10% (to minimise imprecision in proportional Ki67 falls) were included for selection. Data is baseline RNAseq and targeted exome DNA sequencing analysis of POETIC Good/Poor Responders to aromatase inhibitors based on change in Ki67.

Dataset EGAD00001010919

SINGLE CELL ANALYSIS OF IN VITRO ERYTHROPOIESIS

We are developing a protocol to differentiate mouse and human induced pluripotent stem (IPS) and embryonic stem (ES) cells towards the haematopoietic pathway to generate erythrocytes in vitro. This system has many applications such as the study of the role of specific genes and human polymorphisms in infectious diseases such as malaria, as well as haematological diseases such as myelodysplastic syndrome. The nature of the in vitro differentiation process means that a heterogeneous population of cells is generated. In order to understand the types of cells produced with our protocol, we have performed a single cell analysis, which has the power to reveal the different populations of cells and their characteristics. For this, a cDNA library has been made that needs to be sequenced to obtain the gene expression profiles of the different cells. With this information we will be able to assess the quality of the differentiation protocol and improve it in order to produce better cells for the downstream applications.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 EGAD00001000979

Structural rearrangements generate cell-specific, gene-independent CRISPR-Cas9 loss of fitness effects.

CRISPR-Cas9 genome editing is widely used to study gene function, from basic biology to biomedical research. Structural rearrangements are a ubiquitous feature of cancer cells and their impact on the functional consequences of CRISPR-Cas9 gene-editing has not yet been assessed. Utilizing CRISPR-Cas9 knockout screens for 250 cancer cell lines, we demonstrate that targeting structurally rearranged regions, in particular tandem or interspersed amplifications, is highly detrimental to cellular fitness in a gene independent manner. In contrast, amplifications caused by whole chromosomal duplications have little to no impact on fitness. This effect is cell line specific and dependent on the ploidy status. We devise a copy-number ratio metric that substantially improves the detection of gene-independent cell fitness effects in CRISPR-Cas9 screens. Furthermore, we develop a computational tool, called Crispy, to account for these effects on a single sample basis and provide corrected gene fitness effects. Our analysis demonstrates the importance of structural rearrangements in mediating the effect of CRISPR-Cas9-induced DNA damage, with implications for the use of CRISPR-Cas9 gene-editing in cancer cells.

Dataset EGAD00001004124

Autozygosity pilot - Born in Bradford (2014-11-20)

The offspring of first cousin marriages have ~6% of their genome autozygous, i.e. homozygous identical by descent, or even more if there was further consanguinity in their ancestry. In the UK there are large populations with very high first cousin marriage rates of 20-50%. Sequencing the exomes of a sample of these individuals has the potential both to support genetic health programmes in these populations, and to provide genetic research information about rare loss of function mutations. This pilot study based on existing cohort samples from the Born In Bradford study will identify homozygous individuals for almost all variants down to an allele frequency around 1%, plus individuals carrying hundreds of new homozygous rare loss-of-function variants, and will support development of community relations and ethics for a wider study currently being designed. The data deposited in the EGA consist of low coverage whole exome sequencing on these samples.Data Access is controlled by the Wellcome Trust Sanger Institute DAC and the Born In Bradford Executive Group. This dataset contains all the data available for this study on 2014-11-20.

Dataset EGAD00001001079

DAC for Combined MEK and BCL-2/XL inhibition is effective in high-grade serous ovarian cancer patient-derived xenograft models and BIM levels are predictive of responsiveness

Dac EGAC00001001250

cqmuGWAS1

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

Dataset EGAD00010001527

PMBCL IL4R DASL

Expression data for 42 PMBCL patient samples (32 IL4R WT cases and 10 cases with mutations in IL4R)

Dataset EGAD00010001542

Copy number profiling of primary samples and cell lines of retinoblastoma

Illumina SNP-arrays for matching retinoblastoma-blood pairs and retinoblastoma cell lines.

Dataset EGAD00010000908

Serum Proteomics of Aortic Diseases

Comparative proteome-based analysis of different autologous bone entities used for alveolar onlay grafting

Dataset EGAD00010002314

SNParray_PGT_samples_scGBS

PGT samples for haplotyping and copy-number profiling via SNP array

Dataset EGAD00010002169

SNParray_HapMap_samples_scGBS

HapMap samples for haplotyping and copy-number profiling via SNP array

Dataset EGAD00010002168

scRaCH-seq data targeting BTK and SF3B1

Data set for scRaCH-seq manuscript. We base call the fast5 using guppy.

Dataset EGAD50000000235

HMO-microbiome study dataset, 16S sequencing

This dataset contains 16S sequencing of fecal and milk samples for HMO-microbiome Lifelines NEXT study

Dataset EGAD50000000532

Cancer Genomics, ICR, summary data

DAC for summary level data from the Cancer Genomics group, Division of Genetics and Epidemiology at the Institute of Cancer Research, Sutton, UK.

Dac EGAC50000000050

Foundation Medicine Binary Calls

Binary calls for all available samples in 4 Atezo trials.

Dataset EGAD50000000709

Picuris Pueblo Genomic Project (Modern Data)

Data Access Committee for the present-day dataset of the project "Picuris Pueblo oral history and genomics reveal continuity in US Southwest".

Dac EGAC50000000526

DNABR

Genome variation data for 2723 Brazilian individuals from the DNA do Brasil Project.

Dataset EGAD50000001012

DAC for DRS peripheral blood, University Medical Center Johannes Gutenberg University Mainz

This DAC is in relation to the DRS samples of peripheral blood of healthy control

Dac EGAC50000000667

University Clinic Golnik DAC

DAC of the University Clinic Golnik, that will decide regarding the requests for data submitted on behalf of research groups of the Clinic

Dac EGAC50000000820

WGS_TALL_t_14_16_translocation

WGS data for patients with t(14;16) translocated T-ALL.

Dataset EGAD50000001782

NANOPORE_TALL_t_14_16_translocation

Nanopore sequencing for patients with t(14;16) translocated T-ALL

Dataset EGAD50000001785

SF10432 Wildtype Primary GBM Female, 50

Single Cell-RNA Seq of Wildtype Primary GBM for Female, Age 50.

Dataset EGAD00001006019

Fecal WMS HV metadata

Metadata for fecal WMS data from NCT02749630 healthy volunteers.

Dataset EGAD00001008819

RNAseq data for EGAS00001004572

Bulk RNA-seq data of tumours in EGAS00001004572.

Dataset EGAD00001006876

Long-term organoid culture of a small intestinal neuroendocrine tumor

tissue, organoid and normal bam files for whole exome samples

Dataset EGAD00001010162

Direct RNA sequencing of 10 postmortem human brain samples

Contains fast5 data for each of the 10 samples sequenced.

Dataset EGAD00001015347

RNA sequencing BRAF fusion partners

Contains RNAseq data for 14 transduced/non-transduced organoids

Dataset EGAD00001005380

Pulldown DNA methylation study v2

Pilot study to set up sequencing protocols for targeted pulldown methylation profiling

Dataset EGAD00001001242

Reference epigenome IPS01_N_Fibroblast_WGBS data generated from KEP study

A IPS01_N_Fibroblast_WGBS paired end data for iPSC(Oct4)

Dataset EGAD00001003473

Reference epigenome IPS04_X_Fibroblast_WGBS data generated from KEP study

A IPS04_X_Fibroblast_WGBS paired end data for iPSC(Oct4)

Dataset EGAD00001003476

MPB_Bonn

1000Genomes imputed data set of 581 cases and 417 controls for male-pattern baldness

Dataset EGAD00001001456

Bulk RNA-Seq Digital Gene Expression Matrix

Raw read counts and normalized read counts for each gene

Dataset EGAD50000002489

CSER: South-Seq: DNA Sequencing for Newborn Nurseries in the South

SouthSeq, part of the Clinical Sequencing Evidence-Generating Research (CSER2) effort, aims to perform genome sequencing for infants in nurseries at hospitals in Alabama, Mississippi, Louisiana, and Kentucky. SouthSeq aims to sequence 600 newborn babies presenting with congenital anomalies, or other signs of a rare genetic disorder. The primary goals of SouthSeq are to provide early diagnoses to children with rare genetic conditions, to demonstrate the utility of genome sequencing as a frontline test with potential to improve intermediate and long-term clinical outcomes, and to equip non-genetics providers to deliver test results in order to facilitate the conduction of genome sequencing outside of major medical centers where patients may not have access to a genetics provider. SouthSeq also aims to ensure that genomic testing is applied across diverse communities; recruitment efforts focus on patients from rural and minority populations that have been historically under-represented in both medicine and genomics research. Sequencing and analysis for SouthSeq are being conducted at the HudsonAlpha Institute for Biotechnology in Huntsville, Alabama (https://www.hudsonalpha.org/). Enrollment of infants with signs of genetic disease is approved at five clinical sites, including nurseries at the University of Alabama at Birmingham, the University of Mississippi Medical Center, Woman's Hospital in Baton Rouge, Children's Hospital in New Orleans, and the University of Louisville. Both biological parents are enrolled when available so that follow-up testing can be performed to determine inheritance. Further, participant families may choose to opt-into return of secondary findings identified in the proband, which focus on pathogenic and likely pathogenic variation identified within ACMG SFv3.0 genes.

Study phs002307

Rare Diseases Clinical Research Network (RDCRN) Vasculitis Clinical Research Consortium (VCRC) Longitudinal Protocol for Takayasu's Arteritis

The development and validation of accurate biological markers and predictors of disease activity and outcome for Takayasu's arteritis, a form of idiopathic vasculitis, would have a major positive impact on the clinical care of patients with this rare disease, be an important advance in the design of clinical trials and feasibility of new drug development, and provide important insight into the pathogenesis of this condition. Current measures of disease activity, including erythrocyte sedimentation rate (ESR) and C-reactive protein (CRP), while helpful in selected settings, are mostly considered highly flawed and markers whose poor correlation with disease flares makes them inadequate to help guide therapy. Similarly, ESR and CRP are clearly unhelpful as diagnostic tests for vasculitis. The current state of clinical investigation for vasculitis relies heavily upon determination of disease status solely by clinical presentation and investigator opinion. Thus, discovery and validation of more precise markers that quantifiably measure activity, predict disease course, and correlate with an important biological process, would be of great importance. Knowledge to be gained from this study could potentially be highly important. Discovery of effective biomarkers of vasculitis and assessment of long term disease course and complications by the collection of clinical and diagnostic imaging studies proposed could lead to better care and less toxic treatment regimens. Furthermore, insight into the pathophysiology of vasculitis could be gained and this might lead to better treatments and improved outcomes in terms of reduction of vascular and treatment related complications. Ideally, excellent biomarkers reflect a sophisticated understanding of the pathophysiology of a disease.

Study phs000589

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

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

PacBio long-read scRNA-seq

Full-length single-cell cDNA PacBio Kinnex whole-transcriptome sequencing data for adrenal sampleA and sampleB.

Dataset EGAD50000002211

12028 results for "reddit fc 26 coins Visit Buyfc26coins.com for latest FC 26 coins news..FeUn"

in 53.70 milliseconds.