Genomic signatures of evolution in Nautilus-An endangered living fossil

The Rediscovery of a Relict Unlocks the First Global Phylogeny of Whip Spiders (Amblypygi)

Asymmetrical rates of cladogenesis and extinction abound in the tree of life, resulting in numerous minute clades that are dwarfed by larger sister groups. Such taxa are commonly regarded as phylogenetic relicts or "living fossils" when they exhibit an ancient first appearance in the fossil record and prolonged external morphological stasis, particularly in comparison to their more diversified sister groups. Due to their special status, various phylogenetic relicts tend to be well-studied and prioritized for conservation. A notable exception to this trend is found within Amblypygi ("whip spiders"), a visually striking order of functionally hexapodous arachnids that are notable for their antenniform first walking leg pair (the eponymous "whips"). Paleoamblypygi, the putative sister group to the remaining Amblypygi, is known from Late Carboniferous and Eocene deposits but is survived by a single living species, Paracharon caecusHansen (1921), that was last collected in 1899. Due to the absence of genomic sequence-grade tissue for this vital taxon, there is no global molecular phylogeny for Amblypygi to date, nor a fossil-calibrated estimation of divergences within the group. Here, we report a previously unknown species of Paleoamblypygi from a cave site in Colombia. Capitalizing upon this discovery, we generated the first molecular phylogeny of Amblypygi, integrating ultraconserved element sequencing with legacy Sanger datasets and including described extant genera. To quantify the impact of sampling Paleoamblypygi on divergence time estimation, we performed in silico experiments with pruning of Paracharon. We demonstrate that the omission of relicts has a significant impact on the accuracy of node dating approaches that outweighs the impact of excluding ingroup fossils, which bears upon the ancestral range reconstruction for the group. Our results underscore the imperative for biodiversity discovery efforts in elucidating the phylogenetic relationships of "dark taxa," and especially phylogenetic relicts in tropical and subtropical habitats. The lack of reciprocal monophyly for Charontidae and Charinidae leads us to subsume them into one family, Charontidae, new synonymy.

Genomic data reveals habitat partitioning in massive Porites on Guam, Micronesia

Corals in marginal reef habitats generally exhibit less bleaching and associated mortality compared to nearby corals in more pristine reef environments. It is unclear, however, if these differences are due to environmental differences, including turbidity, or genomic differences between the coral hosts in these different environments. One particularly interesting case is in the coral genus Porites, which contains numerous morphologically similar massive Porites species inhabiting a wide range of reef habitats, from turbid river deltas and stagnant back reefs to high-energy fore reefs. Here, we generate ddRAD data for 172 Porites corals from river delta and adjacent (<0.5 km) fore reef populations on Guam to assess the extent of genetic differentiation among massive Porites corals in these two contrasting environments and throughout the island. Phylogenetic and population genomic analyses consistently identify seven different clades of massive Porites, with the two largest clades predominantly inhabiting either river deltas or fore reefs, respectively. No population structure was detected in the two largest clades, and Cladocopium was the dominant symbiont genus in all clades and environments. The perceived bleaching resilience of corals in marginal reefs may therefore be attributed to interspecific differences between morphologically similar species, in addition to potentially mediating environmental differences. Marginal reef environments may therefore not provide a suitable refuge for many reef corals in a heating world, but instead host additional cryptic coral diversity.

Evaluating the Accuracy of Methods for Detecting Correlated Rates of Molecular and Morphological Evolution

Determining the link between genomic and phenotypic change is a fundamental goal in evolutionary biology. Insights into this link can be gained by using a phylogenetic approach to test for correlations between rates of molecular and morphological evolution. However, there has been persistent uncertainty about the relationship between these rates, partly because conflicting results have been obtained using various methods that have not been examined in detail. We carried out a simulation study to evaluate the performance of 5 statistical methods for detecting correlated rates of evolution. Our simulations explored the evolution of molecular sequences and morphological characters under a range of conditions. Of the methods tested, Bayesian relaxed-clock estimation of branch rates was able to detect correlated rates of evolution correctly in the largest number of cases. This was followed by correlations of root-to-tip distances, Bayesian model selection, independent sister-pairs contrasts, and likelihood-based model selection. As expected, the power to detect correlated rates increased with the amount of data, both in terms of tree size and number of morphological characters. Likewise, greater among-lineage rate variation in the data led to improved performance of all 5 methods, particularly for Bayesian relaxed-clock analysis when the rate model was mismatched. We then applied these methods to a data set from flowering plants and did not find evidence of a correlation in evolutionary rates between genomic data and morphological characters. The results of our study have practical implications for phylogenetic analyses of combined molecular and morphological data sets, and highlight the conditions under which the links between genomic and phenotypic rates of evolution can be evaluated quantitatively.

Phylogenomic resolution of Imparidentia (Mollusca: Bivalvia) diversification through mitochondrial genomes

Despite significant advances in the phylogenomics of bivalves over the past decade, the higher-level phylogeny of Imparidentia (a superorder of Heterodonta) remains elusive. Here, a total of five new mitochondrial sequences (Chama asperella, Chama limbula, Chama dunkeri, Barnea manilensis and Ctena divergens) was added to provide resolution in nodes that required additional study. Although the monophyly of Lucinida remains less clear, the results revealed the overall backbone of the Imparidentia tree and the monophyly of Imparidentia. Likewise, most relationships among the five major Imparidentia lineages-Lucinida, Cardiida, Adapedonta, Myida and Venerida-were addressed with a well-supported topology. Basal relationships of Imparidentia recovered Lucinidae as the sister group to all remaining imparidentian taxa. Thyasiridae is a sister group to other imparidentian bivalves (except Lucinidae species) which is split into Cardiida, Adapedonta and the divergent clade of Neoheterodontei. Neoheterodontei was comprised of Venerida and Myida, the former of which now also contains Chamidae as the sister group to all the remaining venerid taxa. Moreover, molecular divergence times were inferred by calibrating nine nodes in the Imparidentia tree of life by extinct taxa. The origin of these major clades ranged from Ordovician to Permian with the diversification through the Palaeozoic to Mesozoic. Overall, the results obtained in this study demonstrate a better-resolved Imparidentia phylogeny based on mitochondrial genomes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00178-x.

Three new species of Nautilus Linnaeus, 1758 (Mollusca, Cephalopoda) from the Coral Sea and South Pacific

Nautiloids are a charismatic group of marine molluscs best known for their rich fossil record, but today they are restricted to a handful of species in the family Nautilidae from around the Coral Triangle. Recent genetic work has shown a disconnect between traditional species, originally defined on shell characters, but now with new findings from genetic structure of various Nautilus populations. Here, three new species of Nautilus from the Coral Sea and South Pacific region are formally named using observations of shell and soft anatomical data augmented by genetic information: N. samoaensissp. nov. (from American Samoa), N. vitiensissp. nov. (from Fiji), and N. vanuatuensissp. nov. (from Vanuatu). The formal naming of these three species is timely considering the new and recently published information on genetic structure, geographic occurrence, and new morphological characters, including color patterns of shell and soft part morphology of hood, and will aid in managing these possibly endangered animals. As recently proposed from genetic analyses, there is a strong geographic component affecting taxonomy, with the new species coming from larger island groups that are separated by at least 200 km of deep water (greater than 800 m) from other Nautilus populations and potential habitats. Nautilid shells implode at depths greater than 800 m and depth therefore acts as a biogeographical barrier separating these species. This isolation, coupled with the unique, endemic species in each locale, are important considerations for the conservation management of the extant Nautilus species and populations.

Genomic diversity and Demographic History of the Dromiciops genus (Marsupialia: Microbiotheriidae)

Three orders represent the South American fauna of marsupials. Of these, Microbiotheria was until recently known as a monotypic genus with the only surviving species Dromiciops gliroides (monito del monte). The recent proposal of a new Dromiciops species (Dromiciops bicinovici), together with new information on the origin and diversification of living microbioterians has changed the prevailing paradigm around the evolutionary history of these emblematic marsupials. Here, we used a RADseq approach to test for evidence of admixture and past or current gene flow among both species of Dromiciops and evaluate the genetic structure within D. gliroides. We analyzed 127 samples of Dromiciops distributed across the known distribution range of both species. We also inferred the joint demographic history of these lineages, thus corroborating the status of D. bozinovici as a distinct species. Demographic history reconstruction indicated that D. bozinovici diverged from D. gliroides around 4my ago and has remained isolated and demographically stable ever since. In contrast, D. gliroides is subdivided into three subclades that experienced recent expansions and moderate gene flow among them (mostly from north to south). Furthermore, genetic distances among populations within D. gliroides were significantly correlated with geographic distances. These results suggest that some of the D. gliroides populations would have survived in glacial refuges, with posterior expansions after ice retreat. Our results have important implications for the systematics of the genus and have profound conservation consequences for the new species, especially considering the fragmentation level of the temperate rainforest.

Significance of the suture line in cephalopod taxonomy revealed by 3D morphometrics in the modern nautilids Nautilus and Allonautilus

Assessing the taxonomic importance of the suture line in shelled cephalopods is a key to better understanding the diversity of this group in Earth history. Because fossils are subject to taphonomic artifacts, an in-depth knowledge of well-preserved modern organisms is needed as an important reference. Here, we examine the suture line morphology of all known species of the modern cephalopods Nautilus and Allonautilus. We applied computed tomography and geometric morphometrics to quantify the suture line morphology as well as the conch geometry and septal spacing. Results reveal that the suture line and conch geometry are useful in distinguishing species, while septal spacing is less useful. We also constructed cluster trees to illustrate the similarity among species. The tree based on conch geometry in middle ontogeny is nearly congruent with those previously reconstructed based on molecular data. In addition, different geographical populations of the same species of Nautilus separate out in this tree. This suggests that genetically distinct (i.e., geographically isolated) populations of Nautilus can also be distinguished using conch geometry. Our results are applicable to closely related fossil cephalopods (nautilids), but may not apply to more distantly related forms (ammonoids).

Genome-wide SNPs redefines species boundaries and conservation units in the freshwater mussel genus Cyprogenia of North America

Detailed information on species delineation and population genetic structure is a prerequisite for designing effective restoration and conservation strategies for imperiled organisms. Phylogenomic and population genomic analyses based on genome-wide double digest restriction-site associated DNA sequencing (ddRAD-Seq) data has identified three allopatric lineages in the North American freshwater mussel genus Cyprogenia. Cyprogenia stegaria is restricted to the Eastern Highlands and displays little genetic structuring within this region. However, two allopatric lineages of C. aberti in the Ozark and Ouachita highlands exhibit substantial levels (mean uncorrected F_ST = 0.368) of genetic differentiation and each warrants recognition as a distinct evolutionary lineage. Lineages of Cyprogenia in the Ouachita and Ozark highlands are further subdivided reflecting structuring at the level of river systems. Species tree inference and species delimitation in a Bayesian framework using single nucleotide polymorphisms (SNP) data supported results from phylogenetic analyses, and supports three species of Cyprogenia over the currently recognized two species. A comparison of SNPs generated from both destructively and non-destructively collected samples revealed no significant difference in the SNP error rate, quality and amount of ddRAD sequence reads, indicating that nondestructive or trace samples can be effectively utilized to generate SNP data for organisms for which destructive sampling is not permitted.

Phylogenomic Resolution of Sea Spider Diversification through Integration of Multiple Data Classes

Despite significant advances in invertebrate phylogenomics over the past decade, the higher-level phylogeny of Pycnogonida (sea spiders) remains elusive. Due to the inaccessibility of some small-bodied lineages, few phylogenetic studies have sampled all sea spider families. Previous efforts based on a handful of genes have yielded unstable tree topologies. Here, we inferred the relationships of 89 sea spider species using targeted capture of the mitochondrial genome, 56 conserved exons, 101 ultraconserved elements, and 3 nuclear ribosomal genes. We inferred molecular divergence times by integrating morphological data for fossil species to calibrate 15 nodes in the arthropod tree of life. This integration of data classes resolved the basal topology of sea spiders with high support. The enigmatic family Austrodecidae was resolved as the sister group to the remaining Pycnogonida and the small-bodied family Rhynchothoracidae as the sister group of the robust-bodied family Pycnogonidae. Molecular divergence time estimation recovered a basal divergence of crown group sea spiders in the Ordovician. Comparison of diversification dynamics with other marine invertebrate taxa that originated in the Paleozoic suggests that sea spiders and some crustacean groups exhibit resilience to mass extinction episodes, relative to mollusk and echinoderm lineages.

Genomic Signatures of Selection along a Climatic Gradient in the Northern Range Margin of the White-Footed Mouse (Peromyscus leucopus)

Identifying genetic variation involved in thermal adaptation is likely to yield insights into how species adapt to different climates. Physiological and behavioral responses associated with overwintering (e.g., torpor) are thought to serve important functions in climate adaptation. In this study, we use 2 isolated Peromyscus leucopus lineages on the northern margin of the species range to identify single nucleotide polymorphisms (SNPs) showing a strong environmental association and test for evidence of parallel evolution. We found signatures of clinal selection in each lineage, but evidence of parallelism was limited, with only 2 SNPs showing parallel allele frequencies across transects. These parallel SNPs map to a gene involved in protection against iron-dependent oxidative stress (Fxn) and to a gene with unknown function but containing a forkhead-associated domain (Fhad1). Furthermore, within transects, we find significant clinal patterns in genes enriched for functions associated with glycogen homeostasis, synaptic function, intracellular Ca2+ balance, H3 histone modification, as well as the G2/M transition of cell division. Our results are consistent with recent literature on the cellular and molecular basis of climate adaptation in small mammals and provide candidate genomic regions for further study.

Genetic differentiation in mountain-dwelling clam shrimp, Paralimnadia (Crustacea : Branchiopoda : Spinicaudata), in eastern Australia

The majority of Australian Spinicaudata Linder, 1945 inhabit the (semi)arid deserts of Australia’s lowlands. However, several closely related species of Paralimnadia Sars, 1896 inhabit small temporary habitats throughout the Great Dividing Range in eastern Australia. By combining analyses of mitochondrial cytochrome c oxidase subunit I (COI) with double-digest restriction-site associated DNA (ddRAD) data, we studied the species diversity and genetic diversity of this group of mountain-dwelling branchiopods. Levels of genetic differentiation in COI are relatively low between putative species (mostly between 1.5 and 6.7%), complicating COI-based species delimitation. Depending on the applied threshold, three to six species are inferred in the studied area, with most putative species being geographically restricted. Particularly notable are the high levels of population differentiation indicated by ddRAD analyses between nearby populations within putative species. This suggests that gene flow is limited, even between populations separated only by a few kilometres. This may lead to fast population differentiation, which in turn might drive speciation. Our data suggest that the species diversity of Paralimnadia in the Great Dividing Range is much higher than currently appreciated.

Resequencing 545 ginkgo genomes across the world reveals the evolutionary history of the living fossil

As Charles Darwin anticipated, living fossils provide excellent opportunities to study evolutionary questions related to extinction, competition, and adaptation. Ginkgo (Ginkgo biloba L.) is one of the oldest living plants and a fascinating example of how people have saved a species from extinction and assisted its resurgence. By resequencing 545 genomes of ginkgo trees sampled from 51 populations across the world, we identify three refugia in China and detect multiple cycles of population expansion and reduction along with glacial admixture between relict populations in the southwestern and southern refugia. We demonstrate multiple anthropogenic introductions of ginkgo from eastern China into different continents. Further analyses reveal bioclimatic variables that have affected the geographic distribution of ginkgo and the role of natural selection in ginkgo’s adaptation and resilience. These investigations provide insights into the evolutionary history of ginkgo trees and valuable genomic resources for further addressing various questions involving living fossil species.

Ginkgo is one of the living fossils from the plant kingdom. Here, authors conduct population genomics analyses to reveal its refugia and demographic history, and provide evidence of multiple anthropogenic introductions of ginkgo from eastern China into different continents.

Evidence that DNA repair genes, a family of tumor suppressor genes, are associated with evolution rate and size of genomes

Adaptive radiation and evolutionary stasis are characterized by very different evolution rates. The main aim of this study was to investigate if any genes have a special role to a high or low evolution rate. The availability of animal genomes permitted comparison of gene content of genomes of 24 vertebrate species that evolved through adaptive radiation (representing high evolutionary rate) and of 20 vertebrate species that are considered as living fossils (representing a slow evolutionary rate or evolutionary stasis). Mammals, birds, reptiles, and bony fishes were included in the analysis. Pathway analysis was performed for genes found to be specific in adaptive radiation or evolutionary stasis respectively. Pathway analysis revealed that DNA repair and cellular response to DNA damage are important (false discovery rate = 8.35 × 10⁻⁵; 7.15 × 10⁻⁶, respectively) for species evolved through adaptive radiation. This was confirmed by further genetic in silico analysis (p = 5.30 × 10⁻³). Nucleotide excision repair and base excision repair were the most significant pathways. Additionally, the number of DNA repair genes was found to be linearly related to the genome size and the protein number (proteome) of the 44 animals analyzed (p < 1.00 × 10⁻⁴), this being compatible with Drake’s rule. This is the first study where radiated and living fossil species have been genetically compared. Evidence has been found that cancer-related genes have a special role in radiated species. Linear association of the number of DNA repair genes with the species genome size has also been revealed. These comparative genetics results can support the idea of punctuated equilibrium evolution.

Postmortem transport in fossil and modern shelled cephalopods

The chambered shells of cephalopod mollusks, such as modern Nautilus and fossil ammonoids, have the potential to float after death, which could result in significant postmortem transport of shells away from living habitats. Such transport would call into question these clades' documented biogeographic distributions and therefore the many (paleo)biological interpretations based on them. It is therefore imperative to better constrain the likelihood and extent of postmortem transport in modern and fossil cephalopods. Here, I combine the results of classic experiments on postmortem buoyancy with datasets on cephalopod shell form to determine that only those shells with relatively high inflation are likely to float for a significant interval after death and therefore potentially experience postmortem transport. Most ammonoid cephalopods have shell forms making postmortem transport unlikely. Data on shell forms and geographic ranges of early Late Cretaceous cephalopod genera demonstrate that even genera with shell forms conducive to postmortem buoyancy do not, in fact, show artificially inflated biogeographic ranges relative to genera with non-buoyant morphologies. Finally, georeferenced locality data for living nautilid specimens and dead drift shells indicate that most species have relatively small geographic ranges and experience limited drift. Nautilus pompilius is the exception, with a broad Indo-Pacific range and drift shells found far from known living populations. Given the similarity of N. pompilius to other nautilids in its morphology and ecology, it seems unlikely that this species would have a significantly different postmortem fate than its close relatives. Rather, it is suggested that drift shells along the east African coast may indicate the existence of modern (or recently extirpated) living populations of nautilus in the western Indian Ocean, which has implications for the conservation of these cephalopods.

Using Genomic Data to Infer Historic Population Dynamics of Nonmodel Organisms

Genome sequence data are now being routinely obtained from many nonmodel organisms. These data contain a wealth of information about the demographic history of the populations from which they originate. Many sophisticated statistical inference procedures have been developed to infer the demographic history of populations from this type of genomic data. In this review, we discuss the different statistical methods available for inference of demography, providing an overview of the underlying theory and logic behind each approach. We also discuss the types of data required and the pros and cons of each method. We then discuss how these methods have been applied to a variety of nonmodel organisms. We conclude by presenting some recommendations for researchers looking to use genomic data to infer demographic history.