Christian Rödelsperger

MutationTaster evaluates disease-causing potential of sequence alterations

(simple_aae) or at alterations causing complex changes in the amino acid sequence (complex_aae). To train the classifier, we generated a dataset with all available and suitable common polymorphisms and known diseasecausing mutations extracted from common databases and the literature. We cross-validated the classifier five times including all three prediction models and obtained an overall accuracy of 91.1 ± 0.1%. We also compared MutationTaster with similar applications (Panther3, Pmut4, PolyPhen and PolyPhen-2 (ref. 5) and ‘screening for non-acceptable polymorphisms’ (SNAP)6) and analyzed the identical 1,000 disease-linked mutations and 1,000 polymorphisms with all programs. For this comparison, we used only alterations causing single amino acid exchanges. MutationTaster performed best in terms of accuracy and speed (Table 1). A description of all training and validation procedures and detailed statistics are available in Supplementary Methods. MutationTaster can be used via an intuitive web interface to analyze single mutations as well as in batch mode. To streamline and to standardize the analysis of NGS data, we provide Perl scripts that can process data from all major platforms (Roche 454, Illumina Genome Analyzer and ABI SOLiD). MutationTaster hence allows the efficient filtering of NGS data for alterations with high disease-causing potential (see Supplementary Methods for an example). Present limitations of the software comprise its inability to analyze insertion-deletions greater than 12 base pairs and alterations spanning an intron-exon border. Also, analysis of non-exonic alterations is restricted to Kozak consensus sequence, splice sites and poly(A) signal. We will add tests for other sequence motifs in the near future. MutationTaster is available at http://www.mutationtaster.org/.

Identity-by-descent filtering of exome sequence data identifies PIGV mutations in hyperphosphatasia mental retardation syndrome

Hyperphosphatasia mental retardation (HPMR) syndrome is an autosomal recessive form of mental retardation with distinct facial features and elevated serum alkaline phosphatase. We performed whole-exome sequencing in three siblings of a nonconsanguineous union with HPMR and performed computational inference of regions identical by descent in all siblings to establish PIGV, encoding a member of the GPI-anchor biosynthesis pathway, as the gene mutated in HPMR. We identified homozygous or compound heterozygous mutations in PIGV in three additional families.

Microindel detection in short-read sequence data

MOTIVATION Several recent studies have demonstrated the effectiveness of resequencing and single nucleotide variant (SNV) detection by deep short-read sequencing platforms. While several reliable algorithms are available for automated SNV detection, the automated detection of microindels in deep short-read data presents a new bioinformatics challenge. RESULTS We systematically analyzed how the short-read mapping tools MAQ, Bowtie, Burrows-Wheeler alignment tool (BWA), Novoalign and RazerS perform on simulated datasets that contain indels and evaluated how indels affect error rates in SNV detection. We implemented a simple algorithm to compute the equivalent indel region eir, which can be used to process the alignments produced by the mapping tools in order to perform indel calling. Using simulated data that contains indels, we demonstrate that indel detection works well on short-read data: the detection rate for microindels (<4 bp) is >90%. Our study provides insights into systematic errors in SNV detection that is based on ungapped short sequence read alignments. Gapped alignments of short sequence reads can be used to reduce this error and to detect microindels in simulated short-read data. A comparison with microindels automatically identified on the ABI Sanger and Roche 454 platform indicates that microindel detection from short sequence reads identifies both overlapping and distinct indels. CONTACT peter.krawitz@googlemail.com; peter.robinson@charite.de SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

A Developmental Switch Coupled to the Evolution of Plasticity Acts through a Sulfatase

Developmental plasticity has been suggested to facilitate phenotypic diversity, but the molecular mechanisms underlying this relationship are little understood. We analyzed a feeding dimorphism in Pristionchus nematodes whereby one of two alternative adult mouth forms is executed after an irreversible developmental decision. By integrating developmental genetics with functional tests in phenotypically divergent populations and species, we identified a regulator of plasticity, eud-1, that acts in a developmental switch. eud-1 mutations eliminate one mouth form, whereas overexpression of eud-1 fixes it. EUD-1 is a sulfatase that acts dosage dependently, is necessary and sufficient to control the sexual dimorphism of feeding forms, and has a conserved function in Pristionchus evolution. It is epistatic to known signaling cascades and results from lineage-specific gene duplications. EUD-1 thus executes a developmental switch for morphological plasticity in the adult stage, showing that regulatory pathways can evolve by terminal addition of new genes.

TCR Repertoire Analysis by Next Generation Sequencing Allows Complex Differential Diagnosis of T Cell–Related Pathology

Clonotype analysis is essential for complete characterization of antigen‐specific T cells. Moreover, knowledge on clonal identity allows tracking of antigen‐specific T cells in whole blood and tissue infiltrates and can provide information on antigenic specificity. Here, we developed a next generation sequencing (NGS)‐based platform for the highly quantitative clonotype characterization of T cells and determined requirements for the unbiased characterization of the input material (DNA, RNA, ex vivo derived or cell culture expanded T cells). Thereafter we performed T cell receptor (TCR) repertoire analysis of various specimens in clinical settings including cytomegalovirus (CMV), polyomavirus BK (BKV) reactivation and acute cellular allograft rejection. Our results revealed dynamic nature of virus‐specific T cell clonotypes; CMV reactivation was linked to appearance of new highly abundant antigen‐specific clonalities. Moreover, analysis of clonotype overlap between BKV‐, alloantigen‐specific T cell–, kidney allograft‐ and urine‐derived lymphocytes provided hints for the differential diagnosis of allograft dysfunction and enabled appropriate therapy adjustment. We believe that the established approach will provide insights into the regulation of virus‐specific/anti‐tumor immunity and has high diagnostic potential in the clinical routine.

Characterization of genetic diversity in the nematode pristionchus pacificus from population-scale resequencing data

The hermaphroditic nematode Pristionchus pacificus is an established model system for comparative studies with Caenorhabditis elegans in developmental biology, ecology, and population genetics. In this study, we present whole-genome sequencing data of 104 P. pacificus strains and the draft assembly of the obligate outcrossing sister species P. exspectatus. We characterize genetic diversity within P. pacificus and investigate the population genetic processes shaping this diversity. P. pacificus is 10 times more diverse than C. elegans and exhibits substantial population structure that allows us to probe its evolution on multiple timescales. Consistent with reduced effective recombination in this self-fertilizing species, we find haplotype blocks that span several megabases. Using the P. exspectatus genome as an outgroup, we polarized variation in P. pacificus and found a site frequency spectrum (SFS) that decays more rapidly than expected in neutral models. The SFS at putatively neutral sites is U shaped, which is a characteristic feature of pervasive linked selection. Based on the additional findings (i) that the majority of nonsynonymous variation is eliminated over timescales on the order of the separation between clades, (ii) that diversity is reduced in gene-rich regions, and (iii) that highly differentiated clades show very similar patterns of diversity, we conclude that purifying selection on many mutations with weak effects is a major force shaping genetic diversity in P. pacificus.

CYNTENATOR: Progressive Gene Order Alignment of 17 Vertebrate Genomes

Whole genome gene order evolution in higher eukaryotes was initially considered as a random process. Gene order conservation or conserved synteny was seen as a feature of common descent and did not imply the existence of functional constraints. This view had to be revised in the light of results from sequencing dozens of vertebrate genomes. It became apparent that other factors exist that constrain gene order in some genomic regions over long evolutionary time periods. Outside of these regions, genomes diverge more rapidly in terms of gene content and order. We have developed CYNTENATOR, a progressive gene order alignment software, to identify genomic regions of conserved synteny over a large set of diverging species. CYNTENATOR does not depend on nucleotide-level alignments and a priori homology assignment. Our software implements an improved scoring function that utilizes the underlying phylogeny. In this manuscript, we report on our progressive gene order alignment approach, a and give a comparison to previous software and an analysis of 17 vertebrate genomes for conservation in gene order. CYNTENATOR has a runtime complexity of and a space complexity of with being the gene number in a genome. CYNTENATOR performs as good as state-of-the-art software on simulated pairwise gene order comparisons, but is the only algorithm that works in practice for aligning dozens of vertebrate-sized gene orders. Lineage-specific characterization of gene order across 17 vertebrate genomes revealed mechanisms for maintaining conserved synteny such as enhancers and coregulation by bidirectional promoters. Genes outside conserved synteny blocks show enrichments for genes involved in responses to external stimuli, stimuli such as immunity and olfactory response in primate genome comparisons. We even see significant gene ontology term enrichments for breakpoint regions of ancestral nodes close to the root of the phylogeny. Additionally, our analysis of transposable elements has revealed a significant accumulation of LINE-1 elements in mammalian breakpoint regions. In summary, CYNTENATOR is a flexible and scalable tool for the identification of conserved gene orders across multiple species over long evolutionary distances.

Neurofibromin (Nf1) is required for skeletal muscle development

Neurofibromatosis type 1 (NF1) is a multi-system disease caused by mutations in the NF1 gene encoding a Ras-GAP protein, neurofibromin, which negatively regulates Ras signaling. Besides neuroectodermal malformations and tumors, the skeletal system is often affected (e.g. scoliosis and long bone dysplasia) demonstrating the importance of neurofibromin for development and maintenance of the musculoskeletal system. Here, we focus on the role of neurofibromin in skeletal muscle development. Nf1 gene inactivation in the early limb bud mesenchyme using Prx1-cre (Nf1(Prx1)) resulted in muscle dystrophy characterized by fibrosis, reduced number of muscle fibers and reduced muscle force. This was caused by an early defect in myogenesis affecting the terminal differentiation of myoblasts between E12.5 and E14.5. In parallel, the muscle connective tissue cells exhibited increased proliferation at E14.5 and an increase in the amount of connective tissue as early as E16.5. These changes were accompanied by excessive mitogen-activated protein kinase pathway activation. Satellite cells isolated from Nf1(Prx1) mice showed normal self-renewal, but their differentiation was impaired as indicated by diminished myotube formation. Our results demonstrate a requirement of neurofibromin for muscle formation and maintenance. This previously unrecognized function of neurofibromin may contribute to the musculoskeletal problems in NF1 patients.

Computational archaeology of the Pristionchus pacificus genome reveals evidence of horizontal gene transfers from insects

BackgroundThe recent sequencing of nematode genomes has laid the basis for comparative genomics approaches to study the impact of horizontal gene transfer (HGT) on the adaptation to new environments and the evolution of parasitism. In the beetle associated nematode Pristionchus pacificus HGT events were found to involve cellulase genes of microbial origin and Diapausin genes that are known from beetles, but not from other nematodes. The insect-to-nematode horizontal transfer is of special interest given that P. pacificus shows a tight association with insects.ResultsIn this study we utilized the observation that horizontally transferred genes often exhibit codon usage patterns more similar to that of the donor than that of the acceptor genome. We introduced GC-normalized relative codon frequencies as a measure to detect characteristic features of P. pacificus orphan genes that show no homology to other nematode genes. We found that atypical codon usage is particularly prevalent in P. pacificus orphans. By comparing codon usage profiles of 71 species, we detected the most significant enrichment in insect-like codon usage profiles. In cross-species comparisons, we identified 509 HGT candidates that show a significantly higher similarity to insect-like profiles than genes with nematode homologs. The most abundant gene family among these genes are non-LTR retrotransposons. Speculating that retrotransposons might have served as carriers of foreign genetic material, we found a significant local clustering tendency of orphan genes in the vicinity of retrotransposons.ConclusionsOur study combined codon usage bias, phylogenetic analysis, and genomic colocalization into a general picture of the computational archaeology of the P. pacificus genome and suggests that a substantial fraction of the gene repertoire is of insect origin. We propose that the Pristionchus-beetle association has facilitated HGT and discuss potential vectors of these events.

Large-scale diversification without genetic isolation in nematode symbionts of figs

A polyphenism discovered to consist of five discrete morphotypes shows that extreme diversification is possible within single species. Diversification is commonly understood to be the divergence of phenotypes accompanying that of lineages. In contrast, alternative phenotypes arising from a single genotype are almost exclusively limited to dimorphism in nature. We report a remarkable case of macroevolutionary-scale diversification without genetic divergence. Upon colonizing the island-like microecosystem of individual figs, symbiotic nematodes of the genus Pristionchus accumulated a polyphenism with up to five discrete adult morphotypes per species. By integrating laboratory and field experiments with extensive genotyping of individuals, including the analysis of 49 genomes from a single species, we show that rapid filling of potential ecological niches is possible without diversifying selection on genotypes. This uncoupling of morphological diversification and speciation in fig-associated nematodes has resulted from a remarkable expansion of discontinuous developmental plasticity.