We describe a new computer program, SnpEff, for rapidly categorizing the effects of variants in genome sequences. Once a genome is sequenced, SnpEff annotates variants based on their genomic locations and predicts coding effects. Annotated genomic locations include intronic, untranslated region, upstream, downstream, splice site, or intergenic regions. Coding effects such as synonymous or non-synonymous amino acid replacement, start codon gains or losses, stop codon gains or losses, or frame shifts can be predicted. Here the use of SnpEff is illustrated by annotating ~356,660 candidate SNPs in ~117 Mb unique sequences, representing a substitution rate of ~1/305 nucleotides, between the Drosophila melanogaster w1118; iso-2; iso-3 strain and the reference y1; cn1 bw1 sp1 strain. We show that ~15,842 SNPs are synonymous and ~4,467 SNPs are non-synonymous (N/S ~0.28). The remaining SNPs are in other categories, such as stop codon gains (38 SNPs), stop codon losses (8 SNPs), and start codon gains (297 SNPs) in...

A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3

Tight control of cell-cell communication is essential for the generation of a normally patterned embryo. A critical mediator of key cell-cell signaling events during embryogenesis is the highly conserved Wnt family of secreted proteins. Recent biochemical and genetic analyses have greatly enriched our understanding of how Wnts signal, and the list of canonical Wnt signaling components has exploded. The data reveal that multiple extracellular, cytoplasmic, and nuclear regulators intricately modulate Wnt signaling levels. In addition, receptor-ligand specificity and feedback loops help to determine Wnt signaling outputs. Wnts are required for adult tissue maintenance, and perturbations in Wnt signaling promote both human degenerative diseases and cancer. The next few years are likely to see novel therapeutic reagents aimed at controlling Wnt signaling in order to alleviate these conditions.

/pdf/the-wnt-signaling-pathway-in-development-and-disease-z4j8s29ao8.pdf

The Wnt signaling pathway in development and disease.

Ribosomal DNA (rDNA) sequences have been aligned and compared in a number of living organisms, and this approach has provided a wealth of information about phylogenetic relationships. Studies of rDNA sequences have been used to infer phylogenetic history across a very broad spectrum, from studies among the basal lineages of life to relationships among closely related species and populations. The reasons for the systematic versatility of rDNA include the numerous rates of evolution among different regions of rDNA (both among and within genes), the presence of many copies of most rDNA sequences per genome, and the pattern of concerted evolution that occurs among repeated copies. These features facilitate the analysis of rDNA by direct RNA sequencing, DNA sequencing (either by cloning or amplification), and restriction enzyme methodologies. Constraints imposed by secondary structure of rRNA and concerted evolution need to be considered in phylogenetic analyses, but these constraints do not appear to impede seriously the usefulness of rDNA. An analysis of aligned sequences of the four nuclear and two mitochondrial rRNA genes identified regions of these genes that are likely to be useful to address phylogenetic problems over a wide range of levels of divergence. In general, the small subunit nuclear sequences appear to be best for elucidating Precambrian divergences, the large subunit nuclear sequences for Paleozoic and Mesozoic divergences, and the organellar sequences of both subunits for Cenozoic divergences. Primer sequences were designed for use in amplifying the entire nuclear rDNA array in 15 sections by use of the polymerase chain reaction; these "universal" primers complement previously described primers for the mitochondrial rRNA genes. Pairs of primers can be selected in conjunction with the analysis of divergence of the rRNA genes to address systematic problems throughout the hierarchy of life.

Ribosomal DNA: molecular evolution and phylogenetic inference.

This revision of the classification of unicellular eukaryotes updates that of Levine et al. (1980) for the protozoa and expands it to include other protists. Whereas the previous revision was primarily to incorporate the results of ultrastructural studies, this revision incorporates results from both ultrastructural research since 1980 and molecular phylogenetic studies. We propose a scheme that is based on nameless ranked systematics. The vocabulary of the taxonomy is updated, particularly to clarify the naming of groups that have been repositioned. We recognize six clusters of eukaryotes that may represent the basic groupings similar to traditional ''kingdoms.'' The multicellular lineages emerged from within monophyletic protist lineages: animals and fungi from Opisthokonta, plants from Archaeplastida, and brown algae from Stramenopiles.

/pdf/the-new-higher-level-classification-of-eukaryotes-with-12im7rq9o7.pdf

The new higher level classification of eukaryotes with emphasis on the taxonomy of protists

The first chordates appear in the fossil record at the time of the Cambrian explosion, nearly 550 million years ago. The modern ascidian tadpole represents a plausible approximation to these ancestral chordates. To illuminate the origins of chordate and vertebrates, we generated a draft of the protein-coding portion of the genome of the most studied ascidian, Ciona intestinalis. The Ciona genome contains approximately 16,000 protein-coding genes, similar to the number in other invertebrates, but only half that found in vertebrates. Vertebrate gene families are typically found in simplified form in Ciona, suggesting that ascidians contain the basic ancestral complement of genes involved in cell signaling and development. The ascidian genome has also acquired a number of lineage-specific innovations, including a group of genes engaged in cellulose metabolism that are related to those in bacteria and fungi.

https://science.sciencemag.org/content/sci/298/5601/2157.full.pdf

The draft genome of Ciona intestinalis : insights into chordate and vertebrate origins

Understanding the early evolution of animal body plans requires knowledge both of metazoan phylogeny and of the genetic and developmental changes involved in the emergence of particular forms. Recent 18S ribosomal RNA phylogenies suggest a three-branched tree for the Bilateria comprising the deuterostomes and two great protostome clades, the lophotrochozoans1 and ecdysozoans2. Here, we show that the complement of Hox genes in critical protostome phyla reflects these phylogenetic relationships and reveals the early evolution of developmental regulatory potential in bilaterians. We have identified Hox genes that are shared by subsets of protostome phyla. These include a diverged pair of posterior (Abdominal-B -like) genes in both a brachiopod and a polychaete annelid, which supports the lophotrochozoan assemblage, and a distinct posterior Hox gene shared by a priapulid, a nematode and the arthropods, which supports the ecdysozoan clade. The ancestors of each of these two major protostome lineages had a minimum of eight to ten Hox genes. The major period of Hox gene expansion and diversification thus occurred before the radiation of each of the three great bilaterian clades.

Hox genes in brachiopods and priapulids and protostome evolution

DNA sequence analysis dictates new interpretation of phylogenic trees. Taxa that were once thought to represent successive grades of complexity at the base of the metazoan tree are being displaced to much higher positions inside the tree. This leaves no evolutionary "intermediates" and forces us to rethink the genesis of bilaterian complexity.

https://www.pnas.org/content/pnas/97/9/4453.full.pdf

The new animal phylogeny: Reliability and implications

A posttranslational modification was detected in the carboxyl-terminal region of axonemal tubulin from Paramecium. Tubulin carboxyl-terminal peptides were isolated and analyzed by Edman degradation sequencing, mass spectrometry, and amino acid analysis. All of the peptides, derived from both alpha and beta tubulin subunits, were modified by polyglycylation, containing up to 34 glycyl units covalently bound to the gamma carboxyl group of glutamyl residues. This modification, present in one of the most stable microtubular systems, may influence microtubule stability or axoneme function, or both.

https://science.sciencemag.org/content/sci/266/5191/1688.full.pdf

Polyglycylation of tubulin: a posttranslational modification in axonemal microtubules

Most of the major invertebrate phyla appear in the fossil record during a relatively short time interval, not exceeding 20 million years (Myr), 540-520 Myr ago. This rapid diversification is known as the `Cambrian explosion9. In the present paper, we ask whether molecular phylogenetic reconstruction provides confirmation for such an evolutionary burst. The expectation is that the molecular phylogenetic trees should take the form of a large unresolved multifurcation of the various animal lineages. Complete 18S rRNA sequences of 69 extant representatives of 15 animal phyla were obtained from data banks. After eliminating a major source of artefact leading to lack of resolution in phylogenetic trees (mutational saturation of sequences), we indeed observe that the major lines of triploblast coelomates (arthropods, molluscs, echinoderms, chordates...) are very poorly resolved i.e. the nodes defining the various clades are not supported by high bootstrap values. Using a previously developed procedure consisting of calculating bootstrap proportions of each node of the tree as a function of increasing amount of nucleotides (Lecointre, G., Philippe, H. Le, H. L. V. and Le Guyader, H. (1994) Mol. Phyl. Evol. , in press) we obtain a more informative indication of the robustness of each node. In addition, this procedure allows us to estimate the number of additional nucleotides that would be required to resolve confidently the currently uncertain nodes; this number turns out to be extremely high and experimentally unfeasible. We then take this approach one step further: using parameters derived from the above analysis, assuming a molecular clock and using palaeontological dates for calibration, we establish a relationship between the number of sites contained in a given data set and the time interval that this data set can confidently resolve (with 95% bootstrap support). Under these assumptions, the presently available 18S rRNA database cannot confidently resolve cladogenetic events separated by less than about 40 Myr. Thus, at the present time, the potential resolution by the palaeontological approach is higher than that by the molecular one.

André Adoutte

Papers

Hox genes in brachiopods and priapulids and protostome evolution

The new animal phylogeny: Reliability and implications

Polyglycylation of tubulin: a posttranslational modification in axonemal microtubules

Can the Cambrian explosion be inferred through molecular phylogeny

Animal evolution. The end of the intermediate taxa