Genome sequences reveal global dispersal routes and suggest convergent genetic adaptations in seahorse evolution

A Fast, Reproducible, High-throughput Variant Calling Workflow for Population Genomics

The increasing availability of genomic resequencing data sets and high-quality reference genomes across the tree of life present exciting opportunities for comparative population genomic studies. However, substantial challenges prevent the simple reuse of data across different studies and species, arising from variability in variant calling pipelines, data quality, and the need for computationally intensive reanalysis. Here, we present snpArcher, a flexible and highly efficient workflow designed for the analysis of genomic resequencing data in nonmodel organisms. snpArcher provides a standardized variant calling pipeline and includes modules for variant quality control, data visualization, variant filtering, and other downstream analyses. Implemented in Snakemake, snpArcher is user-friendly, reproducible, and designed to be compatible with high-performance computing clusters and cloud environments. To demonstrate the flexibility of this pipeline, we applied snpArcher to 26 public resequencing data sets from nonmammalian vertebrates. These variant data sets are hosted publicly to enable future comparative population genomic analyses. With its extensibility and the availability of public data sets, snpArcher will contribute to a broader understanding of genetic variation across species by facilitating the rapid use and reuse of large genomic data sets.

Gene loss and co-option of toll-like receptors facilitate paternal immunological adaptation in the brood pouch of pregnant male seahorses

Male pregnancy in syngnathids (seahorses, pipefishes, and sea dragons) is an evolutionary innovation in the animal kingdom. Paternal immune resistance to the fetus is a critical challenge, particularly in seahorses with fully enclosed brood pouches and sophisticated placentas. In this study, comparative genomic analysis revealed that all syngnathid species lost three vertebrate-conserved Toll-like receptors (TLR1, TLR2, and TLR9), of which all play essential roles in immune protection and immune tolerance in the uterus and placenta. Quantitative real-time PCR (qRT-PCR) analysis showed that the TLR paralog genes including TLR18, TLR25, and TLR21 were highly expressed in the placenta inside the seahorse brood pouch and changed dynamically during the breeding cycle, suggesting the potentially important role of the TLRs during male pregnancy. Furthermore, the immune challenge test in vitro showed a remarkable expression response from all three TLR genes to specific pathogenic antigens, confirming their immune function in seahorse brood pouches. Notably, the altered antigen recognition spectrum of these genes appeared to functionally compensate in part for the lost TLRs, in contrast to that observed in other species. Therefore, we suggest that gene loss and co-option of TLRs may be a typical evolutionary strategy for facilitating paternal immunological adaptation during male pregnancy.

The genetic basis of the leafy seadragon's unique camouflage morphology and avenues for its efficient conservation derived from habitat modeling

The leafy seadragon certainly is among evolution's most "beautiful and wonderful" species aptly named for its extraordinary camouflage mimicking its coastal seaweed habitat. However, limited information is known about the genetic basis of its phenotypes and conspicuous camouflage. Here, we revealed genomic signatures of rapid evolution and positive selection in core genes related to its camouflage, which allowed us to predict population dynamics for this species. Comparative genomic analysis revealed that seadragons have the smallest olfactory repertoires among all ray-finned fishes, suggesting adaptations to the highly specialized habitat. Other positively selected and rapidly evolving genes that serve in bone development and coloration are highly expressed in the leaf-like appendages, supporting a recent adaptive shift in camouflage appendage formation. Knock-out of bmp6 results in dysplastic intermuscular bones with a significantly reduced number in zebrafish, implying its important function in bone formation. Global climate change-induced loss of seagrass beds now severely threatens the continued existence of this enigmatic species. The leafy seadragon has a historically small population size likely due to its specific habitat requirements that further exacerbate its vulnerability to climate change. Therefore, taking climate change-induced range shifts into account while developing future protection strategies.

First record of the Great Seahorse Hippocampus kelloggi Jordan & Snyder, 1901 (Actinopterygii: Syngnathiformes: Syngnathidae) from the northwestern coast of Bay of Bengal

This study reports the first occurrence of the Great Seahorse Hippocampus kelloggi from the state of Odisha in the eastern coast (Bay of Bengal) of India. The seahorse was captured in a ring net during daily fishing activities. The sample was collected from the Ariyapalli fish landing center and identification was carried out based on the morphometric features of the specimen and the seahorse identification guide. The total length of the juvenile seahorse was 12.5 cm (with a tail length of 6.6 cm (52.8%), trunk length of 3.4 cm (27.2%) and head length of 2.5 cm (20%)). The length of the snout was 1 cm. There were 38 tail rings followed by 11 rings on the trunk of the animal. Both eye and cheek spines were present. Northward migration (~1,300 km) of this species can be a response of extensive fishing activities around the southern coast of India. This calls for increased monitoring of the coast coastal ecosystems of India on the east coast for better conservation and management of the remaining seahorse populations.

Independent Evolution of Sex Chromosomes and Male Pregnancy-Related Genes in Two Seahorse Species

Unlike birds and mammals, many teleosts have homomorphic sex chromosomes, and changes in the chromosome carrying the sex-determining locus, termed "turnovers", are common. Recent turnovers allow studies of several interesting questions. One question is whether the new sex-determining regions evolve to become completely non-recombining, and if so, how and why. Another is whether (as predicted) evolutionary changes that benefit one sex accumulate in the newly sex-linked region. To study these questions, we analyzed the genome sequences of two seahorse species of the Syngnathidae, a fish group in which many species evolved a unique structure, the male brood pouch. We find that both seahorse species have XY sex chromosome systems, but their sex chromosome pairs are not homologs, implying that at least one turnover event has occurred. The Y-linked regions occupy 63.9% and 95.1% of the entire sex chromosome of the two species and do not exhibit extensive sequence divergence with their X-linked homologs. We find evidence for occasional recombination between the extant sex chromosomes that may account for their homomorphism. We argue that these Y-linked regions did not evolve by recombination suppression after the turnover, but by the ancestral nature of the low crossover rates in these chromosome regions. With such an ancestral crossover landscape, a turnover can instantly create an extensive Y-linked region. Finally, we test for adaptive evolution of male pouch-related genes after they became Y-linked in the seahorse.

Immunogenetic losses co-occurred with seahorse male pregnancy and mutation in tlx1 accompanied functional asplenia

In the highly derived syngnathid fishes (pipefishes, seadragons & seahorses), the evolution of sex-role reversed brooding behavior culminated in the seahorse lineage's male pregnancy, whose males feature a specialized brood pouch into which females deposit eggs during mating. Then, eggs are intimately engulfed by a placenta-like tissue that facilitates gas and nutrient exchange. As fathers immunologically tolerate allogenic embryos, it was suggested that male pregnancy co-evolved with specific immunological adaptations. Indeed, here we show that a specific amino-acid replacement in the tlx1 transcription factor is associated with seahorses' asplenia (loss of spleen, an organ central in the immune system), as confirmed by a CRISPR-Cas9 experiment using zebrafish. Comparative genomics across the syngnathid phylogeny revealed that the complexity of the immune system gene repertoire decreases as parental care intensity increases. The synchronous evolution of immunogenetic alterations and male pregnancy supports the notion that male pregnancy co-evolved with the immunological tolerance of the embryo.

Leafy and weedy seadragon genomes connect genic and repetitive DNA features to the extravagant biology of syngnathid fishes

Seadragons are widely recognized for their derived traits, which include leaf-like appendages and extreme spinal curvature. Efforts to understand the genetic basis of these unique traits and conserve these species and their relatives have been limited by genomic resource gaps. In this paper we present full, annotated genomes of leafy and weedy seadragons, which we use to uncover surprising features of gene family and genome architecture evolution that likely relate to the extravagant phenotypic traits of seadragons and their pipefish and seahorse relatives. These genomes and their analyses are important advances for the study of elaborate vertebrate traits, leveraging this diverse, morphologically exceptional group of fishes.

Seadragons are a remarkable lineage of teleost fishes in the family Syngnathidae, renowned for having evolved male pregnancy. Comprising three known species, seadragons are widely recognized and admired for their fantastical body forms and coloration, and their specific habitat requirements have made them flagship representatives for marine conservation and natural history interests. Until recently, a gap has been the lack of significant genomic resources for seadragons. We have produced gene-annotated, chromosome-scale genome models for the leafy and weedy seadragon to advance investigations of evolutionary innovation and elaboration of morphological traits in seadragons as well as their pipefish and seahorse relatives. We identified several interesting features specific to seadragon genomes, including divergent noncoding regions near a developmental gene important for integumentary outgrowth, a high genome-wide density of repetitive DNA, and recent expansions of transposable elements and a vesicular trafficking gene family. Surprisingly, comparative analyses leveraging the seadragon genomes and additional syngnathid and outgroup genomes revealed striking, syngnathid-specific losses in the family of fibroblast growth factors (FGFs), which likely involve reorganization of highly conserved gene regulatory networks in ways that have not previously been documented in natural populations. The resources presented here serve as important tools for future evolutionary studies of developmental processes in syngnathids and hold value for conservation of the extravagant seadragons and their relatives.

Phylogenomic analysis of Syngnathidae reveals novel relationships, origins of endemic diversity and variable diversification rates

Background

Seahorses, seadragons, pygmy pipehorses, and pipefishes (Syngnathidae, Syngnathiformes) are among the most recognizable groups of fishes because of their derived morphology, unusual life history, and worldwide distribution. Despite previous phylogenetic studies and recent new species descriptions of syngnathids, the evolutionary relationships among several major groups within this family remain unresolved.

Results

Here, we provide a reconstruction of syngnathid phylogeny based on genome-wide sampling of 1314 ultraconserved elements (UCEs) and expanded taxon sampling to assess the current taxonomy and as a basis for macroevolutionary insights. We sequenced a total of 244 new specimens across 117 species and combined with published UCE data for a total of 183 species of Syngnathidae, about 62% of the described species diversity, to compile the most data-rich phylogeny to date. We estimated divergence times using 14 syngnathiform fossils, including nine fossils with newly proposed phylogenetic affinities, to better characterize current and historical biogeographical patterns, and to reconstruct diversification through time. We present a phylogenetic hypothesis that is well-supported and provides several notable insights into syngnathid evolution. We found nine non-monophyletic genera, evidence for seven cryptic species, five potentially invalid synonyms, and identified a novel sister group to the seahorses, the Indo-Pacific pipefishes Halicampus macrorhynchus and H. punctatus. In addition, the morphologically distinct southwest Pacific seahorse Hippocampus jugumus was recovered as the sister to all other non-pygmy seahorses. As found in many other groups, a high proportion of syngnathid lineages appear to have originated in the Central Indo-Pacific and subsequently dispersed to adjoining regions. Conversely, we also found an unusually high subsequent return of lineages from southern Australasia to the Central Indo-Pacific. Diversification rates rose abruptly during the Middle Miocene Climate Transition and peaked after the closure of the Tethys Sea.

Conclusions

Our results reveal a previously underappreciated diversity of syngnathid lineages. The observed biogeographic patterns suggest a significant role of the southern Australasian region as a source and sink of lineages. Shifts in diversification rates imply possible links to declining global temperatures, the separation of the Atlantic and Pacific faunas, and the environmental changes associated with these events.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01271-w.

A highly contiguous nuclear genome assembly of the mandarinfish Synchiropus splendidus (Syngnathiformes: Callionymidae)

The fish order Syngnathiformes has been referred to as a collection of misfit fishes, comprising commercially important fish such as red mullets as well as the highly diverse seahorses, pipefishes, and seadragons-the well-known family Syngnathidae, with their unique adaptations including male pregnancy. Another ornate member of this order is the species mandarinfish. No less than two types of chromatophores have been discovered in the spectacularly colored mandarinfish: the cyanophore (producing blue color) and the dichromatic cyano-erythrophore (producing blue and red). The phylogenetic position of mandarinfish in Syngnathiformes, and their promise of additional genetic discoveries beyond the chromatophores, made mandarinfish an appealing target for whole-genome sequencing. We used linked sequences to create synthetic long reads, producing a highly contiguous genome assembly for the mandarinfish. The genome assembly comprises 483 Mbp (longest scaffold 29 Mbp), has an N50 of 12 Mbp, and an L50 of 14 scaffolds. The assembly completeness is also high, with 92.6% complete, 4.4% fragmented, and 2.9% missing out of 4584 BUSCO genes found in ray-finned fishes. Outside the family Syngnathidae, the mandarinfish represents one of the most contiguous syngnathiform genome assemblies to date. The mandarinfish genomic resource will likely serve as a high-quality outgroup to syngnathid fish, and furthermore for research on the genomic underpinnings of the evolution of novel pigmentation.

Genome of the estuarine oyster provides insights into climate impact and adaptive plasticity

Understanding the roles of genetic divergence and phenotypic plasticity in adaptation is central to evolutionary biology and important for assessing adaptive potential of species under climate change. Analysis of a chromosome-level assembly and resequencing of individuals across wide latitude distribution in the estuarine oyster (Crassostrea ariakensis) revealed unexpectedly low genomic diversity and population structures shaped by historical glaciation, geological events and oceanographic forces. Strong selection signals were detected in genes responding to temperature and salinity stress, especially of the expanded solute carrier families, highlighting the importance of gene expansion in environmental adaptation. Genes exhibiting high plasticity showed strong selection in upstream regulatory regions that modulate transcription, indicating selection favoring plasticity. Our findings suggest that genomic variation and population structure in marine bivalves are heavily influenced by climate history and physical forces, and gene expansion and selection may enhance phenotypic plasticity that is critical for the adaptation to rapidly changing environments.

Genome and gene evolution of seahorse species revealed by the chromosome-level genome of Hippocampus abdominalis

Seahorses belong to the teleost family Syngnathidae that evolved a distinct body plan and unique male pregnancy compared to other teleosts. As a classic model for studying evolution of viviparity and sexual selection of teleosts, seahorse species still lack a publicly available high‐quality reference genome. Here, we generated the genome assembly of the big‐belly seahorse, Hippocampus abdominalis with long‐read and Hi‐C technologies. We managed to place over 99% of the total length of 444.7 Mb of assembled genome into 21 linkage groups with almost no gaps. We reconstructed a phylogenomic tree with the big‐belly seahorse genome and other representative Syngnathidae and teleost species. We also reconstructed the historical population dynamics of four representative Syngnathidae species. We found the gene families that underwent expansion or contraction in the Syngnathidae ancestor were enriched for immune‐related or ion transporter gene ontology terms. Many of these genes were also reported to show a dynamic expression pattern during the pregnancy stages of H. abdominalis. We also identified putative positively selected genes in the Syngnathidae ancestor or in H. abdominalis, whose mouse mutants are enriched for abnormal craniofacial and limb morphological phenotypes. Overall, our study provides an important genome resource for evolutionary and developmental studies of seahorse species, and candidate genes for future experimental works.

Reinventing the wheel? Reassessing the roles of gene flow, sorting and convergence in repeated evolution

Biologists have long been intrigued by apparently predictable and repetitive evolutionary trajectories inferred across a variety of lineages and systems. In recent years, high-throughput sequencing analyses have started to transform our understanding of such repetitive shifts. While researchers have traditionally categorized such shifts as either "convergent" or "parallel," based on relatedness of the lineages involved, emerging genomic insights provide an opportunity to better describe the actual evolutionary mechanisms at play. A synthesis of recent genomic analyses confirms that convergence is the predominant driver of repetitive evolution among species, whereas repeated sorting of standing variation is the major driver of repeated shifts within species. However, emerging data reveal numerous notable exceptions to these expectations, with recent examples of de novo mutations underpinning convergent shifts among even very closely related lineages, while repetitive sorting processes have occurred among even deeply divergent taxa, sometimes via introgression. A number of very recent analyses have found evidence for both processes occurring on different scales within taxa. We suggest that the relative importance of convergent versus sorting processes depends on the interplay between gene flow among populations, and phylogenetic relatedness of the lineages involved.

Phylogenomics and Historical Biogeography of Seahorses, Dragonets, Goatfishes, and Allies (Teleostei: Syngnatharia): Assessing Factors Driving Uncertainty in Biogeographic Inferences

The charismatic trumpetfishes, goatfishes, dragonets, flying gurnards, seahorses, and pipefishes encompass a recently defined yet extraordinarily diverse clade of percomorph fishes-the series Syngnatharia. This group is widely distributed in tropical and warm-temperate regions, with a great proportion of its extant diversity occurring in the Indo-Pacific. Because most syngnatharians feature long-range dispersal capabilities, tracing their biogeographic origins is challenging. Here, we applied an integrative phylogenomic approach to elucidate the evolutionary biogeography of syngnatharians. We built upon a recently published phylogenomic study that examined ultraconserved elements by adding 62 species (total 169 species) and one family (Draconettidae), to cover ca. 25% of the species diversity and all 10 families in the group. We inferred a set of time-calibrated trees and conducted ancestral range estimations. We also examined the sensitivity of these analyses to phylogenetic uncertainty (estimated from multiple genomic subsets), area delimitation, and biogeographic models that include or exclude the jump-dispersal parameter (j). Of the three factors examined, we found that the j parameter has the strongest effect in ancestral range estimates, followed by number of areas defined, and tree topology and divergence times. After accounting for these uncertainties, our results reveal that syngnatharians originated in the ancient Tethys Sea ca. 87 Ma (84-94 Ma; Late Cretaceous) and subsequently occupied the Indo-Pacific. Throughout syngnatharian history, multiple independent lineages colonized the eastern Pacific (6-8 times) and the Atlantic (6-14 times) from their center of origin, with most events taking place following an east-to-west route prior to the closure of the Tethys Seaway ca. 12-18 Ma. Ultimately, our study highlights the importance of accounting for different factors generating uncertainty in macroevolutionary and biogeographic inferences.