TRAP1 chaperone protein mutations and autoinflammation

We identified a consanguineous kindred, of three affected children with severe autoinflammation, resulting in the death of one sibling and allogeneic stem cell transplantation in the other two. All three were homozygous for MEFV p.S208C mutation; however, their phenotype was more severe than previously reported, prompting consideration of an oligogenic autoinflammation model. Further genetic studies revealed homozygous mutations in TRAP1, encoding the mitochondrial/ER resident chaperone protein tumour necrosis factor receptor associated protein 1 (TRAP1). Identification of a fourth, unrelated patient with autoinflammation and compound heterozygous mutation of TRAP1 alone facilitated further functional studies, confirming the importance of this protein as a chaperone of misfolded proteins with loss of function, which may contribute to autoinflammation. Impaired TRAP1 function leads to cellular stress and elevated levels of serum IL-18. This study emphasizes the importance of considering digenic or oligogenic models of disease in particularly severe phenotypes and suggests that autoinflammatory disease might be enhanced by bi-allelic mutations in TRAP1.

PANNZER2: a rapid functional annotation web server

The unprecedented growth of high-throughput sequencing has led to an ever-widening annotation gap in protein databases. While computational prediction methods are available to make up the shortfall, a majority of public web servers are hindered by practical limitations and poor performance. Here, we introduce PANNZER2 (Protein ANNotation with Z-scoRE), a fast functional annotation web server that provides both Gene Ontology (GO) annotations and free text description predictions. PANNZER2 uses SANSparallel to perform high-performance homology searches, making bulk annotation based on sequence similarity practical. PANNZER2 can output GO annotations from multiple scoring functions, enabling users to see which predictions are robust across predictors. Finally, PANNZER2 predictions scored within the top 10 methods for molecular function and biological process in the CAFA2 NK-full benchmark. The PANNZER2 web server is updated on a monthly schedule and is accessible at http://ekhidna2.biocenter.helsinki.fi/sanspanz/. The source code is available under the GNU Public Licence v3.

HaploForge: a comprehensive pedigree drawing and haplotype visualization web application

Motivation

Haplotype reconstruction is an important tool for understanding the aetiology of human disease. Haplotyping infers the most likely phase of observed genotypes conditional on constraints imposed by the genotypes of other pedigree members. The results of haplotype reconstruction, when visualized appropriately, show which alleles are identical by descent despite the presence of untyped individuals. When used in concert with linkage analysis, haplotyping can help delineate a locus of interest and provide a succinct explanation for the transmission of the trait locus. Unfortunately, the design choices made by existing haplotype visualization programs do not scale to large numbers of markers. Indeed, following haplotypes from generation to generation requires excessive scrolling back and forth. In addition, the most widely used program for haplotype visualization produces inconsistent recombination artefacts for the X chromosome.

Results

To resolve these issues, we developed HaploForge, a novel web application for haplotype visualization and pedigree drawing. HaploForge takes advantage of HTML5 to be fast, portable and avoid the need for local installation. It can accurately visualize autosomal and X-linked haplotypes from both outbred and consanguineous pedigrees. Haplotypes are coloured based on identity by descent using a novel A* search algorithm and we provide a flexible viewing mode to aid visual inspection. HaploForge can currently process haplotype reconstruction output from Allegro, GeneHunter, Merlin and Simwalk.

Availability and implementation

HaploForge is licensed under GPLv3 and is hosted and maintained via GitHub. https://github.com/mtekman/haploforge.

Contact

r.kleta@ucl.ac.uk.

Supplementary information

Supplementary data are available at Bioinformatics online.

TOPAZ: asymmetric suffix array neighbourhood search for massive protein databases

Background

Protein homology search is an important, yet time-consuming, step in everything from protein annotation to metagenomics. Its application, however, has become increasingly challenging, due to the exponential growth of protein databases. In order to perform homology search at the required scale, many methods have been proposed as alternatives to BLAST that make an explicit trade-off between sensitivity and speed. One such method, SANSparallel, uses a parallel implementation of the suffix array neighbourhood search (SANS) technique to achieve high speed and provides several modes to allow for greater sensitivity at the expense of performance.

Results

We present a new approach called asymmetric SANS together with scored seeds and an alternative suffix array ordering scheme called optimal substitution ordering. These techniques dramatically improve both the sensitivity and speed of the SANS approach. Our implementation, TOPAZ, is one of the top performing methods in terms of speed, sensitivity and scalability. In our benchmark, searching UniProtKB for homologous proteins to the Dickeya solani proteome, TOPAZ took less than 3 minutes to achieve a sensitivity of 0.84 compared to BLAST.

Conclusions

Despite the trade-off homology search methods have to make between sensitivity and speed, TOPAZ stands out as one of the most sensitive and highest performance methods currently available.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2290-3) contains supplementary material, which is available to authorized users.

Autoinflammatory periodic fever, immunodeficiency, and thrombocytopenia (PFIT) caused by mutation in actin-regulatory gene WDR1

The importance of actin dynamics in the activation of the inflammasome is becoming increasingly apparent. IL-1β, which is activated by the inflammasome, is known to be central to the pathogenesis of many monogenic autoinflammatory diseases. However, evidence from an autoinflammatory murine model indicates that IL-18, the other cytokine triggered by inflammasome activity, is important in its own right. In this model, autoinflammation was caused by mutation in the actin regulatory gene WDR1 We report a homozygous missense mutation in WDR1 in two siblings causing periodic fevers with immunodeficiency and thrombocytopenia. We found impaired actin dynamics in patient immune cells. Patients had high serum levels of IL-18, without a corresponding increase in IL-18-binding protein or IL-1β, and their cells also secreted more IL-18 but not IL-1β in culture. We found increased caspase-1 cleavage within patient monocytes indicative of increased inflammasome activity. We transfected HEK293T cells with pyrin and wild-type and mutated WDR1 Mutant protein formed aggregates that appeared to accumulate pyrin; this could potentially precipitate inflammasome assembly. We have extended the findings from the mouse model to highlight the importance of WDR1 and actin regulation in the activation of the inflammasome, and in human autoinflammation.

Autoinflammatory periodic fever, immunodeficiency, and thrombocytopenia (PFIT) caused by mutation in actin-regulatory gene WDR1

The importance of actin dynamics in the activation of the inflammasome is becoming increasingly apparent. IL-1β, which is activated by the inflammasome, is known to be central to the pathogenesis of many monogenic autoinflammatory diseases. However, evidence from an autoinflammatory murine model indicates that IL-18, the other cytokine triggered by inflammasome activity, is important in its own right. In this model, autoinflammation was caused by mutation in the actin regulatory gene WDR1 We report a homozygous missense mutation in WDR1 in two siblings causing periodic fevers with immunodeficiency and thrombocytopenia. We found impaired actin dynamics in patient immune cells. Patients had high serum levels of IL-18, without a corresponding increase in IL-18-binding protein or IL-1β, and their cells also secreted more IL-18 but not IL-1β in culture. We found increased caspase-1 cleavage within patient monocytes indicative of increased inflammasome activity. We transfected HEK293T cells with pyrin and wild-type and mutated WDR1 Mutant protein formed aggregates that appeared to accumulate pyrin; this could potentially precipitate inflammasome assembly. We have extended the findings from the mouse model to highlight the importance of WDR1 and actin regulation in the activation of the inflammasome, and in human autoinflammation.

Robust multi-group gene set analysis with few replicates

Background

Competitive gene set analysis is a standard exploratory tool for gene expression data. Permutation-based competitive gene set analysis methods are preferable to parametric ones because the latter make strong statistical assumptions which are not always met. For permutation-based methods, we permute samples, as opposed to genes, as doing so preserves the inter-gene correlation structure. Unfortunately, up until now, sample permutation-based methods have required a minimum of six replicates per sample group.

Results

We propose a new permutation-based competitive gene set analysis method for multi-group gene expression data with as few as three replicates per group. The method is based on advanced sample permutation technique that utilizes all groups within a data set for pairwise comparisons. We present a comprehensive evaluation of different permutation techniques, using multiple data sets and contrast the performance of our method, mGSZm, with other state of the art methods. We show that mGSZm is robust, and that, despite only using less than six replicates, we are able to consistently identify a high proportion of the top ranked gene sets from the analysis of a substantially larger data set. Further, we highlight other methods where performance is highly variable and appears dependent on the underlying data set being analyzed.

Conclusions

Our results demonstrate that robust gene set analysis of multi-group gene expression data is permissible with as few as three replicates. In doing so, we have extended the applicability of such approaches to resource constrained experiments where additional data generation is prohibitively difficult or expensive. An R package implementing the proposed method and supplementary materials are available from the website http://ekhidna.biocenter.helsinki.fi/downloads/pashupati/mGSZm.html.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1403-0) contains supplementary material, which is available to authorized users.

Wasabi: An Integrated Platform for Evolutionary Sequence Analysis and Data Visualization

Wasabi is an open source, web-based environment for evolutionary sequence analysis. Wasabi visualizes sequence data together with a phylogenetic tree within a modern, user-friendly interface: The interface hides extraneous options, supports context sensitive menus, drag-and-drop editing, and displays additional information, such as ancestral sequences, associated with specific tree nodes. The Wasabi environment supports reproducibility by automatically storing intermediate analysis steps and includes built-in functions to share data between users and publish analysis results. For computational analysis, Wasabi supports PRANK and PAGAN for phylogeny-aware alignment and alignment extension, and it can be easily extended with other tools. Along with drag-and-drop import of local files, Wasabi can access remote data through URL and import sequence data, GeneTrees and EPO alignments directly from Ensembl. To demonstrate a typical workflow using Wasabi, we reproduce key findings from recent comparative genomics studies, including a reanalysis of the EGLN1 gene from the tiger genome study: These case studies can be browsed within Wasabi at http://wasabiapp.org:8000?id=usecases. Wasabi runs inside a web browser and does not require any installation. One can start using it at http://wasabiapp.org. All source code is licensed under the AGPLv3.

Séance: reference-based phylogenetic analysis for 18S rRNA studies

Background

Marker gene studies often use short amplicons spanning one or more hypervariable regions from an rRNA gene to interrogate the community structure of uncultured environmental samples. Target regions are chosen for their discriminatory power, but the limited phylogenetic signal of short high-throughput sequencing reads precludes accurate phylogenetic analysis. This is particularly unfortunate in the study of microscopic eukaryotes where horizontal gene flow is limited and the rRNA gene is expected to accurately reflect the species phylogeny. A promising alternative to full phylogenetic analysis is phylogenetic placement, where a reference phylogeny is inferred using the complete marker gene and iteratively extended with the short sequences from a metagenetic sample under study.

Results

Based on the phylogenetic placement approach we built Séance, a community analysis pipeline focused on the analysis of 18S marker gene data. Séance combines the alignment extension and phylogenetic placement capabilities of the Pagan multiple sequence alignment program with a suite of tools to preprocess, cluster and visualise datasets composed of many samples. We showcase Séance by analysing 454 data from a longitudinal study of intestinal parasite communities in wild rufous mouse lemurs (Microcebus rufus) as well as in simulation. We demonstrate both improved OTU picking at higher levels of sequence similarity for 454 data and show the accuracy of phylogenetic placement to be comparable to maximum likelihood methods for lower numbers of taxa.

Conclusions

Séance is an open source community analysis pipeline that provides reference-based phylogenetic analysis for rRNA marker gene studies. Whilst in this article we focus on studying nematodes using the 18S marker gene, the concepts are generic and reference data for alternative marker genes can be easily created. Séance can be downloaded from http://wasabiapp.org/software/seance/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0235-7) contains supplementary material, which is available to authorized users.

SwiftLink: parallel MCMC linkage analysis using multicore CPU and GPU

MOTIVATION

Linkage analysis remains an important tool in elucidating the genetic component of disease and has become even more important with the advent of whole exome sequencing, enabling the user to focus on only those genomic regions co-segregating with Mendelian traits. Unfortunately, methods to perform multipoint linkage analysis scale poorly with either the number of markers or with the size of the pedigree. Large pedigrees with many markers can only be evaluated with Markov chain Monte Carlo (MCMC) methods that are slow to converge and, as no attempts have been made to exploit parallelism, massively underuse available processing power. Here, we describe SWIFTLINK, a novel application that performs MCMC linkage analysis by spreading the computational burden between multiple processor cores and a graphics processing unit (GPU) simultaneously. SWIFTLINK was designed around the concept of explicitly matching the characteristics of an algorithm with the underlying computer architecture to maximize performance.

RESULTS

We implement our approach using existing Gibbs samplers redesigned for parallel hardware. We applied SWIFTLINK to a real-world dataset, performing parametric multipoint linkage analysis on a highly consanguineous pedigree with EAST syndrome, containing 28 members, where a subset of individuals were genotyped with single nucleotide polymorphisms (SNPs). In our experiments with a four core CPU and GPU, SWIFTLINK achieves a 8.5× speed-up over the single-threaded version and a 109× speed-up over the popular linkage analysis program SIMWALK.

AVAILABILITY

SWIFTLINK is available at https://github.com/ajm/swiftlink. All source code is licensed under GPLv3.

Linkage analysis confirms heterogeneity of hereditary gingival fibromatosis

OBJECTIVES

Hereditary Gingival Fibromatosis (HGF) is a rare benign fibrous lesion of the gingival tissues presumably caused by single gene defects. The aim of this study was to identify the genetic defect leading to HGF in an extended pedigree.

MATERIALS AND METHODS

We report the clinical features and genetic analysis of a family affected by HGF. A total of 17 subjects were assessed clinically and had blood samples taken for DNA extraction. Multipoint parametric linkage analysis was performed to identify the possible chromosomal location responsible for HGF in this family.

RESULTS

Presence of severe HGF associated with tooth impaction was confirmed for seven members of this three-generation family. Linkage analysis revealed that loci on chromosomes 7, 10, 13, 15, 16, 17, 19 and 20 were linked to this trait. Previously found mutations in the SOS1 and GINGF loci were therefore excluded by this analysis.

CONCLUSIONS

This study brings further evidence for genetic heterogeneity of HGF and points towards the existence of different, not-yet-identified genes linked to this condition.

Risk HLA-DQA1 and PLA(2)R1 alleles in idiopathic membranous nephropathy

BACKGROUND

Idiopathic membranous nephropathy is a major cause of the nephrotic syndrome in adults, but its etiologic basis is not fully understood. We investigated the genetic basis of biopsy-proven cases of idiopathic membranous nephropathy in a white population.

METHODS

We performed independent genomewide association studies of single-nucleotide polymorphisms (SNPs) in patients with idiopathic membranous nephropathy from three populations of white ancestry (75 French, 146 Dutch, and 335 British patients). The patients were compared with racially matched control subjects; population stratification and quality controls were carried out according to standard criteria. Associations were calculated by means of a chi-square basic allele test; the threshold for significance was adjusted for multiple comparisons (with the Bonferroni method).

RESULTS

In a joint analysis of data from the 556 patients studied (398 men), we identified significant alleles at two genomic loci associated with idiopathic membranous nephropathy. Chromosome 2q24 contains the gene encoding M-type phospholipase A(2) receptor (PLA(2)R1) (SNP rs4664308, P=8.6×10(-29)), previously shown to be the target of an autoimmune response. Chromosome 6p21 contains the gene encoding HLA complex class II HLA-DQ alpha chain 1 (HLA-DQA1) (SNP rs2187668, P=8.0×10(-93)). The association with HLA-DQA1 was significant in all three populations (P=1.8×10(-9), P=5.6×10(-27), and P=5.2×10(-36) in the French, Dutch, and British groups, respectively). The odds ratio for idiopathic membranous nephropathy with homozygosity for both risk alleles was 78.5 (95% confidence interval, 34.6 to 178.2).

CONCLUSIONS

An HLA-DQA1 allele on chromosome 6p21 is most closely associated with idiopathic membranous nephropathy in persons of white ancestry. This allele may facilitate an autoimmune response against targets such as variants of PLA2R1. Our findings suggest a basis for understanding this disease and illuminate how adaptive immunity is regulated by HLA.