CanvasDB: a local database infrastructure for analysis of targeted- and whole genome re-sequencing projects

Gustavson syndrome is caused by an in-frame deletion in RBMX associated with potentially disturbed SH3 domain interactions

RNA binding motif protein X-linked (RBMX) encodes the heterogeneous nuclear ribonucleoprotein G (hnRNP G) that regulates splicing, sister chromatid cohesion and genome stability. RBMX knock down experiments in various model organisms highlight the gene's importance for brain development. Deletion of the RGG/RG motif in hnRNP G has previously been associated with Shashi syndrome, however involvement of other hnRNP G domains in intellectual disability remain unknown. In the current study, we present the underlying genetic and molecular cause of Gustavson syndrome. Gustavson syndrome was first reported in 1993 in a large Swedish five-generation family presented with profound X-linked intellectual disability and an early death. Extensive genomic analyses of the family revealed hemizygosity for a novel in-frame deletion in RBMX in affected individuals (NM_002139.4; c.484_486del, p.(Pro162del)). Carrier females were asymptomatic and presented with skewed X-chromosome inactivation, indicating silencing of the pathogenic allele. Affected individuals presented minor phenotypic overlap with Shashi syndrome, indicating a different disease-causing mechanism. Investigation of the variant effect in a neuronal cell line (SH-SY5Y) revealed differentially expressed genes enriched for transcription factors involved in RNA polymerase II transcription. Prediction tools and a fluorescence polarization assay imply a novel SH3-binding motif of hnRNP G, and potentially a reduced affinity to SH3 domains caused by the deletion. In conclusion, we present a novel in-frame deletion in RBMX segregating with Gustavson syndrome, leading to disturbed RNA polymerase II transcription, and potentially reduced SH3 binding. The results indicate that disruption of different protein domains affects the severity of RBMX-associated intellectual disabilities.

SnpHub: an easy-to-set-up web server framework for exploring large-scale genomic variation data in the post-genomic era with applications in wheat

BACKGROUND

The cost of high-throughput sequencing is rapidly decreasing, allowing researchers to investigate genomic variations across hundreds or even thousands of samples in the post-genomic era. The management and exploration of these large-scale genomic variation data require programming skills. The public genotype querying databases of many species are usually centralized and implemented independently, making them difficult to update with new data over time. Currently, there is a lack of a widely used framework for setting up user-friendly web servers to explore new genomic variation data in diverse species.

RESULTS

Here, we present SnpHub, a Shiny/R-based server framework for retrieving, analysing, and visualizing large-scale genomic variation data that can be easily set up on any Linux server. After a pre-building process based on the provided VCF files and genome annotation files, the local server allows users to interactively access single-nucleotide polymorphisms and small insertions/deletions with annotation information by locus or gene and to define sample sets through a web page. Users can freely analyse and visualize genomic variations in heatmaps, phylogenetic trees, haplotype networks, or geographical maps. Sample-specific sequences can be accessed as replaced by detected sequence variations.

CONCLUSIONS

SnpHub can be applied to any species, and we build up a SnpHub portal website for wheat and its progenitors based on published data in recent studies. SnpHub and its tutorial are available at http://guoweilong.github.io/SnpHub/. The wheat-SnpHub-portal website can be accessed at http://wheat.cau.edu.cn/Wheat_SnpHub_Portal/.

webGQT: A Shiny Server for Genotype Query Tools for Model-Based Variant Filtering

Summary

Genotype Query Tools (GQT) were developed to discover disease-causing variations from billions of genotypes and millions of genomes, processes data at substantially higher speed over other existing methods. While GQT has been available to a wide audience as command-line software, the difficulty of constructing queries among non-IT or non-bioinformatics researchers has limited its applicability. To overcome this limitation, we developed webGQT, an easy-to-use tool with a graphical user interface. With pre-built queries across three modules, webGQT allows for pedigree analysis, case-control studies, and population frequency studies. As a package, webGQT allows researchers with less or no applied bioinformatics/IT experience to mine potential disease-causing variants from billions.

Results

webGQT offers a flexible and easy-to-use interface for model-based candidate variant filtering for Mendelian diseases from thousands to millions of genomes at a reduced computation time. Additionally, webGQT provides adjustable parameters to reduce false positives and rescue missing genotypes across all modules. Using a case study, we demonstrate the applicability of webGQT to query non-human genomes. In addition, we demonstrate the scalability of webGQT on large data sets by implementing complex population-specific queries on the 1000 Genomes Project Phase 3 data set, which includes 8.4 billion variants from 2504 individuals across 26 different populations. Furthermore, webGQT supports filtering single-nucleotide variants, short insertions/deletions, copy number or any other variant genotypes supported by the VCF specification. Our results show that webGQT can be used as an online web service, or deployed on personal computers or local servers within research groups.

Availability

webGQT is made available to the users in three forms: 1) as a webserver available at https://vm1138.kaj.pouta.csc.fi/webgqt/, 2) as an R package to install on personal computers, and 3) as part of the same R package to configure on the user's own servers. The application is available for installation at https://github.com/arumds/webgqt.

Expanding the phenotypic spectrum of osteogenesis imperfecta type V including heterotopic ossification of muscle origins and attachments

Background

Osteogenesis imperfecta (OI) is a clinical and genetic heterogeneous group of connective tissue disorders, characterized by bone fragility and a propensity to fracture.

Methods

In this report we describe the clinical phenotype of two patients, a 28‐year‐old woman and her mother (54 years old), both with a history of short stature and multiple fractures.

Results

Exome sequencing revealed the recurring IFITM5:c.‐14 C>T variant causing OI type V. Both patients had several fractures during childhood. CT‐scan and scintigraphy showed ossification of the origin and attachment of muscles and hypertrophic callus formation.

Conclusion

Ossification of the origin and attachment of muscles seems to be part of the phenotype in patients with OI type V.

A novel ECEL1 mutation expands the phenotype of distal arthrogryposis multiplex congenita type 5D to include pretibial vertical skin creases

Arthrogryposis multiplex congenita (AMC) is a heterogeneous disorder characterized by multiple joint contractures often in association with other congenital abnormalities. Pretibial linear vertical creases are a rare finding associated with arthrogryposis, and the etiology of the specific condition is unknown. We aimed to genetically and clinically characterize a boy from a consanguineous family, presenting with AMC and pretibial vertical linear creases on the shins. Whole exome sequencing and variant analysis revealed homozygous novel missense variants of ECEL1 (c.1163T > C, p.Leu388Pro, NM_004826) and MUSK (c.2572C > T, p.Arg858Cys, NM_005592). Both variants are predicted to have deleterious effects on the protein function, with amino acid positions highly conserved among species. The variants segregated in the family, with healthy mother, father, and sister being heterozygous carriers and the index patient being homozygous for both mutations. We report on a unique patient with a novel ECEL1 homozygous mutation, expanding the phenotypic spectrum of Distal AMC Type 5D to include vertical linear skin creases. The homozygous mutation in MUSK is of unknown clinical significance. MUSK mutations have previously shown to cause congenital myasthenic syndrome, a neuromuscular disorder with defects in the neuromuscular junction.

Benchmarking distributed data warehouse solutions for storing genomic variant information

Database URL

https://github.com/ZSI-Bio/variantsdwh.

SNX10 gene mutation leading to osteopetrosis with dysfunctional osteoclasts

Autosomal recessive osteopetrosis (ARO) is a heterogeneous disorder, characterized by defective osteoclastic resorption of bone that results in increased bone density. We have studied nine individuals with an intermediate form of ARO, from the county of Västerbotten in Northern Sweden. All afflicted individuals had an onset in early infancy with optic atrophy, and in four patients anemia was present at diagnosis. Tonsillar herniation, foramen magnum stenosis, and severe osteomyelitis of the jaw were common clinical features. Whole exome sequencing, verified by Sanger sequencing, identified a splice site mutation c.212 + 1 G > T in the SNX10 gene encoding sorting nexin 10. Sequence analysis of the SNX10 transcript in patients revealed activation of a cryptic splice site in intron 4 resulting in a frame shift and a premature stop (p.S66Nfs * 15). Haplotype analysis showed that all cases originated from a single mutational event, and the age of the mutation was estimated to be approximately 950 years. Functional analysis of osteoclast progenitors isolated from peripheral blood of patients revealed that stimulation with receptor activator of nuclear factor kappa-B ligand (RANKL) resulted in a robust formation of large, multinucleated osteoclasts which generated sealing zones; however these osteoclasts exhibited defective ruffled borders and were unable to resorb bone in vitro.

TheSNPpit-A High Performance Database System for Managing Large Scale SNP Data

The fast development of high throughput genotyping has opened up new possibilities in genetics while at the same time producing considerable data handling issues. TheSNPpit is a database system for managing large amounts of multi panel SNP genotype data from any genotyping platform. With an increasing rate of genotyping in areas like animal and plant breeding as well as human genetics, already now hundreds of thousand of individuals need to be managed. While the common database design with one row per SNP can manage hundreds of samples this approach becomes progressively slower as the size of the data sets increase until it finally fails completely once tens or even hundreds of thousands of individuals need to be managed. TheSNPpit has implemented three ideas to also accomodate such large scale experiments: highly compressed vector storage in a relational database, set based data manipulation, and a very fast export written in C with Perl as the base for the framework and PostgreSQL as the database backend. Its novel subset system allows the creation of named subsets based on the filtering of SNP (based on major allele frequency, no-calls, and chromosomes) and manually applied sample and SNP lists at negligible storage costs, thus avoiding the issue of proliferating file copies. The named subsets are exported for down stream analysis. PLINK ped and map files are processed as in- and outputs. TheSNPpit allows management of different panel sizes in the same population of individuals when higher density panels replace previous lower density versions as it occurs in animal and plant breeding programs. A completely generalized procedure allows storage of phenotypes. TheSNPpit only occupies 2 bits for storing a single SNP implying a capacity of 4 mio SNPs per 1MB of disk storage. To investigate performance scaling, a database with more than 18.5 mio samples has been created with 3.4 trillion SNPs from 12 panels ranging from 1000 through 20 mio SNPs resulting in a database of 850GB. The import and export performance scales linearly with the number of SNPs and is largely independent of panel and database size. Import speed is around 6 mio SNPs/sec, export between 60 and 120 mio SNPs/sec. Being command line based, imports and exports can easily be integrated into pipelines. TheSNPpit is available under the Open Source GNU General Public License (GPL) Version 2.

A novel variant in MYLK causes thoracic aortic dissections: genotypic and phenotypic description

Background

Mutations in MYLK cause non-syndromic familial thoracic aortic aneurysms and dissections (FTAAD). Very little is known about the phenotype of affected families. We sought to characterize the aortic disease and the presence of other vascular abnormalities in FTAAD caused by a deletion in MYLK and to compare thoracic aortic diameter and stiffness in mutation carriers and non-carriers.

Methods

We studied FTAAD in a 5-generation family that included 19 living members. Exome sequencing was performed to identify the underlying gene defect. Aortic elastic properties measured by TTE, MRI and pulse wave velocity were then compared between mutation carriers and non-carriers.

Results

Exome sequencing led to the identification of a 2-bp deletion in MYLK (c3272_3273del, p.Ser1091*) that led to a premature stop codon and nonsense-mediated decay. Eleven people were mutation carriers and eight people were non-carriers. Five aortic ruptures or dissections occurred in this family, with two survivors. There were no differences in aortic diameter or stiffness between carriers and non-carriers of the mutation.

Conclusions

Individuals carrying this deletion in MYLK have a high risk of presenting with an acute aortic dissection or rupture. Aortic events occur over a wide range of ages and are not always preceded by obvious aortic dilatation. Aortic elastic properties do not differ between carriers and non-carriers of this mutation, rendering it uncertain whether and when carriers should undergo elective prophylactic surgery.

FMFilter: A fast model based variant filtering tool

The availability of whole exome and genome sequencing has completely changed the structure of genetic disease studies. It is now possible to solve the disease causing mechanisms within shorter time and budgets. For this reason, mining out the valuable information from the huge amount of data produced by next generation techniques becomes a challenging task. Current tools analyze sequencing data in various methods. However, there is still need for fast, easy to use and efficacious tools. Considering genetic disease studies, there is a lack of publicly available tools which support compound heterozygous and de novo models. Also, existing tools either require advanced IT expertise or are inefficient for handling large variant files. In this work, we provide FMFilter, an efficient sieving tool for next generation sequencing data produced by genetic disease studies. We develop a software which allows to choose the inheritance model (recessive, dominant, compound heterozygous and de novo), the affected and control individuals. The program provides a user friendly Graphical User Interface which eliminates the requirement of advanced computer techniques. It has various filtering options which enable to eliminate the majority of the false alarms. FMFilter requires negligible memory, therefore it can easily handle very large variant files like multiple whole genomes with ordinary computers. We demonstrate the variant reduction capability and effectiveness of the proposed tool with public and in-house data for different inheritance models. We also compare FMFilter with the existing filtering software. We conclude that FMFilter provides an effective and easy to use environment for analyzing next generation sequencing data from Mendelian diseases.

Integrating 400 million variants from 80,000 human samples with extensive annotations: towards a knowledge base to analyze disease cohorts

BACKGROUND

Data from a plethora of high-throughput sequencing studies is readily available to researchers, providing genetic variants detected in a variety of healthy and disease populations. While each individual cohort helps gain insights into polymorphic and disease-associated variants, a joint perspective can be more powerful in identifying polymorphisms, rare variants, disease-associations, genetic burden, somatic variants, and disease mechanisms.

DESCRIPTION

We have set up a Reference Variant Store (RVS) containing variants observed in a number of large-scale sequencing efforts, such as 1000 Genomes, ExAC, Scripps Wellderly, UK10K; various genotyping studies; and disease association databases. RVS holds extensive annotations pertaining to affected genes, functional impacts, disease associations, and population frequencies. RVS currently stores 400 million distinct variants observed in more than 80,000 human samples.

CONCLUSIONS

RVS facilitates cross-study analysis to discover novel genetic risk factors, gene-disease associations, potential disease mechanisms, and actionable variants. Due to its large reference populations, RVS can also be employed for variant filtration and gene prioritization.

AVAILABILITY

A web interface to public datasets and annotations in RVS is available at https://rvs.u.hpc.mssm.edu/.

Mutation in CEP63 co-segregating with developmental dyslexia in a Swedish family

Developmental dyslexia is the most common learning disorder in children. Problems in reading and writing are likely due to a complex interaction of genetic and environmental factors, resulting in reduced power of studies of the genetic factors underlying developmental dyslexia. Our approach in the current study was to perform exome sequencing of affected and unaffected individuals within an extended pedigree with a familial form of developmental dyslexia. We identified a two-base mutation, causing a p.R229L amino acid substitution in the centrosomal protein 63 kDa (CEP63), co-segregating with developmental dyslexia in this pedigree. This mutation is novel, and predicted to be highly damaging for the function of the protein. 3D modelling suggested a distinct conformational change caused by the mutation. CEP63 is localised to the centrosome in eukaryotic cells and is required for maintaining normal centriole duplication and control of cell cycle progression. We found that a common polymorphism in the CEP63 gene had a significant association with brain white matter volume. The brain regions were partly overlapping with the previously reported region influenced by polymorphisms in the dyslexia susceptibility genes DYX1C1 and KIAA0319. We hypothesise that CEP63 is particularly important for brain development and might control the proliferation and migration of cells when those two events need to be highly coordinated.

Phenotypic expansion of visceral myopathy associated with ACTG2 tandem base substitution

Familial visceral myopathy (FVM) is a rare heritable and heterogeneous condition due to impaired smooth muscle function. We identified a family segregating 11 individuals with a spectrum of visceral symptoms involving the small intestine, colon, biliary tract, urinary tract and uterus. Whole-exome sequencing revealed a novel heterozygous tandem base substitution c.806_807delinsAA (p.(Gly269Glu)) in ACTG2, encoding smooth muscle actin γ-2, in affected family members. Variants in ACTG2 were recently identified in FVM with intestinal pseudo-obstruction as well as with the congenital megacystics-microcolon-intestinal hypoperistalsis syndrome. In our family, eight affected members presented with severe complications from the biliary and/or the urinary tracts in addition to gastrointestinal pseudo-obstructions. Furthermore, all affected mothers had a history of assisted deliveries owing to poor progress during labor and weak uterine contractions. The variable involvement of multiple smooth muscle-dependent organs in our family, including the biliary tract and the uterus, add to the phenotypic spectrum associated with ACTG2 missense variants.

MuSK: a new target for lethal fetal akinesia deformation sequence (FADS)

BACKGROUND

Fetal akinesia deformation sequence syndrome (FADS, OMIM 208150) is characterised by decreased fetal movement (fetal akinesia) as well as intrauterine growth restriction, arthrogryposis, and developmental anomalies (eg, cystic hygroma, pulmonary hypoplasia, cleft palate, and cryptorchidism). Mutations in components of the acetylcholine receptor (AChR) pathway have previously been associated with FADS.

METHODS AND RESULTS

We report on a family with recurrent fetal loss, where the parents had five affected fetuses/children with FADS and one healthy child. The fetuses displayed no fetal movements from the gestational age of 17 weeks, extended knee joints, flexed hips and elbows, and clenched hands. Whole exome sequencing of one affected fetus and the parents was performed. A novel homozygous frameshift mutation was identified in muscle, skeletal receptor tyrosine kinase (MuSK), c.40dupA, which segregated with FADS in the family. Haplotype analysis revealed a conserved haplotype block suggesting a founder mutation. MuSK (muscle-specific tyrosine kinase receptor), a component of the AChR pathway, is a main regulator of neuromuscular junction formation and maintenance. Missense mutations in MuSK have previously been reported to cause congenital myasthenic syndrome (CMS) associated with AChR deficiency.

CONCLUSIONS

To our knowledge, this is the first report showing that a mutation in MuSK is associated with FADS. The results support previous findings that CMS and/or FADS are caused by complete or severe functional disruption of components located in the AChR pathway. We propose that whereas milder mutations of MuSK will cause a CMS phenotype, a complete loss is lethal and will cause FADS.