14976 results for "cancer studies using rna-seq"

in 26.37 milliseconds.

• Comparison of sequencing assays for sensitive detection of circulating tumour DNA in stage IA-IV breast cancer

  This study compared different assays for the detection of circulating tumour DNA (ctDNA) in serial plasma from stage IA-IV breast cancer patients, targeting structural variants (SVs), single nucleotide variants (SNVs) and/or somatic copy-number aberrations (SCNAs). SV-multiplex PCR, SNV-/SV-hybrid capture, and different depths of whole-genome sequencing (WGS) were used to evaluate ctDNA levels, demonstrating concordant results. SNV-hybrid capture targeting 1,347-7,491 mutations was the most sensitive assay, detecting 67% (36/54) of samples down to an allele fraction (AF) of 0.00024%. SV-multiplex PCR, targeting 21-47 mutations, detected 63% (34/54) of samples down to 0.00047% AF and has potential as a clinical assay.

  Dataset EGAD00001008589

• Whole genome and Whole exome sequencing of patient-derived xenograft models of endometrial cancer

  Whole genome and Whole exome sequencing of patient-derived xenograft models of endometrial cancer

  Dataset EGAD00001008326

• Chemotherapy induces canalization of cell state in childhood B-cell precursor acute lymphoblastic leukemia

  Comparison of intratumor genetic heterogeneity in cancer at diagnosis and relapse suggests that chemotherapy induces bottleneck selection of subclonal genotypes. However, evolutionary events subsequent to chemotherapy could also explain changes in clonal dominance seen at relapse. We therefore investigated the mechanisms of selection in childhood B-cell precursor acute lymphoblastic leukemia (BCP-ALL) during induction chemotherapy where maximal cytoreduction occurs. To distinguish stochastic versus deterministic events, individual leukemias were transplanted into multiple xenografts and chemotherapy administered. Analyses of the immediate post-treatment leukemic residuum at single-cell resolution revealed that chemotherapy has little impact on genetic heterogeneity. Rather, it acts on extensive, previously unappreciated, transcriptional and epigenetic heterogeneity in BCP-ALL, dramatically reducing the spectrum of cell states represented, leaving a genetically polyclonal but phenotypically uniform population, with hallmark signatures relating to developmental stage, cell cycle and metabolism. Hence, canalization of the cell state accounts for a significant component of bottleneck selection during induction chemotherapy.

  Study EGAS00001004407

• VESPA: Vanderbilt Electronic Systems for Pharmacogenomic Assessment

  The Vanderbilt Electronic Systems for Pharmacogenomic Assessment (VESPA) Project is a large electronic medical record (EMR)- and biobank-based initiative for translational pharmacogenomics discoveries. Key research resources utilized in this effort include the BioVU DNA databank and associated Synthetic Derivative database of clinical information, and software tools developed to identify drugs and clinical events using EHR-derived structured and unstructured ("free text") data. Cohorts include subjects with primarily drug-response phenotypes, and most cases and controls identified include three data types: ICD-9 codes, medication regimens, and medical test results. Advanced informatics methods, such as natural language processing, were also used for some phenotypes. Algorithms used event-sequence analyses to establish temporal relationships between drugs and phenotypes. Both cases and control algorithms excluded records that contained specific comorbidities and were refined to achieve positive predictive value (PPV) > 90%. For automated algorithms failing to meet this threshold, manual review was coupled with algorithms to validate that the included cases were true positives. A total of 11,639 subjects met phenotyping criteria for a least one of the 28 phenotypes investigated. Across all phenotypes cases and controls, 90% were reused as either a case or control for a least one other phenotype. Data are deposited for individual phenotypes in the VESPA substudies. Substudies: VESPA: Genome Wide Association Study of Serum Creatinine during Vancomycin Therapy and Vancomycin Pharmacokinetics - phs000894 VESPA: Genome Wide Study of Clinically Diagnosed Ace Inhibitor Cough - phs000992

  Study phs000991

• Phase II Therapeutic Trial of Mexiletine in Non-Dystrophic Myotonia (IND #77,021)

  Non-dystrophic myotonias are a group of rare muscle disorders caused by abnormalities in different muscle cell membrane proteins. Patients experience delayed muscle relaxation that causes impaired physical activity, stiffness, and pain. Mexiletine has been shown to be beneficial for symptoms of myotonia in case reports of patients with non-dystrophic myotonia, however, current management of these patients primarily depends on individual practitioners' preference. This study is a placebo-controlled, double-blind, cross-over treatment trial using mexiletine in non-dystrophic myotonias. The results, if positive, will permit the identification of mexiletine as an effective therapy for myotonia, making it a first-line therapy for patients with non-dystrophic myotonia. The specific aim of this proposal is to perform a randomized, double-blind, placebo-controlled cross-over study to assess whether mexiletine improves both quantitative and qualitative measures of myotonia in patients with non-dystrophic myotonia. Given the rarity of non-dystrophic myotonia, large randomized controlled treatment trials in these disorders have not been feasible in the past. Primary Outcome Measure: Patient-assessed symptom: stiffness as measured by an Interactive Voice Response Diary (IVR). Secondary Outcome Measures: a) patient-assessed pain, weakness, and fatigue as measured by IVR; b) clinical myotonia assessment; c) quality of life as measured by INQoL, SF36; d) a quantitative measure of grip myotonia; e) measurement of CMAP after short and long exercise; f) grading of myotonia on needle EMG.

  Study phs001311

• The Pioneer 100 Wellness Project (P100)

  108 individuals (age 21-89+ years; 59% males, 41% females; 89% Caucasian; not recruited based on any specific phenotype) participated in an IRB-approved study from April 2014 to January 2015. Each individual had the genome sequenced in full. Blood was collected in clinics every three months. Additionally, participants completed at-home collections of saliva, stool, and first morning void urine every three months. Stool and saliva samples were shipped directly to the vendor by the participant, while urine was given back to the study coordinators for distribution to the proper sample vendor. For each successful participant in each round we carried out 218 clinical laboratory tests, measured up to 643 metabolites and 262 proteins, and measured the abundance of 4616 operational taxonomic units (OTUs) in the gut microbiome using 16S rRNA sequencing. The P100 study had four objectives: 1) Establish cost-efficient procedures for generating, storing and analyzing multiple sources of health data obtained over time from participants and analyzed in combination with genomic data; 2) Develop and deploy analytic tools for integrating these diverse datasets and deriving actionable information from their observed interrelationships; 3) Identify novel patterns within the streams of health data that indicate either wellness, or transitions between wellness and disease; 4) Learn how to best interface with and present longitudinal health information to individuals by studying the reactions and feedback from participants as they are presented actionable information.

  Study phs001363

• NIMH (National Institute of Mental Health) De Novo Mutation Identification in Taiwanese Schizophrenia Trios

  The substantial reproductive impact of schizophrenia, for which affected individuals have fewer than half as many offspring as unaffected individuals do, implies that mutations of largest effect will frequently be de novo mutations. Ascertaining exome sequence variation in father-mother-offspring trios allows such mutations to be identified and distinguished from the far-larger amount of rare variation that is inherited by each individual. The pursuit of this approach in a large, well-powered cohort of trios can also provide lessons that inform the development of such gene discovery strategies more generally in human genetics. Schizophrenia trios from the Taiwanese population are being collected by Dr. Ming Tsuang (PI, UC San Diego, California) and investigators in Taiwan (PI, Dr. Hai Gwo Hwu; both funded by NIMH grant 1R01MH085560; Expanding Rapid Ascertainment Networks of Schizophrenia Families in Taiwan). A total of 3800 trios are anticipated to be collected by May 2013. This represents a highly homogenous national sample from the same ancestral population. DNA samples will be obtained from the NIMH Repository, Rutgers University Cell and DNA Repository (described below) and stored at the Broad Institute. Genetic and data analyses will be performed at the Broad Institute. We propose to sequence the whole exome of trios by hybrid capture and Illumina next generation sequencing and perform targeted genotyping and validation of variants (SNPs, indels and CNVs) using several molecular methods, to include emulsion-based PCR and Sanger sequencing.

  Study phs001196

• Genome and Transcriptome Assembly Reveals SVA-Mediated Aberrant Splicing in X-Linked Dystonia Parkinsonism

  X-linked Dystonia-Parkinsonism (XDP, DYT3) is a long-standing quandary in human disease genetics. XDP is predominantly observed on Panay island in the Philippines. This study is one of the first of its kind to interrogate an unsolved Mendelian disorder by integrating genome and transcriptome assembly methods using Illumina, 10X Genomics, Pacific Biosciences, and Agilent genome targeting technologies. These data provide strong evidence for a pathogenic link between a noncoding SVA retrotransposon and XDP. We demonstrate that this Mendelian disorder is associated with a sine-VNTR-Alu (SVA) retrotransposon that inserted into the TAF1 gene and is shared by all XDP probands, yet never observed in controls from worldwide populations. Transcriptome assembly in iPSC-derived neural stem cells (NSCs) and neurons revealed that this SVA caused aberrant splicing and significant intron retention, which was negatively correlated with TAF1 expression. Remarkably, CRISPR/Cas9 excision of the SVA rescued the aberrant transcriptional signature and normalized expression of TAF1 in patient-derived NSCs. We have also interrogated iPSC-derived microglia to investigate the contribution of glial components to the cellular and molecular deficit, effect of TNF treatment on XDP transcriptional signatures. To further understand what XDP may have in common with other hereditary dystonias (such as DYT6), we are further comparing these signatures to ones associated with dystonia-specific variants in proteins that, like TAF1, are also involved in regulating transcription, i.e. THAP1. 

  Study phs001525

• The preterm infant microbiome: biological, behavioral and health outcomes at 2 and 4 years of age

  This project follows a cohort of 78 Very Low Birth Weight (VLBW) previously enrolled infants in a R21 grant plus additional 25 infants through their Neonatal Intensive Care Unit (NICU) stay until they reach the age of 4 years. The data, gathered over 6 weeks of the NICU stay, includes multiple factors, such as prenatal and postnatal events and illnesses, received human milk amount, weekly means of cytokines, chemokines, growth factors, and secretory Immunoglobulin A in the milk, and weekly levels of fecal calprotectin. These factors could potentially alter the gut microbiome. Microbiome species and diversities will be measured in the laboratory of Dr. Jack Gilbert at Argonne National Laboratory using state of the science deep sequencing and amplification of microbial sRNA genes. The microbiome will again be measured in stool samples from those children at the ages of 2 and 4 years. Relationship between the prenatal and postnatal factors, human milk volume and immunobiology, fecal calprotectin levels, and the very early microbiome will be analyzed. The predictive power of the VLBW infant gut microbiome for determining later childhood microbiomes will be analyzed prospectively. The relationships between microbiomes across time and later growth, development and health will be determined. VLBW infants are at risk for both early and later health effects, and the role of the microbiome in these effects will be measured in this prospective study.

  Study phs001578

• Indexing Genes Impacted by Copy Number Variation in Developmental Disorders

  A challenge in clinical genomics is to predict whether copy number variation (CNV) affecting a gene or multiple genes will manifest as disease. Increasing recognition of gene dosage effects in neurodevelopmental disorders prompted us to develop a computational approach based on critical-exon (highly expressed in brain, highly conserved) examination for potential etiologic effects. Using a large CNV dataset, our updated analyses revealed significant (P < 1.64 x 10-15) enrichment of critical-exons within rare CNVs in cases compared to controls. Separately, we used a weighted gene co-expression network analysis (WGCNA) to construct an unbiased protein module from prenatal and adult tissues and found it significantly enriched for critical exons in prenatal (P < 1.15 x 10-50, OR = 2.11) and adult (P < 6.03 x 10-18, OR = 1.55) development. WGCNA yielded 1,206 proteins for which we prioritized the corresponding genes as likely to have a role in neurodevelopmental disorders. We compared the gene lists obtained from critical-exon and WGCNA analysis and found 438 candidate genes associated with CNVs annotated as pathogenic, or variants of uncertain significance (VOUS), from among 10,619 developmental delay cases. We identified genes containing CNVs previously considered to be VOUS, to be new candidate genes for neurodevelopmental disorders (GIT1, MVB12B and PPP1R9A) demonstrating the utility of this strategy to index the clinical effects of CNVs. Reprinted from Uddin et al, 2016, with permission from Scientific Reports.

  Study phs001154

• Phenotype Risk Scores Identify Patients with Unrecognized Mendelian Disease Patterns

  This study was designed to assess the phenotypic effects of rare variants. Rare variants are difficult to study in a high-throughput manner because most cohorts are underpowered to detect associations. In order to gain power to test rare variants, we use Phenotype Risk Scores (PheRS) based on features of Mendelian diseases. PheRS is calculated using claims data from an EHR that is mapped to clinical features from OMIM's clinical descriptions of Mendelian disease. We calculated PheRS for 1,204 Mendelians diseases in a cohort of 21,701 individuals genotyped on the HumanExome BeadChip. We then tested for association between these diseases and rare variants in causal genes. A phenotype risk score is calculated as the weighted sum of features for a given disease. The features are defined as a set of phecodes (consolidated ICD codes) associated with a particular Mendelian disease as described in OMIM. The weights for each phecode are calculated as the log inverse prevalence of the phecode in our cohort. For an individual, their score equals the sum of the weights for each phecode that is present in their medical record. Our study was based on a cohort of individuals with genotype data linked to de-identified electronic health records (EHR) from Vanderbilt's BioVU resource. We hope that this method will help generate phenotype-genotype correlation data on rare variants which in turn may inform rare variant interpretations use for diagnostics.

  Study phs001516

• Expression of Activation Induced Cytidine Deaminase and Risk of Transformation in Follicular Lymphoma

  It is currently unclear how to identify follicular lymphoma (FL) patients with low disease burden but high risk for early progression. Building on a prior study demonstrating early transformation in FLs with high variant allele frequency (VAF) BCL2 mutations at activation-induced cytidine deaminase (AICDA) sites, we examined 11 AICDA mutational targets, including BCL2, BCL6, PAX5, PIM1, RHOH, SOCS1, and MYC, in a more contemporary cohort of 199 newly diagnosed grade 1 and grade 2 FLs using targeted capture and massively parallel sequencing. BCL2 mutations with VAF >20% occurred in 52% of cases. Among 97 FLs that did not initially receive rituximab-containing therapy, nonsynonymous BCL2 mutations at VAF >20% were associated with shorter, event-free survival (EFS, median 1.4 years with these mutations versus 4.5 years without, p = 0.048), and increased risk of transformation (HR 3.01, 95% CI 1.03-8.77, p = 0.04). Other sequenced genes were less frequently mutated and did not significantly increase the prognostic value of the panel. Across the entire population, BCL2 mutations at VAF >20% were associated with decreased EFS (HR 1.58, 95% CI 1.03-2.40, p = 0.034 after correction for FLIPI and treatment) and decreased overall survival after median 14-year follow-up (HR 1.86, 95% CI 1.07-3.23, p = 0.029). Thus, high VAF BCL2 mutations remain a prognostic factor for FL even in the chemoimmunotherapy era and might need to be considered in future study design.

  Study phs002845

• Diagnostic Genomic Analysis is Prognostic in AYA ALL Patients Treated on a MRD-Stratified Paediatric Protocol

  Background: Over the last decade, next generation sequencing has enabled classification of multiple new recurrent genomic drivers of acute lymphoblastic leukemia (ALL). We provide a comprehensive genomic survey of adolescent and young adult (AYA) patients treated on the Australasian Leukaemia and Lymphoma Group (ALLG) ALL06 study. ALL06 clinical trial described the efficacy and tolerability of AYA ALL patients aged 15-39 years treated on a single minimal residual disease (MRD)-stratified platform. In this study we have taken the ALL06 cohort and retrospectively performed comprehensive genomic analyses assessing response of AYA patients with various genomic subtypes treated with a paediatric approach. Methods: Genomic risk classification of B- and T-ALL patients (n=57) enrolled to the study was based on mRNA-Sequencing, MLPA, immunophenotyping and cytogenetics. Results: Using this approach, 36/40 (90%) B- and 13/17 (76.5%) T-ALL patients were classified according to genomic risk. We have linked outcome data with each patient’s genomic risk and found a strong correlation between adverse genomic risk and MRD at the end of consolidation, translating to inferior overall and relapse free survival. Patients with adverse risk genomics who achieved negative MRD status had improved responses compared to those with persistent MRD. Patients with standard risk genomics had excellent responses regardless of their MRD status. Conclusion: Our data argue strongly for incorporation of a genomic risk classification into future ALL treatment paradigms at the time of diagnosis.

  Study EGAS50000000752

• Ultra-long sequencing for contiguous haplotype resolution of the human immunoglobulin heavy chain locus

  Genetic diversity within the human immunoglobulin heavy chain (IGH) locus influences the expressed antibody repertoire and susceptibility to infectious and autoimmune diseases. However, repetitive sequences and complex structural variation pose significant challenges for large-scale characterization. Here, we introduce a method that combines Oxford Nanopore ultra-long sequencing and adaptive sampling with a bioinformatic pipeline to produce haplotype-resolved, annotated IGH assemblies. Notably, our strategy overcomes prior limitations in phasing resolution, enabling single-contig haplotype assemblies that span the entire IGH locus. We applied this method to four individuals and validated the accuracy of the IGH assemblies using Pacific Biosciences HiFi reads, demonstrating near-complete sequence congruence, detecting only indels with a high degree of confidence. Moreover, when applied to the reference material HG002, our pipeline revealed no base differences and a limited number of indels compared with the Telomere-to-Telomere genome benchmark across the IGH region. Importantly, in the four individuals, our approach uncovered 30 novel alleles and previously uncharacterized large structural variants, including a 120 kb duplication spanning IGHE to IGHA1 and an expanded seven-copy IGHV3-23 gene haplotype. These findings underscore the power of our method to resolve the full complexity of the IGH locus and uncover previously unrecognized variants that may affect immune function and disease susceptibility. Thus, our method provides a strong basis for future immunological research and translational applications.

  Study EGAS50000001042

• Alpha1-Antitrypsin Deficiency Registry (AADR-BioLINCC)

  Data Access NOTE: Please refer to the “Authorized Access” section below for information about how access to the data from this accession differs from many other dbGaP accessions. Objectives: To characterize the clinical and laboratory course of patients with severe alpha 1-antitrypsin deficiency whether or not the patient is undergoing long-term augmentation therapy.Background: A hereditary disorder, patients with low serum levels of alpha-1-antitrypsin are at an increased risk for the early onset of emphysema. The only approved treatment for alpha-1-antitrypsin deficiency is augmentation therapy using a purified preparation of human alpha-1-antitryspin. Sample sizes for a randomized controlled clinical trial of augmentation therapy were determined to be infeasible; therefore, a multi-center registry was initiated in 1988 to explore the natural history of the disease and the relative efficacy of augmentation therapy in patients with a severe deficiency of alpha-1-antitrypsin.Participants: Eligible participants included individuals 18 years of age or greater for whom the Central Laboratory confirmed that the serum alpha 1-antitrypsin level is Conclusions: Participants receiving augmentation therapy had decreased mortality risk during follow-up. FEV1 decline among all participants did not differ by augmentation therapy; however, among participants with FEV1 35-49% predicted, FEV1 decline was significantly slower for participants on augmentation therapy than for those not receiving therapy. (Am J Respir Crit Care Med, 1998; 158:49-59)

  Study phs004187

• DNA methylation at HBV integrants and flanking host genomes

  Integration of DNA viruses into the human genome plays an important role in various types of tumors, including hepatitis B virus (HBV)-related hepatocellular carcinoma. However, the molecular details and clinical impact of HBV integration on either the human or HBV epigenomes are unknown. Here, we show that methylation of the integrated HBV DNA is related to the methylation status of the flanking human genome. We developed a next-generation sequencing-based method for structural methylation analysis of integrated viral genomes (denoted G-NaVI). This method is a novel approach that enables enrichment of viral fragments for sequencing using unique baits based on the sequence of the HBV genome. We detected integrated HBV sequences in the genome of the PLC/PRF/5 cell line and found variable levels of methylation within the integrated HBV genomes. Allele-specific methylation analysis revealed that the HBV genome often became significantly methylated when integrated into highly methylated host sites. After integration into unmethylated human genome regions such as promoters, however, the HBV DNA remains unmethylated and may eventually play an important role in tumorigenesis. The observed dynamic changes in DNA methylation of the host and viral genomes may functionally affect the biological behavior of HBV. These findings may impact public health given that millions of people worldwide are carriers of HBV. We also believe our assay will be a powerful tool to increase our understanding of the various types of DNA virus-associated tumorigenesis.

  Study JGAS000015

• Whole genome study of Hurthle cell thyroid carcinoma

  Oncocytic thyroid carcinoma, also known as Hürthle cell thyroid carcinoma, accounts for only a small percentage of all thyroid cancers. However, this malignancy often presents at an advanced stage and poses unique challenges to patients and clinicians. Surgical resection of the tumor accompanied in some cases by radioactive iodine treatment, radiation and chemotherapy are the established modes of therapy. Knowledge of the perturbed oncogenic pathways can provide better understanding of the mechanism of disease and thus opportunities for more effective clinical management. Initially, two oncocytic thyroid carcinomas and their matched normal tissues were profiled using whole genome sequencing. Subsequently, 72 oncocytic thyroid carcinomas, one cell line and 5 Hürthle cell adenomas were examined by targeted sequencing for the presence of mutations in multiple endocrine neoplasia I (MEN1) gene. We report the identification of MEN1 loss-of-function mutations in approximately 7% of patients diagnosed with oncocytic thyroid carcinoma. Whole genome sequence data also revealed large regions of copy number variation encompassing nearly the entire genomes of these tumours. Menin, a ubiquitously expressed nuclear protein, is a well-characterized tumor suppressor whose loss is the cause of multiple endocrine neoplasia type 1 syndrome. Menin is involved in several major cellular pathways such as regulation of transcription, control of cell cycle, apoptosis and DNA damage repair pathways. Mutations of this gene in a subset of Hürthle cell tumors point to a potential role for this protein and its associated pathways in thyroid tumorigenesis.

  Study EGAS00001000940

• Whole exome sequencing data for paired non-serous endometrial and ovarian carcinomas from 27 patients with concurrent tumours.

  Introduction: Concurrent non-serous endometrial and ovarian tumours are often managed clinically as two separate primary tumours, but almost all exhibit evidence of a genomic relationship. Methodology: This study investigates the extent of relatedness using whole exome sequencing, which was performed on paired non-serous endometrial and ovarian carcinomas from 27 patients with concurrent tumours in a cohort with detailed clinicopathological annotation. Four whole exome sequencing-derived parameters (mutation, mutational burden, mutational signatures and mutant allele tumour heterogeneity scores) were used to develop a novel unsupervised model for the assessment of genomic similarity to infer genomic relatedness of paired tumours. Results: This novel model demonstrated genomic relatedness across all four parameters in all tumour pairs. Mutations in PTEN, ARID1A, CTNNB1, KMT2D and PIK3CA occurred most frequently and 24 of 27 (89 %) tumour pairs shared identical mutations in at least one of these genes, with all pairs sharing mutations in a number of other genes. Ovarian endometriosis, CTNNB1 exon 3 mutation, and progression and death from disease were present across the similarity ranking. Mismatch repair deficiency was associated with less genomically similar pairs. Discussion: Although there was diversity across the cohort, the presence of genomic similarity in all paired tumours supports the hypothesis that concurrent non-serous endometrial and ovarian carcinomas are genomically related and may have arisen from a common origin. Keywords: Co-existing; Concurrent; Endometrial carcinoma; Genomics; Mutation; Ovarian carcinoma; Paired tumours; Tumour mutational burden.

  Study EGAS00001008259

• Spatio-temporal Profiling of a Rhabdoid Tumor Case Study

  Background Immune checkpoint inhibition (ICI) has revolutionized oncology, offering extended survival and long-term remission in previously incurable cancers. While highly effective in tumors with high mutational burden, lowly mutated cancers, including pediatric malignancies, present low response rate and limited predictive biomarkers. Patients and methods We present a framework for the identification and validation of tumor-reactive T cells as a biomarker to quantify ICI efficacy and as candidates for a personalized TCR-T cell therapy. Therefore, we profiled a pediatric malignant rhabdoid tumor patient with complete remission after ICI therapy using deep single-cell T cell receptor (TCR) repertoire sequencing of the tumor microenvironment (TME) and the peripheral blood. Results Tracking T cell dynamics longitudinally from the tumor to cells in circulation over a time course of 12 months revealed a systemic response and durable clonal expansion of tumor-resident and ICI-induced TCR clonotypes. We functionally validated tumor reactivity of TCRs identified from the TME and the blood by co-culturing patient-derived tumor cells with TCR-engineered autologous T cells. Here, we observed unexpectedly high frequencies of tumor-reactive TCR clonotypes in the TME and confirmed T cell dynamics in the blood post-ICI to predict tumor-reactivity. Conclusion These findings strongly support spatio-temporal tracking of T cell activity in response to ICI to inform therapy efficacy and to serve as a source of tumor-reactive TCRs for personalized TCR-T designs.

  Study EGAS00001008174

• Small variants in mtDNA Canary Islands - WGS Illumina (ITER)

  The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. Besides, because of its relevance, we also assess the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short- read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data.

  Study EGAS00001005679

• Small variants in mtDNA Canary Islands - WES Illumina (ITER)

  The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. Besides, because of its relevance, we also assess the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short- read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data.

  Study EGAS00001005678

• Genome-wide analysis of genetic risk factors for rheumatic heart disease in Aboriginal Australians provides support for pathogenic molecular mimicry

  Background. Rheumatic heart disease (RHD) following Group A Streptococcus (GAS) infections is heritable and prevalent in Indigenous populations. Molecular mimicry between human and GAS proteins triggers pro-inflammatory cardiac valve-reactive T-cells. Methods. Genome-wide genetic analysis was undertaken in 1263 Aboriginal Australians (398 RHD cases; 865 controls). Single nucleotide polymorphisms (SNPs) were genotyped using Illumina HumanCoreExome BeadChips. Direct typing and imputation was used to fine-map the human leukocyte antigen (HLA) region. Epitope binding affinities were mapped for human cross-reactive GAS proteins, including M5 and M6. Results. The strongest genetic association was intronic to HLA-DQA1 (rs9272622; P=1.86x10-7). Conditional analyses showed rs9272622 and/or DQA1*AA16 account for the HLA signal. HLA-DQA1*0101_DQB1*0503 (OR 1.44, 95%CI 1.09-1.90, P=9.56x10-3) and HLA-DQA1*0103_DQB1*0601 (OR 1.27, 95%CI 1.07-1.52, P=7.15x10-3) were risk haplotypes; HLA_DQA1*0301-DQB1*0402 (OR 0.30, 95%CI 0.14-0.65, P=2.36x10-3) was protective. Human myosin cross-reactive N-terminal and B repeat epitopes of GAS M5/M6 bind with higher affinity to DQA1/DQB1 alpha/beta dimers for the two risk haplotypes than the protective haplotype. Conclusions. Variation at HLA_DQA1-DQB1 is the major genetic risk factor for RHD in Aboriginal Australians studied here. Cross-reactive epitopes bind with higher affinity to alpha/beta dimers formed by risk haplotypes, supporting molecular mimicry as the key mechanism of RHD pathogenesis.

  Study EGAS00001002678

• GCAT | Genomes for life: cohort study of the genomes of Catalonia

  The GCAT Study have recruited 20 000 participants aged 40–65 years. Participants who agreed to take part in the study completed a self-administered computer-driven questionnaire, and underwent blood pressure, cardiac frequency and anthropometry measurements. For each participant, blood plasma, blood serum and white blood cells are collected at baseline. A total of 5459 genomic profiles have been characterised by comprehensive genotyping. Genome-wide genotypes have been generated using Illumina Infinium SNP-bead array technology. We chose the Multi-Ethnic Global (MEGAEX, V.2) consortium array, a multipurpose, multiethnic genotyping array with two million selected markers (including previously described germline mutations, insertions-deletions (InDels) and SNPs).We have strictly followed the standard manufacturer recommended automated protocol for the Infinium HTS Assay scanned with a HiScan confocal scanner (Illumina, San Diego, California, USA). Genome Studio V.2011.1 has been used for raw data analysis. Genotyping was performed at the Genomics and Bioinformatics Unit of the PMPPC Institute for Health Science Research Germans Trias i Pujol, in Badalona, Spain. Future plans The first follow-up study started in December 2017 and will end by March 2018. Residences of all subjects will be geocoded during the following year. Several genomic analyses are ongoing, and metabolomic and genomic integrations will be performed to identify underlying genetic variants, as well as environmental factors that influence metabolites. http://dx.doi.org/10.1136/bmjopen-2017-018324

  Study EGAS00001003018

• Long-read trio sequencing of unsolved patients with intellectual disability

  Long-read sequencing (LRS) has the potential to comprehensively identify all medically-relevant genome variation, including variation commonly missed by short read sequencing (SRS) approaches. To determine this potential, we performed LRS up to 40x coverage using the Pacific Biosciences Sequel instrument for 5 patient-parent trios with intellectual disability, whose aetiology remained unresolved after SRS exomes and genomes. A first systematic assessment of LRS coverage showed that ~35Mb of the genome was only accessible by LRS and not SRS. Structural variant (SV) calling yielded on average 28,292 SV calls per individual, totalling 12.9 MB of sequence. Trio-based analyses showed concordance for up to 95% of these SV calls across the genome, and 80% of the LRS SV calls were not identified by SRS. De novo mutation analysis did not identify any de novo SVs, confirming these are rare events. Because of high sequence coverage, we were also able to call single nucleotide substitutions. On average, we identified 3 million substitutions per genome, with a Mendelian inheritance concordance of up to 97%. Of these, approximately 100,000 were located in the ~35Mb of the genome that was only captured by LRS. Moreover, these variants affected the coding sequence of 64 genes, including 32 known to cause Mendelian disorders. Although no clinical diagnoses were made in these 5 patients, our data show the potential added value of LRS compared to SRS for identifying medically relevant genome variation.

  Study EGAS00001004319

• Small variants in mtDNA Canary Islands - WGS Oxford Nanopore Technologies (ITER)

  The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. Besides, because of its relevance, we also assess the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short- read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data.

  Study EGAS00001005677