Mudskipper genomes provide insights into the terrestrial adaptation of amphibious fishes

Klumpy: A tool to evaluate the integrity of long-read genome assemblies and illusive sequence motifs

The improvement and decreasing costs of third-generation sequencing technologies has widened the scope of biological questions researchers can address with de novo genome assemblies. With the increasing number of reference genomes, validating their integrity with minimal overhead is vital for establishing confident results in their applications. Here, we present Klumpy, a tool for detecting and visualizing both misassembled regions in a genome assembly and genetic elements (e.g. genes) of interest in a set of sequences. By leveraging the initial raw reads in combination with their respective genome assembly, we illustrate Klumpy's utility by investigating antifreeze glycoprotein (afgp) loci across two icefishes, by searching for a reported absent gene in the northern snakehead fish, and by scanning the reference genomes of a mudskipper and bumblebee for misassembled regions. In the two former cases, we were able to provide support for the noncanonical placement of an afgp locus in the icefishes and locate the missing snakehead gene. Furthermore, our genome scans were able identify an unmappable locus in the mudskipper reference genome and identify a putative repetitive element shared among several species of bees.

Genomic insights into the seawater adaptation in Cyprinidae

BACKGROUND

Cyprinidae, the largest fish family, encompasses approximately 367 genera and 3006 species. While they exhibit remarkable adaptability to diverse aquatic environments, it is exceptionally rare to find them in seawater, with the Far Eastern daces being of few exceptions. Therefore, the Far Eastern daces serve as a valuable model for studying the genetic mechanisms underlying seawater adaptation in Cyprinidae.

RESULTS

Here, we sequenced the chromosome-level genomes of two Far Eastern daces (Pseudaspius brandtii and P. hakonensis), the two known cyprinid fishes found in seawater, and performed comparative genomic analyses to investigate their genetic mechanism of seawater adaptation. Demographic history reconstruction of the two species reveals that their population dynamics are correlated with the glacial-interglacial cycles and sea level changes. Genomic analyses identified Pseudaspius-specific genetic innovations related to seawater adaptation, including positively selected genes, rapidly evolving genes, and conserved non-coding elements (CNEs). Functional assays of Pseudaspius-specific variants of the prolactin (prl) gene showed enhanced cell adaptation to greater osmolarity. Functional assays of Pseudaspius specific CNEs near atg7 and usp45 genes suggest that they exhibit higher promoter activity and significantly induced at high osmolarity.

CONCLUSIONS

Our results reveal the genome-wide evidence for the evolutionary adaptation of cyprinid fishes to seawater, offering valuable insights into the molecular mechanisms supporting the survival of migratory fish in marine environments. These findings are significant as they contribute to our understanding of how cyprinid fishes navigate and thrive in diverse aquatic habitats, providing useful implications for the conservation and management of marine ecosystems.

Genomics comparisons of three chromosome-level mudskipper genome assemblies reveal molecular clues for water-to-land evolution and adaptation

INTRODUCTION

Mudskippers are a large group of amphibious fishes that have developed many morphological and physiological capacities to live on land. Genomics comparisons of chromosome-level genome assemblies of three representative mudskippers, Boleophthalmus pectinirostris (BP), Periophthalmus magnuspinnatus (PM) and P. modestus (PMO), may be able to provide novel insights into the water-to-land evolution and adaptation.

METHODS

Two chromosome-level genome assemblies for BP and PM were respectively sequenced by an integration of PacBio, Nanopore and Hi-C sequencing. A series of standard assembly and annotation pipelines were subsequently performed for both mudskippers. We also re-annotated the PMO genome, downloaded from NCBI, to obtain a redundancy-reduced annotation. Three-way comparative analyses of the three mudskipper genomes in a large scale were carried out to discover detailed genomic differences, such as different gene sizes, and potential chromosomal fission and fusion events. Comparisons of several representative gene families among the three amphibious mudskippers and some other teleosts were also performed to find some molecular clues for terrestrial adaptation.

RESULTS

We obtained two high-quality haplotype genome assemblies with 23 and 25 chromosomes for BP and PM respectively. We also found two specific chromosome fission events in PM. Ancestor chromosome analysis has discovered a common fusion event in mudskipper ancestor. This fusion was then retained in all the three mudskipper species. A loss of some SCPP (secretory calcium-binding phosphoprotein) genes were identified in the three mudskipper genomes, which could lead to reduction of scales for a part-time terrestrial residence. The loss of aanat1a gene, encoding an important enzyme (arylalkylamine N-acetyltransferase 1a, AANAT1a) for dopamine metabolism and melatonin biosynthesis, was confirmed in PM but not in PMO (as previously reported existence in BP), suggesting a better air vision of PM than both PMO and BP. Such a tiny variation within the genus Periophthalmus exemplifies to prove a step-by-step evolution for the mudskippers' water-to-land adaptation.

CONCLUSION

These high-quality mudskipper genome assemblies will become valuable genetic resources for in-depth discovery of genomic evolution for the terrestrial adaptation of amphibious fishes.

Chromosome-level genome assembly and annotation of eel goby (Odontamblyopus rebecca)

The eel gobies fascinate researchers with many important features, including its unique body structure, benthic lifestyle, and degenerated eyes. However, genome assembly and exploration of the unique genomic composition of the eel gobies are still in their infancy. This has severely limited research progress on gobies. In this study, multi-platform sequencing data were generated and used to assemble and annotate the genome of O. rebecca at the chromosome-level. The assembled genome size of O. rebecca is 918.57 Mbp, which is similar to the estimated genome size (903.03 Mbp) using 17-mer. The scaffold N50 is 41.67 Mbp, and 23 chromosomes were assembled using Hi-C technology with a mounting rate of 99.96%. Genome annotation indicates that 53.29% of the genome is repetitive sequences, and 22,999 protein-coding genes are predicted, of which 21,855 have functional annotations. The chromosome-level genome of O. rebecca will not only provide important genomic resources for comparative genomic studies of gobies, but also expand our knowledge of the genetic origin of their unique features fascinating researchers for decades.

Whole-Genome Sequencing Analyses Reveal the Whip-like Tail Formation, Innate Immune Evolution, and DNA Repair Mechanisms of Eupleurogrammus muticus

Smallhead hairtail (Eupleurogrammus muticus) is an important marine economic fish distributed along the northern Indian Ocean and the northwest Pacific coast; however, little is known about the mechanism of its genetic evolution. This study generated the first genome assembly of E. muticus at the chromosomal level using a combination of PacBio SMRT, Illumina Nova-Seq, and Hi-C technologies. The final assembled genome size was 709.27 Mb, with a contig N50 of 25.07 Mb, GC content of 40.81%, heterozygosity rate of 1.18%, and repetitive sequence rate of 35.43%. E. muticus genome contained 21,949 protein-coding genes (97.92% of the genes were functionally annotated) and 24 chromosomes. There were 143 expansion gene families, 708 contraction gene families, and 4888 positively selected genes in the genome. Based on the comparative genomic analyses, we screened several candidate genes and pathways related to whip-like tail formation, innate immunity, and DNA repair in E. muticus. These findings preliminarily reveal some molecular evolutionary mechanisms of E. muticus at the genomic level and provide important reference genomic data for the genetic studies of other trichiurids.

Genomic landscape of introgression from the ghost lineage in a gobiid fish uncovers the generality of forces shaping hybrid genomes

Extinct lineages can leave legacies in the genomes of extant lineages through ancient introgressive hybridization. The patterns of genomic survival of these extinct lineages provide insight into the role of extinct lineages in current biodiversity. However, our understanding on the genomic landscape of introgression from extinct lineages remains limited due to challenges associated with locating the traces of unsampled 'ghost' extinct lineages without ancient genomes. Herein, we conducted population genomic analyses on the East China Sea (ECS) lineage of Chaenogobius annularis, which was suspected to have originated from ghost introgression, with the aim of elucidating its genomic origins and characterizing its landscape of introgression. By combining phylogeographic analysis and demographic modelling, we demonstrated that the ECS lineage originated from ancient hybridization with an extinct ghost lineage. Forward simulations based on the estimated demography indicated that the statistic γ of the HyDe analysis can be used to distinguish the differences in local introgression rates in our data. Consistent with introgression between extant organisms, we found reduced introgression from extinct lineage in regions with low recombination rates and with functional importance, thereby suggesting a role of linked selection that has eliminated the extinct lineage in shaping the hybrid genome. Moreover, we identified enrichment of repetitive elements in regions associated with ghost introgression, which was hitherto little known but was also observed in the re-analysis of published data on introgression between extant organisms. Overall, our findings underscore the unexpected similarities in the characteristics of introgression landscapes across different taxa, even in cases of ghost introgression.

Synergistic modulation of neuroendocrine-inflammation pathway by microRNAs facilitates intertidal adaptation of molluscs

Neuroendocrine-immune system is an evolution-conserved regulatory network in maintaining the homeostasis of animals. While knowledge on the roles of neuroendocrine-immune system in the disease and stress responses of organisms is growing, the ecological roles of neuroendocrine-immune system, especially how it shapes the unique lifestyle of organisms remain insufficiently investigated. As an endemic and dominant mollusc in intertidal region, oysters have evolved with a primitive neuroendocrine-immune system and with a sessile lifestyle. Recently, a novel neuroendocrine-immune pathway, Ca2+/calmodulin (CaM)-nitrite oxide synthase (NOS)/nitrite oxide (NO)-tumor necrosis factor (TNF) pathway, is identified in oysters and found altered dynamically during aerial exposure, one common but challenging stresses for intertidal organisms and a decisive factor shaping their habitat. Since the pathway proves fatal in prolonged aerial exposure, we hypothesized that the activation/deactivation of pathway could be strictly modulated in adaptation to the sessile lifestyle of oysters. Here, a synergistic modulation on the Ca2+/CaM-NOS/NO-TNF pathway by four members of miR-92 family and two oyster-specific miRNAs was identified, which further hallmarks the resilience and survival strategy of oysters to aerial exposure. Briefly, these six miRNAs were down-regulating CgCaM24243 post-transcriptionally and deactivating the pathway during the early-stage of stress. However, a robust recession of these miRNAs occurred at the late-stage of stress, resulting in the reactivation of pathway and overwhelming accumulation of cytokines. These results demonstrated a complicated interaction between miRNAs and ancient neuroendocrine-immune system, which facilitates the environmental adaptation of intertidal oysters and provides novel insight on the function and evolution of neuroendocrine-immune system in ecological context.

The translation initiation factor homolog eif4e1c regulates cardiomyocyte metabolism and proliferation during heart regeneration

The eIF4E family of translation initiation factors bind 5' methylated caps and act as the limiting step for mRNA translation. The canonical eIF4E1A is required for cell viability, yet other related eIF4E families exist and are utilized in specific contexts or tissues. Here, we describe a family called Eif4e1c, for which we find roles during heart development and regeneration in zebrafish. The Eif4e1c family is present in all aquatic vertebrates but is lost in all terrestrial species. A core group of amino acids shared over 500 million years of evolution forms an interface along the protein surface, suggesting that Eif4e1c functions in a novel pathway. Deletion of eif4e1c in zebrafish caused growth deficits and impaired survival in juveniles. Mutants surviving to adulthood had fewer cardiomyocytes and reduced proliferative responses to cardiac injury. Ribosome profiling of mutant hearts demonstrated changes in translation efficiency of mRNA for genes known to regulate cardiomyocyte proliferation. Although eif4e1c is broadly expressed, its disruption had most notable impact on the heart and at juvenile stages. Our findings reveal context-dependent requirements for translation initiation regulators during heart regeneration.

Ancient fish lineages illuminate toll-like receptor diversification in early vertebrate evolution

Since its initial discovery over 50 years ago, understanding the evolution of the vertebrate RAG- mediated adaptive immune response has been a major area of research focus for comparative geneticists. However, how the evolutionary novelty of an adaptive immune response impacted the diversity of receptors associated with the innate immune response has received considerably less attention until recently. Here, we investigate the diversification of vertebrate toll-like receptors (TLRs), one of the most ancient and well conserved innate immune receptor families found across the Tree of Life, integrating genomic data that represent all major vertebrate lineages with new transcriptomic data from Polypteriformes, the earliest diverging ray-finned fish lineage. Our analyses reveal TLR sequences that reflect the 6 major TLR subfamilies, TLR1, TLR3, TLR4, TLR5, TLR7, and TLR11, and also currently unnamed, yet phylogenetically distinct TLR clades. We additionally recover evidence for a pulse of gene gain coincident with the rise of the RAG-mediated adaptive immune response in jawed vertebrates, followed by a period of rapid gene loss during the Cretaceous. These gene losses are primarily concentrated in marine teleost fish and synchronous with the mid Cretaceous anoxic event, a period of rapid extinction for marine species. Finally, we reveal a mismatch between phylogenetic placement and gene nomenclature for up to 50% of TLRs found in clades such as ray-finned fishes, cyclostomes, amphibians, and elasmobranchs. Collectively, these results provide an unparalleled perspective of TLR diversity and offer a ready framework for testing gene annotations in non-model species.

Genome Sequencing Provides Novel Insights into Mudflat Burrowing Adaptations in Eel Goby Taenioides sp. (Teleost: Amblyopinae)

Amblyopinae is one of the lineage of bony fish that preserves amphibious traits living in tidal mudflat habitats. In contrast to other active amphibious fish, Amblyopinae species adopt a seemly more passive lifestyle by living in deep burrows of mudflat to circumvent the typical negative effects associated with terrestriality. However, little is known about the genetic origin of these mudflat deep-burrowing adaptations in Amblyopinae. Here we sequenced the first genome of Amblyopinae species, Taenioides sp., to elucidate their mudflat deep-burrowing adaptations. Our results revealed an assembled genome size of 774.06 Mb with 23 pseudochromosomes anchored, which predicted 22,399 protein-coding genes. Phylogenetic analyses indicated that Taenioides sp. diverged from the active amphibious fish of mudskipper approximately 28.3 Ma ago. In addition, 185 and 977 putative gene families were identified to be under expansion, contraction and 172 genes were undergone positive selection in Taenioides sp., respectively. Enrichment categories of top candidate genes under significant expansion and selection were mainly associated with hematopoiesis or angiogenesis, DNA repairs and the immune response, possibly suggesting their involvement in the adaptation to the hypoxia and diverse pathogens typically observed in mudflat burrowing environments. Some carbohydrate/lipid metabolism, and insulin signaling genes were also remarkably alterated, illustrating physiological remolding associated with nutrient-limited subterranean environments. Interestingly, several genes related to visual perception (e.g., crystallins) have undergone apparent gene losses, pointing to their role in the small vestigial eyes development in Taenioides sp. Our work provide valuable resources for understanding the molecular mechanisms underlying mudflat deep-burrowing adaptations in Amblyopinae, as well as in other tidal burrowing teleosts.

Skeletal anatomy of the pectoral fin in mudskipper species from terrestrial and aquatic habitats

Mudskippers are a group of amphibious fishes in the family Oxudercidae, whose species inhabit a range of habitats from mostly aquatic to mostly terrestrial. Most of our understanding about habitat preference comes from natural history observations, particularly where they are collected (i.e., low intertidal vs. high intertidal regions). Mudskippers have undergone several morphological changes to accommodate a terrestrial life, including major changes to the pectoral and pelvic girdles. These changes result in a novel crutching gait, which mudskippers use to move over land. Though the appendicular morphology and crutching gait of mudskippers have been described in some species, few studies have compared skeletal structures across the family. In our study, we use microcomputed tomography (µCT) scans to compare the skeletal anatomy of 16 species of aquatic and terrestrial mudskippers. Linear discriminant analysis is used to analyze measurements obtained through geometric morphometrics (landmarks). We found bone structures of the pectoral region in the terrestrial group were significantly longer and wider than those in the aquatic group. Furthermore, a significant difference in anatomy is shown between terrestrial and aquatic genera with both axial and appendicular elements contributing to the separation between groups. This work describes the differences in skeletal morphology associated with terrestriality in mudskippers and provides valuable insights into specific anatomical characteristics contributing to their adaptation to novel environments.

Molecular evolution of the hemoglobin gene family across vertebrates

Adaptation to various altitudes and oxygen levels is a major aspect of vertebrate evolution. Hemoglobin is an erythrocyte protein belonging to the globin superfamily, and the α-, β-globin genes of jawed vertebrates encode tetrameric ((α₂β₂) hemoglobin, which contributes to aerobic metabolism by delivering oxygen from the respiratory exchange surfaces into cells. However, there are various gaps in knowledge regarding hemoglobin gene evolution, including patterns in cartilaginous fish and the roles of gene conversion in various taxa. Hence, we evaluated the evolutionary history of the vertebrate hemoglobin gene family by analyses of 97 species representing all classes of vertebrates. By genome-wide analyses, we extracted 879 hemoglobin sequences. Members of the hemoglobin gene family were conserved in birds and reptiles but variable in mammals, amphibians, and teleosts. Gene motifs, structures, and synteny were relatively well-conserved among vertebrates. Our results revealed that purifying selection contributed substantially to the evolution of all vertebrate hemoglobin genes, with mean d_N/d_S (ω) values ranging from 0.057 in teleosts to 0.359 in reptiles. In general, after the fish-specific genome duplication, the teleost hemoglobin genes showed variation in rates of evolution, and the β-globin genes showed relatively high ω values after a gene transposition event in amniotes. We also observed that the frequency of gene conversion was high in amniotes, with fewer hemoglobin genes and higher rates of evolution. Collectively, our findings provide detail insight into complex evolutionary processes shaping the vertebrate hemoglobin gene family, involving gene duplication, gene loss, purifying selection, and gene conversion.

A progestin regulates the prostaglandin pathway in the neuroendocrine system in female mudskipper Boleophthalmus pectinirostris

Sex hormones regulate the reproductive cycle through brain-pituitary axis, but the molecular mechanism is still enigmatic. In the reproductive season, the mudskipper Boleophthalmus pectinirostris possesses a semilunar periodicity spawning rhythm, which coincides with the semilunar periodicity variations in 17α-hydroxyprogesterone, the precursor of 17α,20β-dihydroxy-4-pregnen-3-one (DHP), a sexual progestin in teleosts. In the present study, we investigated in vitro the brain transcriptional differences between DHP-treated tissues and control groups using RNA-seq. Differential expression analysis revealed that 2700 genes significantly differentially expressed, including 1532 up-regulated and 1168 down-regulated genes. The majority of prostaglandin pathway-related genes were dramatically up-regulated, especially the prostaglandin receptor 6 (ptger6). Tissue distribution analysis revealed that ptger6 gene was ubiquitously expressed. In situ hybridization results showed that ptger6, nuclear progestin receptor (pgr), and DHP-induced c-fos mRNA were co-expressed in the ventral telencephalic area, the ventral nucleus of ventral telencephalic area, the anterior part of parvocellular preoptic nucleus, the magnocellular part of magnocellular preoptic nucleus, the ventral zone of periventricular hypothalamus, the anterior tubercular nucleus, the periventricular nucleus of posterior tuberculum, and the torus longitudinalis. DHP significantly enhanced promoter activities of ptger6 via Pgr. Together, this study suggested that DHP regulates the prostaglandin pathway in the neuroendocrine system of teleost fish.

Amblyopinae Mitogenomes Provide Novel Insights into the Paraphyletic Origin of Their Adaptation to Mudflat Habitats

The water-to-land transition is one of the most important events in evolutionary history of vertebrates. However, the genetic basis underlying many of the adaptations during this transition remains unclear. Mud-dwelling gobies in the subfamily Amblyopinae are one of the teleosts lineages that show terrestriality and provide a useful system for clarifying the genetic changes underlying adaptations to terrestrial life. Here, we sequenced the mitogenome of six species in the subfamily Amblyopinae. Our results revealed a paraphyletic origin of Amblyopinae with respect to Oxudercinae, which are the most terrestrial fishes and lead an amphibious life in mudflats. This partly explains the terrestriality of Amblyopinae. We also detected unique tandemly repeated sequences in the mitochondrial control region in Amblyopinae, as well as in Oxudercinae, which mitigate oxidative DNA damage stemming from terrestrial environmental stress. Several genes, such as ND2, ND4, ND6 and COIII, have experienced positive selection, suggesting their important roles in enhancing the efficiency of ATP production to cope with the increased energy requirements for life in terrestrial environments. These results strongly suggest that the adaptive evolution of mitochondrial genes has played a key role in terrestrial adaptions in Amblyopinae, as well as in Oxudercinae, and provide new insights into the molecular mechanisms underlying the water-to-land transition in vertebrates.

The gastric proton pump in gobiid and mudskipper fishes. Evidence of stomach loss?

Stomach loss has occurred independently multiple times during gnathostome evolution with notable frequency within the Teleostei. Significantly, this loss of acid-peptic digestion has been found to correlate with the secondary genomic loss of the gastric proton pump subunits (atp4a, atp4b) and pepsinogens/pepsins (pga, pgc). Gastric glands produce gastric juice containing the acid and pepsin and thus their presence is a hallmark feature of a digestive system capable of acid-peptic digestion. However, in gobiid fishes although oesogaster and gastric glands have been identified histologically, their functional significance has been questioned. In the present study we address whether the gastric proton pump is present and expressed in gastric glands of the goby Neogobius species (Gobiidae) and in members of the family Oxudercidae, a group of amphibious gobiid fishes commonly known as mudskippers (genera: Periophthalmus, Boleophthalmus, Periophthalmodon and Scartelaos). We confirmed the presence of gastric glands and have immunohistochemically localized gastric proton pump expression to these glands in Neogobius fluviatilis and Periophthalmus novemradiatus, Periophthalmus barbarus and Boleophthalmus boddarti. Genome analysis in Neogobius melanostomus, Periophthalmus magnuspinnatus, Scartelaos histophorus, Boleophthalmus pectinirostris, and Periophthalmodon schlosseri revealed the presence of both atp4a and atp4b subunit orthologues in all species in a conserved genomic loci organization. Moreover, it was possible to deduce that the complete open reading frame and the key functional amino acid residues are present. The conserved expression of the gastric proton pump provides clear evidence of the potential for gastric acid secretion indicating that acid digestion is retained in these gobiid fishes and not lost.

One size does not fit all: Variation in anatomical traits associated with emersion behavior in mudskippers (Gobiidae: Oxudercinae)

Life histories involving transitions between differing habitats (i.e., aquatic to terrestrial or marine to freshwater) require numerous anatomical, physiological, and behavioral changes. Often, the traits associated with these changes are thought to come in suites, but all traits thought to be associated with particular life histories may not be required. While some traits are found in all species with a particular habitat transition, a grab bag approach may apply to other traits in that any trait may be sufficient for successful habitat transitions. We examine patterns of morphological traits associated with prolonged emersion in mudskipper, an amphibious fishes clade, where prolonged emersion appears twice. We test the evolutionary history of multiple characteristics associated with cutaneous respiration. We find most traits thought to be key for prolonged emersion show no phylogenetic signal and no tight correlation with prolonged emersion. Such traits appear in species with prolonged emersion but also non-emerging species. Only capillary density, which, when increased, allows for increased oxygen absorption, shows strong phylogenetic signal and correlation with prolonged emersion. Further experimental, functional genomics, and observational studies are needed to fully understand the mechanisms associated with each of these traits. With respect to traits associated with other particular behaviors, a comparative framework can be helpful in identifying evolutionary correlates.

Expression of Antimicrobic Peptide Piscidin1 in Gills Mast Cells of Giant Mudskipper Periophthalmodon schlosseri (Pallas, 1770)

The amphibious teleost Giant mudskipper (Periophthalmodon schlosseri, Pallas 1770) inhabit muddy plains and Asian mangrove forests. It spends more than 90% of its life outside of the water, using its skin, gills, and buccal-pharyngeal cavity mucosa to breathe in oxygen from the surrounding air. All vertebrates have been found to have mast cells (MCs), which are part of the innate immune system. These cells are mostly found in the mucous membranes of the organs that come in contact with the outside environment. According to their morphology, MCs have distinctive cytoplasmic granules that are released during the degranulation process. Additionally, these cells have antimicrobial peptides (AMPs) that fight a variety of infections. Piscidins, hepcidins, defensins, cathelicidins, and histonic peptides are examples of fish AMPs. Confocal microscopy was used in this study to assess Piscidin1 expression in Giant Mudskipper branchial MCs. Our results demonstrated the presence of MCs in the gills is highly positive for Piscidin1. Additionally, colocalized MCs labeled with TLR2/5-HT and Piscidin1/5-HT supported our data. The expression of Piscidin1 in giant mudskipper MCs highlights the involvement of this peptide in the orchestration of teleost immunity, advancing the knowledge of the defense system of this fish.

Whole genome sequencing of the fast-swimming Southern bluefin tuna (Thunnus maccoyii)

The economically important Southern bluefin tuna (Thunnus maccoyii) is a world-famous fast-swimming fish, but its genomic information is limited. Here, we performed whole genome sequencing and assembled a draft genome for Southern bluefin tuna, aiming to generate useful genetic data for comparative functional prediction. The final genome assembly is 806.54 Mb, with scaffold and contig N50 values of 3.31 Mb and 67.38 kb, respectively. Genome completeness was evaluated to be 95.8%. The assembled genome contained 23,403 protein-coding genes and 236.1 Mb of repeat sequences (accounting for 29.27% of the entire assembly). Comparative genomics analyses of this fast-swimming tuna revealed that it had more than twice as many hemoglobin genes (18) as other relatively slow-moving fishes (such as seahorse, sunfish, and tongue sole). These hemoglobin genes are mainly localized in two big clusters (termed as “MNˮ and “LAˮ respectively), which is consistent with other reported fishes. However, Thr39 of beta-hemoglobin in the MN cluster, conserved in other fishes, was mutated as cysteine in tunas including the Southern bluefin tuna. Since hemoglobins are reported to transport oxygen efficiently for aerobic respiration, our genomic data suggest that both high copy numbers of hemoglobin genes and an adjusted function of the beta-hemoglobin may support the fast-swimming activity of tunas. In summary, we produced a primary genome assembly and predicted hemoglobin-related roles for the fast-swimming Southern bluefin tuna.

Whole-Genome Survey Analyses Provide a New Perspective for the Evolutionary Biology of Shimofuri Goby, Tridentiger bifasciatus

The shimofuri goby (Tridentiger bifasciatus) is a small and highly adaptable goby, distributed along the coasts of China, the Sea of Japan, and the west coastal and estuarine areas of the Northwest Pacific. Next-generation sequencing was used to generate genome-wide survey data to provide essential characterization of the shimofuri goby genome and for the further mining of genomic information. The genome size of the shimofuri goby was estimated to be approximately 887.60 Mb through K-mer analysis, with a heterozygosity ratio and repeat sequence ratio of 0.47% and 32.60%, respectively. The assembled genome was used to identify microsatellite motifs (Simple Sequence Repeats, SSRs), extract single-copy homologous genes and assemble the mitochondrial genome. A total of 288,730 SSRs were identified. The most frequent SSRs were dinucleotide repeats (with a frequency of 61.15%), followed by trinucleotide (29.87%), tetranucleotide (6.19%), pentanucleotide (1.13%), and hexanucleotide repeats (1.66%). The results of the phylogenetic analysis based on single-copy homologous genes showed that the shimofuri goby and Rhinogobius similis can be clustered into one branch. The shimofuri goby was originally thought to be the same as the chameleon goby (Tridentiger trigonocephalus) due to their close morphological resemblance. However, a complete mitochondrial genome was assembled and the results of the phylogenetic analysis support the inclusion of the shimofuri goby as a separate species. PSMC analysis indicated that the shimofuri goby experienced a bottleneck event during the Pleistocene Glacial Epoch, in which its population size decreased massively, and then it began to recover gradually after the Last Glacial Maximum. This study provides a reference for the further assembly of the complete genome map of the shimofuri goby, and is a valuable genomic resource for the study of its evolutionary biology.

Coastal ecological impacts from pumice rafts

An explosive volcanic eruption occurred in the Ogasawara Islands on 13-15 August 2021, bringing unprecedented amounts of floating pumice to the coast of Okinawa Island in the Ryukyu Archipelago, 1300 km west of the volcano, approximately 2 months later. The coast of Okinawa Island, especially along the northern part, is home to many typical subtropical seascapes, including coral reefs and mangrove forests, so the possible impact of the large amount of pumice is attracting attention. Here, we report early evidence of ecosystem changes as a result of large-scale pumice stranding on coastal beaches, in estuaries and mangrove forests and passage across fringing coral reefs. Massive pumice drifts are major obstacles to fishing activities and ship traffic, but short and long-term changes in coastal ecosystems can also occur. The phenomena observed on Okinawa Island can be a preview of coastal impacts for the Kyushu, Shikoku, Honshu Islands, where pumice has subsequently washed ashore.

A chromosome-level genome assembly of the jade perch (Scortum barcoo)

Endemic to Australia, jade perch (Scortum barcoo) is a highly profitable freshwater bass species. It has extraordinarily high levels of omega-3 polyunsaturated fatty acids (PUFAs), which detailed genes involved in are largely unclear. Meanwhile, there were four chromosome-level bass species have been previous sequenced, while the bass ancestor genome karyotypes have not been estimated. Therefore, we sequenced, assembled and annotated a genome of jade perch to characterize the detailed genes for biosynthesis of omega-3 PUFAs and to deduce the bass ancestor genome karyotypes. We constructed a chromosome-level genome assembly with 24 pairs of chromosomes, 657.7 Mb in total length, and the contig and the scaffold N50 of 4.8 Mb and 28.6 Mb respectively. We also identified repetitive elements (accounting for 19.7% of the genome assembly) and predicted 26,905 protein-coding genes. Meanwhile, we performed genome-wide localization and characterization of several important genes encoding some key enzymes in the biosynthesis pathway of PUFAs. These genes may contribute to the high concentration of omega-3 in jade perch. Moreover, we conducted a series of comparative genomic analyses among four representative bass species at a chromosome level, resulting in a series of sequences of a deductive bass ancestor genome.

Measurement(s)	whole genome sequencing
Technology Type(s)	Illumina Sequencing • Oxford Nanopore Sequencing

Genome-wide identification and expression patterns of opsin genes during larval development in Chinese perch (Siniperca chuatsi)

Vision is important for fish to forage food and fishes express opsin genes to receive visual signals. Chinese perch (Siniperca chuatsi) larvae prey on other fish species larvae at firstfeeding but donoteat any zooplankton, the expression of opsin genes in S. chuatsilarvae is unknown. In this study, we conducted a whole-genome analysis and demonstrated that S. chuatsihave5cone opsin genes (sws1, sws2Aα, sws2Aβ, rh2and lws)and 2 rod opsin genes (rh1and rh1-exorh). The syntenicanalysisshowedthe flanking genes ofall opsin genes were conserved during fish evolution, but the ancestorof S. chuatsimightlost some opsin gene copies duringtheevolution.The phylogeneticanalysisshowed sws1of S. chuatsiwas closest to those of Lates calcariferwhich had a truncated sws1gene; the sws2Aα, sws2Aβ,lws,rh2,rh1 andrh1-exorh of S. chuatsihad a closer relationship with those of Percomorpha fishes.Importantly, results of in situhybridization showed the sws1 opsingene,which is related to forage zooplankton,had extremely low levelexpression in retinaat early stages.Surprisingly, the rh2 opsin gene had a high level expression at firstfeeding stage. The sws2Aα, sws2Aβand lwshad a little expression at early stages but the lwsshowed a increasing trend with larval development, rh1 opsin gene expression appeared at15 dph. In thisstudy, we found a specialpattern of visual opsin genes expression in S. chuatsi, it might influence the larval first feeding and feeding habit.

Evolutionary Genomics Reveals Multiple Functions of Arylalkylamine N-Acetyltransferase in Fish

As an important hormone, melatonin participates in endocrine regulation of diverse functions in vertebrates. Its biosynthesis is catalyzed by four cascaded enzymes, among them, arylalkylamine N-acetyltransferase (AANAT) is the most critical one. Although only single aanat gene has been identified in most groups of vertebrates, researchers including us have determined that fish have the most diverse of aanat genes (aanat1a, aanat1b, and aanat2), playing various potential roles such as seasonal migration, amphibious aerial vision, and cave or deep-sea adaptation. With the rapid development of genome and transcriptome sequencing, more and more putative sequences of fish aanat genes are going to be available. Related phylogeny and functional investigations will enrich our understanding of AANAT functions in various fish species.

New Insight Into Visual Adaptation in the Mudskipper Cornea: From Morphology to the Cornea-Related COL8A2 Gene

Much research has focused on visual system evolution in bony fishes. The capacity of visual systems to perceive and respond to external signals is integral to evolutionary success. However, integrated research on the mechanisms of adaptive evolution based on corneal structure and related genes remains limited. In this study, scanning electron microscopy (SEM) was used to assess the microstructure and adaptation of corneal epithelial cells. Then, the evolution of the cornea-related COL8A2 gene was investigated. We found various projections (microridges, microplicae, microholes, and microvilli) on the corneal epithelial cells of amphibious mudskippers. Compared with those of fully aquatic fishes, these microstructures were considered adaptations to the variable environments experienced by amphibious mudskippers, as they can resist dryness in terrestrial environments and infection in aquatic environments. Moreover, strong purifying selection was detected for COL8A2. In addition, some specific amino acid substitution sites were also identified in the COL8A2 sequence in mudskippers. Interestingly, the evolutionary rate of the COL8A2 gene was significantly and positively correlated with maximum diving depth in our dataset. Specifically, with increasing diving depth, the evolutionary rate of the COL8A2 gene seemed to gradually accelerate. The results indicated that the cornea of bony fishes has evolved through adaptation to cope with the different diving depths encountered during the evolutionary process, with the corneal evolution of the amphibious mudskipper group showing a unique pattern.

Phylogeny, diversification, and biogeography of a hemiclonal hybrid system of native Australian freshwater fishes (Gobiiformes: Gobioidei: Eleotridae: Hypseleotris)

Background

Carp gudgeons (genus Hypseleotris) are a prominent part of the Australian freshwater fish fauna, with species distributed around the western, northern, and eastern reaches of the continent. We infer a calibrated phylogeny of the genus based on nuclear ultraconserved element (UCE) sequences and using Bayesian estimation of divergence times, and use this phylogeny to investigate geographic patterns of diversification with GeoSSE. The southeastern species have hybridized to form hemiclonal lineages, and we also resolve relationships of hemiclones and compare their phylogenetic placement in the UCE phylogeny with a hypothesis based on complete mitochondrial genomes. We then use phased SNPs extracted from the UCE sequences for population structure analysis among the southeastern species and hemiclones.

Results

Hypseleotris cyprinoides, a widespread euryhaline species known from throughout the Indo-Pacific, is resolved outside the remainder of the species. Two Australian radiations comprise the bulk of Hypseleotris, one primarily in the northwestern coastal rivers and a second inhabiting the southeastern region including the Murray–Darling, Bulloo-Bancannia and Lake Eyre basins, plus coastal rivers east of the Great Dividing Range. Our phylogenetic results reveal cytonuclear discordance between the UCE and mitochondrial hypotheses, place hemiclone hybrids among their parental taxa, and indicate that the genus Kimberleyeleotris is nested within the northwestern Hypseleotris radiation along with three undescribed species. We infer a crown age for Hypseleotris of 17.3 Ma, date the radiation of Australian species at roughly 10.1 Ma, and recover the crown ages of the northwestern (excluding H. compressa) and southeastern radiations at 5.9 and 7.2 Ma, respectively. Range-dependent diversification analyses using GeoSSE indicate that speciation and extinction rates have been steady between the northwestern and southeastern Australian radiations and between smaller radiations of species in the Kimberley region and the Arnhem Plateau. Analysis of phased SNPs confirms inheritance patterns and reveals high levels of heterozygosity among the hemiclones.

Conclusions

The northwestern species have restricted ranges and likely speciated in allopatry, while the southeastern species are known from much larger areas, consistent with peripatric speciation or allopatric speciation followed by secondary contact. Species in the northwestern Kimberley region differ in shape from those in the southeast, with the Kimberley species notably more elongate and slender than the stocky southeastern species, likely due to the different topographies and flow regimes of the rivers they inhabit.

Epigenomic and Proteomic Changes in Fetal Spleens Persistently Infected with Bovine Viral Diarrhea Virus: Repercussions for the Developing Immune System, Bone, Brain, and Heart

Bovine viral diarrhea virus (BVDV) infection during early gestation results in persistently infected (PI) immunotolerant calves that are the primary reservoirs of the virus. Pathologies observed in PI cattle include congenital defects of the brain, heart, and bone as well as marked functional defects in their immune system. It was hypothesized that fetal BVDV infection alters T cell activation and signaling genes by epigenetic mechanisms. To test this, PI and control fetal splenic tissues were collected on day 245 of gestation, 170 days post maternal infection. DNA was isolated for reduced representation bisulfite sequencing, protein was isolated for proteomics, both were analyzed with appropriate bioinformatic methods. Within set parameters, 1951 hypermethylated and 691 hypomethylated DNA regions were identified in PI compared to control fetuses. Pathways associated with immune system, neural, cardiac, and bone development were associated with heavily methylated DNA. The proteomic analysis revealed 12 differentially expressed proteins in PI vs. control animals. Upregulated proteins were associated with protein processing, whereas downregulated proteins were associated with lymphocyte migration and development in PI compared to control fetal spleens. The epigenetic changes in DNA may explain the immune dysfunctions, abnormal bone formation, and brain and heart defects observed in PI animals.

Phylogenetics And The Evolution Of Terrestriality In Mudskippers (Gobiidae: Oxudercinae)

The initial vertebrate conquest of land by stegocephalians (Sarcopterygia) allowed access to new resources and exploitation of untapped niches precipitating a major phylogenetic diversification. However, a paucity of fossils has left considerable uncertainties about phylogenetic relationships and the eco-morphological stages in this key transition in Earth history. Among extant actinopterygians, three genera of mudskippers (Gobiidae: Oxudercinae), Boleophthalmus, Periophthalmus and Periophthalmodon are the most terrestrialized, with vertebral, appendicular, locomotory, respiratory, and epithelial specializations enabling overland excursions up to 14 hours. Unlike early stegocephalians, the ecologies and morphologies of the 45 species of oxudercines are well known, making them viable analogs for the initial vertebrate conquest of land. Nevertheless, they have received little phylogenetic attention. We compiled the largest molecular dataset to date, with 29 oxudercine species, and 5 nuclear and mitochondrial loci. Phylogenetic and comparative analyses revealed strong support for two independent terrestrial transitions, and a complex suit of ecomorphological forms in estuarine environments. Furthermore, neither Oxudercinae nor their presumed sister-group the eel gobies (Amblyopinae, a group of elongated gobies) were monophyletic with respect to each other, requiring a merging of these two subfamilies and revealing an expansion of phenotypic variation within the "mudskipper" clade. We did not find support for the expected linear model of ecomorphological and locomotory transition from fully aquatic, to mudswimming, to pectoral-aided mudswimming, to lobe-finned terrestrial locomotion proposed by earlier morphological studies. This high degree of convergent or parallel transitions to terrestriality, and apparent divergent directions of estuarine adaptation, promises even greater potential for this clade to illuminate the conquest of land. Future work should focus on these less-studied species with "transitional" and other mud-habitat specializations to fully resolve the dynamics of this diversification.

Chromosome-level genome assembly of the shuttles hoppfish, Periophthalmus modestus

Background

The shuttles hoppfish (mudskipper), Periophthalmus modestus, is one of the mudskippers, which are the largest group of amphibious teleost fishes, which are uniquely adapted to live on mudflats. Because mudskippers can survive on land for extended periods by breathing through their skin and through the lining of the mouth and throat, they were evaluated as a model for the evolutionary sea-land transition of Devonian protoamphibians, ancestors of all present tetrapods.

Results

A total of 39.6, 80.2, 52.9, and 33.3 Gb of Illumina, Pacific Biosciences, 10X linked, and Hi-C data, respectively, was assembled into 1,419 scaffolds with an N50 length of 33 Mb and BUSCO score of 96.6%. The assembly covered 117% of the estimated genome size (729 Mb) and included 23 pseudo-chromosomes anchored by a Hi-C contact map, which corresponded to the top 23 longest scaffolds above 20 Mb and close to the estimated one. Of the genome, 43.8% were various repetitive elements such as DNAs, tandem repeats, long interspersed nuclear elements, and simple repeats. Ab initio and homology-based gene prediction identified 30,505 genes, of which 94% had homology to the 14 Actinopterygii transcriptomes and 89% and 85% to Pfam familes and InterPro domains, respectively. Comparative genomics with 15 Actinopterygii species identified 59,448 gene families of which 12% were only in P. modestus.

Conclusions

We present the high quality of the first genome assembly and gene annotation of the shuttles hoppfish. It will provide a valuable resource for further studies on sea-land transition, bimodal respiration, nitrogen excretion, osmoregulation, thermoregulation, vision, and mechanoreception.

Expressions of melanopsins in telencephalon imply their function in synchronizing semilunar spawning rhythm in the mudskipper Boleophthalmus pectinirostris

The mudskipper Boleophthalmus pectinirostris inhabits intertidal mudflats, exhibiting semilunar reproductive rhythms. To investigate whether melanopsin is possibly involved in the synchronization of the semilunar spawning rhythm in the female mudskipper, we first cloned all four melanopsin subtypes (opn4m1, opn4m3, opn4x1, opn4x2) in B. pectinirostris. Results from RTq-PCR showed that significantly higher transcription levels of all four melanopsin subtypes were observed in the eyes rather than other tissues. In brain, all four melanopsin subtypes were also detectable in different regions, including the telencephalon, in which the expression of melanopsin has not been reported in other teleosts. The transcription levels of opn4m3 and opn4x1 in the telencephalon exhibited a daily fluctuation pattern. When females entered the spawning season, opn4m1 and opn4x1 transcript levels increased significantly in the telencephalon. During the spawning season, the transcript levels of opn4m3 and opn4x1 in the telencephalon appeared to have a cyclic pattern associated with semilunar periodicity, exhibiting two cycles with a peak around the first or the last lunar quarters. Results from ISH showed that, opn4x1 mRNA was localized in the medial of dorsal telencephalic area, dorsal nucleus of ventral telencephalic area (Vd), ventral nucleus of ventral telencephalic area (Vv), anterior part of parvocellular preoptic nucleus, magnocellular part of the magnocellular preoptic nucleus (PMmc), habenular and ventral zone of hypothalamus. Intriguingly, gnrh3 mRNA was also located in Vd, Vv and PMmc. Taken together, our results suggested that melanopsins, e.g. opn4x1, expressed in the telencephalon might mediate semilunar spawning activity in the female mudskipper.

Temporal dynamics of teleost populations during the Pleistocene: a report from publicly available genome data

Background

Global climate oscillation, as a selection dynamic, is an ecologically important element resulting in global biodiversity. During the glacial geological periods, most organisms suffered detrimental selection pressures (such as food shortage and habitat loss) and went through population declines. However, during the mild interglacial periods, many species re-flourished. These temporal dynamics of effective population sizes (N_e) provide essential information for understanding and predicting evolutionary outcomes during historical and ongoing global climate changes.

Results

Using high-quality genome assemblies and corresponding sequencing data, we applied the Pairwise Sequentially Markovian Coalescent (PSMC) method to quantify N_e changes of twelve representative teleost species from approximately 10 million years ago (mya) to 10 thousand years ago (kya). These results revealed multiple rounds of population contraction and expansion in most of the examined teleost species during the Neogene and the Quaternary periods. We observed that 83% (10/12) of the examined teleosts had experienced a drastic decline in N_e before the last glacial period (LGP, 110–12 kya), slightly earlier than the reported pattern of N_e changes in 38 avian species. In comparison with the peaks, almost all of the examined teleosts maintained long-term lower N_e values during the last few million years. This is consistent with increasingly dramatic glaciation during this period.

Conclusion

In summary, these findings provide a more comprehensive understanding of the historical N_e changes in teleosts. Results presented here could lead to the development of appropriate strategies to protect species in light of ongoing global climate changes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07816-7.

Balancing selection maintains ancient polymorphisms at conserved enhancers for the olfactory receptor genes of a Chinese marine fish

The study of balancing selection, as a selective force maintaining adaptive genetic variation in gene pools longer than expected by drift, is currently experiencing renewed interest due to the increased availability of new data, methods of analysis, and case studies. In this investigation, evidence of balancing selection operating on conserved enhancers of the olfactory receptor (OR) genes is presented for the Chinese sleeper (Bostrychus sinensis), a coastal marine fish that is emerging as a model species for evolutionary studies in the Northwest Pacific marginal seas. Coupled with tests for Gene Ontology enrichment and transcription factor binding, population genomic data allow for the identification of an OR cluster in the sleeper with a downstream flanking region containing three enhancers that are conserved with human and other fish species. Phylogenetic and population genetic analyses indicate that the enhancers are under balancing selection as evidenced by their translineage polymorphisms, excess common alleles, and increased within-group diversities. Age comparisons between the translineage polymorphisms and most recent common ancestors of neutral genealogies substantiate that the former are old, and thus, due to ancient balancing selection. The survival and reproduction of vertebrates depend on their sense of smell, and thereby, on their ORs. In addition to locus duplication and allelic variation of structural genes, this study highlights a third mechanism by which receptor diversity can be achieved for detecting and responding to the huge variety of environmental odorants (i.e., by balancing selection acting on OR gene expression through their enhancer variability).

A Chromosome-Level Genome Assembly of the Mandarin Fish (Siniperca chuatsi)

The mandarin fish, Siniperca chuatsi, is an economically important perciform species with widespread aquaculture practices in China. Its special feeding habit, acceptance of only live prey fishes, contributes to its delicious meat. However, little is currently known about related genetic mechanisms. Here, we performed whole-genome sequencing and assembled a 758.78 Mb genome assembly of the mandarin fish, with the scaffold and contig N50 values reaching 2.64 Mb and 46.11 Kb, respectively. Approximately 92.8% of the scaffolds were ordered onto 24 chromosomes (Chrs) with the assistance of a previously established genetic linkage map. The chromosome-level genome contained 19,904 protein-coding genes, of which 19,059 (95.75%) genes were functionally annotated. The special feeding behavior of mandarin fish could be attributable to the interaction of a variety of sense organs (such as vision, smell, and endocrine organs). Through comparative genomics analysis, some interesting results were found. For example, olfactory receptor (OR) genes (especially the beta and delta types) underwent a significant expansion, and endocrinology/vision related npy, spexin, and opsin genes presented various functional mutations. These may contribute to the special feeding habit of the mandarin fish by strengthening the olfactory and visual systems. Meanwhile, previously identified sex-related genes and quantitative trait locis (QTLs) were localized on the Chr14 and Chr17, respectively. 155 toxin proteins were predicted from mandarin fish genome. In summary, the high-quality genome assembly of the mandarin fish provides novel insights into the feeding habit of live prey and offers a valuable genetic resource for the quality improvement of this freshwater fish.

Large-scale genome sampling reveals unique immunity and metabolic adaptations in bats

Comprising more than 1,400 species, bats possess adaptations unique among mammals including powered flight, unexpected longevity, and extraordinary immunity. Some of the molecular mechanisms underlying these unique adaptations includes DNA repair, metabolism and immunity. However, analyses have been limited to a few divergent lineages, reducing the scope of inferences on gene family evolution across the Order Chiroptera. We conducted an exhaustive comparative genomic study of 37 bat species, one generated in this study, encompassing a large number of lineages, with a particular emphasis on multi-gene family evolution across immune and metabolic genes. In agreement with previous analyses, we found lineage-specific expansions of the APOBEC3 and MHC-I gene families, and loss of the proinflammatory PYHIN gene family. We inferred more than 1,000 gene losses unique to bats, including genes involved in the regulation of inflammasome pathways such as epithelial defence receptors, the natural killer gene complex and the interferon-gamma induced pathway. Gene set enrichment analyses revealed genes lost in bats are involved in defence response against pathogen-associated molecular patterns and damage-associated molecular patterns. Gene family evolution and selection analyses indicate bats have evolved fundamental functional differences compared to other mammals in both innate and adaptive immune system, with the potential to enhance antiviral immune response while dampening inflammatory signalling. In addition, metabolic genes have experienced repeated expansions related to convergent shifts to plant-based diets. Our analyses support the hypothesis that, in tandem with flight, ancestral bats had evolved a unique set of immune adaptations whose functional implications remain to be explored.

Chromosome-level genome assembly of the endangered humphead wrasse Cheilinus undulatus: Insight into the expansion of opsin genes in fishes

Wrasses are dominant components of major coral reef systems. Among wrasses, Cheilinus undulatus is an endangered species with high economic and ecological value that exhibits sex reversal of females to males, while sexual selection occurs in breeding aggregations. However, the molecular-associated mechanism underlying this remains unclear. Opsin gene diversification is regarded as a potent force in sexual selection. Here we present a genome assembly of C. undulatus, using Illumina, Nanopore and Hi-C sequencing. The 1.17 Gb genome was generated from 328 contigs with an N50 length of 16.5 Mb and anchored to 24 chromosomes. In total, 22,218 genes were functionally annotated, and 96.36% of BUSCO genes were fully represented. Transcriptomic analyses showed that 96.79% of the predicted genes were expressed. Transposons were most abundant, accounting for 39.88% of the genome, with low divergence, owing to their evolution with close species ~60.53 million years ago. In total, 567/1,826 gene families were expanded and contracted in the reconstructed phylogeny, respectively. Forty-six genes were under positive selection. Comparative genomic analyses with other fish revealed expansion of opsin SWS2B, LWS1 and Rh2. The elevated duplicates of SWS2B were generated by gene conversions via transposition of transposons followed by nonallelic homologous recombination. Amino acid substitutions of opsin paralogues occurred at key tuning sites, causing a spectral shift in maximal absorbance of visual pigment to capture functional changes. Among these opsin genes, SWS2B-3 and 4 and Rh1 are expressed in the retina. The genome sequence of C. undulatus provides valuable resources for future investigation of the conservation, evolution and behaviour of fishes.

Evolutionary divergence of a Hoxa2b hindbrain enhancer in syngnathids mimics results of functional assays

Hoxa2 genes provide critical patterning signals during development, and their regulation and function have been extensively studied. We report a previously uncharacterized significant sequence divergence of a highly conserved hindbrain hoxa2b enhancer element in the family syngnathidae (pipefishes, seahorses, pipehorses, seadragons). We compared the hox cis-regulatory element variation in the Gulf pipefish and two species of seahorse against eight other species of fish, as well as human and mouse. We annotated the hoxa2b enhancer element binding sites across three species of seahorse, four species of pipefish, and one species of ghost pipefish. Finally, we performed in situ hybridization analysis of hoxa2b expression in Gulf pipefish embryos. We found that all syngnathid fish examined share a modified rhombomere 4 hoxa2b enhancer element, despite the fact that this element has been found to be highly conserved across all vertebrates examined previously. Binding element sequence motifs and spacing between binding elements have been modified for the hoxa2b enhancer in several species of pipefish and seahorse, and that the loss of the Prep/Meis binding site and further space shortening happened after ghost pipefish split from the rest of the syngnathid clade. We showed that expression of this gene in rhombomere 4 is lower relative to the surrounding rhombomeres in developing Gulf pipefish embryos, reflecting previously published functional tests for this enhancer. Our findings highlight the benefits of studying highly derived, diverse taxa for understanding of gene regulatory evolution and support the hypothesis that natural mutations can occur in deeply conserved pathways in ways potentially related to phenotypic diversity.

Genome-wide identification and characterization of 14-3-3 genes in fishes

The 14-3-3 family genes are highly conserved regulatory factors in eukaryotes with involvement in multiple important cellular processes. However, detailed investigations of this family in fishes are very limited. Here, a comparative genomic and transcriptomic survey were performed to investigate the 14-3-3 family in fishes. We confirmed that the numbers of 14-3-3 genes ranged from 5 to 7 in non-teleost fishes, as well as additional 14-3-3 genes (9 to 11) in teleost fishes. In addition, some special teleost fishes possess 17 to 25 14-3-3s, which undergone the fourth whole-genome duplication (WGD). We also found that six pairs of fish 14-3-3 genes were clustered with mammalian ε, γ, ς, η, τand β isotypes, respectively, while σ was absent with a potential specificity within mammals, on the basis of their phylogenetic and synteny analyses. According to our results, we inferred that the diversity of 14-3-3 genes in fishes seems to be generated from a combination of WGD and gene loss. Comparative transcriptomic analysis revealed that there are differences in tissue distribution, and we speculated that 14-3-3 genes may contribute to terrestrial adaptations in mudskippers. In addition, protein sequence alignments of 14-3-3s supported their differential roles in fishes. In summary, our present comparative genomic and transcriptomic survey will benefit for further functional investigations of these fish genes.

Whole genome sequencing of a snailfish from the Yap Trench (~7,000 m) clarifies the molecular mechanisms underlying adaptation to the deep sea

Hadal environments (depths below 6,000 m) are characterized by extremely high hydrostatic pressures, low temperatures, a scarce food supply, and little light. The evolutionary adaptations that allow vertebrates to survive in this extreme environment are poorly understood. Here, we constructed a high-quality reference genome for Yap hadal snailfish (YHS), which was captured at a depth of ~7,000 m in the Yap Trench. The final YHS genome assembly was 731.75 Mb, with a contig N50 of 0.75 Mb and a scaffold N50 of 1.26 Mb. We predicted 24,329 protein-coding genes in the YHS genome, and 24,265 of these genes were successfully functionally annotated. Phylogenetic analyses suggested that YHS diverged from a Mariana Trench snailfish approximately 0.92 million years ago. Many genes associated with DNA repair show evidence of positive selection and have expanded copy numbers in the YHS genome, possibly helping to maintain the integrity of DNA under increased hydrostatic pressure. The levels of trimethylamine N-oxide (TMAO), a potent protein stabilizer, are much higher in the muscles of YHS than in those of shallow-water fish. This difference is perhaps due to the five copies of the TMAO-generating enzyme flavin-containing monooxygenase-3 gene (fmo3) in the YHS genome and the abundance of trimethylamine (TMA)-generating bacteria in the YHS gut. Thus, the high TMAO content might help YHS adapt to high hydrostatic pressure by improving protein stability. Additionally, the evolutionary features of the YHS genes encoding sensory-related proteins are consistent with the scarce food supply and darkness in the hadal environments. These results clarify the molecular mechanisms underlying the adaptation of hadal organisms to the deep-sea environment and provide valuable genomic resources for in-depth investigations of hadal biology.

Hadal environments (depths below 6,000 m) are characterized by extremely high hydrostatic pressures, low temperatures, a scarce food supply, and little light. Fish are the only vertebrates inhabiting the hadal zone, and hadal snailfishes have been found in at least five geographically separated marine trenches. However, the genetic mechanisms that allow vertebrates to live in such extreme conditions are not well understood. Here, we constructed a high-quality reference genome for Yap hadal snailfish (YHS) captured at a depth of ~7,000 m in the Yap Trench, using long reads obtained by Pacific Biosciences Sequel sequencing. Comparative genomic analyses revealed that many genes associated with DNA repair show evidence of positive selection and have expanded copy numbers in the YHS genome, which potentially reflect the difficulty of maintaining DNA integrity under high hydrostatic pressure. Moreover, the five copies of the trimethylamine N-oxide (TMAO)-generating enzyme flavin-containing monooxygenase-3 gene (fmo3) and the abundance of trimethylamine (TMA)-generating bacteria in the YHS gut could provide enough TMAO to improve protein stability under hadal conditions. In addition, characteristics of the YHS sensory system genes were consistent with the scarce food supply and darkness in the hadal zone. Our results provide new insights into the molecular mechanisms underlying the adaptation of hadal organisms to the deep-sea environment and valuable genomic resources that will help further clarify hadal adaptations.

Genome of the butterfly hillstream loach provides insights into adaptations to torrential mountain stream life

Butterfly hillstream loach (Beaufortia kweichowensis), a benthic fish in the torrential mountain streams, possesses a totally flat ventrum, flattened craniofacial and body skeletons, and enlarged paired fins covered by substantially small keratinous structures. However, little is known about the genetic basis of these specialized morphological adaptations. Here we present a 448.52-Mb genome assembly with contig N50 length of 5.53 Mb by integrating Illumina short-read sequencing, Nanopore long-read sequencing and HiC-based chromatin map. Demographic history reconstruction of the butterfly hillstream loach reveals that the population dynamics is correlated with the different stages of uplifting of the Tibetan Plateau. Comparative genomic analysis finds evidence of six keratin genes in butterfly hillstream loach evolving under positive selection. Within these genes, two keratin genes exhibit species-specific and divergent amino acid changes, suggesting a role in the formation of the unculi. Additionally, a series of positively selected genes, rapid evolving genes, specific variant genes and expanded gene families are found, including genes related to Hedgehog, Notch and BMP pathways, which may be involved in craniofacial development. These findings may have important implications for understanding the genetic basis of phenotypic adaptation to torrential mountain stream life.

Genomic characteristics and profile of microsatellite primers for Acanthogobius ommaturus by genome survey sequencing

Acanthogobius ommaturus is one of the suitable species to study the genetic mechanism of adaptive evolution, but there are few reports on its genetics. In the present study, the genomic survey sequencing method was used to analyze the genome characters of A. ommaturus. A total of 50.50 G high-quality sequence data were obtained in the present study. From the 19-mer distribution frequency, the estimated genome size was 928.01 Mb. The calculated sequence repeat rate was about 38.31%, the heterozygosity was approximately 0.17%, and the GC% content was approximately 40.88%. Moreover, 475,724 simple sequence repeats (SSRs) were identified. Among them, dinucleotide repeats were the most (53.70% of the total SSRs), followed by tri- (35.36%), hexa- (4.59%), tetra- (4.57%) and penta- (1.77%) nucleotide repeats type. This is the first genome-wide feature of this species to be reported.

Genetic Adaptations in Mudskipper and Tetrapod Give Insights into Their Convergent Water-to-Land Transition

Water-to-land transition has been independently evolved in multiple vertebrate lineages including the most recent common ancestor of tetrapod and multiple fish clades, and among them, mudskippers uniquely adapted to the mudflat. Even though physiological and morphological adaptation of mudskippers is thought to resemble that of the ancestral tetrapod, it is unclear if they share genome-wide evolutionary signatures. To detect potential signatures of positive selection in mudskipper and tetrapods, we analyzed 4118 singleton orthologues of terrestrial tetrapods, coelacanth, mudskipper, and fully aquatic fishes. Among positively selected genes identified in mudskipper and tetrapod lineages, genes involved in immune responses, mitochondrial oxidative phosphorylation, and kidney development were detected. On the other hand, tetrapod-specific and mudskipper-specific positively selected genes were functionally enriched for DNA repair processes, which could be associated with higher exposure to UV light. We also performed gene family analysis and discovered convergent contraction of eight gene families, including βγ-crystallin coding genes in both tetrapod and mudskipper lineages. Findings of this study suggest the similar genetic adaptation against environmental constraints between the ancient tetrapod and mudskippers for their land adaptation.

Genomic and physiological mechanisms underlying skin plasticity during water to air transition in an amphibious fish

The terrestrial radiation of vertebrates required changes in skin that resolved the dual demands of maintaining a mechanical and physiological barrier while also facilitating ion and gas transport. Using the amphibious killifish Kryptolebias marmoratus, we found that transcriptional regulation of skin morphogenesis was quickly activated upon air exposure (1 h). Rapid regulation of cell-cell adhesion complexes and pathways that regulate stratum corneum formation was consistent with barrier function and mechanical reinforcement. Unique blood vessel architecture and regulation of angiogenesis likely supported cutaneous respiration. Differences in ionoregulatory transcripts and ionocyte morphology were correlated with differences in salinity acclimation and resilience to air exposure. Evolutionary analyses reinforced the adaptive importance of these mechanisms. We conclude that rapid plasticity of barrier, respiratory and ionoregulatory functions in skin evolved to support the amphibious lifestyle of K. marmoratus; similar processes may have facilitated the terrestrial radiation of other contemporary and ancient fishes.

The genome of walking catfish Clarias magur (Hamilton, 1822) unveils the genetic basis that may have facilitated the development of environmental and terrestrial adaptation systems in air-breathing catfishes

The walking catfish Clarias magur (Hamilton, 1822) (magur) is an important catfish species inhabiting the Indian subcontinent. It is considered as a highly nutritious food fish and has the capability to walk to some distance, and survive a considerable period without water. Assembly, scaffolding and several rounds of iterations resulted in 3,484 scaffolds covering ∼94% of estimated genome with 9.88 Mb largest scaffold, and N50 1.31 Mb. The genome possessed 23,748 predicted protein encoding genes with annotation of 19,279 orthologous genes. A total of 166 orthologous groups represented by 222 genes were found to be unique for this species. The Computational Analysis of gene Family Evolution (CAFE) analysis revealed expansion of 207 gene families and 100 gene families have rapidly evolved. Genes specific to important environmental and terrestrial adaptation, viz. urea cycle, vision, locomotion, olfactory and vomeronasal receptors, immune system, anti-microbial properties, mucus, thermoregulation, osmoregulation, air-breathing, detoxification, etc. were identified and critically analysed. The analysis clearly indicated that C. magur genome possessed several unique and duplicate genes similar to that of terrestrial or amphibians’ counterparts in comparison to other teleostean species. The genome information will be useful in conservation genetics, not only for this species but will also be very helpful in such studies in other catfishes.

Trained innate immunity

The innate immune system acts rapidly in an identical and nonspecific way every time the body is exposed to pathogens. As such, it cannot build and maintain immunological memory to help prevent reinfection. Researchers contend that trained immunity is influenced by intracellular metabolic pathways and epigenetic remodeling. The purpose of this review was to explore the topic of trained innate immunity based on the results of relevant previous studies. This systematic review entailed identifying articles related to trained innate immunity. The sources were obtained from PubMed using different search terms that included "trained innate immunity," "trained immunity," "trained," "innate," "immunity," and "immune system." Boolean operators were used to combine terms and phrases. A review of previous study results revealed that little is currently known about the molecular and cellular processes that mediate or induce a trained immune response in animals. However, it is believed that alterations in the phenotypes of cell populations and the numbers of specific cells may play a critical role in mediating the trained immune response. Increasing evidence shows that the protective processes and actions that occur during a secondary infection are not entirely linked to the adaptive immune system. Instead, these events also involve heightened activation of innate immune cells. While trained innate immune cells may have a shorter memory, they assist in the fight against pathogens and provide cross-protection. Identification of the mechanisms and molecules that underlie trained innate immunity has highlighted important features of the human immune response. Such advances continue to open doors for future research on how the body responds to disease-causing pathogens.

Whole-genome sequencing reveals sex determination and liver high-fat storage mechanisms of yellowstripe goby (Mugilogobius chulae)

As a promising novel marine fish model for future research on marine ecotoxicology as well as an animal model of human disease, the genome information of yellowstripe goby (Mugilogobius chulae) remains unknown. Here we report the first annotated chromosome-level reference genome assembly for yellowstripe goby. A 20.67-cM sex determination region was discovered on chromosome 5 and seven potential sex-determining genes were identified. Based on combined genome and transcriptome data, we identified three key lipid metabolic pathways for high-fat accumulation in the liver of yellowstripe goby. The changes in the expression patterns of MGLL and CPT1 at different development stage of the liver, and the expansion of the ABCA1 gene, innate immune gene TLR23, and TRIM family genes may help in balancing high-fat storage in hepatocytes and steatohepatitis. These results may provide insights into understanding the molecular mechanisms of sex determination and high-fat storage in the liver of marine fishes.

Tides: A key environmental driver of osteichthyan evolution and the fish-tetrapod transition?

Tides are a major component of the interaction between the marine and terrestrial environments, and thus play an important part in shaping the environmental context for the evolution of shallow marine and coastal organisms. Here, we use a dedicated tidal model and palaeogeographic reconstructions from the Late Silurian to early Late Devonian (420 Ma, 400 Ma and 380 Ma, Ma = millions of years ago) to explore the potential significance of tides for the evolution of osteichthyans (bony fish) and tetrapods (land vertebrates). The earliest members of the osteichthyan crown-group date to the Late Silurian, approximately 425 Ma, while the earliest evidence for tetrapods is provided by trackways from the Middle Devonian, dated to approximately 393 Ma, and the oldest tetrapod body fossils are Late Devonian, approximately 373 Ma. Large tidal ranges could have fostered both the evolution of air-breathing organs in osteichthyans to facilitate breathing in oxygen-depleted tidal pools, and the development of weight-bearing tetrapod limbs to aid navigation within the intertidal zones. We find that tidal ranges over 4 m were present around areas of evolutionary significance for the origin of osteichthyans and the fish-tetrapod transition, highlighting the possible importance of tidal dynamics as a driver for these evolutionary processes.

Acclimation to prolonged aquatic hypercarbia or air enhances hemoglobin‑oxygen affinity in an amphibious fish

When the amphibious mangrove rivulus (Kryptolebias marmoratus) leaves water for extended periods, hemoglobin-O₂ binding affinity increases. We tested the hypothesis that the change in affinity was a consequence of hemoglobin isoform switching driven by exposure to environments associated with increased internal CO₂ levels. We exposed K. marmoratus to either water (control, pH 8.1), air, aquatic hypercarbia (5.1 kPa CO₂, pH 6.6-6.8), or aquatic acid (isocarbic control, pH 6.6-6.8), for 7 days, and measured hemoglobin-O₂ affinity spectrophotometrically. We found that mangrove rivulus compensated for elevated CO₂ and aquatic acid exposure by shifting hemoglobin-O₂ affinity back to aquatic (control) levels when measured at an ecologically-relevant high CO₂ level that would be experienced in vivo. Using proteomics, we found that the hemoglobin subunits present in the blood did not change between treatments, but air and aquatic acid exposure altered the abundance of cathodic hemoglobin subunits. We therefore conclude that hemoglobin isoform switching is not a primary strategy used by mangrove rivulus to adjust P₅₀ under these conditions. Abundances of other RBC proteins also differed between treatment groups relative to control fish (e.g. Rhesus protein type A, band 3 anion exchanger). Overall, our data indicate that both aquatic hypercarbia and aquatic acidosis create similar changes in hemoglobin-O₂ affinity as air exposure. However, the protein-level consequences differ between these groups, indicating that the red blood cell response of mangrove rivulus can be modulated depending on the environmental cue received.

To the Land and Beyond: Crab Microbiomes as a Paradigm for the Evolution of Terrestrialization

The transition to terrestrial environments by formerly aquatic species has occurred repeatedly in many animal phyla and lead to the vast diversity of extant terrestrial species. The differences between aquatic and terrestrial habitats are enormous and involved remarkable morphological and physiological changes. Convergent evolution of various traits is evident among phylogenetically distant taxa, but almost no information is available about the role of symbiotic microbiota in such transition. Here, we suggest that intertidal and terrestrial brachyuran crabs are a perfect model to study the evolutionary pathways and the ecological role of animal-microbiome symbioses, since their transition to land is happening right now, through a number of independent lineages. The microorganisms colonizing the gut of intertidal and terrestrial crabs are expected to play a major role to conquer the land, by reducing water losses and permitting the utilization of novel food sources. Indeed, it has been shown that the microbiomes hosted in the digestive system of terrestrial isopods has been critical to digest plant items, but nothing is known about the microbiomes present in the gut of truly terrestrial crabs. Other important physiological regulations that could be facilitated by microbiomes are nitrogen excretion and osmoregulation in the new environment. We also advocate for advances in comparative and functional genomics to uncover physiological aspects of these ongoing evolutionary processes. We think that the multidisciplinary study of microorganisms associated with terrestrial crabs will shed a completely new light on the biological and physiological processes involved in the sea-land transition.

Mudskipper interleukin-34 modulates the functions of monocytes/macrophages via the colony-stimulating factor-1 receptor 1

Interleukin-34 (IL-34) is a novel cytokine that plays an important role in innate immunity and inflammatory processes by binding to the colony-stimulating factor-1 receptor (CSF-1R). However, information on the function of IL-34 in fish remains limited. In the present study, we identified an IL-34 homolog from mudskippers (Boleophthalmus pectinirostris). In silico analysis showed that the mudskipper IL-34 (BpIL-34) was similar to other known IL-34 variants in sequence and structure and was most closely related to an orange-spotted grouper (Epinephelus coioides) homolog. BpIL-34 transcripts were constitutively expressed in various tissues, with the highest level of expression found in the brain. Edwardsiella tarda infection significantly up-regulated the mRNA expression of BpIL-34 in the mudskipper tissues. The recombinant mature BpIL-34 peptide (rBpIL-34) was purified and used to produce anti-rBpIL-34 IgG. Western blot analysis combined with PNGase F digestion revealed that native BpIL-34 in monocytes/macrophages (MOs/MФs) was N-glycosylated. In vitro, rBpIL-34 treatment enhanced the phagocytotic and bactericidal activity of mudskipper MOs/MФs, as well as the mRNA expression of pro-inflammatory cytokines like tumor necrosis factor α (BpTNF-α) and BpIL-1β in these cells. Furthermore, the knockdown of mudskipper CSF-1R1 (BpCSF-1R1), but not mudskipper BpCSF-1R2, significantly inhibited the rBpIL-34-mediated enhanced effect on MO/MФ function. In conclusion, our results indicate that mudskipper BpIL-34 modulates the functions of MOs/MФs via BpCSF-1R1.

Functional characterisation of fatty acyl desaturase, Fads2, and elongase, Elovl5, in the Boddart's goggle-eyed goby Boleophthalmus boddarti (Gobiidae) suggests an incapacity for long-chain polyunsaturated fatty acid biosynthesis

The biosynthesis of long-chain polyunsaturated fatty acids (LC-PUFA), a process to convert C18 polyunsaturated fatty acids into eicosapentaenoic acid (EPA), docosahexaenoic acid (DHA) or arachidonic acid (ARA), requires the concerted activities of two enzymes, the fatty acyl desaturase (Fads) and elongase (Elovl). This study highlights the cloning, functional characterisation and tissue expression pattern of a Fads and an Elovl from the Boddart's goggle-eyed goby (Boleophthalmus boddarti), a mudskipper species widely distributed in the Indo-Pacific region. Phylogenetic analysis revealed that the cloned fads and elovl are clustered with other teleost orthologs, respectively. The investigation of the genome of several mudskipper species, namely Boleophthalmus pectinirostris, Periophthalmus schlosseri and Periophthalmus magnuspinnatus, revealed a single Fads2 and two elongases, Elovl5 and Elovl4 for each respective species. A heterologous yeast assay indicated that the B. boddarti Fads2 possessed low desaturation activity on C18 PUFA and no desaturation on C20 and C22 PUFA substrates. In comparison, the Elovl5 showed a wide range of substrate specificity, with a capacity to elongate C18, C20 and C22 PUFA substrates. An amino acid residue that affects the capacity to elongate C22:5n-3 was identified in the B. boddarti Elovl5. Both genes are highly expressed in brain tissue. Among all tissues, DHA is highly concentrated in neuron-rich tissues, whereas EPA is highly deposited in gills. Taken together, the results showed that due to the inability to perform desaturation steps, B. boddarti is unable to biosynthesise LC-PUFA, relying on dietary intake to acquire these nutrients.