22929 results for "ega"

in 35.07 milliseconds.

• Unifying recovery dynamics in heterogeneous diseases exemplified by COVID-19

  Studying the recovery dynamics in different diseases poses many difficulties, as patients often display high heterogeneity in their recovery courses. Moreover, most attempts to study disease dynamics focus on disease progression rather than disease recovery mechanisms. To model recovery dynamics, using severe COVID-19 as the example, we align heterogeneous recovery trajectories via a novel computational scheme applied to longitudinally sampled blood transcriptomes. We thus generate pseudotime trajectories, which we then link to cellular and molecular mechanisms based on cell deconvolution analysis and molecular pathway prediction, thus presenting a unique framework for studying recovery processes over time. Specifically, mature neutrophils displayed a gradual decrease during recovery, allowing superior useability for outcome prediction compared to currently used clinical markers. Further, we discovered a recovery-related regulatory change in gene programs resulting in immune rebalancing between interferon and NFkB activity and the restoration of cell homeostasis. We thus propose regulatory mechanisms governing COVID-19 recovery and suggest mature neutrophils and additional gene markers, as novel clinical biomarkers for disease outcome. Overall, we present a novel clinically relevant computational framework for modeling disease recovery, paving the way for future studies of the recovery dynamics in other diseases and tissues.

  Study EGAS00001005735

• Pancreatic tropism of metastatic renal cell carcinoma

  Renal cell carcinoma (RCC) is characterized by a particularly broad metastatic swath, and enigmatically, when the pancreas is a destination, the disease is associated with improved survival. Intrigued by this observation, we sought to characterize the clinical behavior, therapeutic implications, and underlying biology. While pancreatic metastases (PM) are infrequent, we identified 31 patients across two institutional cohorts and show that improved survival is independent of established prognostic variables, that these tumors are exquisitely sensitive to anti-angiogenic agents and resistant to immune checkpoint inhibitors (ICI), and that they are characterized by a distinctive biology. Primary tumors of patients with PM exhibited frequent PBRM1 mutations, 3p loss, and 5q amplification, along with a lower frequency of aggressive features such as BAP1 mutations and loss of 9p, 14q, and 4q. Gene expression analyses revealed constrained evolution with remarkable uniformity, reduced T-effector gene signatures and increased angiogenesis. Similar findings were observed histopathologically. Thus, RCC metastatic to the pancreas is characterized by indolent biology, heightened angiogenesis, and an uninflamed stroma, likely underlying its good prognosis, sensitivity to anti-angiogenic therapies, and refractoriness to ICI. These data suggest that metastatic organotropism may be an indicator of a particular biology with prognostic and treatment implications for patients.

  Study EGAS00001004208

• CancerLocator: Non-Invasive Cancer Diagnosis and Tissue-of-Origin Prediction Using Methylation Profiles of Cell-Free DNA

  Background: The detection and characterization of cell-free DNA in plasma is one of the most promising new areas in cancer diagnosis. Liquid biopsy, unlike traditional tissue biopsy, has the potential to diagnose a variety of different malignancies. Results: Here we propose a probabilistic method, CancerLocator, which exploits the diagnostic potential of cell-free DNA by determining not only the presence but also the location of tumors. CancerLocator simultaneously infers the proportions and the tissue-of-origin of tumor-derived cell-free DNA in a blood sample using genome-wide DNA methylation data. We comprehensively evaluate CancerLocator with simulations and real data, and compare its performance with that of two established multi-class classification methods. We show that the predicted tumor burdens are highly consistent with the true values. In addition, when the proportion of tumor-derived DNAs in the cell-free DNAs is low, the two popular machine learning methods completely fail for cancer diagnosis, while CancerLocator successfully overcomes the challenge. CancerLocator also achieves promising results on patient plasma samples, despite the fact that the DNA methylation data from these samples has very low sequencing coverage. Conclusions: CfDNA methylation may be developed as an important approach for non-invasive early cancer diagnosis.

  Study EGAS00001002211

• Host pathogen interaction long read transcriptome

  Our incomplete knowledge of the human transcriptome impairs the detection of disease-causing variants, in particular in transcripts only expressed under certain conditions. These transcripts are often lacking from reference transcript sets, such as Ensembl/GENCODE and RefSeq, and could be relevant for establishing genetic diagnoses. We present SUsPECT (Solving Unsolved Patient Exomes/gEnomes using Custom Transcriptomes), a pipeline based on the Ensembl Variant Effect Predictor (VEP) to predict variant impact on custom transcript sets, such as those generated by long-read RNA-sequencing, for downstream prioritization. Our pipeline predicts the functional consequence and likely deleteriousness scores for missense variants in the context of novel open reading frames predicted from any transcriptome. We demonstrate the utility of SUsPECT by uncovering potential mutational mechanisms of pathogenic variants in ClinVar that are predicted to be benign using the reference transcript annotation. In further support of SUsPECT’s utility, we identified an enrichment of immune-related variants predicted to have a more severe molecular consequence when annotating with a newly generated transcriptome from stimulated immune cells instead of the reference transcriptome. Our pipeline outputs crucial information for further prioritization of potentially disease-causing variants for any disease and will become increasingly useful as more long-read RNA sequencing datasets become available.

  Study EGAS00001006779

• Direct Detection of Early-stage Cancers Using Circulating Tumor DNA

  Early detection and intervention are likely to be the most effective means for reducing morbidity and mortality of human cancer. However, development of methods for noninvasive detection of early-stage tumors has remained a challenge. We have developed an approach called targeted error correction sequencing (TEC-Seq) that allows ultrasensitive direct evaluation of sequence changes in circulating cell-free DNA using massively parallel sequencing. We have used this approach to examine 58 cancer-related genes encompassing 81 kb. Analysis of plasma from 44 healthy individuals identified genomic changes related to clonal hematopoiesis in 16% of asymptomatic individuals but no alterations in driver genes related to solid cancers. Evaluation of 200 patients with colorectal, breast, lung, or ovarian cancer detected somatic mutations in the plasma of 71, 59, 59, and 68%, respectively, of patients with stage I or II disease. Analyses of mutations in the circulation revealed high concordance with alterations in the tumors of these patients. In patients with resectable colorectal cancers, higher amounts of preoperative circulating tumor DNA were associated with disease recurrence and decreased overall survival. These analyses provide a broadly applicable approach for noninvasive detection of early-stage tumors that may be useful for screening and management of patients with cancer.

  Study EGAS00001002577

• Spatiotemporal Genomic Profiling of Intestinal Metaplasia Reveals Clonal Dynamics of Gastric Cancer Progression

  Intestinal metaplasia (IM) is a pre-malignant condition of the gastric mucosa associated with increased gastric cancer (GC) risk. We analyzed 1256 gastric samples (1152 IMs) from 692 subjects from a 10-year prospective study. We identified 26 IM driver genes in diverse pathways including chromatin regulation (ARID1A) and intestinal homeostasis (SOX9), largely occurring as subclonal events. Analysis of clonal dynamics between and within subjects, and also longitudinally across time, revealed that IM clones are likely transient but increase in size upon progression to dysplasia, with eventual transmission of genetic events to paired GCs. Single-cell and spatial profiling highlighted changes in tissue ecology and lineage heterogeneity in IM, including an intestinal stem-cell dominant cellular compartment linked to early malignancy. Expanded transcriptome profiling revealed expression-based molecular subtypes of IM, including a body-resident “pseudoantralized” subtype associated with incomplete histology, antral/intestinal cell types, ARID1A mutations, inflammation, and microbial communities normally associated with the healthy oral tract. We demonstrate that combined clinical-genomic models outperform clinical-only models in predicting IMs likely to progress. Our results raise opportunities for GC precision prevention and interception by highlighting strategies for accurately identifying IM patients at high GC risk and a role for microbial dysbiosis in IM progression.

  Study EGAS00001007067

• Long-term organoid culture of a small intestinal neuroendocrine tumor

  Gastroenteropancreatic neuroendocrine tumors (GEP-NETs) are rare and highly heterogeneous neoplasms whose incidence has markedly increased over the last decades. A grading system based on the tumor cells' proliferation index predicts high-risk for G3 NETs. However, low-to-intermediate grade (G1/G2) NETs have an unpredictable clinical course that varies from indolent to highly malignant. Cultures of human cancer cells enable to perform functional perturbation analyses that are instrumental to enhance our understanding of cancer biology. To date, no tractable and reliable long-term culture of G1/G2 NET has been reported to permit disease modeling and pharmacological screens. Here, we report of the first long-term culture of a G2 metastatic small intestinal NET that preserves the main genetic drivers of the tumor and retains expression patterns of the endocrine cell lineage. Replicating the tissue, this long-term culture showed a low proliferation index, and yet it could be propagated continuously without dramatic changes in the karyotype. The model was readily available for pharmacological screens using targeted agents and as expected, showed low tumorigenic capacity in vivo. Overall, this is the first long-term culture of NETs to faithfully recapitulate many aspects of the original neuroendocrine tumor.

  Study EGAS00001007093

• High Grade Serous Ovarian Carcinomas Originate in the Fallopian Tube

  High-grade serous ovarian carcinoma (HGSOC) is the most frequent type of ovarian cancer and has a poor outcome. It has been proposed that fallopian tube cancers may be precursors of HGSOC but evolutionary evidence for this hypothesis has been limited. To provide insights into the origins of HGSOC, we We performed whole-exome sequence and copy number analyses of laser-capture microdissected fallopian tube lesions (p53 signatures, serous tubal intraepithelial carcinomas (STICs), and fallopian tube carcinomas), ovarian cancers, and metastases from nine patients with HGSOC. The majority of tumor-specific sequence and structural alterations in ovarian cancers were present in STICs, including those affecting TP53, BRCA1, BRCA2 or PTEN genes. An eEvolutionary analyseis revealed that p53 signatures and STICs were are the precursors of ovarian carcinoma which in turn gave rise to metastatic lesionsand identify a. In one patient we identified a second STIC as a metastasis in the fallopian tube opposite from the affected ovary. These analyses revealed a window of seven years between the development of a STIC and the initiation of ovarian carcinomas, with development of metastases following rapidly thereafter. Our results provide insights into the etiology of ovarian cancer and have implications for the prevention, early detection and therapeutic intervention of this disease.

  Study EGAS00001002589

• Mapping_gene_environment_interactions_in_macrophages

  Immune cells sense and respond to external stimuli such as pathogens and initiate an inflammatory response. Genetic variants altering the behaviour and responses of immune cells may inform the results of disease association studies and provide an insight into the regulation of immune response. Studies mapping quantitative trait loci (QTL) are usually performed in quiescent or naïve cells and therefore any variants with a condition-specific manner (e.g. those with an effect only detectable in a specific time-point after cell stimulation) are likely to remain undiscovered. So far, only a handful of studies have attempted to study QTL in multiple states and even on those studies the number of states has remained relatively low (<5). Only very recently Horst and colleagues published the largest study to date investigating the effect of host variable (e.g. age, sex, contraceptive use etc.) on protein levels in macrophages stimulated with 19 different conditions (Horst et al, 2016). This project aims to develop a high-throughput experimental platform to detect QTL in a large number of states using iPSC derived macrophages. 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 EGAS00001002268

• Single cell transcriptomic landscape of pediatric B-cell acute lymphoblastic leukemia: dissection of transcriptional heterogeneity and B-cell developmental state

  Introduction: Genomic analyses of bulk ALL samples have improved our understanding of the genetic basis and risk stratification of B-ALL, but do not directly examine intratumor heterogeneity or enable inference of leukemia developmental state and cell of origin. Methods: We profiled 89 B-ALL samples by single-cell RNA-seq (scRNA-seq) (10X Genomics 5’v2) and compared them to a scRNA-seq reference map of normal human B-cell development. Results: Analysis of heterogeneity of inferred DNA copy number alterations at single cell level showed that aneuploid ALL with near haploidy and hyperdiploid harbor chromosomal losses or gains in all blasts, consistent with early, synchronous origin rather than sequential accumulation. Intra-sample heterogeneity was driven by five gene expression signatures: cell cycle, metabolism, differentiation, inflammation, refining leukemic subtyping. An atlas of B cell development was established to map B-ALL samples to specific developmental states. Hyperdiploid samples resembled pro-B cells, while MEF2D-r and TCF3::PBX1 showed high pre-B abundance. Patients with high early lymphoid abundance had poorer outcomes. Conclusions: Understanding variation in transcriptional programs and developmental states of B-ALL blasts by sc transcriptome refines existing clinical and genomic classifications and provides novel prognostic markers.

  Study EGAS00001007512

• Ewing and Ewing-like-sarcoma tumoroid biobank reveals distinct drug sensitivities: translocation makes the difference

  Ewing sarcoma (ES) and Ewing-like sarcomas (ELS) are highly aggressive tumors found in soft tissues and bone of mostly children and young adults. Despite being different in many aspects, including genetics, possible cell-of-origin, and pathology, patients with any of these entities all receive the same therapeutic regimen. Currently, there is a lack of preclinical cell models that can predict drug response of these entities in vitro. Here, we describe a pediatric ES and ELS tumor organoid (tumoroid) biobank containing tumoroids with different translocations, including EWSR1::FLI1, EWSR1::ERG, CIC::DUX4, and BCOR-rearrangements. Using histology, whole genome sequencing and RNA sequencing, we demonstrate that these tumoroids retain key properties of their matching patient tumors. In addition, we traced tumor subclones in the tumoroids across patient-matched longitudinal samples, which showed that cellular heterogeneity was maintained. We subsequently performed drug screening on the tumoroid models, which unveiled entity-specific drug sensitivity to various cytotoxic compounds and targeted compounds, including MCL-1 inhibitors for CIC::DUX4 sarcomas. Taken together, this newly established ES and ELS patient-derived tumoroid biobank generates a unique source of material for future basic cancer research and drug-screening.

  Study EGAS00001007941

• Integrated genomic analysis for HCC

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

  Study EGAS00001007957

• Detection of cancers three years prior to diagnosis using plasma cell-free DNA

  A major goal of early cancer detection is to identify subclinical disease when the tumor burden is low, so that treatments are more effective. But how early can cancers be detected prior to clinical signs or symptoms? This question can be answered only through the evaluation of participants whose clinical course has not been altered by the study itself. We here describe such an evaluation, performed on prospectively collected plasma samples from the Atherosclerosis Risk in Communities (ARIC) study, including 26 participants diagnosed with cancer and 26 matched controls. At the index time point, eight of these 52 participants scored positively with a multicancer early detection (MCED) blood test. All eight of these participants were diagnosed with cancer within 4 months after blood collection. In six of these 8 participants, we were able to assess an earlier plasma sample collected 3.1 to 3.5 years prior to clinical diagnosis. In four of these six participants, the same mutations detected by the MCED test could be identified, but at 8.6 to 79-fold lower levels. These results demonstrate that it is possible to detect circulating tumor DNA (ctDNA) more than three years prior to clinical diagnosis, and provide benchmark sensitivities required for the success of ctDNA-based tests for this purpose.

  Study EGAS00001008068

• Tumor-derived exosomes modulate PD-L1 expression in monocytes

  In chronic lymphocytic leukemia (CLL), monocytes and macrophages are skewed toward protumorigenic phenotypes, including the release of tumor-supportive cytokines and the expression of immunosuppressive molecules such as programmed cell death 1 ligand 1 (PD-L1). To understand the mechanism driving protumorigenic skewing in CLL, we evaluated the role of tumor cell–derived exosomes in the cross-talk with monocytes. We carried out RNA sequencing and proteome analyses of CLL-derived exosomes and identified noncoding Y RNA hY4 as a highly abundant RNA species that is enriched in exosomes from plasma of CLL patients compared with healthy donor samples. Transfer of CLL-derived exosomes or hY4 alone to monocytes resulted in key CLL-associated phenotypes, including the release of cytokines, such as C-C motif chemokine ligand 2 (CCL2), CCL4, and interleukin-6, and the expression of PD-L1. These responses were abolished in Toll-like receptor 7 (TLR7)–deficient monocytes, suggesting exosomal hY4 as a driver of TLR7 signaling. Pharmacologic inhibition of endosomal TLR signaling resulted in a substantially reduced activation of monocytes in vitro and attenuated CLL development in vivo. Our results indicate that exosome-mediated transfer of noncoding RNAs to monocytes contributes to cancer-related inflammation and concurrent immune escape via PD-L1 expression.

  Study EGAS00001002377

• Pooled scRNA-seq of iPSC-derived neural stem cells from ADHD and control individuals

  This dataset comprises paired-end FASTQ files from single-cell RNA sequencing (scRNA-seq) of neural stem cells (NSCs) derived from induced pluripotent stem cells (iPSCs). The iPSCs were reprogrammed from somatic cells (keratinocytes, PBMCs, or lymphoblastoid cell lines) collected from individuals diagnosed with attention-deficit hyperactivity disorder (ADHD) and age-matched healthy controls (ages 6–18 years). Reprogramming was carried out using the non-integrative Sendai virus system (CytoTune-iPS 2.0, Thermo Fisher Scientific), introducing OCT3/4, SOX2, KLF4, and cMYC. iPSCs were cultured in Essential 8 medium on vitronectin and underwent rigorous quality control including immunocytochemistry, qRT-PCR, mycoplasma testing, and SNP genotyping (Illumina Global Screening Array). NSCs were generated using a standardized neural induction protocol and validated via qRT-PCR and immunostaining for markers including PAX6, SOX2, NESTIN, TUJ1, and FOXG1. Two independent scRNA-seq multiplexed experiments were performed: Experiment 1 (ADHD: n = 7; Control: n = 6) and Experiment 2 (ADHD: n = 9; Control: n = 7). Samples were pooled and demultiplexed using the Ensemblex algorithm, incorporating prior genotype data. After quality control filtering (>500 unique RNA molecules per cell), 30,433 single-cell transcriptomes were retained. Differential gene expression analysis was conducted to identify transcriptional signatures associated with ADHD compared to controls.

  Study EGAS00001008169

• Patient-derived pediatric brain tumor organoids faithfully recapitulate primary tumors

  Extensive molecular analysis by many groups has revealed a heterogeneous landscape of embryonal brain tumors, including medulloblastoma (MB) and atypical teratoid/rhabdoid tumor (ATRT). MB, the most common malignant pediatric brain tumor, is comprised of four molecularly distinct subgroups: Wingless (WNT), Sonic Hedgehog (SHH), Group 3 (G3), and Group 4 (G4). Among these, MB-SHH harboring MYCN amplification and TP53 mutation and MB-G3 with MYC overexpression or amplification constitute the most aggressive tumors, have the poorest clinical outcomes, and lack effective therapeutic options for recurrent disease. Here we describe the establishment of tumor organoid (TO) models from patient-derived orthotopic xenografts (PDOX) of MYCN-amplified and TP53-mutant MB-SHH, MYC-amplified MB-G3 and two ATRT which all grew in mice when reimplanted into the cortices of naïve recipient animals (TOX). Methylation, bulk and single-cell RNA sequencing, and whole genome sequencing demonstrated that TOs and TOXs faithfully recapitulated the molecular landscape and heterogeneity of PDOXs and primary patient tumors compared to established cell lines. The development of TOs will facilitate in vitro high-throughput drug or CRISPR screens without the need for use tumor-bearing mice. This will also accelerate the identification and validation of new vulnerabilities and in vivo preclinical testing toward clinical trials.

  Study EGAS00001008305

• Botswana 15 autosomal unlinked microsatellites

  Population structure was investigated in 990 Botswana individuals according to ethno-linguistics, Bantu and Khoisan, and geography (the nine administrative districts) using the AmpFLSTR Identifiler autosomal microsatellite markers. Genetic diversity and forensic parameters were calculated for the overall population, and according to ethno-linguistics and geography. The overall combined power of exclusion (CPE) was 0.9999965412 and the combined match probability 6,28 x 10-19. CPE was highest for the Khoisan Tuu ethnolinguistic group and the Northeast District at 0.9999582029 and 0.9999922652 respectively. CMP ranged from 6.28 x 10-19 (Khoisan Tuu) to 1,02 x 10-18 (Northwest district). Using pairwise genetic distances (FST), analysis of molecular variance (AMOVA), factorial correspondence analysis (FCA), and the unsupervised Bayesian clustering method found in STRUCTURE and landscape genetics software TESS, ethno-linguistics were found to have a greater influence on population structure than geography. FCA showed clustering between Bantu and Khoisan, and within the Bantu. This Bantu sub-structuring was not seen with STRUCTURE and TESS, which detected clustering only between Bantu and Khoisan. The patterns of population structure revealed highlight the need for regional reference databases that include ethno-linguistic and geographic location information. These markers have important potential for bio-anthropological studies as well as for forensic applications.

  Study EGAS00001002380

• Inhibiton of the GABPB1L-containing GABP tetramer is sufficient to reverse replicative immortality in TERT promoter mutant glioblastoma cells.

  Over 50 types of cancer acquire TERT promoter mutations. These single point mutations reactivate telomerase, allowing for indefinite maintenance of telomere length and enabling cellular immortalization. The transcription factor binding site created by the point mutations specifically recruit the ETS factor GABP, a multimeric transcription factor composed of the GABPα and GABPβ subunits. GABP can form two functionally independent transcription factor species – a dimer or a tetramer – depending on which of the structurally distinct GABPβ isoforms is incorporated into the complex. In this study, we show that genetic disruption of GABPβ1L, a tetramer forming isoform of GABPβ that is dispensable in normal development, results in TERT silencing in a TERT promoter mutation dependent manner. Failure to activate TERT expression by GABPβ1L culminates in telomere dysfunction, DNA damage, and mitotic cell death exclusively in TERT promoter mutant cells. Furthermore, exogenous expression of TERT is sufficient to prevent telomere degradation and loss of cell viability in GABPβ1L-reduced lines bearing the mutant TERT promoter. Orthotopic injection of tetramer-deficient mutant TERT promoter GBM cells rendered lower tumor burden and prolonged the overall survival of the tumor-bearing mice. These results highlight the potentially widespread role of GABPB1L in enabling replicative immortality of TERT promoter cancers.

  Study EGAS00001002582

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

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

  Study EGAS00001004290

• Allele Balance Bias Identifies Systematic Genotyping Errors and False Disease Associations

  In recent years, Next Generation Sequencing (NGS) has become a cornerstone of clinical genetics and diagnostics. Many clinical applications require high precision, especially if rare events such as somatic mutations in cancer or genetic variants causing rare diseases need to be identified. Although random sequencing errors can be modeled statistically and deep sequencing minimizes their impact, systematic errors remain a problem even at high depth of coverage. Understanding their source is crucial to increase precision of clinical NGS applications. In this work, we studied the relation between recurrent biases in allele balance (AB), systematic errors and false positive variant calls across a large cohort of human samples analyzed by whole exome sequencing (WES). We have modeled the allele balance distribution for biallelic genotypes in 987 WES samples in order to identify positions recurrently deviating significantly from the expectation, a phenomenon we termed allele balance bias (ABB). Furthermore, we have developed a genotype callability score based on ABB for all positions of the human exome, which detects false positive variant calls that passed state-of-the-art filters. Finally, we demonstrate the use of ABB for detection of false associations proposed by rare variant association studies (RVAS).

  Study EGAS00001003027