Molecular phylogeny and divergence times of deuterostome animals

Comparison of oocyte mRNA localization patterns in sterlet Acipenser ruthenus and African clawed frog Xenopus laevis

In oocytes, RNA localization has critical implications, as assembly of proteins in particular subcellular domains is crucial to embryo development. The distribution of mRNA molecules can identify and characterize localized transcripts. The goal of this study was to clarify the origin of primordial germ cells in the oocyte body plan and to reveal the generation of cell lineages by localized RNAs. The distribution of 12 selected mRNAs in sterlet Acipenser ruthenus oocytes was investigated by qPCR tomography and compared with known patterns of mRNA localization in Xenopus laevis. We investigated the distribution of two gene clusters in the ooplasm along the animal-vegetal axis of the sturgeon oocyte, both of which showed clearly defined intracellular gradient pattern remarkably similar to their distribution in the frog oocyte. We elucidated the localization of sturgeon egg germplasm markers belonging to the vegetal group of mRNAs. The mRNAs coding otx1, wnt11, and veg1 found to be localized in the sturgeon animal hemisphere are, in contrast, distributed in the vegetal hemisphere in amphibian. Actinopterygii and Sarcopterygii, two major lineages of osteichthyan vertebrates, split about 476 Ma (Blair & Hedges, ), albeit basal lineages share conserved biological features. Acipenseriformes is one the most basal living lineages of Actinopterygii, having evolved about 200 Ma (Bemis, Birstein, & Waldman, ), contemporaneous with modern amphibians (Roelants et al., ).

A death in the family: Sea lamprey (Petromyzon marinus) avoidance of confamilial alarm cues diminishes with phylogenetic distance

Alarm signals released after predator attack function as reliable public information revealing areas of high risk. The utility of this information can extend beyond species boundaries, benefiting heterospecifics capable of recognizing and responding appropriately to the signal. Nonmutually exclusive hypotheses explaining the acquisition of heterospecific reactivity to cues suggest it could be conserved phylogenetically following its evolution in a common ancestor (a species‐level effect) and/or learned during periods of shared risk (a population‐level effect; e.g., shared predators). Using a laboratory‐based space‐use behavioral assay, we tested the response of sea lamprey (Petromyzon marinus) to the damage‐released alarm cues of five confamilial (sympatric and allopatric) species and two distantly related out‐groups: a sympatric teleost (white sucker Catostomus commersonii) and an allopatric agnathan (Atlantic hagfish Myxine glutinosa). We found that sea lamprey differed in their response to conspecific and heterospecific odors; exhibiting progressively weaker avoidance of cues derived from more phylogenetically distant confamilials regardless of current overlap in distribution. Odors from out‐groups elicited no response. These findings suggest that a damage‐released alarm cue is at least partially conserved within the Petromyzontidae and that sea lamprey perceives predator attacks directed to closely related taxa. These findings are consistent with similar observations from gastropod, amphibian and bony fish taxa, and we discuss this in an eco‐evo context to provide a plausible explanation for the acquisition and maintenance of the response in sea lamprey.

Cytotoxic activity of Holothuria tubulosa (Echinodermata) coelomocytes

The immune system of marine invertebrates, in particular that of holothurians, still requires further study. Our research showed that coelomocyte cells contained in the coelomic fluid of the sea cucumber, Holothuria tubulosa, are able to lyse, in vitro, red blood cells in rabbits and sheep. A plaque-forming assay showed spherule cells to be the effector cells, able to release cytotoxic molecules after xenogenic cell contact. The coelomocyte lysate supernatant, analysed by polyacrylamide gel electrophoresis overlay technique, using rabbit and sheep erythrocytes, showed two different haemolytic protein patterns: one calcium dependent and the other calcium independent. The fractions of each pattern were resolved on a polyacrylamide gel and calcium-dependent and independent coelomocyte lysate patterns were compared.

The sea cucumber genome provides insights into morphological evolution and visceral regeneration

Apart from sharing common ancestry with chordates, sea cucumbers exhibit a unique morphology and exceptional regenerative capacity. Here we present the complete genome sequence of an economically important sea cucumber, A. japonicus, generated using Illumina and PacBio platforms, to achieve an assembly of approximately 805 Mb (contig N50 of 190 Kb and scaffold N50 of 486 Kb), with 30,350 protein-coding genes and high continuity. We used this resource to explore key genetic mechanisms behind the unique biological characters of sea cucumbers. Phylogenetic and comparative genomic analyses revealed the presence of marker genes associated with notochord and gill slits, suggesting that these chordate features were present in ancestral echinoderms. The unique shape and weak mineralization of the sea cucumber adult body were also preliminarily explained by the contraction of biomineralization genes. Genome, transcriptome, and proteome analyses of organ regrowth after induced evisceration provided insight into the molecular underpinnings of visceral regeneration, including a specific tandem-duplicated prostatic secretory protein of 94 amino acids (PSP94)-like gene family and a significantly expanded fibrinogen-related protein (FREP) gene family. This high-quality genome resource will provide a useful framework for future research into biological processes and evolution in deuterostomes, including remarkable regenerative abilities that could have medical applications. Moreover, the multiomics data will be of prime value for commercial sea cucumber breeding programs.

Echinoderms, ubiquitous in the marine environment, are important from evolutionary, ecological, and socioeconomic perspectives. Together with chordates and hemichordates, they form the deuterostome clade, making them a crucial node in the study of chordate ancestry. Within echinoderms, class Holothuroidea is unique; its members (the sea cucumbers) display remarkable regenerative abilities and play key roles as sediment bioturbators and symbiotic hosts, and many are prized in the seafood and pharmaceutical industries. The sea cucumber genome therefore has the potential to significantly contribute to our understanding of important evolutionary and biological processes and help enhance aquaculture programs. Here we present a high-quality genome sequence for the economically important species Apostichopus japonicus. Through comparative analyses, we identified 763 echinoderm-specific gene families enriched in genes encoding membrane proteins, ion channels, and signal transduction proteins. Marker genes associated with the notochord and gill slits were also found, providing valuable insight into the origin of chordates. The reduced number and low expression levels of biomineralization genes reflect the skeletal degeneration seen in sea cucumbers. Importantly, 2 gene families appeared to be expanded in A. japonicus and may play crucial roles in its heightened regenerative potential. Together, findings from the sea cucumber genome provide important and novel insights into echinoderm and deuterostome biology.

An evolutionary perspective on the systems of adaptive immunity

We propose an evolutionary perspective to classify and characterize the diverse systems of adaptive immunity that have been discovered across all major domains of life. We put forward a new function-based classification according to the way information is acquired by the immune systems: Darwinian immunity (currently known from, but not necessarily limited to, vertebrates) relies on the Darwinian process of clonal selection to 'learn' by cumulative trial-and-error feedback; Lamarckian immunity uses templated targeting (guided adaptation) to internalize heritable information on potential threats; finally, shotgun immunity operates through somatic mechanisms of variable targeting without feedback. We argue that the origin of Darwinian (but not Lamarckian or shotgun) immunity represents a radical innovation in the evolution of individuality and complexity, and propose to add it to the list of major evolutionary transitions. While transitions to higher-level units entail the suppression of selection at lower levels, Darwinian immunity re-opens cell-level selection within the multicellular organism, under the control of mechanisms that direct, rather than suppress, cell-level evolution for the benefit of the individual. From a conceptual point of view, the origin of Darwinian immunity can be regarded as the most radical transition in the history of life, in which evolution by natural selection has literally re-invented itself. Furthermore, the combination of clonal selection and somatic receptor diversity enabled a transition from limited to practically unlimited capacity to store information about the antigenic environment. The origin of Darwinian immunity therefore comprises both a transition in individuality and the emergence of a new information system - the two hallmarks of major evolutionary transitions. Finally, we present an evolutionary scenario for the origin of Darwinian immunity in vertebrates. We propose a revival of the concept of the 'Big Bang' of vertebrate immunity, arguing that its origin involved a 'difficult' (i.e. low-probability) evolutionary transition that might have occurred only once, in a common ancestor of all vertebrates. In contrast to the original concept, we argue that the limiting innovation was not the generation of somatic diversity, but the regulatory circuitry needed for the safe operation of amplifiable immune responses with somatically acquired targeting. Regulatory complexity increased abruptly by genomic duplications at the root of the vertebrate lineage, creating a rare opportunity to establish such circuitry. We discuss the selection forces that might have acted at the origin of the transition, and in the subsequent stepwise evolution leading to the modern immune systems of extant vertebrates.

Compartmentalized expression patterns of pancreatic- and gastric-related genes in the alimentary canal of the ascidian Ciona intestinalis: evolutionary insights into the functional regionality of the gastrointestinal tract in Olfactores

Many heterotrophic animals have a one-way alimentary canal that is essential for their nutrition and sequential steps of the digestive system, namely ingestion, digestion, absorption and elimination, are widely shared among bilaterians. Morphological, functional and molecular knowledge of the alimentary canal has been obtained in particular from mammalian research but the shared features and evolution of these aspects of the highly diverged alimentary canal in the animal kingdom are still unclear. We therefore investigate spatial gene expression patterns of pancreatic- and gastric-related molecules of ascidians (a sister group of vertebrates) with special reference to the functional regionality of the gastrointestinal tract. Genome-wide surveys of ascidian homologs to mammalian exocrine digestive enzyme genes revealed that pancreatic enzymes, namely alpha-amylase, lipase, phospholipase A2, trypsin, chymotrypsin and carboxypeptidase, exist in the ascidian genome. However, an ascidian homolog of the mammalian gastric enzyme pepsin has not been identified, although molecules resembling cathepsin D, a pepsin relative, are indeed present. Spatial expression analyses in the ascidian Ciona intestinalis, by means of whole-mount in situ hybridization, have elucidated that the expression of Ciona homologs of pancreatic- and gastric-related exocrine enzyme genes and of their transcriptional regulator genes is restricted to the Ciona stomach. Furthermore, the expression of these genes is localized to specific regions of the stomach epithelium according to their regionality in the vertebrate digestive system. The compartmentalized expression patterns of Ciona homologs imply primitive and/or ancestral aspects of molecular, functional and morphological bases among Olfactores.

The Tetraspanin-Associated Uroplakins Family (UPK2/3) Is Evolutionarily Related to PTPRQ, a Phosphotyrosine Phosphatase Receptor

Uroplakins are a widespread group of vertebrate integral membrane proteins that belong to two different families: UPK1a and UPK1b belong to the large tetraspanin (TSPAN) gene family, and UPK3a, UPK3b, UPK3c, UPK3d, UPK2a and UPK2b form a family of their own, the UPK2/3 tetraspanin-associated family. In a previous study, we reported that uroplakins first appeared in vertebrates, and that uroplakin tetraspanins (UPK1a and UPK1b) should have originated by duplication of an ancestor tetraspanin gene. However, the evolutionary origin of the UPK2/3 family remains unclear. In this study, we provide evidence that the UPK2/3 family originated by gene duplication and domain loss from a protoPTPRQ-like basal deuterostome gene. PTPRQs are members of the subtype R3 tyrosine phosphatase receptor (R3 PTPR) family, which are characterized by having a unique modular composition of extracellular fibronectin (FN3) repeats, a transmembrane helix, and a single intra-cytoplasmic phosphotyrosine phophatase (PTP) domain. Our assumption of a deuterostome protoPTPRQ-like gene as an ancestor of the UPK2/3 family by gene duplication and loss of its PTP and fibronectin (FN3) domains, excluding the one closest to the transmembrane helix, is based on the following: (i) phylogenetic analyses, (ii) the existence of an identical intron/exon gene pattern between UPK2/3 and the corresponding genetic region in R3 PTPRs, (iii) the conservation of cysteine patterns and protein motifs between UPK2/3 and PTPRQ proteins and, (iv) the existence in tunicates, the closest organisms to vertebrates, of two sequences related to PTPRQ; one with the full subtype R3 modular characteristic and another without the PTP domain but with a short cytoplasmic tail with some sequence similarity to that of UPK3a. This finding will facilitate further studies on the structure and function of these important proteins with implications in human diseases.

Molecular evolution of troponin I and a role of its N-terminal extension in nematode locomotion

The troponin complex, composed of troponin T (TnT), troponin I (TnI), and troponin C (TnC), is the major calcium-dependent regulator of muscle contraction, which is present widely in both vertebrates and invertebrates. Little is known about evolutionary aspects of troponin in the animal kingdom. Using a combination of data mining and functional analysis of TnI, we report evidence that an N-terminal extension of TnI is present in most of bilaterian animals as a functionally important domain. Troponin components have been reported in species in most of representative bilaterian phyla. Comparison of TnI sequences shows that the core domains are conserved in all examined TnIs, and that N- and C-terminal extensions are variable among isoforms and species. In particular, N-terminal extensions are present in all protostome TnIs and chordate cardiac TnIs but lost in a subset of chordate TnIs including vertebrate skeletal-muscle isoforms. Transgenic rescue experiments in Caenorhabditis elegans striated muscle show that the N-terminal extension of TnI (UNC-27) is required for coordinated worm locomotion but not in sarcomere assembly and single muscle-contractility kinetics. These results suggest that N-terminal extensions of TnIs are retained from a TnI ancestor as a functional domain.

Structural Biology and Evolution of the TGF-β Family

We review the evolution and structure of members of the transforming growth factor β (TGF-β) family, antagonistic or agonistic modulators, and receptors that regulate TGF-β signaling in extracellular environments. The growth factor (GF) domain common to all family members and many of their antagonists evolved from a common cystine knot growth factor (CKGF) domain. The CKGF superfamily comprises six distinct families in primitive metazoans, including the TGF-β and Dan families. Compared with Wnt/Frizzled and Notch/Delta families that also specify body axes, cell fate, tissues, and other families that contain CKGF domains that evolved in parallel, the TGF-β family was the most fruitful in evolution. Complexes between the prodomains and GFs of the TGF-β family suggest a new paradigm for regulating GF release by conversion from closed- to open-arm procomplex conformations. Ternary complexes of the final step in extracellular signaling show how TGF-β GF dimers bind type I and type II receptors on the cell surface, and enable understanding of much of the specificity and promiscuity in extracellular signaling. However, structures suggest that when GFs bind repulsive guidance molecule (RGM) family coreceptors, type I receptors do not bind until reaching an intracellular, membrane-enveloped compartment, blurring the line between extra- and intracellular signaling. Modulator protein structures show how structurally diverse antagonists including follistatins, noggin, and members of the chordin family bind GFs to regulate signaling; complexes with the Dan family remain elusive. Much work is needed to understand how these molecular components assemble to form signaling hubs in extracellular environments in vivo.

Evolution of the Rax family of developmental transcription factors in vertebrates

Rax proteins comprise a small family of paired-type, homeodomain-containing transcription factors with essential functions in eye and forebrain development. While invertebrates possess only one Rax gene, vertebrates can have several Rax paralogue genes, but the evolutionary history of the members of the family has not been studied in detail. Here, we present a thorough analysis of the evolutionary relationships between vertebrate Rax genes and proteins available in diverse genomic databases. Phylogenetic and synteny analyses indicate that Rax genes went through a duplication in an ancestor of all jawed vertebrates (Gnathostomata), giving rise to the ancestral vertebrate Rax1 and Rax2 genes. This duplication event is likely related to the proposed polyploidisations that occurred during early vertebrate evolution. Subsequent genome-wide duplications in the lineage of ray-finned fish (Actinopterygii) originated new Rax2 paralogues in the genomes of teleosts. In the lobe-finned fish lineage (Sarcopterygii), the N-terminal octapeptide domain of Rax2 was lost in a common ancestor of tetrapods, giving rise to a shorter version of Rax2 in this lineage. Within placental mammals, the Rax2 gene was lost altogether in an ancestor of rodents and lagomorphs (Glires). Finally, we discuss the scientific literature in the light of Rax gene evolution and propose new avenues of research on the function of this important family of transcriptional regulators.

Physiology and Evolution of Voltage-Gated Calcium Channels in Early Diverging Animal Phyla: Cnidaria, Placozoa, Porifera and Ctenophora

Voltage-gated calcium (Ca_v) channels serve dual roles in the cell, where they can both depolarize the membrane potential for electrical excitability, and activate transient cytoplasmic Ca²⁺ signals. In animals, Ca_v channels play crucial roles including driving muscle contraction (excitation-contraction coupling), gene expression (excitation-transcription coupling), pre-synaptic and neuroendocrine exocytosis (excitation-secretion coupling), regulation of flagellar/ciliary beating, and regulation of cellular excitability, either directly or through modulation of other Ca²⁺-sensitive ion channels. In recent years, genome sequencing has provided significant insights into the molecular evolution of Ca_v channels. Furthermore, expanded gene datasets have permitted improved inference of the species phylogeny at the base of Metazoa, providing clearer insights into the evolution of complex animal traits which involve Ca_v channels, including the nervous system. For the various types of metazoan Ca_v channels, key properties that determine their cellular contribution include: Ion selectivity, pore gating, and, importantly, cytoplasmic protein-protein interactions that direct sub-cellular localization and functional complexing. It is unclear when these defining features, many of which are essential for nervous system function, evolved. In this review, we highlight some experimental observations that implicate Ca_v channels in the physiology and behavior of the most early-diverging animals from the phyla Cnidaria, Placozoa, Porifera, and Ctenophora. Given our limited understanding of the molecular biology of Ca_v channels in these basal animal lineages, we infer insights from better-studied vertebrate and invertebrate animals. We also highlight some apparently conserved cellular functions of Ca_v channels, which might have emerged very early on during metazoan evolution, or perhaps predated it.

The Evolution of Antiviral Defense Systems

Self-replicating genetic material presumably provided the architecture necessary for generating the last universal ancestor of all nucleic-acid-based life. As biological complexity increased in the billions of years that followed, the same genetic material also morphed into a wide spectrum of viruses and other parasitic genetic elements. The resulting struggle for existence drove the evolution of host defenses, giving rise to a perpetual arms race. This Perspective summarizes the antiviral mechanisms evident across the tree of life, discussing each in their evolutionary context to postulate how the coevolution of host and pathogen shaped the cellular antiviral defenses we know today.

Molecular Evolution of Human Coronavirus Genomes

Human coronaviruses (HCoVs), including SARS-CoV and MERS-CoV, are zoonotic pathogens that originated in wild animals. HCoVs have large genomes that encode a fixed array of structural and nonstructural components, as well as a variety of accessory proteins that differ in number and sequence even among closely related CoVs. Thus, in addition to recombination and mutation, HCoV genomes evolve through gene gains and losses. In this review we summarize recent findings on the molecular evolution of HCoV genomes, with special attention to recombination and adaptive events that generated new viral species and contributed to host shifts and to HCoV emergence. VIDEO ABSTRACT.

Enamel ribbons, surface nodules, and octacalcium phosphate in C57BL/6 Amelx-/- mice and Amelx+/- lyonization

BACKGROUND

Amelogenin is required for normal enamel formation and is the most abundant protein in developing enamel.

METHODS

Amelx^+/+, Amelx^+/- , and Amelx^-/- molars and incisors from C57BL/6 mice were characterized using RT-PCR, Western blotting, dissecting and light microscopy, immunohistochemistry (IHC), transmission electron microscopy (TEM), scanning electron microscopy (SEM), backscattered SEM (bSEM), nanohardness testing, and X-ray diffraction.

RESULTS

No amelogenin protein was detected by Western blot analyses of enamel extracts from Amelx^-/- mice. Amelx^-/- incisor enamel averaged 20.3 ± 3.3 μm in thickness, or only 1/6th that of the wild type (122.3 ± 7.9 μm). Amelx^-/- incisor enamel nanohardness was 1.6 Gpa, less than half that of wild-type enamel (3.6 Gpa). Amelx^+/- incisors and molars showed vertical banding patterns unique to each tooth. IHC detected no amelogenin in Amelx^-/- enamel and varied levels of amelogenin in Amelx^+/- incisors, which correlated positively with enamel thickness, strongly supporting lyonization as the cause of the variations in enamel thickness. TEM analyses showed characteristic mineral ribbons in Amelx^+/+ and Amelx^-/- enamel extending from mineralized dentin collagen to the ameloblast. The Amelx^-/- enamel ribbons were not well separated by matrix and appeared to fuse together, forming plates. X-ray diffraction determined that the predominant mineral in Amelx^-/- enamel is octacalcium phosphate (not calcium hydroxyapatite). Amelx^-/- ameloblasts were similar to wild-type ameloblasts except no Tomes' processes extended into the thin enamel. Amelx^-/- and Amelx^+/- molars both showed calcified nodules on their occlusal surfaces. Histology of D5 and D11 developing molars showed nodules forming during the maturation stage.

CONCLUSION

Amelogenin forms a resorbable matrix that separates and supports, but does not shape early secretory-stage enamel ribbons. Amelogenin may facilitate the conversion of enamel ribbons into hydroxyapatite by inhibiting the formation of octacalcium phosphate. Amelogenin is necessary for thickening the enamel layer, which helps maintain ribbon organization and development and maintenance of the Tomes' process.

Novel pentapeptide, PALAL, derived from a bony fish elicits contraction of the muscle in starfish Patiria pectinifera

A bioactive peptide mimicking peptide-signaling molecules has been isolated from the skin extract of fish Channa argus which caused contraction of the apical muscle of a starfish Patiria pectinifera, a deuterostomian invertebrate. The primary structure of the isolated pentapeptide comprises amino acid sequence of H-Pro-Ala-Leu-Ala-Leu-OH (PALAL) with a molecular mass of 483.7 Da. Pharmacological activity of PALAL, dosage ranging from 10^-9 to 10^-5 M, revealed concentration-dependent contraction of the apical muscles of P. pectinifera and Asterias amurensis. However, PALAL was not active on the intestinal smooth muscle of the goldfish Carassius auratus and has presumably other physiological roles in fish skin. Investigation of structure-activity relationship using truncated and substituted analogs of PALAL demonstrated that H-Ala-Leu-Ala-Leu-OH was necessary and should be sufficient to constrict apical muscle of P. pectinifera. Furthermore, the second alanine residue was required to display the activity, and the fifth leucine residue was responsible for its potency. Comparison with PALAL's primary structure with those of other known bioactive peptides from fish and starfish revealed that PALAL does not have any significant homology. Consequently, PALAL is a bioactive peptide that elicits a muscle contraction in starfish, and the isolation of PALAL may lead to develop other bioactive peptides sharing its similar sequence and/or activity. Copyright © 2016 European Peptide Society and John Wiley & Sons, Ltd.

Origin of a major infectious disease in vertebrates: The timing of Cryptosporidium evolution and its hosts

Protozoan parasites of the genus Cryptosporidium infect all vertebrate groups and display some host specificity in their infections. It is therefore possible to assume that Cryptosporidium parasites evolved intimately aside with vertebrate lineages. Here we propose a scenario of Cryptosporidium–Vertebrata coevolution testing the hypothesis that the origin of Cryptosporidium parasites follows that of the origin of modern vertebrates. We use calibrated molecular clocks and cophylogeny analyses to provide and compare age estimates and patterns of association between these clades. Our study provides strong support for the evolution of parasitism of Cryptosporidium with the rise of the vertebrates about 600 million years ago (Mya). Interestingly, periods of increased diversification in Cryptosporidium coincides with diversification of crown mammalian and avian orders after the Cretaceous-Palaeogene (K-Pg) boundary, suggesting that adaptive radiation to new mammalian and avian hosts triggered the diversification of this parasite lineage. Despite evidence for ongoing host shifts we also found significant correlation between protozoan parasites and vertebrate hosts trees in the cophylogenetic analysis. These results help us to understand the underlying macroevolutionary mechanisms driving evolution in Cryptosporidium and may have important implications for the ecology, dynamics and epidemiology of cryptosporidiosis disease in humans and other animals.

Alternative haplotypes of antigen processing genes in zebrafish diverged early in vertebrate evolution

Antigen processing and presentation genes found within the MHC are among the most highly polymorphic genes of vertebrate genomes, providing populations with diverse immune responses to a wide array of pathogens. Here, we describe transcriptome, exome, and whole-genome sequencing of clonal zebrafish, uncovering the most extensive diversity within the antigen processing and presentation genes of any species yet examined. Our CG2 clonal zebrafish assembly provides genomic context within a remarkably divergent haplotype of the core MHC region on chromosome 19 for six expressed genes not found in the zebrafish reference genome: mhc1uga, proteasome-β 9b (psmb9b), psmb8f, and previously unknown genes psmb13b, tap2d, and tap2e We identify ancient lineages for Psmb13 within a proteasome branch previously thought to be monomorphic and provide evidence of substantial lineage diversity within each of three major trifurcations of catalytic-type proteasome subunits in vertebrates: Psmb5/Psmb8/Psmb11, Psmb6/Psmb9/Psmb12, and Psmb7/Psmb10/Psmb13. Strikingly, nearby tap2 and MHC class I genes also retain ancient sequence lineages, indicating that alternative lineages may have been preserved throughout the entire MHC pathway since early diversification of the adaptive immune system ∼500 Mya. Furthermore, polymorphisms within the three MHC pathway steps (antigen cleavage, transport, and presentation) are each predicted to alter peptide specificity. Lastly, comparative analysis shows that antigen processing gene diversity is far more extensive than previously realized (with ancient coelacanth psmb8 lineages, shark psmb13, and tap2t and psmb10 outside the teleost MHC), implying distinct immune functions and conserved roles in shaping MHC pathway evolution throughout vertebrates.

The Identification of the Closest Living Relative(s) of Tetrapods: Phylogenomic Lessons for Resolving Short Ancient Internodes

Identifying the closest living relative(s) of tetrapods is an important, yet still contested question in vertebrate phylogenetics. Three hypotheses are possible and ruling out alternatives has proven difficult even with large molecular data sets due to weak phylogenetic signal coupled nonphylogenetic noise resulting from relatively rapid speciation events that occurred a long time ago ([Formula: see text]400 Ma). Here, we revisit the identity of the closest living relative of land vertebrates from a phylogenomic perspective and include new genomic data for all extant lungfish genera. RNA-seq proves to be a great alternative to genomic sequencing, which currently is technically not feasible in lungfishes due to their huge (50-130 Gb) and repetitive genomes. We examined the most important sources of systematic error, namely long-branch attraction (LBA), compositional heterogeneity and distribution of missing data and applied different correction techniques. A multispecies coalescent approach is used to account for deep coalescence that might come from the short and deep internodes separating early sarcopterygian splits. Concatenation methods favored lungfishes as the closest living relatives of tetrapods with strong statistical support. Amino acid profile mixture models can unambiguously resolve this difficult internode thanks to their ability to avoid systematic error. We assessed the performance of different site-heterogeneous models and data partitioning and compared the ability of different strategies designed to overcome LBA, including taxon manipulation, reduction of among-lineage rate heterogeneity and removal of fast-evolving or compositionally heterogeneous positions. The identification of lungfish as sister group of tetrapods is robust regarding the effects of nonstationary composition and distribution of missing data. The multispecies coalescent method reconstructed strongly supported topologies that were congruent with concatenation, despite pervasive gene tree heterogeneity. We reject alternative topologies for early sarcopterygian relationships by increasing the signal-to-noise ratio in our alignments. The analytical pipeline outlined here combines probabilistic phylogenomic inference with methods for evaluating data quality, model adequacy, and assessing systematic error, and thus is likely to help resolve similarly difficult internodes in the tree of life. [Coalescence; coelacanth; compositional heterogeneity; gene tree; long-branch attraction; lungfish; missing data; model misspecification; phylogenomic; species tree; systematic error.].

Significant variations in alternative splicing patterns and expression profiles between human-mouse orthologs in early embryos

Human and mouse orthologs are expected to have similar biological functions; however, many discrepancies have also been reported. We systematically compared human and mouse orthologs in terms of alternative splicing patterns and expression profiles. Human-mouse orthologs are divergent in alternative splicing, as human orthologs could generally encode more isoforms than their mouse orthologs. In early embryos, exon skipping is far more common with human orthologs, whereas constitutive exons are more prevalent with mouse orthologs. This may correlate with divergence in expression of splicing regulators. Orthologous expression similarities are different in distinct embryonic stages, with the highest in morula. Expression differences for orthologous transcription factor genes could play an important role in orthologous expression discordance. We further detected largely orthologous divergence in differential expression between distinct embryonic stages. Collectively, our study uncovers significant orthologous divergence from multiple aspects, which may result in functional differences and dynamics between human-mouse orthologs during embryonic development.

Explaining large-scale patterns of vertebrate diversity

The major clades of vertebrates differ dramatically in their current species richness, from 2 to more than 32,000 species each, but the causes of this variation remain poorly understood. For example, a previous study noted that vertebrate clades differ in their diversification rates, but did not explain why they differ. Using a time-calibrated phylogeny and phylogenetic comparative methods, I show that most variation in diversification rates among 12 major vertebrate clades has a simple ecological explanation: predominantly terrestrial clades (i.e. birds, mammals, and lizards and snakes) have higher net diversification rates than predominantly aquatic clades (i.e. amphibians, crocodilians, turtles and all fish clades). These differences in diversification rates are then strongly related to patterns of species richness. Habitat may be more important than other potential explanations for richness patterns in vertebrates (such as climate and metabolic rates) and may also help explain patterns of species richness in many other groups of organisms.

The Ancient Evolutionary History of Polyomaviruses

Polyomaviruses are a family of DNA tumor viruses that are known to infect mammals and birds. To investigate the deeper evolutionary history of the family, we used a combination of viral metagenomics, bioinformatics, and structural modeling approaches to identify and characterize polyomavirus sequences associated with fish and arthropods. Analyses drawing upon the divergent new sequences indicate that polyomaviruses have been gradually co-evolving with their animal hosts for at least half a billion years. Phylogenetic analyses of individual polyomavirus genes suggest that some modern polyomavirus species arose after ancient recombination events involving distantly related polyomavirus lineages. The improved evolutionary model provides a useful platform for developing a more accurate taxonomic classification system for the viral family Polyomaviridae.

Polyomaviruses are a family of DNA-based viruses that are known to infect various terrestrial vertebrates, including humans. In this report, we describe our discovery of highly divergent polyomaviruses associated with various marine fish. Searches of public deep sequencing databases unexpectedly revealed the existence of polyomavirus-like sequences in scorpion and spider datasets. Our analysis of these new sequences suggests that polyomaviruses have slowly co-evolved with individual host animal lineages through an established mechanism known as intrahost divergence. The proposed model is similar to the mechanisms through with other DNA viruses, such as papillomaviruses, are thought to have evolved. Our analysis also suggests that distantly related polyomaviruses sometimes recombine to produce new chimeric lineages. We propose a possible taxonomic scheme that can account for these inferred ancient recombination events.

Molecular Evolution and Functional Divergence of Trace Amine-Associated Receptors

Trace amine-associated receptors (TAARs) are a member of the G-protein-coupled receptor superfamily and are known to be expressed in olfactory sensory neurons. A limited number of molecular evolutionary studies have been done for TAARs so far. To elucidate how lineage-specific evolution contributed to their functional divergence, we examined 30 metazoan genomes. In total, 493 TAAR gene candidates (including 84 pseudogenes) were identified from 26 vertebrate genomes. TAARs were not identified from non-vertebrate genomes. An ancestral-type TAAR-like gene appeared to have emerged in lamprey. We found four therian-specific TAAR subfamilies (one eutherian-specific and three metatherian-specific) in addition to previously known nine subfamilies. Many species-specific TAAR gene duplications and losses contributed to a large variation of TAAR gene numbers among mammals, ranging from 0 in dolphin to 26 in flying fox. TAARs are classified into two groups based on binding preferences for primary or tertiary amines as well as their sequence similarities. Primary amine-detecting TAARs (TAAR1-4) have emerged earlier, generally have single-copy orthologs (very few duplication or loss), and have evolved under strong functional constraints. In contrast, tertiary amine-detecting TAARs (TAAR5-9) have emerged more recently and the majority of them experienced higher rates of gene duplications. Protein members that belong to the tertiary amine-detecting TAAR group also showed the patterns of positive selection especially in the area surrounding the ligand-binding pocket, which could have affected ligand-binding activities and specificities. Expansions of the tertiary amine-detecting TAAR gene family may have played important roles in terrestrial adaptations of therian mammals. Molecular evolution of the TAAR gene family appears to be governed by a complex, species-specific, interplay between environmental and evolutionary factors.

Structure and vascularization of the ventricular myocardium in Holocephali: their evolutionary significance

It was generally assumed that the ventricle of the primitive vertebrate heart was composed of trabeculated, or spongy, myocardium, supplied by oxygen-poor luminal blood. In addition, it was presumed that the mixed ventricular myocardium, consisting of a compacta and a spongiosa, and its supply through coronary arteries appeared several times throughout fish evolution. Recent work has suggested, however, that a fully vascularized, mixed myocardium may be the primitive condition in gnathostomes. The present study of the heart ventricles of four holocephalan species aimed to clarify this controversy. Our observations showed that the ventricular myocardium of Chimaera monstrosa and Harriotta raleighana consists of a very thin compacta overlying a widespread spongiosa. The ventricle of Hydrolagus affinis is composed exclusively of trabeculated myocardium. In these three species there is a well-developed coronary artery system. The main coronary artery trunks run along the outflow tract, giving off subepicardial ventricular arteries. The trabeculae of the spongiosa are irrigated by branches of the subepicardial arteries and by penetrating arterial vessels arising directly from the main coronary trunks at the level of the conoventricular junction. The ventricle of Rhinochimaera atlantica has only spongy myocardium supplied by luminal blood. Small coronary arterial vessels are present in the subepicardium, but they do not enter the myocardial trabeculae. The present findings show for the first time that in a wild living vertebrate species, specifically H. affinis, an extensive coronary artery system supplying the whole cardiac ventricle exists in the absence of a well-developed compact ventricular myocardium. This is consistent with the notion derived from experimental work that myocardial cell proliferation and coronary vascular growth rely on distinct developmental programs. Our observations, together with data in the literature on elasmobranchs, support the view that the mixed ventricular myocardium is primitive for chondrichthyans. The reduction or even lack of compacta in holocephali has to be regarded as a derived anatomical trait. Our findings also fit in with the view that the mixed myocardium was the primitive condition in gnathostomes, and that the absence of compact ventricular myocardium in different actinopterygian groups is the result of a repeated loss of such type of cardiac muscle during fish evolution.

Evolutionary Analysis and Classification of OATs, OCTs, OCTNs, and Other SLC22 Transporters: Structure-Function Implications and Analysis of Sequence Motifs

The SLC22 family includes organic anion transporters (OATs), organic cation transporters (OCTs) and organic carnitine and zwitterion transporters (OCTNs). These are often referred to as drug transporters even though they interact with many endogenous metabolites and signaling molecules (Nigam, S.K., Nature Reviews Drug Discovery, 14:29-44, 2015). Phylogenetic analysis of SLC22 supports the view that these transporters may have evolved over 450 million years ago. Many OAT members were found to appear after a major expansion of the SLC22 family in mammals, suggesting a physiological and/or toxicological role during the mammalian radiation. Putative SLC22 orthologs exist in worms, sea urchins, flies, and ciona. At least six groups of SLC22 exist. OATs and OCTs form two Major clades of SLC22, within which (apart from Oat and Oct subclades), there are also clear Oat-like, Octn, and Oct-related subclades, as well as a distantly related group we term "Oat-related" (which may have different functions). Based on available data, it is arguable whether SLC22A18, which is related to bacterial drug-proton antiporters, should be assigned to SLC22. Disease-causing mutations, single nucleotide polymorphisms (SNPs) and other functionally analyzed mutations in OAT1, OAT3, URAT1, OCT1, OCT2, OCTN1, and OCTN2 map to the first extracellular domain, the large central intracellular domain, and transmembrane domains 9 and 10. These regions are highly conserved within subclades, but not between subclades, and may be necessary for SLC22 transporter function and functional diversification. Our results not only link function to evolutionarily conserved motifs but indicate the need for a revised sub-classification of SLC22.

Molecular signatures that are distinctive characteristics of the vertebrates and chordates and supporting a grouping of vertebrates with the tunicates

Members of the phylum Chordata and the subphylum Vertebrata are presently distinguished solely on the basis of morphological characteristics. The relationship of the vertebrates to the two non-vertebrate chordate subphyla is also a subject of debate. Analyses of protein sequences have identified multiple conserved signature indels (CSIs) that are specific for Chordata or for Vertebrata. Five CSIs in 4 important proteins are specific for the Vertebrata, whereas two other CSIs are uniquely found in all sequenced chordate species including Ciona intestinalis and Oikapleura dioica (Tunicates) as well as Branchiostoma floridae (Cephalochordates). The shared presence of these molecular signatures by all vertebrates/chordate species, but in no other animal taxa, strongly indicates that the genetic changes represented by the identified CSIs diagnose monophyletic groups. Two other discovered CSIs are uniquely shared by different vertebrate species and by either one (Ciona intestinalis) or both tunicate (Ciona and Oikapleura) species, but they are not found in Branchiostoma or other animal species. Specific presence of these CSIs in different vertebrates and either one or both tunicate species provides strong independent evidence that the vertebrate species are more closely related to the urochordates (tunicates) than to the cephalochordates.

Is there convergence in the molecular pathways underlying the repeated evolution of sociality in African cichlids?

Despite wide variation in the complexity of social interactions across taxa, the basic behavioral components of sociality appear to be modulated by conserved hormone pathways. Specifically, the nonapeptide hormones oxytocin and vasopressin and their receptors have been implicated in regulating diverse social behaviors across vertebrates. Here, we took advantage of the repeated evolution of cooperative breeding in African cichlids to investigate whether there are consistent brain gene expression patterns of isotocin and arginine vasotocin (teleost homologues of oxytocin and vasopressin), as well as their receptors, between four closely related pairs of social (cooperative) and non-social (non-cooperative) species. We first found that the coding sequences for the five genes studied were highly conserved across the eight species. This is the first study to examine the expression of both isotocin receptors, and so we performed a phylogenetic analysis that suggests that these two isotocin receptors are paralogues that arose during the teleost genome duplication. When we then examined brain gene expression patterns relative to social system, we found that there were whole-brain gene expression differences between the social and non-social species in many of the species pairs. However, these relationships varied in both the direction and magnitude among the four species pairs. In conclusion, our results suggest high sequence conservation and species-specific gene expression patterns relative to social behavior for these candidate hormone pathways in the cichlid fishes.

Up in Arms: Immune and Nervous System Response to Sea Star Wasting Disease

Echinoderms, positioned taxonomically at the base of deuterostomes, provide an important system for the study of the evolution of the immune system. However, there is little known about the cellular components and genes associated with echinoderm immunity. The 2013-2014 sea star wasting disease outbreak is an emergent, rapidly spreading disease, which has led to large population declines of asteroids in the North American Pacific. While evidence suggests that the signs of this disease, twisting arms and lesions, may be attributed to a viral infection, the host response to infection is still poorly understood. In order to examine transcriptional responses of the sea star Pycnopodia helianthoides to sea star wasting disease, we injected a viral sized fraction (0.2 μm) homogenate prepared from symptomatic P. helianthoides into apparently healthy stars. Nine days following injection, when all stars were displaying signs of the disease, specimens were sacrificed and coelomocytes were extracted for RNA-seq analyses. A number of immune genes, including those involved in Toll signaling pathways, complement cascade, melanization response, and arachidonic acid metabolism, were differentially expressed. Furthermore, genes involved in nervous system processes and tissue remodeling were also differentially expressed, pointing to transcriptional changes underlying the signs of sea star wasting disease. The genomic resources presented here not only increase understanding of host response to sea star wasting disease, but also provide greater insight into the mechanisms underlying immune function in echinoderms.

Molecular and neuronal homology between the olfactory systems of zebrafish and mouse

Studies of the two major olfactory organs of rodents, the olfactory mucosa (OM) and the vomeronasal organ (VNO), unraveled the molecular basis of smell in vertebrates. However, some vertebrates lack a VNO. Here we generated and analyzed the olfactory transcriptome of the zebrafish and compared it to the olfactory transcriptomes of mouse to investigate the evolutionary and molecular relationship between single and dual olfactory systems. Our analyses revealed a high degree of molecular conservation, with orthologs of mouse olfactory cell-specific markers and all but one of their chemosensory receptor classes expressed in the single zebrafish olfactory organ. Zebrafish chemosensory receptor genes are expressed across a large dynamic range and their RNA abundance correlates positively with the number of neurons expressing that RNA. Thus we estimate the relative proportions of neuronal sub-types expressing different chemosensory receptors. Receptor repertoire size drives the absolute abundance of different classes of neurons, but we find similar underlying patterns in both species. Finally, we identified novel marker genes that characterize rare neuronal populations in both mouse and zebrafish. In sum, we find that the molecular and cellular mechanisms underpinning olfaction in teleosts and mammals are similar despite 430 million years of evolutionary divergence.

The deuterostome context of chordate origins

Our understanding of vertebrate origins is powerfully informed by comparative morphology, embryology and genomics of chordates, hemichordates and echinoderms, which together make up the deuterostome clade. Striking body-plan differences among these phyla have historically hindered the identification of ancestral morphological features, but recent progress in molecular genetics and embryology has revealed deep similarities in body-axis formation and organization across deuterostomes, at stages before morphological differences develop. These developmental genetic features, along with robust support of pharyngeal gill slits as a shared deuterostome character, provide the foundation for the emergence of chordates.

Identification and response to metals of metallothionein in two ancient fishes: white sturgeon (Acipenser transmontanus) and lake sturgeon (Acipenser fulvescens)

White sturgeon (Acipenser transmontanus) are among the most sensitive species of fishes to Cu, Cd, and Zn, but there is no information about sensitivity of lake sturgeon (Acipenser fulvescens). To begin to elucidate molecular mechanism(s) of sensitivity of sturgeons to metals a cDNA encoding metallothionein (MT) was amplified from livers of white sturgeon (WS-MT) and lake sturgeon (LS-MT), and expression in response to Cu, Cd, or Zn was characterized in liver explants from each species. The primary structure of WS-MT and LS-MT contained 20 cysteine residues, which is the same as MTs of teleost fishes. However, the primary structure of WS-MT and LS-MT contained 63 amino acids, which is longer than any MT identified in teleost fishes. Abundance of transcripts of WS-MT in explants exposed to 0.3, 3, 30, or 100 μg/L of Cu was 1.7-, 1.7-, 2.1-, and 2.6-fold less than in controls, respectively. In contrast, abundances of transcripts of WS-MT were 3.3- and 2.4-fold greater in explants exposed to 30 μg/L of Cd and 1000 μg/L of Zn, respectively. Abundance of transcripts of LS-MT was not significantly different at any concentration of Cu, Cd, or Zn. MT is hypothesized to represent a critical mechanism for detoxification of metals. Therefore, results of this study suggest that sensitivity of sturgeons to exposure to Cu, Cd, or Zn might be a result of the relatively lesser maximal response of MT to metals. The study also suggestslake sturgeon might be more sensitive than white sturgeon to metals.

Comparative analysis of cystatin superfamily in platyhelminths

The cystatin superfamily is comprised of cysteine proteinase inhibitors and encompasses at least 3 subfamilies: stefins, cystatins and kininogens. In this study, the platyhelminth cystatin superfamily was identified and grouped into stefin and cystatin subfamilies. The conserved domain of stefins (G, QxVxG) was observed in all members of platyhelminth stefins. The three characteristics of cystatins, the cystatin-like domain (G, QxVxG, PW), a signal peptide, and one or two conserved disulfide bonds, were observed in platyhelminths, with the exception of cestodes, which lacked the conserved disulfide bond. However, it is noteworthy that cestode cystatins had two tandem repeated domains, although the second tandem repeated domain did not contain a cystatin-like domain, which has not been previously reported. Tertiary structure analysis of Taenia solium cystatin, one of the cestode cystatins, demonstrated that the N-terminus of T. solium cystatin formed a five turn α-helix, a five stranded β-pleated sheet and a hydrophobic edge, similar to the structure of chicken cystatin. Although no conserved disulfide bond was found in T. solium cystatin, the models of T. solium cystatin and chicken cystatin corresponded at the site of the first disulfide bridge of the chicken cystatin. However, the two models were not similar regarding the location of the second disulfide bridge of chicken cystatin. These results showed that T. solium cystatin and chicken cystatin had similarities and differences, suggesting that the biochemistry of T. solium cystatin could be similar to chicken cystatin in its inhibitory function and that it may have further functional roles. The same results were obtained for other cestode cystatins. Phylogenetic analysis showed that cestode cystatins constituted an independent clade and implied that cestode cystatins should be considered to have formed a new clade during evolution.

A Protocol for Diagnosing the Effect of Calibration Priors on Posterior Time Estimates: A Case Study for the Cambrian Explosion of Animal Phyla

We present a procedure to test the effect of calibration priors on estimated times, which applies a recently developed calibration-free approach (RelTime) method that produces relative divergence times for all nodes in the tree. We illustrate this protocol by applying it to a timetree of metazoan diversification (Erwin DH, Laflamme M, Tweedt SM, Sperling EA, Pisani D, Peterson KJ. 2011. The Cambrian conundrum: early divergence and later ecological success in the early history of animals. Science 334:1091-1097.), which placed the divergence of animal phyla close to the time of the Cambrian explosion inferred from the fossil record. These analyses revealed that the two maximum-only calibration priors in the pre-Cambrian are the primary determinants of the young divergence times among animal phyla in this study. In fact, these two maximum-only calibrations produce divergence times that severely violate minimum boundaries of almost all of the other 22 calibration constraints. The use of these 22 calibrations produces dates for metazoan divergences that are hundreds of millions of years earlier in the Proterozoic. Our results encourage the use of calibration-free approaches to identify most influential calibration constraints and to evaluate their impact in order to achieve biologically robust interpretations.

Nuclear genomic signals of the 'microturbellarian' roots of platyhelminth evolutionary innovation

Flatworms number among the most diverse invertebrate phyla and represent the most biomedically significant branch of the major bilaterian clade Spiralia, but to date, deep evolutionary relationships within this group have been studied using only a single locus (the rRNA operon), leaving the origins of many key clades unclear. In this study, using a survey of genomes and transcriptomes representing all free-living flatworm orders, we provide resolution of platyhelminth interrelationships based on hundreds of nuclear protein-coding genes, exploring phylogenetic signal through concatenation as well as recently developed consensus approaches. These analyses robustly support a modern hypothesis of flatworm phylogeny, one which emphasizes the primacy of the often-overlooked 'microturbellarian' groups in understanding the major evolutionary transitions within Platyhelminthes: perhaps most notably, we propose a novel scenario for the interrelationships between free-living and vertebrate-parasitic flatworms, providing new opportunities to shed light on the origins and biological consequences of parasitism in these iconic invertebrates.

Evolution of basal deuterostome nervous systems

Understanding the evolution of deuterostome nervous systems has been complicated by the ambiguous phylogenetic position of the Xenocoelomorpha (Xenoturbellids, acoel flat worms, nemertodermatids), which has been placed either as basal bilaterians, basal deuterostomes or as a sister group to the hemichordate/echinoderm clade (Ambulacraria), which is a sister group of the Chordata. None of these groups has a single longitudinal nerve cord and a brain. A further complication is that echinoderm nerve cords are not likely to be evolutionarily related to the chordate central nervous system. For hemichordates, opinion is divided as to whether either one or none of the two nerve cords is homologous to the chordate nerve cord. In chordates, opposition by two secreted signaling proteins, bone morphogenetic protein (BMP) and Nodal, regulates partitioning of the ectoderm into central and peripheral nervous systems. Similarly, in echinoderm larvae, opposition between BMP and Nodal positions the ciliary band and regulates its extent. The apparent loss of this opposition in hemichordates is, therefore, compatible with the scenario, suggested by Dawydoff over 65 years ago, that a true centralized nervous system was lost in hemichordates.

Did homeobox gene duplications contribute to the Cambrian explosion?

The Cambrian explosion describes an apparently rapid increase in the diversity of bilaterian animals around 540-515 million years ago. Bilaterian animals explore the world in three-dimensions deploying forward-facing sense organs, a brain, and an anterior mouth; they possess muscle blocks enabling efficient crawling and burrowing in sediments, and they typically have an efficient 'through-gut' with separate mouth and anus to process bulk food and eject waste, even when burrowing in sediment. A variety of ecological, environmental, genetic, and developmental factors have been proposed as possible triggers and correlates of the Cambrian explosion, and it is likely that a combination of factors were involved. Here, I focus on a set of developmental genetic changes and propose these are part of the mix of permissive factors. I describe how ANTP-class homeobox genes, which encode transcription factors involved in body patterning, increased in number in the bilaterian stem lineage and earlier. These gene duplications generated a large array of ANTP class genes, including three distinct gene clusters called NK, Hox, and ParaHox. Comparative data supports the idea that NK genes were deployed primarily to pattern the bilaterian mesoderm, Hox genes coded position along the central nervous system, and ParaHox genes most likely originally specified the mouth, midgut, and anus of the newly evolved through-gut. It is proposed that diversification of ANTP class genes played a role in the Cambrian explosion by contributing to the patterning systems used to build animal bodies capable of high-energy directed locomotion, including active burrowing.

Decelerated genome evolution in modern vertebrates revealed by analysis of multiple lancelet genomes

Vertebrates diverged from other chordates ~500 Myr ago and experienced successful innovations and adaptations, but the genomic basis underlying vertebrate origins are not fully understood. Here we suggest, through comparison with multiple lancelet (amphioxus) genomes, that ancient vertebrates experienced high rates of protein evolution, genome rearrangement and domain shuffling and that these rates greatly slowed down after the divergence of jawed and jawless vertebrates. Compared with lancelets, modern vertebrates retain, at least relatively, less protein diversity, fewer nucleotide polymorphisms, domain combinations and conserved non-coding elements (CNE). Modern vertebrates also lost substantial transposable element (TE) diversity, whereas lancelets preserve high TE diversity that includes even the long-sought RAG transposon. Lancelets also exhibit rapid gene turnover, pervasive transcription, fastest exon shuffling in metazoans and substantial TE methylation not observed in other invertebrates. These new lancelet genome sequences provide new insights into the chordate ancestral state and the vertebrate evolution.

The lancelet, or amphioxus, is an extant basal chordate that diverged from other chordate lineages about 550 million years ago. Here the authors sequence and assemble the diploid genome of a male adult of the Chinese lancelet, B. belcheri, and highlight genomic features that may have played an important role in the origin and evolution of vertebrates.

Phylostratigraphic profiles in zebrafish uncover chordate origins of the vertebrate brain

An elaborated tripartite brain is considered one of the important innovations of vertebrates. Other extant chordate groups have a more basic brain organization. For instance, cephalochordates possess a relatively simple brain possibly homologous to the vertebrate forebrain and hindbrain, whereas tunicates display the tripartite organization, but without the specialized brain centers. The difference in anatomical complexity is even more pronounced if one compares chordates with other deuterostomes that have only a diffuse nerve net or alternatively a rather simple central nervous system. To gain a new perspective on the evolutionary roots of the complex vertebrate brain, we made here a phylostratigraphic analysis of gene expression patterns in the developing zebrafish (Danio rerio). The recovered adaptive landscape revealed three important periods in the evolutionary history of the zebrafish brain. The oldest period corresponds to preadaptive events in the first metazoans and the emergence of the nervous system at the metazoan-eumetazoan transition. The origin of chordates marks the next phase, where we found the overall strongest adaptive imprint in almost all analyzed brain regions. This finding supports the idea that the vertebrate brain evolved independently of the brains within the protostome lineage. Finally, at the origin of vertebrates we detected a pronounced signal coming from the dorsal telencephalon, in agreement with classical theories that consider this part of the cerebrum a genuine vertebrate innovation. Taken together, these results reveal a stepwise adaptive history of the vertebrate brain where most of its extant organization was already present in the chordate ancestor.

Multiple thyrotropin β-subunit and thyrotropin receptor-related genes arose during vertebrate evolution

Thyroid-stimulating hormone (TSH) is composed of a specific β subunit and an α subunit that is shared with the two pituitary gonadotropins. The three β subunits derive from a common ancestral gene through two genome duplications (1R and 2R) that took place before the radiation of vertebrates. Analysis of genomic data from phylogenetically relevant species allowed us to identify an additional Tshβ subunit-related gene that was generated through 2R. This gene, named Tshβ2, present in cartilaginous fish, little skate and elephant shark, and in early lobe-finned fish, coelacanth and lungfish, was lost in ray-finned fish and tetrapods. The absence of a second type of TSH receptor (Tshr) gene in these species suggests that both TSHs act through the same receptor. A novel Tshβ sister gene, named Tshβ3, was generated through the third genomic duplication (3R) that occurred early in the teleost lineage. Tshβ3 is present in most teleost groups but was lostin tedraodontiforms. The 3R also generated a second Tshr, named Tshrb. Interestingly, the new Tshrb was translocated from its original chromosomic position after the emergence of eels and was then maintained in its new position. Tshrb was lost in tetraodontiforms and in ostariophysians including zebrafish although the latter species have two TSHs, suggesting that TSHRb may be dispensable. The tissue distribution of duplicated Tshβs and Tshrs was studied in the European eel. The endocrine thyrotropic function in the eel would be essentially mediated by the classical Tshβ and Tshra, which are mainly expressed in the pituitary and thyroid, respectively. Tshβ3 and Tshrb showed a similar distribution pattern in the brain, pituitary, ovary and adipose tissue, suggesting a possible paracrine/autocrine mode of action in these non-thyroidal tissues. Further studies will be needed to determine the binding specificity of the two receptors and how these two TSH systems are interrelated.

On the age of eukaryotes: evaluating evidence from fossils and molecular clocks

Our understanding of the phylogenetic relationships among eukaryotic lineages has improved dramatically over the few past decades thanks to the development of sophisticated phylogenetic methods and models of evolution, in combination with the increasing availability of sequence data for a variety of eukaryotic lineages. Concurrently, efforts have been made to infer the age of major evolutionary events along the tree of eukaryotes using fossil-calibrated molecular clock-based methods. Here, we review the progress and pitfalls in estimating the age of the last eukaryotic common ancestor (LECA) and major lineages. After reviewing previous attempts to date deep eukaryote divergences, we present the results of a Bayesian relaxed-molecular clock analysis of a large dataset (159 proteins, 85 taxa) using 19 fossil calibrations. We show that for major eukaryote groups estimated dates of divergence, as well as their credible intervals, are heavily influenced by the relaxed molecular clock models and methods used, and by the nature and treatment of fossil calibrations. Whereas the estimated age of LECA varied widely, ranging from 1007 (943-1102) Ma to 1898 (1655-2094) Ma, all analyses suggested that the eukaryotic supergroups subsequently diverged rapidly (i.e., within 300 Ma of LECA). The extreme variability of these and previously published analyses preclude definitive conclusions regarding the age of major eukaryote clades at this time. As more reliable fossil data on eukaryotes from the Proterozoic become available and improvements are made in relaxed molecular clock modeling, we may be able to date the age of extant eukaryotes more precisely.

Comparative evolutionary genomics of the STAT family of transcription factors

The STAT signaling pathway is one of the seven common pathways that govern cell fate decisions during animal development. Comparative genomics revealed multiple incidences of stat gene duplications throughout metazoan evolutionary history. While pseudogenization is a frequent fate of duplicated genes, many of these STAT duplications evolved into novel genes through rapid sequence diversification and neofunctionalization. Additionally, the core of STAT gene regulatory networks, comprising stat1 through 4, stat5 and stat6, arose early in vertebrate evolution, probably through the two whole genome duplication events that occurred after the split of Cephalochordates but before the rise of Chondrichthyes. While another complete genome duplication event took place during the evolution of bony fish after their separation from the tetrapods about 450 million years ago (Mya), modern fish have only one set of these core stats, suggesting the rapid loss of most duplicated stat genes. The two stat5 genes in mammals likely arose from a duplication event in early Eutherian evolution, a period from about 310 Mya at the avian-mammal divergence to the separation of marsupials from other mammals about 130 Mya. These analyses indicate that whole genome duplications and gene duplications by unequal chromosomal crossing over were likely the major mechanisms underlying the evolution of STATs.

The structural analysis of shark IgNAR antibodies reveals evolutionary principles of immunoglobulins

Sharks and other cartilaginous fish are the phylogenetically oldest living organisms that rely on antibodies as part of their adaptive immune system. They produce the immunoglobulin new antigen receptor (IgNAR), a homodimeric heavy chain-only antibody, as a major part of their humoral adaptive immune response. Here, we report the atomic resolution structure of the IgNAR constant domains and a structural model of this heavy chain-only antibody. We find that despite low sequence conservation, the basic Ig fold of modern antibodies is already present in the evolutionary ancient shark IgNAR domains, highlighting key structural determinants of the ubiquitous Ig fold. In contrast, structural differences between human and shark antibody domains explain the high stability of several IgNAR domains and allowed us to engineer human antibodies for increased stability and secretion efficiency. We identified two constant domains, C1 and C3, that act as dimerization modules within IgNAR. Together with the individual domain structures and small-angle X-ray scattering, this allowed us to develop a structural model of the complete IgNAR molecule. Its constant region exhibits an elongated shape with flexibility and a characteristic kink in the middle. Despite the lack of a canonical hinge region, the variable domains are spaced appropriately wide for binding to multiple antigens. Thus, the shark IgNAR domains already display the well-known Ig fold, but apart from that, this heavy chain-only antibody employs unique ways for dimerization and positioning of functional modules.

Identification and expression of aryl hydrocarbon receptors (AhR1 and AhR2) provide insight in an evolutionary context regarding sensitivity of white sturgeon (Acipenser transmontanus) to dioxin-like compounds

Sturgeons are ancient fishes, which are endangered in many parts of the world. Due to their benthic nature and longevity, sturgeon are at great risk of exposure to bioaccumulative contaminants such as dioxin-like compounds (DLCs). Despite their endangered status, little research has been conducted to characterize the relative sensitivity of sturgeons to DLCs. Proper assessment of risk of DLCs posed to these fishes therefore, requires a better understanding of this sensitivity and the factors that are driving it. Adverse effects associated with exposure to DLCs are mediated by the aryl hydrocarbon receptor (AhR). This study identified and characterized two distinct AhRs, AhR1 and AhR2, in white sturgeon (Acipenser transmontanus) for the first time as a first step in studying the relative sensitivities of sturgeons to DLCs. Furthermore, tissue-specific expression of both AhRs under basal conditions and in response to exposure to the model DLC, β-naphthoflavone (βNF), was determined. The sequence of amino acids of AhR1 of white sturgeon had greater similarity to AhRs of tetrapods, including amphibians, birds, and mammals, than to AhR1s of other fishes. The sequence of amino acids in the ligand binding domain of the AhR1 had greater than 80% similarity to AhRs known to bind DLCs and was less similar to AhRs not known to bind DLCs. AhR2 of white sturgeon had greatest similarity to AhR2 of other fishes. Profiles of expression of AhR1 and AhR2 in white sturgeon were distinct from those known in other fishes and appear more similar to profiles observed in birds. Expressions of both AhR1 and AhR2 of white sturgeon were greatest in liver and heart, which are target organs for DLCs. Furthermore, abundances of transcripts of AhR1 and AhR2 in all tissues from white sturgeon were greater than controls (up to 35-fold) following exposure to βNF. Based upon both AhRs having similar abundances of transcript in target organs of DLC toxicity, both AhRs being up-regulated following exposure to βNF, and both AhRs having greatest similarity to AhRs known to bind DLCs, it is hypothesized that both AhR1 and AhR2 of white sturgeon might mediate effects of DLCs in this species. Since current risk assessments are based on data derived largely from highly divergent fishes within the Salmonidae, presence of two functional AhRs in white sturgeon, one of which has greatest similarity to AhRs of birds, might have significant implications for the sensitivity of sturgeons to DLCs compared to other fishes.

The characters of Palaeozoic jawed vertebrates

Newly discovered fossils from the Silurian and Devonian periods are beginning to challenge embedded perceptions about the origin and early diversification of jawed vertebrates (gnathostomes). Nevertheless, an explicit cladistic framework for the relationships of these fossils relative to the principal crown lineages of the jawed vertebrates (osteichthyans: bony fishes and tetrapods; chondrichthyans: sharks, batoids, and chimaeras) remains elusive. We critically review the systematics and character distributions of early gnathostomes and provide a clearly stated hierarchy of synapomorphies covering the jaw-bearing stem gnathostomes and osteichthyan and chondrichthyan stem groups. We show that character lists, designed to support the monophyly of putative groups, tend to overstate their strength and lack cladistic corroboration. By contrast, synapomorphic hierarchies are more open to refutation and must explicitly confront conflicting evidence. Our proposed synapomorphy scheme is used to evaluate the status of the problematic fossil groups Acanthodii and Placodermi, and suggest profitable avenues for future research. We interpret placoderms as a paraphyletic array of stem-group gnathostomes, and suggest what we regard as two equally plausible placements of acanthodians: exclusively on the chondrichthyan stem, or distributed on both the chondrichthyan and osteichthyan stems.