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

Phylogenomics resolves the puzzling phylogeny of banded newts (genus Ommatotriton)

Resolving the order of speciation events that occurred in rapid succession is inherently hard and typically requires a phylogenomic approach. A case in point concerns the previously unresolved phylogeny of the three species of banded newt (genus Ommatotriton). We obtain c. 7k nuclear DNA markers using target enrichment by sequence capture and analyze the dataset using maximum likelihood inference of concatenated data with RAxML, summary multi-species coalescent analysis with ASTRAL and Bayesian species tree inference using a diffusion model with SNAPPER, and use TreeMix and PhyloNet to test for interspecific gene flow. All analyses recover three distinct species with no evidence of interspecific gene flow. All analyses retrieved the topology (O. nesterovi, (O. ophryticus, O. vittatus)), with high support. SNAPPER did show the tendency to get stuck in a local optimum, resulting in a different but still highly supported topology. Furthermore, we notice that fewer SNAPPER runs get stuck in a local optimum when we include an outgroup. Therefore, we recommend the exploration of multiple independent runs and the use of an outgroup with this approach. The banded newt radiation illustrates the use of genome-wide data to tackle formerly unresolved phylogenies.

Untangling Zebrafish Genetic Annotation: Addressing Complexities and Nomenclature Issues in Orthologous Evaluation of TCOF1 and NOLC1

Treacher Collins syndrome (TCS) is a genetic disorder affecting facial development, primarily caused by mutations in the TCOF1 gene. TCOF1, along with NOLC1, play important roles in ribosomal RNA transcription and processing. Previously, a zebrafish model of TCS successfully recapitulated the main characteristics of the syndrome by knocking down the expression of a gene on chromosome 13 (coding for Uniprot ID B8JIY2), which was identified as the TCOF1 orthologue. However, database updates renamed this gene as nolc1 and the zebrafish database (ZFIN) identified a different gene on chromosome 14 as the TCOF1 orthologue (coding for Uniprot ID E7F9D9). NOLC1 and TCOF1 are large proteins with unstructured regions and repetitive sequences that complicate alignments and comparisons. Also, the additional whole genome duplication of teleosts sets further difficulty. In this study, we present evidence that endorses that NOLC1 and TCOF1 are paralogs, and that the zebrafish gene on chromosome 14 is a low-complexity LisH domain-containing factor that displays homology to NOLC1 but lacks essential sequence features to accomplish TCOF1 nucleolar functions. Our analysis also supports the idea that zebrafish, as has been suggested for other non-tetrapod vertebrates, lack the TCOF1 gene that is associated with tripartite nucleolus. Using BLAST searches in a group of teleost genomes, we identified fish-specific sequences similar to E7F9D9 zebrafish protein. We propose naming them "LisH-containing Low Complexity Proteins" (LLCP). Interestingly, the gene on chromosome 13 (nolc1) displays the sequence features, developmental expression patterns, and phenotypic impact of depletion that are characteristic of TCOF1 functions. These findings suggest that in teleost fish, the nucleolar functions described for both NOLC1 and TCOF1 mediated by their repeated motifs, are carried out by a single gene, nolc1. Our study, which is mainly based on computational tools available as free web-based algorithms, could help to solve similar conflicts regarding gene orthology in zebrafish.

The genomes of all lungfish inform on genome expansion and tetrapod evolution

The genomes of living lungfishes can inform on the molecular-developmental basis of the Devonian sarcopterygian fish-tetrapod transition. We de novo sequenced the genomes of the African (Protopterus annectens) and South American lungfishes (Lepidosiren paradoxa). The Lepidosiren genome (about 91 Gb, roughly 30 times the human genome) is the largest animal genome sequenced so far and more than twice the size of the Australian (Neoceratodus forsteri)1 and African2 lungfishes owing to enlarged intergenic regions and introns with high repeat content (about 90%). All lungfish genomes continue to expand as some transposable elements (TEs) are still active today. In particular, Lepidosiren's genome grew extremely fast during the past 100 million years (Myr), adding the equivalent of one human genome every 10 Myr. This massive genome expansion seems to be related to a reduction of PIWI-interacting RNAs and C2H2 zinc-finger and Krüppel-associated box (KRAB)-domain protein genes that suppress TE expansions. Although TE abundance facilitates chromosomal rearrangements, lungfish chromosomes still conservatively reflect the ur-tetrapod karyotype. Neoceratodus' limb-like fins still resemble those of their extinct relatives and remained phenotypically static for about 100 Myr. We show that the secondary loss of limb-like appendages in the Lepidosiren-Protopterus ancestor was probably due to loss of sonic hedgehog limb-specific enhancers.

Large carrion and burying beetles evolved from Staphylinidae (Coleoptera, Staphylinidae, Silphinae): a review of the evidence

Large carrion beetles (Silphidae) are the focus of ongoing behavioral ecology, forensic, ecological, conservation, evolutionary, systematic, and other research, and were recently reclassified as a subfamily of Staphylinidae. Twenty-three analyses in 21 publications spanning the years 1927-2023 that are relevant to the question of the evolutionary origin and taxonomic classification of Silphidae are reviewed. Most of these analyses (20) found Silphidae nested inside Staphylinidae (an average of 4.38 branches deep), two found Silphidae in an ambiguous position, and one found Silphidae outside Staphylinidae, as sister to Hydrophilidae. There is strong evidence supporting the hypothesis that large carrion beetles evolved from within Staphylinidae and good justification for their classification as the subfamily Silphinae of the megadiverse, and apparently now monophyletic, Staphylinidae. Considerable uncertainty remains regarding the interrelationships and monophyly of many staphylinid subfamilies. Nonetheless, the subfamily Tachyporinae was found to be the sister of Silphinae in more analyses (7) than any other subfamily.

A comprehensive phylogenomic study unveils evolutionary patterns and challenges in the mitochondrial genomes of Carcharhiniformes: A focus on Triakidae

The complex evolutionary patterns in the mitochondrial genome (mitogenome) of the most species-rich shark order, the Carcharhiniformes (ground sharks) has led to challenges in the phylogenomic reconstruction of the families and genera belonging to the order, particularly the family Triakidae (houndsharks). The current state of Triakidae phylogeny remains controversial, with arguments for both monophyly and paraphyly within the family. We hypothesize that this variability is triggered by the selection of different a priori partitioning schemes to account for site and gene heterogeneity within the mitogenome. Here we used an extensive statistical framework to select the a priori partitioning scheme for inference of the mitochondrial phylogenomic relationships within Carcharhiniformes, tested site heterogeneous CAT + GTR + G4 models and incorporated the multi-species coalescent model (MSCM) into our analyses to account for the influence of gene tree discordance on species tree inference. We included five newly assembled houndshark mitogenomes to increase resolution of Triakidae. During the assembly procedure, we uncovered a 714 bp-duplication in the mitogenome of Galeorhinus galeus. Phylogenetic reconstruction confirmed monophyly within Triakidae and the existence of two distinct clades of the expanded Mustelus genus. The latter alludes to potential evolutionary reversal of reproductive mode from placental to aplacental, suggesting that reproductive mode has played a role in the trajectory of adaptive divergence. These new sequences have the potential to contribute to population genomic investigations, species phylogeography delineation, environmental DNA metabarcoding databases and, ultimately, improved conservation strategies for these ecologically and economically important species.

A Practical Guide to Design and Assess a Phylogenomic Study

Over the last decade, molecular systematics has undergone a change of paradigm as high-throughput sequencing now makes it possible to reconstruct evolutionary relationships using genome-scale datasets. The advent of "big data" molecular phylogenetics provided a battery of new tools for biologists but simultaneously brought new methodological challenges. The increase in analytical complexity comes at the price of highly specific training in computational biology and molecular phylogenetics, resulting very often in a polarized accumulation of knowledge (technical on one side and biological on the other). Interpreting the robustness of genome-scale phylogenetic studies is not straightforward, particularly as new methodological developments have consistently shown that the general belief of "more genes, more robustness" often does not apply, and because there is a range of systematic errors that plague phylogenomic investigations. This is particularly problematic because phylogenomic studies are highly heterogeneous in their methodology, and best practices are often not clearly defined. The main aim of this article is to present what I consider as the ten most important points to take into consideration when planning a well-thought-out phylogenomic study and while evaluating the quality of published papers. The goal is to provide a practical step-by-step guide that can be easily followed by nonexperts and phylogenomic novices in order to assess the technical robustness of phylogenomic studies or improve the experimental design of a project.

A phylogenomic backbone for gastropod molluscs

Gastropods have survived several mass extinctions during their evolutionary history resulting in extraordinary diversity in morphology, ecology, and developmental modes, which complicate the reconstruction of a robust phylogeny. Currently, gastropods are divided into six subclasses: Caenogastropoda, Heterobranchia, Neomphaliones, Neritimorpha, Patellogastropoda, and Vetigastropoda. Phylogenetic relationships among these taxa historically lack consensus, despite numerous efforts using morphological and molecular information. We generated sequence data for transcriptomes derived from twelve taxa belonging to clades with little or no prior representation in previous studies in order to infer the deeper cladogenetic events within Gastropoda and, for the first time, infer the position of the deep-sea Neomphaliones using a phylogenomic approach. We explored the impact of missing data, homoplasy, and compositional heterogeneity on the inferred phylogenetic hypotheses. We recovered a highly supported backbone for gastropod relationships that is congruent with morphological and mitogenomic evidence, in which Patellogastropoda, true limpets, are the sister lineage to all other gastropods (Orthogastropoda) which are divided into two main clades (i) Vetigastropoda s.l. (including Pleurotomariida + Neomphaliones) and (ii) Neritimorpha + (Caenogastropoda + Heterobranchia). As such, our results support the recognition of five subclasses (or infraclasses) in Gastropoda: Patellogastropoda, Vetigastropoda, Neritimorpha, Caenogastropoda and Heterobranchia.

OUP accepted manuscript

Whole genome sequences are beginning to revolutionize our understanding of phylogenetic relationships. Yet, even whole genome sequences can fail to resolve the evolutionary history of the most rapidly radiating lineages, where incomplete lineage sorting, standing genetic variation, introgression, and other factors obscure the phylogenetic history of the group. To overcome such challenges, one emerging strategy is to integrate results across different methods. Most such approaches have been implemented on reduced representation genomic data sets, but whole genomes should provide the maximum possible evidence approach. Here, we test the ability of single nucleotide polymorphisms extracted from whole genome resequencing data, implemented in an integrative genomic approach, to resolve key nodes in the phylogeny of the mbuna, rock-dwelling cichlid fishes of Lake Malaŵi, which epitomize the phylogenetic intractability that often accompanies explosive lineage diversification. This monophyletic radiation has diversified at an unparalleled rate into several hundred species in less than 2 million years. Using an array of phylogenomic methods, we consistently recovered four major clades of mbuna, but a large basal polytomy among them. Although introgression between clades apparently contributed to the challenge of phylogenetic reconstruction, reduction of the data set to nonintrogressed sites still did not help to resolve the basal polytomy. On the other hand, relationships among six congeneric species pairs were resolved without ambiguity, even in one case where existing data led us to predict that resolution would be difficult. We conclude that the bursts of diversification at the earliest stages of the mbuna radiation may be phylogenetically unresolvable, but other regions of the tree are phylogenetically clearly supported. Integration of multiple phylogenomic approaches will continue to increase confidence in relationships inferred from these and other whole-genome data sets. [Incomplete lineage sorting; introgression; linkage disequilibrium; multispecies coalescence; rapid radiation; soft polytomy.]

A Morphological and Histological Investigation of Imperfect Lungfish Fin Regeneration

Regeneration, the replacement of body parts in a living animal, has excited scientists for centuries and our knowledge of vertebrate appendage regeneration has increased significantly over the past decades. While the ability of amniotes to regenerate body parts is very limited, members of other vertebrate clades have been shown to have rather high regenerative capacities. Among tetrapods (four-limbed vertebrates), only salamanders show unparalleled capacities of epimorphic tissue regeneration including replacement of organ and body parts in an apparently perfect fashion. The closest living relatives of Tetrapoda, the lungfish, show regenerative abilities that are comparable to those of salamanders and recent studies suggest that these high regenerative capacities may indeed be ancestral for bony fish (osteichthyans) including tetrapods. While great progress has been made in recent years in understanding the cellular and molecular mechanisms deployed during appendage regeneration, comparatively few studies have investigated gross morphological and histological features of regenerated fins and limbs. Likewise, rather little is known about how fin regeneration compares morphologically to salamander limb regeneration. In this study, we investigated the morphology and histology of regenerated fins in all three modern lungfish families. Data from histological serial sections, 3D reconstructions, and x-ray microtomography scans were analyzed to assess morphological features, quality and pathologies in lungfish fin regenerates. We found several anomalies resulting from imperfect regeneration in regenerated fins in all investigated lungfish species, including fusion of skeletal elements, additional or fewer elements, and distal branching. The similarity of patterns in regeneration abnormalities compared to salamander limb regeneration lends further support to the hypothesis that high regenerative capacities are plesiomorphic for sarcopterygians.

Phylogenomics and evolutionary history of Oreobates (Anura: Craugastoridae) Neotropical frogs along elevational gradients

Mountain ranges offer opportunities for understanding how species evolved and diversified across different environmental conditions. Neotropical frogs of the genus Oreobates (Anura: Craugastoridae) are adapted to highland and lowland habitats along the Andes, but many aspects of their evolution remain unknown. We studied their evolutionary history using ~18,000 exons enriched by targeted sequence-capture. Since capture success was very variable across samples, we evaluated to what degree differing data filtering produced robust inferences. The inferred evolutionary framework evidenced phylogenetic discordances among lowland species that can be explained by taxonomic misidentification or admixture of ancestral lineages. Highland species showed smaller effective populations than lowland frogs, probably due to greater habitat fragmentation in montane environments. Stronger genetic drift likely decreased the power of purifying selection and led to an increased proportion of nonsynonymous mutations in highland populations that could play an important role in their adaptation. Overall, our work sheds light on the evolutionary history and diversification of this group of Neotropical frogs along elevational gradients in the Andes as well as on their patterns of intraspecific diversity.

Giant lungfish genome elucidates the conquest of land by vertebrates

Lungfishes belong to lobe-fined fish (Sarcopterygii) that, in the Devonian period, 'conquered' the land and ultimately gave rise to all land vertebrates, including humans1-3. Here we determine the chromosome-quality genome of the Australian lungfish (Neoceratodus forsteri), which is known to have the largest genome of any animal. The vast size of this genome, which is about 14× larger than that of humans, is attributable mostly to huge intergenic regions and introns with high repeat content (around 90%), the components of which resemble those of tetrapods (comprising mainly long interspersed nuclear elements) more than they do those of ray-finned fish. The lungfish genome continues to expand independently (its transposable elements are still active), through mechanisms different to those of the enormous genomes of salamanders. The 17 fully assembled lungfish macrochromosomes maintain synteny to other vertebrate chromosomes, and all microchromosomes maintain conserved ancient homology with the ancestral vertebrate karyotype. Our phylogenomic analyses confirm previous reports that lungfish occupy a key evolutionary position as the closest living relatives to tetrapods4,5, underscoring the importance of lungfish for understanding innovations associated with terrestrialization. Lungfish preadaptations to living on land include the gain of limb-like expression in developmental genes such as hoxc13 and sall1 in their lobed fins. Increased rates of evolution and the duplication of genes associated with obligate air-breathing, such as lung surfactants and the expansion of odorant receptor gene families (which encode proteins involved in detecting airborne odours), contribute to the tetrapod-like biology of lungfishes. These findings advance our understanding of this major transition during vertebrate evolution.

Automated Removal of Non-homologous Sequence Stretches with PREQUAL

Large-scale multigene datasets used in phylogenomics and comparative genomics often contain sequence errors inherited from source genomes and transcriptomes. These errors typically manifest as stretches of non-homologous characters and derive from sequencing, assembly, and/or annotation errors. The lack of automatic tools to detect and remove sequence errors leads to the propagation of these errors in large-scale datasets. PREQUAL is a command line tool that identifies and masks regions with non-homologous adjacent characters in sets of unaligned homologous sequences. PREQUAL uses a full probabilistic approach based on pair hidden Markov models. On the front end, PREQUAL is user-friendly and simple to use while also allowing full customization to adjust filtering sensitivity. It is primarily aimed at amino acid sequences but can handle protein-coding nucleotide sequences. PREQUAL is computationally efficient and shows high sensitivity and accuracy. In this chapter, we briefly introduce the motivation for PREQUAL and its underlying methodology, followed by a description of basic and advanced usage, and conclude with some notes and recommendations. PREQUAL fills an important gap in the current bioinformatics tool kit for phylogenomics, contributing toward increased accuracy and reproducibility in future studies.

Sarcopterygian fin ontogeny elucidates the origin of hands with digits

How the hand and digits originated from fish fins during the Devonian fin-to-limb transition remains unsolved. Controversy in this conundrum stems from the scarcity of ontogenetic data from extant lobe-finned fishes. We report the patterning of an autopod-like domain by hoxa13 during fin development of the Australian lungfish, the most closely related extant fish relative of tetrapods. Differences from tetrapod limbs include the absence of digit-specific expansion of hoxd13 and hand2 and distal limitation of alx4 and pax9, which potentially evolved through an enhanced response to shh signaling in limbs. These developmental patterns indicate that the digit program originated in postaxial fin radials and later expanded anteriorly inside of a preexisting autopod-like domain during the evolution of limbs. Our findings provide a genetic framework for the transition of fins into limbs that supports the significance of classical models proposing a bending of the tetrapod metapterygial axis.

A mitogenomic phylogeny of chitons (Mollusca: Polyplacophora)

BACKGROUND

Polyplacophora, or chitons, have long fascinated malacologists for their distinct and rather conserved morphology and lifestyle compared to other mollusk classes. However, key aspects of their phylogeny and evolution remain unclear due to the few morphological, molecular, or combined phylogenetic analyses, particularly those addressing the relationships among the major chiton lineages.

RESULTS

Here, we present a mitogenomic phylogeny of chitons based on 13 newly sequenced mitochondrial genomes along with eight available ones and RNAseq-derived mitochondrial sequences from four additional species. Reconstructed phylogenies largely agreed with the latest advances in chiton systematics and integrative taxonomy but we identified some conflicts that call for taxonomic revisions. Despite an overall conserved gene order in chiton mitogenomes, we described three new rearrangements that might have taxonomic utility and reconstructed the most likely scenario of gene order change in this group. Our phylogeny was time-calibrated using various fossils and relaxed molecular clocks, and the robustness of these analyses was assessed with several sensitivity analyses. The inferred ages largely agreed with previous molecular clock estimates and the fossil record, but we also noted that the ambiguities inherent to the chiton fossil record might confound molecular clock analyses.

CONCLUSIONS

In light of the reconstructed time-calibrated framework, we discuss the evolution of key morphological features and call for a continued effort towards clarifying the phylogeny and evolution of chitons.

Pattern of nitrergic cells and fibers organization in the central nervous system of the Australian lungfish, Neoceratodus forsteri (Sarcopterygii: Dipnoi)

The Australian lungfish Neoceratodus forsteri is the only extant species of the order Ceratodontiformes, which retained most of the primitive features of ancient lobe finned-fishes. Lungfishes are the closest living relatives of land vertebrates and their study is important for deducing the neural traits that were conserved, modified, or lost with the transition from fishes to land vertebrates. We have investigated the nitrergic system with neural nitric oxide synthase (NOS) immunohistochemistry and NADPH-diaphorase (NADPH-d) histochemistry, which yielded almost identical results except for the primary olfactory projections and the terminal and preoptic nerve fibers labeled only for NADPH-d. Combined immunohistochemistry was used for simultaneous detection of NOS with catecholaminergic, cholinergic, and serotonergic structures, aiming to establish accurately the localization of the nitrergic elements and to assess possible interactions between these neurotransmitter systems. The results demonstrated abundant nitrergic cells in the basal ganglia, amygdaloid complex, preoptic area, basal hypothalamus, mesencephalic tectum and tegmentum, laterodorsal tegmental nucleus, reticular formation, spinal cord, and retina. In addition, low numbers of nitrergic cells were observed in the olfactory bulb, all pallial divisions, lateral septum, suprachiasmatic nucleus, prethalamic and thalamic areas, posterior tubercle, pretectum, torus semicircularis, cerebellar nucleus, interpeduncular nucleus, the medial octavolateral nucleus, nucleus of the solitary tract, and the dorsal column nucleus. Colocalization of NOS and tyrosine hydroxylase was observed in numerous cells of the ventral tegmental area/substantia nigra complex. Comparison with other vertebrates, using a neuromeric analysis, reveals that the nitrergic system of Neoceratodus shares many neuroanatomical features with tetrapods and particularly with amphibians.

Pax6 expression highlights regional organization in the adult brain of lungfishes, the closest living relatives of land vertebrates

The Pax6 gene encodes a regulatory transcription factor that is key in brain development. The molecular structure of Pax6, the roles it plays and its patterns of expression in the brain have been highly conserved during vertebrate evolution. As neurodevelopment proceeds, the Pax6 expression changes from the mitotic germinal zone in the ventricular zone to become distributed in cell groups in the adult brain. Studies in various vertebrates, from fish to mammals, found that the Pax6 expression is maintained in adults in most regions that express it during development. Specifically, in amphibians, Pax6 is widely expressed in the adult brain and its distribution pattern serves to highlight regional organization of the brain. In the present study, we analyzed the detailed distribution of Pax6 cells in the adult central nervous system of lungfishes, the closest living relatives of all tetrapods. Immunohistochemistry performed using double labeling techniques with several neuronal markers of known distribution patterns served to evaluate the actual location of Pax6 cells. Our results show that the Pax6 expression is maintained in the adult brain of lungfishes, in distinct regions of the telencephalon (pallium and subpallium), diencephalon, mesencephalon, hindbrain, spinal cord, and retina. The pattern of Pax6 expression is largely shared with amphibians and helps to understand the primitive condition that would have characterized the common ancestors to all sarcopterygians (lobe-finned fishes and tetrapods), in which Pax6 would be needed to maintain specific entities of subpopulations of neurons.

Globin E is a myoglobin-related, respiratory protein highly expressed in lungfish oocytes

Globins are a classical model system for the studies of protein evolution and function. Recent studies have shown that – besides the well-known haemoglobin and myoglobin – additional globin-types occur in vertebrates that serve different functions. Globin E (GbE) was originally identified as an eye-specific protein of birds that is distantly related to myoglobin. GbE is also present in turtles and the coelacanth but appeared to have been lost in other vertebrates. Here, we show that GbE additionally occurs in lungfish, the closest living relatives of the tetrapods. Each lungfish species harbours multiple (≥5) GbE gene copies. Surprisingly, GbE is exclusively and highly expressed in oocytes, with mRNA levels that exceed that of myoglobin in the heart. Thus, GbE is the first known oocyte-specific globin in vertebrates. No GbE transcripts were found in the ovary or egg transcriptomes of other vertebrates, suggesting a lungfish-specific function. Spectroscopic analysis and kinetic studies of recombinant GbE1 of the South American lungfish Lepidosiren paradoxa revealed a typical pentacoordinate globin with myoglobin-like O₂-binding kinetics, indicating similar functions. Our findings suggest that the multiple copies of GbE evolved to enhance O₂-supply in the developing embryo of lungfish, analogous to the embryonic and fetal haemoglobins of other vertebrates. In evolution, GbE must have changed its expression site from oocytes to eyes, or vice versa.

New patellogastropod mitogenomes help counteracting long-branch attraction in the deep phylogeny of gastropod mollusks

Long-branch attraction (LBA) is a well-known artifact in phylogenetic reconstruction. Sparse taxon sampling and extreme heterogeneity of evolutionary rates among lineages generate propitious situations for LBA, even defying probabilistic methods of phylogenetic inference. A clear example illustrating LBA challenges is the difficulty of reconstructing the deep gastropod phylogeny, particularly using mitochondrial (mt) genomes. Previous studies consistently obtained unorthodox phylogenetic relationships due to the LBA between the mitogenomes of patellogastropods (true limpets, represented only by Lottia digitalis), heterobranchs, and outgroup taxa. Here, we use the reconstruction of the gastropod mitogenomic phylogeny as a case exercise to test the effect of key methodological approaches proposed to counteract LBA, including the selection of slow-evolving representatives, the use of different outgroups, the application of site-heterogeneous evolutionary models, and the removal of fast-evolving sites. In this regard, we sequenced three new patellogastropod mt genomes, which displayed shorter branches than the one of Lottia as well as gene organizations more similar to that of the hypothetical gastropod ancestor. Phylogenetic analyses incorporating the mt genomes of Patella ferruginea, Patella vulgata, and Cellana radiata allowed eliminating the artificial clustering of Patellogastropoda and Heterobranchia that had prevailed in previous studies. Furthermore, the use of site-heterogeneous models with certain combinations of lineages within the outgroup allowed eliminating also the LBA between Heterobranchia and the outgroup, and recovering Apogastropoda (i.e., Caenogastropoda + Heterobranchia). Hence, for the first time, we were able to obtain a mitogenomic phylogeny of gastropods that is congruent with both morphological and nuclear datasets.

Global Rate Variation in Bony Vertebrates

This study investigated long-term substitution rate differences using three calibration points, divergences between lobe-finned vertebrates and ray-finned fish, between mammals and sauropsids, and between holosteans (gar and bowfin) and teleost fish with amino acid sequence data of 625 genes for 25 bony vertebrates. The result showed that the substitution rate was two to three times higher in the stem branches of lobe-finned vertebrates before the mammal-sauropsid divergence than in amniotes. The rate in the stem branch of ray-finned fish before the holostean-teleost fish divergence was also a few times higher than the holostean rate, whereas it was similar to or somewhat slower than the teleost fish rate. The phylogenetic relationship of coelacanth and lungfish with tetrapod was difficult to determine because of the short interval of the divergences. Considering the high rate in the stem branches, the divergences of coelacanth and lungfish from the stem branch were estimated as 408–427 Ma and 399–414 Ma, respectively, with the interval of 9–13 Myr. With the external calibration of the mammal-sauropsid split, the estimated times for ordinal divergences within eutherian mammals tend to be smaller than those in previous studies that used the calibration points within the lineage, with deeper divergences before the Cretaceous–Paleogene boundary and shallower ones after the boundary. In contrast the estimated times within birds were larger than those of previous studies, with the divergence between Galliformes and Anseriformes ∼80 Ma and that between Galloanserae and Neoaves 110 Ma.

A Comparative View on Sex Differentiation and Gametogenesis Genes in Lungfish and Coelacanths

Gonadal sex differentiation and reproduction are the keys to the perpetuation of favorable gene combinations and positively selected traits. In vertebrates, several gonad development features that differentiate tetrapods and fishes are likely to be, at least in part, related to the water-to-land transition. The collection of information from basal sarcopterygians, coelacanths, and lungfishes, is crucial to improve our understanding of the molecular evolution of pathways involved in reproductive functions, since these organisms are generally regarded as “living fossils” and as the direct ancestors of tetrapods. Here, we report for the first time the characterization of >50 genes related to sex differentiation and gametogenesis in Latimeria menadoensis and Protopterus annectens. Although the expression profiles of most genes is consistent with the intermediate position of basal sarcopterygians between actinopterygian fish and tetrapods, their phylogenetic placement and presence/absence patterns often reveal a closer affinity to the tetrapod orthologs. On the other hand, particular genes, for example, the male gonad factor gsdf (Gonadal Soma-Derived Factor), provide examples of ancestral traits shared with actinopterygians, which disappeared in the tetrapod lineage.

Comparative genomics and the nature of placozoan species

Placozoans are a phylum of nonbilaterian marine animals currently represented by a single described species, Trichoplax adhaerens, Schulze 1883. Placozoans arguably show the simplest animal morphology, which is identical among isolates collected worldwide, despite an apparently sizeable genetic diversity within the phylum. Here, we use a comparative genomics approach for a deeper appreciation of the structure and causes of the deeply diverging lineages in the Placozoa. We generated a high-quality draft genome of the genetic lineage H13 isolated from Hong Kong and compared it to the distantly related T. adhaerens. We uncovered substantial structural differences between the two genomes that point to a deep genomic separation and provide support that adaptation by gene duplication is likely a crucial mechanism in placozoan speciation. We further provide genetic evidence for reproductively isolated species and suggest a genus-level difference of H13 to T. adhaerens, justifying the designation of H13 as a new species, Hoilungia hongkongensis nov. gen., nov. spec., now the second described placozoan species and the first in a new genus. Our multilevel comparative genomics approach is, therefore, likely to prove valuable for species distinctions in other cryptic microscopic animal groups that lack diagnostic morphological characters, such as some nematodes, copepods, rotifers, or mites.

Regional chemoarchitecture of the brain of lungfishes based on calbindin D-28K and calretinin immunohistochemistry

Lungfishes are the closest living relatives of land vertebrates, and their neuroanatomical organization is particularly relevant for deducing the neural traits that have been conserved, modified, or lost with the transition from fishes to land vertebrates. The immunohistochemical localization of calbindin (CB) and calretinin (CR) provides a powerful method for discerning segregated neuronal populations, fiber tracts, and neuropils and is here applied to the brains of Neoceratodus and Protopterus, representing the two extant orders of lungfishes. The results showed abundant cells containing these proteins in pallial and subpallial telencephalic regions, with particular distinct distribution in the basal ganglia, amygdaloid complex, and septum. Similarly, the distribution of CB and CR containing cells supports the division of the hypothalamus of lungfishes into neuromeric regions, as in tetrapods. The dense concentrations of CB and CR positive cells and fibers highlight the extent of the thalamus. As in other vertebrates, the optic tectum is characterized by numerous CB positive cells and fibers and smaller numbers of CR cells. The so-called cerebellar nucleus contains abundant CB and CR cells with long ascending axons, which raises the possibility that it could be homologized to the secondary gustatory nucleus of other vertebrates. The corpus of the cerebellum is devoid of CB and CR and cells positive for both proteins are found in the cerebellar auricles and the octavolateralis nuclei. Comparison with other vertebrates reveals that lungfishes share most of their features of calcium binding protein distribution with amphibians, particularly with salamanders.

Back to the roots: Reducing evolutionary rate heterogeneity among sequences gives support for the early morphological hypothesis of the root of Siluriformes (Teleostei: Ostariophysi)

Catfishes (Teleostei: Siluriformes) are a highly diverse order within Ostariophysi that is distributed worldwide. At the base of this clade emerge three lineages with well-defined monophylies: Diplomystidae, Loricarioidei, and Siluroidei. Morphological phylogeny studies place the Diplomystidae as the earliest branching of these three lineages, but studies based on molecular phylogenetics consistently find the fast-evolving Loricarioidei instead. The high lineage evolutionary rate heterogeneity in this order and the fact that the lineage placed closest to the root in the molecular phylogenies is fast evolving, including many long branches, raises the possibility that the discrepancy between morphological and molecular phylogenies may be the result of a long branch attraction inference artifact. We test this hypothesis by using a 10-gene dataset to evaluate the arrangement of the three main siluriform lineages, and apply the LS³ and LS⁴ taxon sequence subsampling methods to reduce evolutionary rate heterogeneity among lineages. The initial and complete dataset supports the basal branching of Loricarioidei as in all previous molecular phylogenies, but once lineage rate heterogeneity is reduced with LS³ or LS⁴ through the removal of sequences disrupting homogeneity, the phylogeny shows Diplomystidae as the earliest branching group, with high supports, as proposed by morphological phylogeny. The result obtained with LS³, however, introduces the misplacement of one of the species with the highest amount of missing data, Scoloplax sp. Because the sequence sub-selection criterion of LS⁴ has been optimized to reduce data removal, the phylogeny resulting from the LS⁴-processed data is in agreement with the known intra-lineage relationships in addition to supporting the morphologically-based rooting hypothesis. Our results are the first instance in which a consensus between molecular and morphological phylogeny is reached concerning the root of this order.

Comparative physiology of the ductus arteriosus among vertebrates

The ductus arteriosus is typically viewed as a mammalian fetal blood vessel providing a right-to-left shunt of right ventricular outflow away from the lungs and to the systemic circuit, that must close at birth. This review provides a wider comparative examination of the ductus arteriosus in lungfish, reptiles, birds, and mammals. The ductus arteriosus evolved with the lung in the ancestors of the lungfish as a connection between the pulmonary arteries and dorsal aorta. During embryonic development, reptiles, birds, and mammals all possess either one or two paired ductus arteriosi that provide a fetal shunt of blood away from the lungs. Differences in the fetal circulatory arrangement are seen between these groups and this influences the importance of the ductus arteriosus as an embryonic shunt. The ductus arteriosus from lungfish and tetrapod vertebrates is an oxygen sensitive blood vessel, with shared conserved pathways involved in oxygen sensing. By expanding studies into more comparative models such as lungfish or developing birds a better understanding of the physiology of the ductus arteriosus can be developed.

Comparative Analysis of Nkx2.1 and Islet-1 Expression in Urodele Amphibians and Lungfishes Highlights the Pattern of Forebrain Organization in Early Tetrapods

Expression patterns of Nkx2.1 and Islet-1 (Isl1), which encode transcription factors that are key in the regionalization of the forebrain, were analyzed by combined immunohistochemical methods in young adult specimens of two lungfishes (Neoceratodus forsteri and Protopterus dolloi) and a urodele amphibian (Pleurodeles waltl). We aimed to get insights into the possible organization of the forebrain in the common ancestor of all tetrapods because of the pivotal phylogenetic significance of these two groups, being lungfishes the closest living relatives of tetrapods, and representing urodeles a model of simple brain organization with most shared features with amniotes. These transcription factors display regionally restricted expression domains in adult (juvenile) brains that are best interpreted according to the current prosomeric model. The regional patterns observed serve to identify regions and compare between the three species studied, and with previous data reported mainly for amniotes. We corroborate that Nkx2.1 and Isl1 expressions have very similar topologies in the forebrain. Common features in all sarcopterygians (lungfishes and tetrapods) have been observed, such as the Isl1 expression in most striatal neurons, whereas Nkx2.1 is restricted to migrated interneurons that reach the ventral pallium (VP). In the pallidal derivatives, the combination of both markers allows the identification of the boundaries between the ventral septum, the bed nucleus of the stria terminalis (BST) and the preoptic commissural region. In addition, the high Isl1 expression in the central amygdala (CeA), its boundary with the lateral amygdala (LA), and the scattered Nkx2.1 expression in the medial amygdala (MeA) are also shared features. The alar and basal hypothalamic territories, and the prethalamus and posterior tubercle (TP) in the diencephalon, have maintained a common pattern of expression. This regional distribution of Isl1 and Nkx2.1 observed in the forebrain of urodeles and lungfishes contributes further to our understanding of the first terrestrial vertebrates and their ancestors.

Phylogenomic Resolution of the Phylogeny of Laurasiatherian Mammals: Exploring Phylogenetic Signals within Coding and Noncoding Sequences

The interordinal relationships of Laurasiatherian mammals are currently one of the most controversial questions in mammalian phylogenetics. Previous studies mainly relied on coding sequences (CDS) and seldom used noncoding sequences. Here, by data mining public genome data, we compiled an intron data set of 3,638 genes (all introns from a protein-coding gene are considered as a gene) (19,055,073 bp) and a CDS data set of 10,259 genes (20,994,285 bp), covering all major lineages of Laurasiatheria (except Pholidota). We found that the intron data contained stronger and more congruent phylogenetic signals than the CDS data. In agreement with this observation, concatenation and species-tree analyses of the intron data set yielded well-resolved and identical phylogenies, whereas the CDS data set produced weakly supported and incongruent results. Further analyses showed that the phylogeny inferred from the intron data is highly robust to data subsampling and change in outgroup, but the CDS data produced unstable results under the same conditions. Interestingly, gene tree statistical results showed that the most frequently observed gene tree topologies for the CDS and intron data are identical, suggesting that the major phylogenetic signal within the CDS data is actually congruent with that within the intron data. Our final result of Laurasiatheria phylogeny is (Eulipotyphla,((Chiroptera, Perissodactyla),(Carnivora, Cetartiodactyla))), favoring a close relationship between Chiroptera and Perissodactyla. Our study 1) provides a well-supported phylogenetic framework for Laurasiatheria, representing a step towards ending the long-standing "hard" polytomy and 2) argues that intron within genome data is a promising data resource for resolving rapid radiation events across the tree of life.

Resolving the systematic positions of enigmatic taxa: Manipulating the chloroplast genome data of Saxifragales

Accurately resolving the phylogeny of enigmatic taxa is always a challenge in phylogenetic inference. Such uncertainties could be due to systematic errors or model violations. Here, we provide an example demonstrating how these factors affect the positioning of Paeoniaceae within Saxifragales based on chloroplast genome data. We newly assembled 14 chloroplast genomes from Saxifragales, and by combining these genomes with those of 63 other angiosperms, three datasets were assembled to test different hypotheses proposed by recent studies. These datasets were subjected to maximum parsimony, maximum likelihood and Bayesian analyses with site-homogeneous/heterogeneous models, different data partitioning strategies, and the inclusion/exclusion of weak phylogenetic signals. Three datasets exhibited remarkable heterogeneity among sites and among taxa of Saxifragales. Phylogenetic analyses under homogeneous models or maximum parsimony showed a closer relationship of Paeoniaceae with herbaceous families in the order. Data partitioning strategies did not change the general tree topology. However, PhyloBayes analysis under the CAT+GTR model resulted in a relationship closer to woody families. We conclude that although genomic data significantly increase the phylogenetic resolution of enigmatic taxa with high support, the phylogenetic results inferred from such data might be analysis or signal dependent. The analytical pipeline outlined here combines phylogenomic inference methods with evaluation of lineage-specific rates of substitution, model selection, and assessment of systematic error. These methods would be applicable to resolve similar difficult questions in the tree of life.

Phylotranscriptomic consolidation of the jawed vertebrate timetree

Phylogenomics is extremely powerful but introduces new challenges as no agreement exists on "standards" for data selection, curation and tree inference. We use jawed vertebrates (Gnathostomata) as model to address these issues. Despite considerable efforts in resolving their evolutionary history and macroevolution, few studies have included a full phylogenetic diversity of gnathostomes and some relationships remain controversial. We tested a novel bioinformatic pipeline to assemble large and accurate phylogenomic datasets from RNA sequencing and find this phylotranscriptomic approach successful and highly cost-effective. Increased sequencing effort up to ca. 10Gbp allows recovering more genes, but shallower sequencing (1.5Gbp) is sufficient to obtain thousands of full-length orthologous transcripts. We reconstruct a robust and strongly supported timetree of jawed vertebrates using 7,189 nuclear genes from 100 taxa, including 23 new transcriptomes from previously unsampled key species. Gene jackknifing of genomic data corroborates the robustness of our tree and allows calculating genome-wide divergence times by overcoming gene sampling bias. Mitochondrial genomes prove insufficient to resolve the deepest relationships because of limited signal and among-lineage rate heterogeneity. Our analyses emphasize the importance of large curated nuclear datasets to increase the accuracy of phylogenomics and provide a reference framework for the evolutionary history of jawed vertebrates.

Phylogenetic classification of bony fishes

BACKGROUND

Fish classifications, as those of most other taxonomic groups, are being transformed drastically as new molecular phylogenies provide support for natural groups that were unanticipated by previous studies. A brief review of the main criteria used by ichthyologists to define their classifications during the last 50 years, however, reveals slow progress towards using an explicit phylogenetic framework. Instead, the trend has been to rely, in varying degrees, on deep-rooted anatomical concepts and authority, often mixing taxa with explicit phylogenetic support with arbitrary groupings. Two leading sources in ichthyology frequently used for fish classifications (JS Nelson's volumes of Fishes of the World and W. Eschmeyer's Catalog of Fishes) fail to adopt a global phylogenetic framework despite much recent progress made towards the resolution of the fish Tree of Life. The first explicit phylogenetic classification of bony fishes was published in 2013, based on a comprehensive molecular phylogeny ( www.deepfin.org ). We here update the first version of that classification by incorporating the most recent phylogenetic results.

RESULTS

The updated classification presented here is based on phylogenies inferred using molecular and genomic data for nearly 2000 fishes. A total of 72 orders (and 79 suborders) are recognized in this version, compared with 66 orders in version 1. The phylogeny resolves placement of 410 families, or ~80% of the total of 514 families of bony fishes currently recognized. The ordinal status of 30 percomorph families included in this study, however, remains uncertain (incertae sedis in the series Carangaria, Ovalentaria, or Eupercaria). Comments to support taxonomic decisions and comparisons with conflicting taxonomic groups proposed by others are presented. We also highlight cases were morphological support exist for the groups being classified.

CONCLUSIONS

This version of the phylogenetic classification of bony fishes is substantially improved, providing resolution for more taxa than previous versions, based on more densely sampled phylogenetic trees. The classification presented in this study represents, unlike any other, the most up-to-date hypothesis of the Tree of Life of fishes.