Multiple global radiations in tadpole shrimps challenge the concept of 'living fossils'

Mechanisms of Rapid Karyotype Evolution in Mammals

Chromosome reshuffling events are often a foundational mechanism by which speciation can occur, giving rise to highly derivative karyotypes even amongst closely related species. Yet, the features that distinguish lineages prone to such rapid chromosome evolution from those that maintain stable karyotypes across evolutionary time are still to be defined. In this review, we summarize lineages prone to rapid karyotypic evolution in the context of Simpson's rates of evolution-tachytelic, horotelic, and bradytelic-and outline the mechanisms proposed to contribute to chromosome rearrangements, their fixation, and their potential impact on speciation events. Furthermore, we discuss relevant genomic features that underpin chromosome variation, including patterns of fusions/fissions, centromere positioning, and epigenetic marks such as DNA methylation. Finally, in the era of telomere-to-telomere genomics, we discuss the value of gapless genome resources to the future of research focused on the plasticity of highly rearranged karyotypes.

Comparative genomics of Hox and ParaHox genes among major lineages of Branchiopoda with emphasis on tadpole shrimps

Hox and ParaHox genes (HPHGs) are key developmental genes that pattern regional identity along the anterior–posterior body axis of most animals. Here, we identified HPHGs in tadpole shrimps (Pancrustacea, Branchiopoda, Notostraca), an iconic example of the so-called “living fossils” and performed a comparative genomics analysis of HPHGs and the Hox cluster among major branchiopod lineages. Notostraca possess the entire Hox complement, and the Hox cluster seems to be split into two different subclusters, although we were not able to support this finding with chromosome-level assemblies. However, the genomic structure of Hox genes in Notostraca appears more derived than that of Daphnia spp., which instead retains the plesiomorphic condition of a single compact cluster. Spinicaudata and Artemia franciscana show instead a Hox cluster subdivided across two or more genomic scaffolds with some orthologs either duplicated or missing. Yet, branchiopod HPHGs are similar among the various clades in terms of both intron length and number, as well as in their pattern of molecular evolution. Sequence substitution rates are in fact generally similar for most of the branchiopod Hox genes and the few differences we found cannot be traced back to natural selection, as they are not associated with any signals of diversifying selection or substantial switches in selective modes. Altogether, these findings do not support a significant stasis in the Notostraca Hox cluster and further confirm how morphological evolution is not tightly associated with genome dynamics.

Comparative genomics of tadpole shrimps (Crustacea, Branchiopoda, Notostraca): Dynamic genome evolution against the backdrop of morphological stasis

This analysis presents five genome assemblies of four Notostraca taxa. Notostraca origin dates to the Permian/Upper Devonian and the extant forms show a striking morphological similarity to fossil taxa. The comparison of sequenced genomes with other Branchiopoda genomes shows that, despite the morphological stasis, Notostraca share a dynamic genome evolution with high turnover for gene families' expansion/contraction and a transposable elements content comparable to other branchiopods. While Notostraca substitutions rate appears similar or lower in comparison to other branchiopods, a subset of genes shows a faster evolutionary pace, highlighting the difficulty of generalizing about genomic stasis versus dynamism. Moreover, we found that the variation of Triops cancriformis transposable elements content appeared linked to reproductive strategies, in line with theoretical expectations. Overall, besides providing new genomic resources for the study of these organisms, which appear relevant for their ecology and evolution, we also confirmed the decoupling of morphological and molecular evolution.

Phylogeny and Biogeography of Spinicaudata (Crustacea: Branchiopoda)

Spinicaudata (spiny clam shrimp) is a taxon of Branchiopoda occurring since the Devonian and today it occurs nearly globally in temporary water bodies. We present the most species-rich phylogenetic analyses of this taxon based on four molecular loci: COI, 16S rRNA, EF1α and 28S rRNA. Our results support previous findings that Cyzicidae sensu lato is paraphyletic. To render Cyzicidae monophyletic we establish a fourth extant spinicaudatan family to accommodate Eocyzicus. Within Cyzicidae, none of the genera Cyzicus, Caenestheria or Caenestheriella are monophyletic, and the morphological characters used to define these genera (condyle length and rostrum shape) are not associated with well-delimited clades within Cyzicidae. There is insufficient resolution to elucidate the relationships within Leptestheriidae. However, there is sufficient evidence to show that the leptestheriid genera Eoleptestheria and Leptestheria are non-monophyletic, and there is no support for the genus Leptestheriella. Molecular clock analyses suggest that the wide geographic distribution of many spinicaudatan taxa across multiple continents is largely based on vicariance associated with the break-up of Pangea and Gondwana. Trans-oceanic dispersal has occurred in some taxa (e.g., Eulimnadia and within Leptestheriidae) but has been relatively rare. Our results highlight the need to revise the taxonomy of Cyzicidae and Leptestheriidae and provide evidence that the global spinicaudatan diversity may be underestimated due to the presence of numerous cryptic species. We establish Eocyzicidae fam. nov. to accommodate the genus Eocyzicus. Consequently, Cyzicidae comprises only two genera -Cyzicus and Ozestheria. Ozestheria occurs also in Africa and Asia and Ozestheria pilosa new comb. is assigned to this genus.

Characterization of the complete mitochondrial genome of Chinese Triops granarius and implications for species delimitation

The mt genome of T. granarius (15.121 bp), containing 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and one control region, present the typical pattern of phyllopod mt genomes from the perspective of gene content and order. All PCGs initiate with standard start codon of ATN and share the complete stop codon of TAA or TAG, whereas nd1 initiates with GTG, nd2 terminates with AGA and nd4 terminates with a single T. The secondary structures of rrnL and rrnS comprise of 48 helices and 32 helices, respectively. The patterns of evolutionary rate heterogeneity were observed at different levels in Triops. Phylogenetic analyses were performed based on the nucleotide sequences and amino acid sequences showed that Asian T. granarius was a sister group to (T. longicaudatus + T. australiensis), and T. cancriformis was the basal clade of Triops. Furthermore, we investigated the genetic diversity and phylogenetic relationships of different T. granarius s.l. populations. Ten phylogenetic lineages are separated on species level by mean maximum likelihood distance ≥ 12.2% in the cox1 gene. While the six populations occurring in Inner Mongolia of China, Mongolia and Russia, investigated for haplotypes belonged to the same haplotype networks.

Evolution of the Torso activation cassette, a pathway required for terminal patterning and moulting

Embryonic terminal patterning and moulting are critical developmental processes in insects. In Drosophila and Tribolium both of these processes are regulated by the Torso-activation cassette (TAC). The TAC consists of a common receptor, Torso, ligands Trunk and prothoracicotropic hormone (PTTH), and the spatially restricted protein Torso-like, with combinations of these elements acting mechanistically to activate the receptor in different developmental contexts. In order to trace the evolutionary history of the TAC we determined the presence or absence of TAC components in the genomes of arthropods. Our analyses reveal that Torso, Trunk and PTTH are evolutionarily labile components of the TAC with multiple individual or combined losses occurring in the arthropod lineages leading to and within the insects. These losses are often correlated, with both ligands and receptor missing from the genome of the same species. We determine that the PTTH gene evolved in the common ancestor of Hemiptera and Holometabola, and is missing from the genomes of a number of species with experimentally demonstrated PTTH activity, implying another molecule may be involved in ecdysis in these species. In contrast, the torso-like gene is a common component of pancrustacean genomes.

Mitochondrial genome diversity and evolution in Branchiopoda (Crustacea)

BACKGROUND

The crustacean class Branchiopoda includes fairy shrimps, clam shrimps, tadpole shrimps, and water fleas. Branchiopods, which are well known for their great variety of reproductive strategies, date back to the Cambrian and extant taxa can be mainly found in freshwater habitats, also including ephemeral ponds. Mitochondrial genomes of the notostracan taxa Lepidurus apus lubbocki (Italy), L. arcticus (Iceland) and Triops cancriformis (an Italian and a Spanish population) are here characterized for the first time and analyzed together with available branchiopod mitogenomes.

RESULTS

Overall, branchiopod mitogenomes share the basic structure congruent with the ancestral Pancrustacea model. On the other hand, rearrangements involving tRNAs and the control region are observed among analyzed taxa. Remarkably, an unassigned region in the L. apus lubbocki mitogenome showed a chimeric structure, likely resulting from a non-homologous recombination event between the two flanking trnC and trnY genes. Notably, Anostraca and Onychocaudata mitogenomes showed increased GC content compared to both Notostraca and the common ancestor, and a significantly higher substitution rate, which does not correlate with selective pressures, as suggested by dN/dS values.

CONCLUSIONS

Branchiopod mitogenomes appear rather well-conserved, although gene rearrangements have occurred. For the first time, it is reported a putative non-homologous recombination event involving a mitogenome, which produced a pseudogenic tRNA sequence. In addition, in line with data in the literature, we explain the higher substitution rate of Anostraca and Onychocaudata with the inferred GC substitution bias that occurred during their evolution.

Rethinking Living Fossils

Biologists would be mistaken if they relegated living fossils to paleontological inquiry or assumed that the concept is dead. It is now used to describe entities ranging from viruses to higher taxa, despite recent warnings of misleading inferences. Current work on character evolution illustrates how analyzing living fossils and stasis in terms of parts (characters) and wholes (e.g., organisms and lineages) advances our understanding of prolonged stasis at many hierarchical levels. Instead of viewing the concept's task as categorizing living fossils, we show how its primary role is to mark out what is in need of explanation, accounting for the persistence of both molecular and morphological traits. Rethinking different conceptions of living fossils as specific hypotheses reveals novel avenues for research that integrate phylogenetics, ecological and evolutionary modeling, and evo-devo to produce a more unified theoretical outlook.

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.

A molecular palaeobiological exploration of arthropod terrestrialization

Understanding animal terrestrialization, the process through which animals colonized the land, is crucial to clarify extant biodiversity and biological adaptation. Arthropoda (insects, spiders, centipedes and their allies) represent the largest majority of terrestrial biodiversity. Here we implemented a molecular palaeobiological approach, merging molecular and fossil evidence, to elucidate the deepest history of the terrestrial arthropods. We focused on the three independent, Palaeozoic arthropod terrestrialization events (those of Myriapoda, Hexapoda and Arachnida) and showed that a marine route to the colonization of land is the most likely scenario. Molecular clock analyses confirmed an origin for the three terrestrial lineages bracketed between the Cambrian and the Silurian. While molecular divergence times for Arachnida are consistent with the fossil record, Myriapoda are inferred to have colonized land earlier, substantially predating trace or body fossil evidence. An estimated origin of myriapods by the Early Cambrian precedes the appearance of embryophytes and perhaps even terrestrial fungi, raising the possibility that terrestrialization had independent origins in crown-group myriapod lineages, consistent with morphological arguments for convergence in tracheal systems.This article is part of the themed issue 'Dating species divergences using rocks and clocks'.

Transition in sexual system and sex chromosome evolution in the tadpole shrimp Triops cancriformis

Transitions in sexual system and reproductive mode may affect the course of sex chromosome evolution, for instance by altering the strength of sexually antagonistic selection. However, there have been few studies of sex chromosomes in systems where such transitions have been documented. The European tadpole shrimp, Triops cancriformis, has undergone a transition from dioecy to androdioecy (a sexual system where hermaphrodites and males coexist), offering an excellent opportunity to test the impact of this transition on the evolution of sex chromosomes. To identify sex-linked markers, to understand mechanisms of sex determination and to investigate differences between sexual systems, we carried out a genome-wide association study using restriction site-associated DNA sequencing (RAD-seq) of 47 males, females and hermaphrodites from one dioecious and one androdioecious population. We analysed 22.9 Gb of paired-end sequences and identified and scored >3000 high coverage novel genomic RAD markers. Presence–absence of markers, single-nucleotide polymorphism association and read depth identified 52 candidate sex-linked markers. We show that sex is genetically determined in T. cancriformis, with a ZW system conserved across dioecious and androdioecious populations and that hermaphrodites have likely evolved from females. We also show that the structure of the sex chromosomes differs strikingly, with a larger sex-linked region in the dioecious population compared with the androdioecious population.

Identification, expression, and molecular evolution of microRNAs in the "living fossil" Triops cancriformis (tadpole shrimp)

MicroRNAs have been identified and analyzed in various model species, but an investigation of miRNAs in nonmodel species is required for a more complete understanding of miRNA evolution. In this study, we investigated the miRNAs of the nonmodel species Triops cancriformis (tadpole shrimp), a "living fossil," whose morphological form has not changed in almost 200 million years. Dramatic ontogenetic changes occur during its development. To clarify the evolution of miRNAs, we comparatively analyzed its miRNAs and the components of its RNAi machinery. We used deep sequencing to analyze small RNA libraries from the six different developmental stages of T. cancriformis (egg, first-fourth instars, and adult), and also analyzed its genomic DNA with deep sequencing. We identified 180 miRNAs (87 conserved miRNAs and 93 novel candidate miRNAs), and deduced the components of its RNAi machinery: the DICER1, AGO1-3, PIWI, and AUB proteins. A comparative miRNA analysis of T. cancriformis and Drosophila melanogaster showed inconsistencies in the expression patterns of four conserved miRNAs. This suggests that although the miRNA sequences of the two species are very similar, their roles differ across the species. An miRNA conservation analysis revealed that most of the conserved T. cancriformis miRNAs share sequence similarities with those of arthropods, although T. cancriformis is called a "living fossil." However, we found that let-7 and DICER1 of T. cancriformis are more similar to those of the vertebrates than to those of the arthropods. These results suggest that miRNA systems of T. cancriformis have evolved in a unique fashion.

Boom and bust: ancient and recent diversification in bichirs (Polypteridae: Actinopterygii), a relictual lineage of ray-finned fishes

Understanding the history that underlies patterns of species richness across the Tree of Life requires an investigation of the mechanisms that not only generate young species-rich clades, but also those that maintain species-poor lineages over long stretches of evolutionary time. However, diversification dynamics that underlie ancient species-poor lineages are often hidden due to a lack of fossil evidence. Using information from the fossil record and time calibrated molecular phylogenies, we investigate the history of lineage diversification in Polypteridae, which is the sister lineage of all other ray-finned fishes (Actinopterygii). Despite originating at least 390 million years (Myr) ago, molecular timetrees support a Neogene origin for the living polypterid species. Our analyses demonstrate polypterids are exceptionally species depauperate with a stem lineage duration that exceeds 380 million years (Ma) and is significantly longer than the stem lineage durations observed in other ray-finned fish lineages. Analyses of the fossil record show an early Late Cretaceous (100.5-83.6 Ma) peak in polypterid genus richness, followed by 60 Ma of low richness. The Neogene species radiation and evidence for high-diversity intervals in the geological past suggest a "boom and bust" pattern of diversification that contrasts with common perceptions of relative evolutionary stasis in so-called "living fossils."

Molecular phylogeny of the Notostraca

We used a combined analysis of one nuclear (28S rDNA) and three mitochondrial markers (COI, 12S rDNA, 16S rDNA) to infer the molecular phylogeny of the Notostraca, represented by samples from the six continents that are inhabited by this group of branchiopod crustaceans. Our results confirm the monophyly of both extant notostracan genera Triops and Lepidurus with good support in model based and maximum parsimony analyses. We used branchiopod fossils as a calibration to infer divergence times among notostracan lineages and accounted for rate heterogeneity among lineages by applying relaxed-clock models. Our divergence date estimates indicate an initial diversification into the genera Triops and Lepidurus in the Mesozoic, most likely at a minimum age of 152.3-233.5 Ma, i.e., in the Triassic or Jurassic. Implications for the interpretation of fossils and the evolution of notostracan morphology are discussed. We further use the divergence date estimates to formulate a biogeographic hypothesis that explains distributions of extant lineages predominantly by overland dispersal routes. We identified an additional hitherto unrecognised highly diverged lineage within Lepidurus apus lubbocki and three additional previously unknown major lineages within Triops. Within T. granarius we found deep differentiation, with representatives distributed among three major phylogenetic lineages. One of these major lineages comprises T. cancriformis, the T. mauritanicus species group and two hitherto unrecognised T. granarius lineages. Samples that were morphologically identified as T. granarius diverged from the most basal nodes within this major lineage, and divergence dates suggested an approximate age of 23.7-49.6 Ma for T. cancriformis, indicating the need for a taxonomic revision of Triassic and Permian fossils that are currently attributed to the extant T. cancriformis. We thus elevate T. cancriformis minor to full species status as Triops minorTrusheim, 1938 and include in this species the additional Upper Triassic samples that were attributed to T. cancriformis. We further elevate T. cancriformis permiensis to full species status as Triops permiensisGand et al., 1997.

Stable phylogenetic patterns in scutigeromorph centipedes (Myriapoda : Chilopoda : Scutigeromorpha): dating the diversification of an ancient lineage of terrestrial arthropods

Although stable and well-supported relationships are in place for the three main clades (families) of Scutigeromorpha, the interrelationships of particular taxa within the most diverse family, Scutigeridae, are less clearly resolved. Novel molecular data for taxa from Mesoamerica, the Caribbean, southern Africa, New Guinea and previously unsampled parts of the Pacific are incorporated into phylogenetic analyses. Relationships across the tree are stable under variable analytical conditions, whether these are homology-based (multiple sequence alignment versus implied alignment; untrimmed versus trimmed datasets) or method-based (parsimony versus maximum likelihood). Hypervariable regions, contrary to common belief, add phylogenetic structure to the data, as measured by the increased support for many nodes when compared with the same alignments trimmed with Gblocks. Our analyses show that a Yule-3-rate model best explained the diversification of Scutigeromorpha during their 400 million years of history. More complete molecular data for the New Guinea genus Ballonema stabilise its position as sister group to Thereuoneminae. To reconcile scutigeromorph systematics with the phylogeny, the monotypic genus Madagassophora Verhoeff, 1936, is placed in synonymy with Scutigerina Silvestri, 1901 (n. syn.), its type species M. hova becoming Scutigerina hova (de Saussure & Zehntner, 1902) new comb. (from Scutigera), and Lassophora Verhoeff, 1905, is re-established for an Afro-Malagasy clade containing Lassophora nossibei (de Saussure & Zehntner, 1902) new comb. (from Scutigera) and a newly sequenced species from Mozambique that diverged at the base of the lineage to Thereuoneminae. The dated phylogeny of Scutigeromorpha is more consistent with ancient vicariant splits between Madagascar–southern Africa and Australia–New Caledonia than with younger dispersal scenarios, though some geologically young Pacific islands that harbour lineages dating to the Cretaceous demonstrate the potential for trans-oceanic dispersal.