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]

FJRW-rings and Mirror Symmetry

We verify the Landau-Ginzburg Mirror Symmetry Conjecture for Arnol'd's list of unimodal and bimodal quasi-homogeneous singularities with G the maximal diagonal symmetry group, and include a discussion of eight axioms which facilitate the computation of FJRW-rings.Comment: 29 pages, 9 page introduction followed by computations. This replacement corrects a small mistake on page 25 in the W-degree tabl

Integrative analysis of key candidate genes and signaling pathways in autoimmune thyroid dysfunction related to anti-CTLA-4 therapy by bioinformatics

Cytotoxic T lymphocyte-associated antigen-4 (CTLA-4), the first immune checkpoint to be targeted clinically, has provided an effective treatment option for various malignancies. However, the clinical advantages associated with CTLA-4 inhibitors can be offset by the potentially severe immune-related adverse events (IRAEs), including autoimmune thyroid dysfunction. To investigate the candidate genes and signaling pathways involving in autoimmune thyroid dysfunction related to anti-CTLA-4 therapy, integrated differentially expressed genes (DEGs) were extracted from the intersection of genes from Gene Expression Omnibus (GEO) datasets and text mining. The functional enrichment was performed by gene ontology (GO) annotation and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis. Protein-protein interaction (PPI) network, module enrichment, and hub gene identification were constructed and visualized by the online Search Tool for the Retrieval of Interacting Genes (STRING) and Cytoscape software. A total of 22 and 17 integrated human DEGs in hypothyroidism and hyperthyroidism group related to anti-CTLA-4 therapy were identified, respectively. Functional enrichment analysis revealed 24 GO terms and 1 KEGG pathways in the hypothyroid group and 21 GO terms and 2 KEGG pathways in the hyperthyroid group. After PPI network construction, the top five hub genes associated with hypothyroidism were extracted, including ALB, MAPK1, SPP1, PPARG, and MIF, whereas those associated with hyperthyroidism were ALB, FCGR2B, CD44, LCN2, and CD74. The identification of the candidate key genes and enriched signaling pathways provides potential biomarkers for autoimmune thyroid dysfunction related to anti-CTLA-4 therapy and might contribute to the future diagnosis and management of IRAEs for cancer patients

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

Whole-exome re-sequencing in a family quartet identifies POP1 mutations as the cause of a novel skeletal dysplasia

Recent advances in DNA sequencing have enabled mapping of genes for monogenic traits in families with small pedigrees and even in unrelated cases. We report the identification of disease-causing mutations in a rare, severe, skeletal dysplasia, studying a family of two healthy unrelated parents and two affected children using whole-exome sequencing. The two affected daughters have clinical and radiographic features suggestive of anauxetic dysplasia (OMIM 607095), a rare form of dwarfism caused by mutations of RMRP. However, mutations of RMRP were excluded in this family by direct sequencing. Our studies identified two novel compound heterozygous loss-of-function mutations in POP1, which encodes a core component of the RNase mitochondrial RNA processing (RNase MRP) complex that directly interacts with the RMRP RNA domains that are affected in anauxetic dysplasia. We demonstrate that these mutations impair the integrity and activity of this complex and that they impair cell proliferation, providing likely molecular and cellular mechanisms by which POP1 mutations cause this severe skeletal dysplasia