Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences

Odorant transport in a hagfish

Odorant transport is of fundamental and applied importance. Using computational simulations, we studied odorant transport in an anatomically accurate model of the nasal passage of a hagfish (probably Eptatretus stoutii). We found that ambient water is sampled widely, with a significant ventral element. Additionally, there is a bilateral element to olfactory flow, which enters the single nostril in two narrow, laminar streams that are then split prior to the nasal chamber by the anterior edge of the central olfactory lamella. An appendage on this lamella directs a small portion (10-14%) of the overall nasal flow to the olfactory sensory channels. Much of the remaining flow is diverted away from the sensory channels by two peripheral channels. The anterior edge of the central olfactory lamella, together with a jet-impingement mechanism, disperses flow over the olfactory surfaces. Diffusion of odorant from bulk water to the olfactory surfaces is facilitated by the large surface area:volume ratio of the sensory channels, and by a resistance-based hydrodynamic mechanism that leads to long residence times (up to 4.5 s) in the sensory channels. With increasing volumetric flow rate, the rate of odorant transfer to the olfactory surfaces increases, but the efficiency of odorant uptake decreases, falling in the range 2-6%. Odorant flux decreases caudally across the olfactory surfaces, suggesting in vivo a preponderance of olfactory sensory neurons on the anterior part of each olfactory surface. We conclude that the hagfish has a subtle anatomy for locating and capturing odorant molecules.

Molecular Dating of the Teleost Whole Genome Duplication (3R) Is Compatible With the Expectations of Delayed Rediploidization

Vertebrate evolution has been punctuated by three whole genome duplication events that have been implicated causally in phenotypic evolution, from the origin of phenotypic novelties to explosive diversification. Arguably, the most dramatic of these is the 3R whole genome duplication event associated with the origin of teleost fishes which comprise more than half of all living vertebrate species. However, tests of a causal relationship between whole genome duplication and teleost diversification have proven difficult due to the challenge of establishing the timing of these phenomena. Here we show, based on molecular clock dating of concatenated gene alignments, that the 3R whole genome duplication event occurred in the early-middle Permian (286.18 to 267.20 million years ago; Ma), 52.02 to 12.84 million years (Myr) before the divergence of crown-teleosts in the latest Permian-earliest Late Triassic (254.36 to 234.16 Ma) and long before the major pulses of teleost diversification in Ostariophysi and Percomorpha (56.37 to 100.17 Myr and at least 139.24 to 183.29 Myr later, respectively). The extent of this temporal gap between putative cause and effect precludes 3R as a deterministic driver of teleost diversification. However, these age constraints remain compatible with the expectations of a prolonged rediploidization process following whole genome duplication which, through the effects of chromosome rearrangement and gene loss, remains a viable mechanism to explain the evolution of teleost novelties and diversification.

ETV2 induces endothelial, but not hematopoietic, lineage specification in birds

Cardiovascular system develops from the lateral plate mesoderm. Its three primary cell lineages (hematopoietic, endothelial, and muscular) are specified by the sequential actions of conserved transcriptional factors. ETV2, a master regulator of mammalian hemangioblast development, however, is absent in the chicken genome and acts downstream of NPAS4L in zebrafish. Here, we investigated the epistatic relationship between NPAS4L and ETV2 in avian hemangioblast development. We showed that ETV2 is deleted in all 363 avian genomes analyzed. Mouse ETV2 induced LMO2, but not NPAS4L or SCL, expression in chicken mesoderm. Squamate (lizards, geckos, and snakes) genomes contain both NPAS4L and ETV2 In Madagascar ground gecko, both genes were expressed in developing hemangioblasts. Gecko ETV2 induced only LMO2 in chicken mesoderm. We propose that both NPAS4L and ETV2 were present in ancestral amniote, with ETV2 acting downstream of NPAS4L in endothelial lineage specification. ETV2 may have acted as a pioneer factor by promoting chromatin accessibility of endothelial-specific genes and, in parallel with NPAS4L loss in ancestral mammals, has gained similar function in regulating blood-specific genes.

slc26a12-A novel member of the slc26 family, is located in tandem with slc26a2 in coelacanths, amphibians, reptiles, and birds

Solute carrier family 26 (Slc26) is a family of anion exchangers with 11 members in mammals (named Slc26a1-a11). Here, we identified a novel member of the slc26 family, slc26a12, located in tandem with slc26a2 in the genomes of several vertebrate lineages. BLAST and synteny analyses of various jawed vertebrate genome databases revealed that slc26a12 is present in coelacanths, amphibians, reptiles, and birds but not in cartilaginous fishes, lungfish, mammals, or ray-finned fishes. In some avian and reptilian lineages such as owls, penguins, egrets, and ducks, and most turtles examined, slc26a12 was lost or pseudogenized. Phylogenetic analysis showed that Slc26a12 formed an independent branch with the other Slc26 members and Slc26a12, Slc26a1 and Slc26a2 formed a single branch, suggesting that these three members formed a subfamily in Slc26. In jawless fish, hagfish have two genes homologous to slc26a2 and slc26a12, whereas lamprey has a single gene homologous to slc26a2. African clawed frogs express slc26a12 in larval gills, skin, and fins. These results show that slc26a12 was present at least before the separation of lobe-finned fish and tetrapods; the name slc26a12 is appropriate because the gene duplication occurred in the distant past.

Dollo Parsimony Overestimates Ancestral Gene Content Reconstructions

Ancestral reconstruction is a widely used technique that has been applied to understand the evolutionary history of gain and loss of gene families. Ancestral gene content can be reconstructed via different phylogenetic methods, but many current and previous studies employ Dollo parsimony. We hypothesize that Dollo parsimony is not appropriate for ancestral gene content reconstruction inferences based on sequence homology, as Dollo parsimony is derived from the assumption that a complex character cannot be regained. This premise does not accurately model molecular sequence evolution, in which false orthology can result from sequence convergence or lateral gene transfer. The aim of this study is to test Dollo parsimony's suitability for ancestral gene content reconstruction and to compare its inferences with a maximum likelihood-based approach that allows a gene family to be gained more than once within a tree. We first compared the performance of the two approaches on a series of artificial data sets each of 5,000 genes that were simulated according to a spectrum of evolutionary rates without gene gain or loss, so that inferred deviations from the true gene count would arise only from errors in orthology inference and ancestral reconstruction. Next, we reconstructed protein domain evolution on a phylogeny representing known eukaryotic diversity. We observed that Dollo parsimony produced numerous ancestral gene content overestimations, especially at nodes closer to the root of the tree. These observations led us to the conclusion that, confirming our hypothesis, Dollo parsimony is not an appropriate method for ancestral reconstruction studies based on sequence homology.

Extensive DNA methylome rearrangement during early lamprey embryogenesis

DNA methylation (5mC) is a repressive gene regulatory mark widespread in vertebrate genomes, yet the developmental dynamics in which 5mC patterns are established vary across species. While mammals undergo two rounds of global 5mC erasure, teleosts, for example, exhibit localized maternal-to-paternal 5mC remodeling. Here, we studied 5mC dynamics during the embryonic development of sea lamprey, a jawless vertebrate which occupies a critical phylogenetic position as the sister group of the jawed vertebrates. We employed 5mC quantification in lamprey embryos and tissues, and discovered large-scale maternal-to-paternal epigenome remodeling that affects ~30% of the embryonic genome and is predominantly associated with partially methylated domains. We further demonstrate that sequences eliminated during programmed genome rearrangement (PGR), are hypermethylated in sperm prior to the onset of PGR. Our study thus unveils important insights into the evolutionary origins of vertebrate 5mC reprogramming, and how this process might participate in diverse developmental strategies.

Extended Cleavage Specificity of two Hematopoietic Serine Proteases from a Ray-Finned Fish, the Spotted Gar (Lepisosteus oculatus)

The extended cleavage specificities of two hematopoietic serine proteases originating from the ray-finned fish, the spotted gar (Lepisosteus oculatus), have been characterized using substrate phage display. The preference for particular amino acids at and surrounding the cleavage site was further validated using a panel of recombinant substrates. For one of the enzymes, the gar granzyme G, a strict preference for the aromatic amino acid Tyr was observed at the cleavable P1 position. Using a set of recombinant substrates showed that the gar granzyme G had a high selectivity for Tyr but a lower activity for cleaving after Phe but not after Trp. Instead, the second enzyme, gar DDN1, showed a high preference for Leu in the P1 position of substrates. This latter enzyme also showed a high preference for Pro in the P2 position and Arg in both P4 and P5 positions. The selectivity for the two Arg residues in positions P4 and P5 suggests a highly specific substrate selectivity of this enzyme. The screening of the gar proteome with the consensus sequences obtained by substrate phage display for these two proteases resulted in a very diverse set of potential targets. Due to this diversity, a clear candidate for a specific immune function of these two enzymes cannot yet be identified. Antisera developed against the recombinant gar enzymes were used to study their tissue distribution. Tissue sections from juvenile fish showed the expression of both proteases in cells in Peyer's patch-like structures in the intestinal region, indicating they may be expressed in T or NK cells. However, due to the lack of antibodies to specific surface markers in the gar, it has not been possible to specify the exact cellular origin. A marked difference in abundance was observed for the two proteases where gar DDN1 was expressed at higher levels than gar granzyme G. However, both appear to be expressed in the same or similar cells, having a lymphocyte-like appearance.