A likelihood ratio based method to predict exact pedigrees for complex families from next-generation sequencing data

MOTIVATION: Next generation sequencing (NGS) technology considerably changed the way we screen for pathogenic mutations in rare Mendelian disorders. However, the identification of the disease- causing mutation amongst thousands of variants of partly unknown relevance is still challenging and efficient techniques that reduce the genomic search space play a decisive role. Often segregation- or linkage analysis are used to prioritize candidates, however, these approaches require correct information about the degree of relationship among the sequenced samples. For quality assurance an automated control of pedigree structures and sample assignment is therefore highly desirable in order to detect label mix-ups that might otherwise corrupt downstream analysis. RESULTS: We developed an algorithm based on likelihood ratios that discriminates between different classes of relationship for an arbitrary number of genotyped samples. By identifying the most likely class we are able to reconstruct entire pedigrees iteratively, even for highly consanguineous families. We tested our approach on exome data of different sequencing studies and achieved high precision for all pedigree predictions. By analyzing the precision for varying degrees of relatedness or inbreeding we could show that a prediction is robust down to magnitudes of a few hundred loci. AVAILABILITY: A java standalone application that computes the relationships between multiple samples as well as a Rscript that visualizes the pedigree information is available for download as well as a web service at www.gene-talk.de CONTACT: [email protected]

Communication efficient multi-processor FFT

Computing the Fast Fourier Transform on a distributed memory architecture by a direct pipelined radix-2 algorithm, a bi-section or multi-section algorithm, all yield the same communications requirement, if communication for all FFT stages can be performed concurrently, the input data is in normal order, and the data allocation consecutive. With a cyclic data allocation, or bit-reversed input data and a consecutive allocation, multi-sectioning offers a reduced communications requirement by approximately a factor of two. For a consecutive data allocation, normal input order, a decimation-in-time FFT requires that (P/N) + d - 2 twiddle factors be stored for P elements distributed evenly over N processors, and the axis subject to transformation distributed over 2d processors. No communication of twiddle factors is required. The same storage requirements hold for a decimation-in-frequency FFT, bit-reversed input order, and consecutive data allocation. The opposite combination of FFT type and data ordering requires a factor of log2 N more storage for N processors. The peak performance for a Connection Machine system CM-200 implementation is 12.9 Gflops/s in 32-bit precision, and 10.7 Gflops/s in 64-bit precision for unordered transforms local to each processor. The corresponding execution rates for ordered transforms are 11.1 Gflops/s and 8.5 Gflops/s, respectively. For distributed one- and two-dimensional transforms the peak performance for unordered transforms exceeds 5 Gflops/s in 32-bit precision, and 3 Gflops/s in 64-bit precision. Three-dimensional transforms executes at a slightly lower rate. Distributed ordered transforms executes at a rate of about 1/2 to 2/3 of the unordered transforms.Engineering and Applied Science

Identity-by-descent filtering of exome sequence data for disease–gene identification in autosomal recessive disorders

Motivation: Next-generation sequencing and exome-capture technologies are currently revolutionizing the way geneticists screen for disease-causing mutations in rare Mendelian disorders. However, the identification of causal mutations is challenging due to the sheer number of variants that are identified in individual exomes. Although databases such as dbSNP or HapMap can be used to reduce the plethora of candidate genes by filtering out common variants, the remaining set of genes still remains on the order of dozens

Somatic neurofibromatosis type 1 (NF1) inactivation events in cutaneous neurofibromas of a single NF1 patient

Neurofibromatosis type 1 (NF1) (MIM#162200) is a relatively frequent genetic condition that predisposes to tumor formation. The main types of tumors occurring in NF1 patients are cutaneous and subcutaneous neurofibromas, plexiform neurofibromas, optic pathway gliomas, and malignant peripheral nerve sheath tumors. To search for somatic mutations in cutaneous (dermal) neurofibromas, whole-exome sequencing (WES) was performed on seven spatially separated tumors and two reference tissues (blood and unaffected skin) from a single NF1 patient. Validation of WES findings was done using routine Sanger sequencing or Sequenom IPlex SNP genotyping. Exome sequencing confirmed the existence of a known familial splice-site mutation NM_000267.3:c.3113+1G>A in exon 23 of NF1 gene (HGMD ID CS951480) in blood, unaffected skin, and all tumor samples. In five out of seven analyzed tumors, we additionally detected second-hit mutations in the NF1 gene. Four of them were novel and one was previously observed. Each mutation was distinct, demonstrating the independent origin of each tumor. Only in two of seven tumors we detected an additional somatic mutation that was not associated with NF1. Our study demonstrated that somatic mutations of NF1 are likely the main drivers of cutaneous tumor formation. The study provides evidence for the rareness of single base pair level alterations in the exomes of benign NF1 cutaneous tumors.European Journal of Human Genetics advance online publication, 8 October 2014; doi:10.1038/ejhg.2014.210

Carbon Nanotubes by a CVD Method. Part I: Synthesis and Characterization of the (Mg, Fe)O Catalysts

The controlled synthesis of carbon nanotubes by chemical vapor deposition requires tailored and wellcharacterized catalyst materials. We attempted to synthesize Mg1-xFexO oxide solid solutions by the combustion route, with the aim of performing a detailed investigation of the influence of the synthesis conditions (nitrate/urea ratio and the iron content) on the valency and distribution of the iron ions and phases. Notably, characterization of the catalyst materials is performed using 57Fe Mo¨ssbauer spectroscopy, X-ray diffraction, and electron microscopy. Several iron species are detected including Fe2+ ions substituting for Mg2+ in the MgO lattice, Fe3+ ions dispersed in the octahedral sites of MgO, different clusters of Fe3+ ions, and MgFe2O4-like nanoparticles. The dispersion of these species and the microstructure of the oxides are discussed. Powders markedly different from one another that may serve as model systems for further study are identified. The formation of carbon nanotubes upon reduction in a H2/CH4 gas atmosphere of the selected powders is reported in a companion paper

Combining callers improves the detection of copy number variants from whole-genome sequencing

Copy Number Variants (CNVs) are deletions, duplications or insertions larger than 50 base pairs. They account for a large percentage of the normal genome variation and play major roles in human pathology. While array-based approaches have long been used to detect them in clinical practice, whole-genome sequencing (WGS) bears the promise to allow concomitant exploration of CNVs and smaller variants. However, accurately calling CNVs from WGS remains a difficult computational task, for which a consensus is still lacking. In this paper, we explore practical calling options to reach the best compromise between sensitivity and sensibility. We show that callers based on different signal (paired-end reads, split reads, coverage depth) yield complementary results. We suggest approaches combining four selected callers (Manta, Delly, ERDS, CNVnator) and a regenotyping tool (SV2), and show that this is applicable in everyday practice in terms of computation time and further interpretation. We demonstrate the superiority of these approaches over array-based Comparative Genomic Hybridization (aCGH), specifically regarding the lack of resolution in breakpoint definition and the detection of potentially relevant CNVs. Finally, we confirm our results on the NA12878 benchmark genome, as well as one clinically validated sample. In conclusion, we suggest that WGS constitutes a timely and economically valid alternative to the combination of aCGH and whole-exome sequencing

A CRISPR-Cas9-engineered mouse model for GPI anchor deficiency mirrors human phenotype and shows hippocampal synaptic dysfunctions

Pathogenic germline mutations in PIGV lead to glycosylphosphatidylinositol biosynthesis deficiency. Individuals with pathogenic biallelic mutations in genes of the glycosylphosphatidylinositol anchor pathway show cognitive impairments, a motor delay and in many cases epilepsy. Thus far, the pathophysiology underlying the disease remains unclear and suitable rodent models that mirror human pathophysiology have not been available. We therefore generated a mouse model using CRISPR-Cas9 to introduce the most prevalent hypomorphic missense mutation in European patients, at a site that is also conserved in mice, Pigv:c.1022C>A (p.A341E). Reflecting the human pathology mutant Pigv(341E) mice showed deficits in motor coordination and cognitive impairment with poorer long-term spatial memory than wild-type mice, as well as alterations in sociability and sleep patterns. Furthermore, immunohistochemistry showed decreased synaptophysin-immunoreactivity and electrophysiology recordings demonstrated reduced hippocampal synaptic transmission in Pigv(341E) mice that may underlie impaired memory formation. To gain a deeper and broader molecular understanding of the consequences of glycosylphosphatidylinositol anchor deficiency, we performed single-cell RNA sequencing on acutely isolated hippocampal cells of Pigv(341E) and wild-type mice. We found that hippocampal cells from adult Pigv(341E) mice exhibited changes in gene expression, most prominently in a subtype of microglia and subicular neurons. A significant reduction of Abl1 transcripts in several cell clusters suggests a link to the signaling pathway of glycosylphosphatidylinositol-anchored ephrins. We also observed increased levels of Hdc that might affect histamine metabolism with consequences in circadian rhythm. In summary, we present here the first mouse model with a patient-specific hypomorphic mutation that mirrors the human phenotype and shows a hippocampal synaptic defect. This new mouse model will not only open the doors for further investigation into the pathophysiology of glycosylphosphatidylinositol biosynthesis deficiency in future studies, but will also deepen our understanding in the role of glycosylphosphatidylinositol-anchor related pathways in brain development

Identity-by-descent filtering as a tool for the identification of disease alleles in exome sequence data from distant relatives

Large-scale, deep resequencing may be the next logical step in the genetic investigation of common complex diseases. Because each individual is likely to carry many thousands of variants, the identification of causal alleles requires an efficient strategy to reduce the number of candidate variants. Under many genetic models, causal alleles can be expected to reside within identity-by-descent (IBD) regions shared by affected relatives. In distant relatives, IBD regions constitute a small portion of the genome and can thus greatly reduce the search space for causal alleles. However, the effectiveness of this strategy is unknown. We test the simulated mini-exome data set in extended pedigrees provided by Genetic Analysis Workshop 17. At the fourth- and fifth-degree level of relatedness, case-case pairs shared between 1% and 9% of the genome identical by descent. As expected, no genes were shared identical by descent by all case subjects, but 43 genes were shared by many case subjects across at least 50 replicates. We filtered variants in these genes based on population frequency, function, informativeness, and evidence of association using the family-based association test. This analysis highlighted five genes previously implicated in triglyceride, lipid, and cholesterol metabolism. Comparison with the list of true risk alleles revealed that strict IBD filtering followed by association testing of the rarest alleles was the most sensitive strategy. IBD filtering may be a useful strategy for narrowing down the list of candidate variants in exome data, but the optimal degree of relatedness of affected pairs will depend on the genetic architecture of the disease under study

Inhibitory activity of a standardized elderberry liquid extract against clinically-relevant human respiratory bacterial pathogens and influenza A and B viruses

<p>Abstract</p> <p>Background</p> <p>Black elderberries (<it>Sambucus nigra </it>L.) are well known as supportive agents against common cold and influenza. It is further known that bacterial super-infection during an influenza virus (IV) infection can lead to severe pneumonia. We have analyzed a standardized elderberry extract (Rubini, BerryPharma AG) for its antimicrobial and antiviral activity using the microtitre broth micro-dilution assay against three Gram-positive bacteria and one Gram-negative bacteria responsible for infections of the upper respiratory tract, as well as cell culture experiments for two different strains of influenza virus.</p> <p>Methods</p> <p>The antimicrobial activity of the elderberry extract was determined by bacterial growth experiments in liquid cultures using the extract at concentrations of 5%, 10%, 15% and 20%. The inhibitory effects were determined by plating the bacteria on agar plates. In addition, the inhibitory potential of the extract on the propagation of human pathogenic H5N1-type influenza A virus isolated from a patient and an influenza B virus strain was investigated using MTT and focus assays.</p> <p>Results</p> <p>For the first time, it was shown that a standardized elderberry liquid extract possesses antimicrobial activity against both Gram-positive bacteria of <it>Streptococcus pyogenes </it>and group C and G <it>Streptococci</it>, and the Gram-negative bacterium <it>Branhamella catarrhalis </it>in liquid cultures. The liquid extract also displays an inhibitory effect on the propagation of human pathogenic influenza viruses.</p> <p>Conclusion</p> <p>Rubini elderberry liquid extract is active against human pathogenic bacteria as well as influenza viruses. The activities shown suggest that additional and alternative approaches to combat infections might be provided by this natural product.</p