Agile parallel bioinformatics workflow management using Pwrake

Aberrant hypomethylation at imprinted differentially methylated regions is involved in biparental placental mesenchymal dysplasia

BACKGROUND

Placental mesenchymal dysplasia (PMD) is a morphological abnormality resembling partial hydatidiform moles. It is often associated with androgenetic/biparental mosaicism (ABM) and complicated by Beckwith-Wiedemann syndrome (BWS), an imprinting disorder. These phenomena suggest an association between PMD and aberrant genomic imprinting, particularly of CDKN1C and IGF2. The existence of another type of PMD containing the biparental genome has been reported. However, the frequency and etiology of biparental PMD are not yet fully understood.

RESULTS

We examined 44 placental specimens from 26 patients with PMD: 19 of these were macroscopically normal and 25 exhibited macroscopic PMD. Genotyping by DNA microarray or short tandem repeat analysis revealed that approximately 35% of the macroscopic PMD specimens could be classified as biparental, while the remainder were ABM. We performed a DNA methylation analysis using bisulfite pyrosequencing of 15 placenta-specific imprinted differentially methylated regions (DMRs) and 36 ubiquitous imprinted DMRs. As expected, most DMRs in the macroscopic PMD specimens with ABM exhibited the paternal epigenotype. Importantly, the biparental macroscopic PMD specimens exhibited frequent aberrant hypomethylation at seven of the placenta-specific DMRs. Allelic expression analysis using single-nucleotide polymorphisms revealed that five imprinted genes associated with these aberrantly hypomethylated DMRs were biallelically expressed. Frequent aberrant hypomethylation was observed at five ubiquitous DMRs, including GRB10 but not ICR2 or ICR1, which regulate the expression of CDKN1C and IGF2, respectively. Whole-exome sequencing performed on four biparental macroscopic PMD specimens did not reveal any pathological genetic abnormalities. Clinical and molecular analyses of babies born from pregnancies with PMD revealed four cases with BWS, each exhibiting different molecular characteristics, and those between BWS and PMD specimens were not always the same.

CONCLUSION

These data clarify the prevalence of biparental PMD and ABM-PMD and strongly implicate hypomethylation of DMRs in the pathogenesis of biparental PMD, particularly placenta-specific DMRs and the ubiquitous GRB10, but not ICR2 or ICR1. Aberrant hypomethylation of DMRs was partial, indicating that it occurs after fertilization. PMD is an imprinting disorder, and it may be a missing link between imprinting disorders and placental disorders incompatible with life, such as complete hydatidiform moles and partial hydatidiform moles.

Heterozygous missense variant of the proteasome subunit β-type 9 causes neonatal-onset autoinflammation and immunodeficiency

Impaired proteasome activity due to genetic variants of certain subunits might lead to proteasome-associated autoinflammatory syndromes (PRAAS). Here we report a de novo heterozygous missense variant of the PSMB9 proteasome subunit gene in two unrelated Japanese infants resulting in amino acid substitution of the glycine (G) by aspartic acid (D) at position 156 of the encoded protein β1i. In addition to PRAAS-like manifestations, these individuals suffer from pulmonary hypertension and immunodeficiency, which are distinct from typical PRAAS symptoms. The missense variant results in impaired immunoproteasome maturation and activity, yet ubiquitin accumulation is hardly detectable in the patients. A mouse model of the heterozygous human genetic variant (Psmb9^G156D/+) recapitulates the proteasome defects and the immunodeficiency phenotype of patients. Structurally, PSMB9 G156D interferes with the β-ring-βring interaction of the wild type protein that is necessary for 20S proteasome formation. We propose the term, proteasome-associated autoinflammatory syndrome with immunodeficiency (PRAAS-ID), to indicate a separate category of autoinflammatory diseases, similar to, but distinct from PRAAS, that describes the patients in this study.

Genetic variants of proteasome subunit genes have been shown to associate with perturbed immune function. Here authors show that a heterozygous missense variant of the immunoproteasome subunit β-type 9 causes an autoinflammatory/immune deficiency syndrome in humans and in a mouse model.

A case of ezetimibe-effective hypercholesterolemia with a novel heterozygous variant in ABCG5

Sitosterolemia is caused by homozygous or compound heterozygous gene mutations in either ATP-binding cassette subfamily G member 5 (ABCG5) or 8 (ABCG8). Since ABCG5 and ABCG8 play pivotal roles in the excretion of neutral sterols into feces and bile, patients with sitosterolemia present elevated levels of serum plant sterols and in some cases also hypercholesterolemia. A 48-year-old woman was referred to our hospital for hypercholesterolemia. She had been misdiagnosed with familial hypercholesterolemia at the age of 20 and her serum low-density lipoprotein cholesterol (LDL-C) levels had remained about 200-300 mg/dL at the former clinic. Although the treatment of hydroxymethylglutaryl-CoA (HMG-CoA) reductase inhibitors was ineffective, her serum LDL-C levels were normalized by ezetimibe, a cholesterol transporter inhibitor. We noticed that her serum sitosterol and campesterol levels were relatively high. Targeted analysis sequencing identified a novel heterozygous ABCG5 variant (c.203A>T; p.Ile68Asn) in the patient, whereas no mutations were found in low-density lipoprotein receptor (LDLR), proprotein convertase subtilisin/kexin type 9 (PCSK9), or Niemann-Pick C1-like intracellular cholesterol transporter 1 (NPC1L1). While sitosterolemia is a rare disease, a recent study has reported that the incidence of loss-of-function mutation in the ABCG5 or ABCG8 gene is higher than we thought at 1 in 220 individuals. The present case suggests that serum plant sterol levels should be examined and ezetimibe treatment should be considered in patients with hypercholesterolemia who are resistant to HMG-CoA reductase inhibitors.

Open Agile text mining for bioinformatics: the PubAnnotation ecosystem

MOTIVATION

Most currently available text mining tools share two characteristics that make them less than optimal for use by biomedical researchers: they require extensive specialist skills in natural language processing and they were built on the assumption that they should optimize global performance metrics on representative datasets. This is a problem because most end-users are not natural language processing specialists and because biomedical researchers often care less about global metrics like F-measure or representative datasets than they do about more granular metrics such as precision and recall on their own specialized datasets. Thus, there are fundamental mismatches between the assumptions of much text mining work and the preferences of potential end-users.

RESULTS

This article introduces the concept of Agile text mining, and presents the PubAnnotation ecosystem as an example implementation. The system approaches the problems from two perspectives: it allows the reformulation of text mining by biomedical researchers from the task of assembling a complete system to the task of retrieving warehoused annotations, and it makes it possible to do very targeted customization of the pre-existing system to address specific end-user requirements. Two use cases are presented: assisted curation of the GlycoEpitope database, and assessing coverage in the literature of pre-eclampsia-associated genes.

AVAILABILITY AND IMPLEMENTATION

The three tools that make up the ecosystem, PubAnnotation, PubDictionaries and TextAE are publicly available as web services, and also as open source projects. The dictionaries and the annotation datasets associated with the use cases are all publicly available through PubDictionaries and PubAnnotation, respectively.

HaTSPiL: A modular pipeline for high-throughput sequencing data analysis

Background

Next generation sequencing methods are widely adopted for a large amount of scientific purposes, from pure research to health-related studies. The decreasing costs per analysis led to big amounts of generated data and to the subsequent improvement of software for the respective analyses. As a consequence, many approaches have been developed to chain different software in order to obtain reliable and reproducible workflows. However, the large range of applications for NGS approaches entails the challenge to manage many different workflows without losing reliability.

Methods

We here present a high-throughput sequencing pipeline (HaTSPiL), a Python-powered CLI tool designed to handle different approaches for data analysis with a high level of reliability. The software relies on the barcoding of filenames using a human readable naming convention that contains any information regarding the sample needed by the software to automatically choose different workflows and parameters. HaTSPiL is highly modular and customisable, allowing the users to extend its features for any specific need.

Conclusions

HaTSPiL is licensed as Free Software under the MIT license and it is available at https://github.com/dodomorandi/hatspil.

Nonsense mutation in CFAP43 causes normal-pressure hydrocephalus with ciliary abnormalities

Objective

To identify genes related to normal-pressure hydrocephalus (NPH) in one Japanese family with several members with NPH.

Methods

We performed whole-exome sequencing (WES) on a Japanese family with multiple individuals with NPH and identified a candidate gene. Then we generated knockout mouse using CRISPR/Cas9 to confirm the effect of the candidate gene on the pathogenesis of hydrocephalus.

Results

In WES, we identified a loss-of-function variant in CFAP43 that segregated with the disease. CFAP43 encoding cilia- and flagella-associated protein is preferentially expressed in the testis. Recent studies have revealed that mutations in this gene cause male infertility owing to morphologic abnormalities of sperm flagella. We knocked out mouse ortholog Cfap43 using CRISPR/Cas9 technology, resulting in Cfap43-deficient mice that exhibited a hydrocephalus phenotype with morphologic abnormality of motile cilia.

Conclusion

Our results strongly suggest that CFAP43 is responsible for morphologic or movement abnormalities of cilia in the brain that result in NPH.

Identification of a homozygous frameshift variant in RFLNA in a patient with a typical phenotype of spondylocarpotarsal synostosis syndrome

Spondylocarpotarsal synostosis syndrome, a rare syndromic skeletal disorder characterized by disrupted vertebral segmentation with vertebral fusion, scoliosis, short stature, and carpal/tarsal synostosis, has been associated with biallelic truncating mutations in the filamin B gene or monoallelic mutations in the myosin heavy chain 3 gene. We herein report the case of a patient with a typical phenotype of spondylocarpotarsal synostosis syndrome who had a homozygous frameshift mutation in the refilin A gene (RFLNA) [c.241delC, p.(Leu81Cysfs*111)], which encodes one of the filamin-binding proteins. Refilins, filamins, and myosins play critical roles in forming perinuclear actin caps, which change the nuclear morphology during cell migration and differentiation. The present study implies that RFLNA is an additional causative gene for spondylocarpotarsal synostosis syndrome in humans and a defect in forming actin bundles and perinuclear actin caps may be a critical mechanism for the development of spondylocarpotarsal synostosis syndrome.

Patients with SATB2-associated syndrome exhibiting multiple odontomas

Special AT-rich sequence-binding protein 2 (SATB2)-associated syndrome (SAS) is characterized by alterations of SATB2. Its clinical features include intellectual disability and craniofacial abnormalities, such as cleft palate, dysmorphic features, and dental abnormalities. Here, we describe three previously undiagnosed, unrelated patients with SAS who exhibited dental abnormalities, including multiple odontomas. Although isolated odontomas are common, multiple odontomas are rare. Individuals in families 1 and 3 underwent whole-exome sequencing. Patient 2 and parents underwent targeted amplicon sequencing. On the basis of the hg19/GRCh37 reference and the RefSeq mRNA NM_001172517, respective heterozygous mutations were found and validated in Patients 1, 2, and 3: a splice-site mutation (chr2:g.200137396C > T, c.1741-1G > A), a nonsense mutation (chr2:g.200213750G > A, c.847C > T, p.R283*), and a frame-shift mutations (chr2:g.200188589_200188590del, c.1478_1479del, p.Q493Rfs*19). All mutations occurred de novo. The mutations in Patients 1 and 3 were novel; the mutation in Patient 2 has been described previously. Tooth mesenchymal cells derived from Patient 2 showed diminished SATB2 expression. Multiple odontomas were evident in the patients in this report; however, this has not been recognized previously as a SAS-associated phenotype. We propose that multiple odontomas be considered as an occasional manifestation of SAS.

Whole-exome sequencing and gene-based rare variant association tests suggest that PLA2G4E might be a risk gene for panic disorder

Panic disorder (PD) is characterized by recurrent and unexpected panic attacks, subsequent anticipatory anxiety, and phobic avoidance. Recent epidemiological and genetic studies have revealed that genetic factors contribute to the pathogenesis of PD. We performed whole-exome sequencing on one Japanese family, including multiple patients with panic disorder, which identified seven rare protein-altering variants. We then screened these genes in a Japanese PD case-control group (384 sporadic PD patients and 571 controls), resulting in the detection of three novel single nucleotide variants as potential candidates for PD (chr15: 42631993, T>C in GANC; chr15: 42342861, G>T in PLA2G4E; chr20: 3641457, G>C in GFRA4). Statistical analyses of these three genes showed that PLA2G4E yielded the lowest p value in gene-based rare variant association tests by Efficient and Parallelizable Association Container Toolbox algorithms; however, the p value did not reach the significance threshold in the Japanese. Likewise, in a German case-control study (96 sporadic PD patients and 96 controls), PLA2G4E showed the lowest p value but again did not reach the significance threshold. In conclusion, we failed to find any significant variants or genes responsible for the development of PD. Nonetheless, our results still leave open the possibility that rare protein-altering variants in PLA2G4E contribute to the risk of PD, considering the function of this gene.

Deep sequencing reveals variations in somatic cell mosaic mutations between monozygotic twins with discordant psychiatric disease

Monozygotic (MZ) twins have been thought to be genetically identical. However, recent studies have shown discordant variants between them. We performed whole-exome sequencing (WES) in five MZ twin pairs with discordant neurodevelopmental disorders and one healthy control MZ twin to detect discordant variants. We identified three discordant variants confirmed by deep sequencing after analysis by personalized next-generation sequencing (NGS). Three mutations in FBXO38 (chr5:147774428;T>G), SMOC2 (chr6:169051385;A>G) and TDRP (chr8:442616;A>G), were detected with low allele frequency of mutant alleles on deep sequencing, suggesting that these loci are mosaic due to somatic mutations in a developmental stage. Our results suggest that deep sequencing analysis would be an adequate method to detect discordant mutations in candidate genes responsible for heritable diseases.

Clinical and histological findings of autosomal dominant renal-limited disease with LMX1B mutation

AIM

Mutations of LMX1B cause nail-patella syndrome, a rare autosomal dominant disorder. Recently, LMX1B R246Q heterozygous mutations were recognised in nephropathy without extrarenal manifestation. The aim of this study was to clarify characteristics of nephropathy caused by R246Q mutation.

METHODS

Whole exome sequencing was performed on a large family with nonsyndromic autosomal dominant nephropathy without extrarenal manifestation. Clinical and histological findings of patients with LMX1B mutation were investigated.

RESULTS

LMX1B R246Q heterozygous mutation was identified in five patients over three generations. Proteinuria or haematoproteinuria was recognized by urinary screening from all patients in childhood. Proteinuria gradually increased to nephrotic levels and renal function decreased in adolescence. Two patients progressed to end-stage renal disease in adulthood. Renal histology demonstrated minimal change in childhood and focal segmental glomerulosclerosis in adulthood. Using electron microscopy, focal collagen deposition could be detected in glomeruli even when a "moth-eaten appearance" was not apparent in the glomerular basement membrane. In addition, podocin expression in glomerular podocytes was significantly decreased, even in the early stages of disease progression.

CONCLUSION

Comprehensive genetic analyses and collagen or tannic acid staining may be useful for diagnosis of LMX1B-associated nephropathy. While renal prognosis of R246Q may be worse than that of typical NPS nephropathy, signs of podocytopathy can be detected during the infantile period; thus, childhood urinary screening may facilitate early detection.

De novo mutations in SMCHD1 cause Bosma arhinia microphthalmia syndrome and abrogate nasal development

Bosma arhinia microphthalmia syndrome (BAMS) is an extremely rare and striking condition characterized by complete absence of the nose with or without ocular defects. We report here that missense mutations in the epigenetic regulator SMCHD1 mapping to the extended ATPase domain of the encoded protein cause BAMS in all 14 cases studied. All mutations were de novo where parental DNA was available. Biochemical tests and in vivo assays in Xenopus laevis embryos suggest that these mutations may behave as gain-of-function alleles. This finding is in contrast to the loss-of-function mutations in SMCHD1 that have been associated with facioscapulohumeral muscular dystrophy (FSHD) type 2. Our results establish SMCHD1 as a key player in nasal development and provide biochemical insight into its enzymatic function that may be exploited for development of therapeutics for FSHD.

Novel compound heterozygous DNA ligase IV mutations in an adolescent with a slowly-progressing radiosensitive-severe combined immunodeficiency

We herein describe a case of a 17-year-old boy with intractable common warts, short stature, microcephaly and slowly-progressing pancytopenia. Simultaneous quantification of T-cell receptor recombination excision circles (TREC) and immunoglobulin κ-deleting recombination excision circles (KREC) suggested very poor generation of both T-cells and B-cells. By whole exome sequencing, novel compound heterozygous mutations were identified in the patient's DNA ligase IV (LIG4) gene. The diagnosis of LIG4 syndrome was confirmed by delayed DNA double-strand break repair kinetics in γ-irradiated fibroblasts from the patient and their restoration by an introduction of wild-type LIG4. Although the patient received allogeneic hematopoietic stem cell transplantation from his haploidentical mother, he unfortunately expired due to an insufficiently reconstructed immune system. An earlier definitive diagnosis using TREC/KREC quantification and whole exome sequencing would thereby allow earlier intervention, which would be essential for improving long-term survival in similar cases with slowly-progressing LIG4 syndrome masked in adolescents.

DRAW+SneakPeek: analysis workflow and quality metric management for DNA-seq experiments

Summary: We report our new DRAW+SneakPeek software for DNA-seq analysis. DNA resequencing analysis workflow (DRAW) automates the workflow of processing raw sequence reads including quality control, read alignment and variant calling on high-performance computing facilities such as Amazon elastic compute cloud. SneakPeek provides an effective interface for reviewing dozens of quality metrics reported by DRAW, so users can assess the quality of data and diagnose problems in their sequencing procedures. Both DRAW and SneakPeek are freely available under the MIT license, and are available as Amazon machine images to be used directly on Amazon cloud with minimal installation.

Availability: DRAW+SneakPeek is released under the MIT license and is available for academic and nonprofit use for free. The information about source code, Amazon machine images and instructions on how to install and run DRAW+SneakPeek locally and on Amazon elastic compute cloud is available at the National Institute on Aging Genetics of Alzheimer’s Disease Data Storage Site (http://www.niagads.org/) and Wang lab Web site (http://wanglab.pcbi.upenn.edu/).

Contact:gerardsc@mail.med.upenn.edu or lswang@mail.med.upenn.edu

The Ruby UCSC API: accessing the UCSC genome database using Ruby

Background

The University of California, Santa Cruz (UCSC) genome database is among the most used sources of genomic annotation in human and other organisms. The database offers an excellent web-based graphical user interface (the UCSC genome browser) and several means for programmatic queries. A simple application programming interface (API) in a scripting language aimed at the biologist was however not yet available. Here, we present the Ruby UCSC API, a library to access the UCSC genome database using Ruby.

Results

The API is designed as a BioRuby plug-in and built on the ActiveRecord 3 framework for the object-relational mapping, making writing SQL statements unnecessary. The current version of the API supports databases of all organisms in the UCSC genome database including human, mammals, vertebrates, deuterostomes, insects, nematodes, and yeast.

The API uses the bin index—if available—when querying for genomic intervals. The API also supports genomic sequence queries using locally downloaded *.2bit files that are not stored in the official MySQL database. The API is implemented in pure Ruby and is therefore available in different environments and with different Ruby interpreters (including JRuby).

Conclusions

Assisted by the straightforward object-oriented design of Ruby and ActiveRecord, the Ruby UCSC API will facilitate biologists to query the UCSC genome database programmatically. The API is available through the RubyGem system. Source code and documentation are available at https://github.com/misshie/bioruby-ucsc-api/ under the Ruby license. Feedback and help is provided via the website at http://rubyucscapi.userecho.com/.