Systematic evaluation of NIPT aneuploidy detection software tools with clinically validated NIPT samples.

Non-invasive prenatal testing (NIPT) is a powerful screening method for fetal aneuploidy detection, relying on laboratory and computational analysis of cell-free DNA. Although several published computational NIPT analysis tools are available, no prior comprehensive, head-to-head accuracy comparison of the various tools has been published. Here, we compared the outcome accuracies obtained for clinically validated samples with five commonly used computational NIPT aneuploidy analysis tools (WisecondorX, NIPTeR, NIPTmer, RAPIDR, and GIPseq) across various sequencing depths (coverage) and fetal DNA fractions. The sample set included cases of fetal trisomy 21 (Down syndrome), trisomy 18 (Edwards syndrome), and trisomy 13 (Patau syndrome). We determined that all of the compared tools were considerably affected by lower sequencing depths, such that increasing proportions of undetected trisomy cases (false negatives) were observed as the sequencing depth decreased. We summarised our benchmarking results and highlighted the advantages and disadvantages of each computational NIPT software. To conclude, trisomy detection for lower coverage NIPT samples (e.g. 2.5M reads per sample) is technically possible but can, with some NIPT tools, produce troubling rates of inaccurate trisomy detection, especially in low-FF samples.

pBRIT: gene prioritization by correlating functional and phenotypic annotations through integrative data fusion.

Motivation: Computational gene prioritization can aid in disease gene identification. Here, we propose pBRIT (prioritization using Bayesian Ridge regression and Information Theoretic model), a novel adaptive and scalable prioritization tool, integrating Pubmed abstracts, Gene Ontology, Sequence similarities, Mammalian and Human Phenotype Ontology, Pathway, Interactions, Disease Ontology, Gene Association database and Human Genome Epidemiology database, into the prediction model. We explore and address effects of sparsity and inter-feature dependencies within annotation sources, and the impact of bias towards specific annotations. Results: pBRIT models feature dependencies and sparsity by an Information-Theoretic (data driven) approach and applies intermediate integration based data fusion. Following the hypothesis that genes underlying similar diseases will share functional and phenotype characteristics, it incorporates Bayesian Ridge regression to learn a linear mapping between functional and phenotype annotations. Genes are prioritized on phenotypic concordance to the training genes. We evaluated pBRIT against nine existing methods, and on over 2000 HPO-gene associations retrieved after construction of pBRIT data sources. We achieve maximum AUC scores ranging from 0.92 to 0.96 against benchmark datasets and of 0.80 against the time-stamped HPO entries, indicating good performance with high sensitivity and specificity. Our model shows stable performance with regard to changes in the underlying annotation data, is fast and scalable for implementation in routine pipelines. Availability and implementation: http://biomina.be/apps/pbrit/; https://bitbucket.org/medgenua/pbrit. Supplementary information: Supplementary data are available at Bioinformatics online.

Towards practical privacy-preserving genome-wide association study.

BACKGROUND: The deployment of Genome-wide association studies (GWASs) requires genomic information of a large population to produce reliable results. This raises significant privacy concerns, making people hesitate to contribute their genetic information to such studies. RESULTS: We propose two provably secure solutions to address this challenge: (1) a somewhat homomorphic encryption (HE) approach, and (2) a secure multiparty computation (MPC) approach. Unlike previous work, our approach does not rely on adding noise to the input data, nor does it reveal any information about the patients. Our protocols aim to prevent data breaches by calculating the χ2 statistic in a privacy-preserving manner, without revealing any information other than whether the statistic is significant or not. Specifically, our protocols compute the χ2 statistic, but only return a yes/no answer, indicating significance. By not revealing the statistic value itself but only the significance, our approach thwarts attacks exploiting statistic values. We significantly increased the efficiency of our HE protocols by introducing a new masking technique to perform the secure comparison that is necessary for determining significance. CONCLUSIONS: We show that full-scale privacy-preserving GWAS is practical, as long as the statistics can be computed by low degree polynomials. Our implementations demonstrated that both approaches are efficient. The secure multiparty computation technique completes its execution in approximately 2 ms for data contributed by one million subjects.

Beegle: from literature mining to disease-gene discovery.

Disease-gene identification is a challenging process that has multiple applications within functional genomics and personalized medicine. Typically, this process involves both finding genes known to be associated with the disease (through literature search) and carrying out preliminary experiments or screens (e.g. linkage or association studies, copy number analyses, expression profiling) to determine a set of promising candidates for experimental validation. This requires extensive time and monetary resources. We describe Beegle, an online search and discovery engine that attempts to simplify this process by automating the typical approaches. It starts by mining the literature to quickly extract a set of genes known to be linked with a given query, then it integrates the learning methodology of Endeavour (a gene prioritization tool) to train a genomic model and rank a set of candidate genes to generate novel hypotheses. In a realistic evaluation setup, Beegle has an average recall of 84% in the top 100 returned genes as a search engine, which improves the discovery engine by 12.6% in the top 5% prioritized genes. Beegle is publicly available at http://beegle.esat.kuleuven.be/.

Candidate gene prioritization with Endeavour.

Genomic studies and high-throughput experiments often produce large lists of candidate genes among which only a small fraction are truly relevant to the disease, phenotype or biological process of interest. Gene prioritization tackles this problem by ranking candidate genes by profiling candidates across multiple genomic data sources and integrating this heterogeneous information into a global ranking. We describe an extended version of our gene prioritization method, Endeavour, now available for six species and integrating 75 data sources. The performance (Area Under the Curve) of Endeavour on cross-validation benchmarks using 'gold standard' gene sets varies from 88% (for human phenotypes) to 95% (for worm gene function). In addition, we have also validated our approach using a time-stamped benchmark derived from the Human Phenotype Ontology, which provides a setting close to prospective validation. With this benchmark, using 3854 novel gene-phenotype associations, we observe a performance of 82%. Altogether, our results indicate that this extended version of Endeavour efficiently prioritizes candidate genes. The Endeavour web server is freely available at https://endeavour.esat.kuleuven.be/.

Galahad: a web server for drug effect analysis from gene expression.

Galahad (https://galahad.esat.kuleuven.be) is a web-based application for analysis of drug effects. It provides an intuitive interface to be used by anybody interested in leveraging microarray data to gain insights into the pharmacological effects of a drug, mainly identification of candidate targets, elucidation of mode of action and understanding of off-target effects. The core of Galahad is a network-based analysis method of gene expression. As an input, Galahad takes raw Affymetrix human microarray data from treatment versus control experiments and provides quality control and data exploration tools, as well as computation of differential expression. Alternatively, differential expression values can be uploaded directly. Using these differential expression values, drug target prioritization and both pathway and disease enrichment can be calculated and visualized. Drug target prioritization is based on the integration of the gene expression data with a functional protein association network. The web site is free and open to all and there is no login requirement.

eXtasy: variant prioritization by genomic data fusion.

Massively parallel sequencing greatly facilitates the discovery of novel disease genes causing Mendelian and oligogenic disorders. However, many mutations are present in any individual genome, and identifying which ones are disease causing remains a largely open problem. We introduce eXtasy, an approach to prioritize nonsynonymous single-nucleotide variants (nSNVs) that substantially improves prediction of disease-causing variants in exome sequencing data by integrating variant impact prediction, haploinsufficiency prediction and phenotype-specific gene prioritization.

Identification of a Novel Homozygous Mutation in PRDM12 Gene in a Patient with Hereditary Sensory and Autonomic Neuropathy Type VIII.

Hereditary sensory autonomic neuropathy type VIII (HSAN-VIII) is a rare genetic disease that occurs due to mutations in the PRDM12 gene. Here, we describe a novel homozygous mutation c.826_840dupTGCAACCGCCGCTTC (p.Cys276_Phe280dup) on exon 5 in the PRDM12 gene identified by WES and confirmed using Sanger sequencing method.

Towards solving the genetic diagnosis odyssey in Iranian patients with congenital anomalies.

Understanding the underlying causes of congenital anomalies (CAs) can be a complex diagnostic journey. We aimed to assess the efficiency of exome sequencing (ES) and chromosomal microarray analysis (CMA) in patients with CAs among a population with a high fraction of consanguineous marriage. Depending on the patient's symptoms and family history, karyotype/Quantitative Fluorescence- Polymerase Chain Reaction (QF-PCR) (n = 84), CMA (n = 81), ES (n = 79) or combined CMA and ES (n = 24) were performed on 168 probands (66 prenatal and 102 postnatal) with CAs. Twelve (14.28%) probands were diagnosed by karyotype/QF-PCR and seven (8.64%) others were diagnosed by CMA. ES findings were conclusive in 39 (49.36%) families, and 61.90% of them were novel variants. Also, 64.28% of these variants were identified in genes that follow recessive inheritance in CAs. The diagnostic rate (DR) of ES was significantly higher than that of CMA in children from consanguineous families (P = 0·0001). The highest DR by CMA was obtained in the non-consanguineous postnatal subgroup and by ES in the consanguineous prenatal subgroup. In a population that is highly consanguineous, our results suggest that ES may have a higher diagnostic yield than CMA and should be considered as the first-tier test in the evaluation of patients with congenital anomalies.

Clinical whole-exome sequencing analysis reveals a novel missense COL11A1 mutation resulting in an 18-week Iranian male aborted fetus with Fibrochondrogenesis 1: A case report.

Fibrochondrogenesis 1, an autosomal recessive syndrome, is a rare disease that causes short-limbed skeletal dysplasia. Mutations in the gene encoding the α1 chain of type XI collagen (COL11A1) are seen to be the main cause of this disease. We present an 18-week Iranian male aborted fetus with Fibrochondrogenesis 1 from consanguineous parents. Whole-exome sequencing revealed a novel missense variant from G to A in exon 45 of 68 in the COL11A1 gene (NM_080629.2: c.3440G > A, [p.G1147E, g.103404625]). The mutation was confirmed by Sanger sequencing and further, MutationTaster predicted this variant to be disease-causing. Bioinformatic analysis suggests that this variant is highly conserved in both nucleotide and protein levels, suggesting that it has an important function in the proper role of COL11A1 protein. In silico analysis suggests that this mutation alters the COL11A1 protein structure through a Glycine to Glutamic acid substitution.

Recommendations for whole genome sequencing in diagnostics for rare diseases.

In 2016, guidelines for diagnostic Next Generation Sequencing (NGS) have been published by EuroGentest in order to assist laboratories in the implementation and accreditation of NGS in a diagnostic setting. These guidelines mainly focused on Whole Exome Sequencing (WES) and targeted (gene panels) sequencing detecting small germline variants (Single Nucleotide Variants (SNVs) and insertions/deletions (indels)). Since then, Whole Genome Sequencing (WGS) has been increasingly introduced in the diagnosis of rare diseases as WGS allows the simultaneous detection of SNVs, Structural Variants (SVs) and other types of variants such as repeat expansions. The use of WGS in diagnostics warrants the re-evaluation and update of previously published guidelines. This work was jointly initiated by EuroGentest and the Horizon2020 project Solve-RD. Statements from the 2016 guidelines have been reviewed in the context of WGS and updated where necessary. The aim of these recommendations is primarily to list the points to consider for clinical (laboratory) geneticists, bioinformaticians, and (non-)geneticists, to provide technical advice, aid clinical decision-making and the reporting of the results.

Multiple Sclerosis Data Alliance - A global multi-stakeholder collaboration to scale-up real world data research.

The Multiple Sclerosis Data Alliance (MSDA), a global multi-stakeholder collaboration, is working to accelerate research insights for innovative care and treatment for people with multiple sclerosis (MS) through better use of real-world data (RWD). Despite the increasing reliance on RWD, challenges and limitations complicate the generation, collection, and use of these data. MSDA aims to tackle sociological and technical challenges arising with scaling up RWD, specifically focused on MS data. MSDA envisions a patient-centred data ecosystem in which all stakeholders contribute and use big data to co-create the innovations needed to advance timely treatment and care of people with MS.

NGS-Logistics: federated analysis of NGS sequence variants across multiple locations.

As many personal genomes are being sequenced, collaborative analysis of those genomes has become essential. However, analysis of personal genomic data raises important privacy and confidentiality issues. We propose a methodology for federated analysis of sequence variants from personal genomes. Specific base-pair positions and/or regions are queried for samples to which the user has access but also for the whole population. The statistics results do not breach data confidentiality but allow further exploration of the data; researchers can negotiate access to relevant samples through pseudonymous identifiers. This approach minimizes the impact on data confidentiality while enabling powerful data analysis by gaining access to important rare samples. Our methodology is implemented in an open source tool called NGS-Logistics, freely available at https://ngsl.esat.kuleuven.be.

Effective diagnosis of genetic disease by computational phenotype analysis of the disease-associated genome.

Less than half of patients with suspected genetic disease receive a molecular diagnosis. We have therefore integrated next-generation sequencing (NGS), bioinformatics, and clinical data into an effective diagnostic workflow. We used variants in the 2741 established Mendelian disease genes [the disease-associated genome (DAG)] to develop a targeted enrichment DAG panel (7.1 Mb), which achieves a coverage of 20-fold or better for 98% of bases. Furthermore, we established a computational method [Phenotypic Interpretation of eXomes (PhenIX)] that evaluated and ranked variants based on pathogenicity and semantic similarity of patients' phenotype described by Human Phenotype Ontology (HPO) terms to those of 3991 Mendelian diseases. In computer simulations, ranking genes based on the variant score put the true gene in first place less than 5% of the time; PhenIX placed the correct gene in first place more than 86% of the time. In a retrospective test of PhenIX on 52 patients with previously identified mutations and known diagnoses, the correct gene achieved a mean rank of 2.1. In a prospective study on 40 individuals without a diagnosis, PhenIX analysis enabled a diagnosis in 11 cases (28%, at a mean rank of 2.4). Thus, the NGS of the DAG followed by phenotype-driven bioinformatic analysis allows quick and effective differential diagnostics in medical genetics.