Large-scale sequencing of flatfish genomes provides insights into the polyphyletic origin of their specialized body plan

Evolution of rhodopsin in flatfishes (Pleuronectiformes) is associated with depth and migratory behavior

Visual signals are involved in many fitness-related tasks and are therefore essential for survival in many species. Aquatic organisms are ideal systems to study visual evolution, as the high diversity of spectral properties in aquatic environments generates great potential for adaptation to different light conditions. Flatfishes are an economically important group, with over 800 described species distributed globally, including halibut, flounder, sole, and turbot. The diversity of flatfish species and wide array of environments they occupy provides an excellent opportunity to understand how this variation translates to molecular adaptation of vision genes. Using models of molecular evolution, we investigated how the light environments inhabited by different flatfish lineages have shaped evolution in the rhodopsin gene, which is responsible for mediating dim-light visual transduction. We found strong evidence for positive selection in rhodopsin, and this was correlated with both migratory behavior and several fundamental aspects of habitat, including depth and freshwater/marine evolutionary transitions. We also identified several mutations that likely affect the wavelength of peak absorbance of rhodopsin, and outline how these shifts in absorbance correlate with the response to the light spectrum present in different habitats. This is the first study of rhodopsin evolution in flatfishes that considers their extensive diversity, and our results highlight how ecologically-driven molecular adaptation has occurred across this group in response to transitions to novel light environments.

Genome sequencing and analysis of black flounder (Paralichthys orbignyanus) reveals new insights into Pleuronectiformes genomic size and structure

Black flounder (Paralichthys orbignyanus, Pleuronectiformes) is a commercially significant marine fish with promising aquaculture potential in Argentina. Despite extensive studies on Black flounder aquaculture, its limited genetic information available hampers the crucial role genetics plays in the development of this activity. In this study, we first employed Illumina sequencing technology to sequence the entire genome of Black flounder. Utilizing two independent libraries-one from a female and another from a male-with 150 bp paired-end reads, a mean insert length of 350 bp, and over 35 X-fold coverage, we achieved assemblies resulting in a genome size of ~ 538 Mbp. Analysis of the assemblies revealed that more than 98% of the core genes were present, with more than 78% of them having more than 50% coverage. This indicates a somehow complete and accurate genome at the coding sequence level. This genome contains 25,231 protein-coding genes, 445 tRNAs, 3 rRNAs, and more than 1,500 non-coding RNAs of other types. Black flounder, along with pufferfishes, seahorses, pipefishes, and anabantid fish, displays a smaller genome compared to most other teleost groups. In vertebrates, the number of transposable elements (TEs) is often correlated with genome size. However, it remains unclear whether the sizes of introns and exons also play a role in determining genome size. Hence, to elucidate the potential factors contributing to this reduced genome size, we conducted a comparative genomic analysis between Black flounder and other teleost orders to determine if the small genomic size could be explained by repetitive elements or gene features, including the whole genome genes and introns sizes. We show that the smaller genome size of flounders can be attributed to several factors, including changes in the number of repetitive elements, and decreased gene size, particularly due to lower amount of very large and small introns. Thus, these components appear to be involved in the genome reduction in Black flounder. Despite these insights, the full implications and potential benefits of genome reduction in Black flounder for reproduction and aquaculture remain incompletely understood, necessitating further research.

Diverse origins of near-identical antifreeze proteins in unrelated fish lineages provide insights into evolutionary mechanisms of new gene birth and protein sequence convergence

Determining the origins of novel genes and the genetic mechanisms underlying the emergence of new functions is challenging yet crucial for understanding evolutionary innovations. The novel fish antifreeze proteins, exemplifying convergent evolution, represent excellent opportunities to investigate the evolutionary origins and pathways of new genes. Particularly notable is the near-identical type I antifreeze proteins (AFPI) in four phylogenetically divergent fish taxa. This study tested the hypothesis of protein sequence convergence beyond functional convergence in three unrelated AFPI-bearing fish lineages, revealing different paths by which a similar protein arose from diverse genomic resources. Comprehensive comparative analyses of de novo sequenced genome of the winter flounder and grubby sculpin, available high-quality genome of the cunner, and those of 14 other relevant species found that the near-identical AFPI originated from a distinct genetic precursor in each lineage, and independently evolved coding regions for the novel ice-binding protein while retaining sequence identity in the regulatory regions with their respective ancestor. The deduced evolutionary processes and molecular mechanisms is consistent with the Innovation-Amplification-Divergence (IAD) model applicable to AFPI formation in all three lineages, a new Duplication-Degeneration-Divergence (DDD) model we propose for the sculpin lineage, and a DDD model with gene fission for the cunner lineage. This investigation illustrates the multiple ways by which a novel functional gene with sequence convergence at the protein level could evolve across divergent species, advancing our understanding of the mechanistic intricacies in new gene formation.

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.

Single-cell transcriptomics reveals the brain evolution of web-building spiders

Spiders are renowned for their efficient capture of flying insects using intricate aerial webs. How the spider nervous systems evolved to cope with this specialized hunting strategy and various environmental clues in an aerial space remains unknown. Here we report a brain-cell atlas of >30,000 single-cell transcriptomes from a web-building spider (Hylyphantes graminicola). Our analysis revealed the preservation of ancestral neuron types in spiders, including the potential coexistence of noradrenergic and octopaminergic neurons, and many peptidergic neuronal types that are lost in insects. By comparing the genome of two newly sequenced plesiomorphic burrowing spiders with three aerial web-building spiders, we found that the positively selected genes in the ancestral branch of web-building spiders were preferentially expressed (42%) in the brain, especially in the three mushroom body-like neuronal types. By gene enrichment analysis and RNAi experiments, these genes were suggested to be involved in the learning and memory pathway and may influence the spiders' web-building and hunting behaviour. Our results provide key sources for understanding the evolution of behaviour in spiders and reveal how molecular evolution drives neuron innovation and the diversification of associated complex behaviours.

Insights into the Mitochondrial Genetic Makeup and Miocene Colonization of Primitive Flatfishes (Pleuronectiformes: Psettodidae) in the East Atlantic and Indo-West Pacific Ocean

The mitogenomic evolution of the Psettodes flatfishes is still poorly known from their range distribution in eastern Atlantic and Indo-West Pacific Oceans. The study delves into the matrilineal evolutionary pathway of these primitive flatfishes, with a specific focus on the complete mitogenome of the Psettodes belcheri species, as determined through next-generation sequencing. The mitogenome in question spans a length of 16,747 base pairs and comprises a total of 37 genes, including 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a control region. Notably, the mitogenome of P. belcheri exhibits a bias towards AT base pairs, with a composition of 54.15%, mirroring a similar bias observed in its close relative, Psettodes erumei, which showcases percentages of 53.07% and 53.61%. Most of the protein-coding genes commence with an ATG initiation codon, except for Cytochrome c oxidase I (COI), which initiates with a GTG codon. Additionally, four protein-coding genes commence with a TAA termination codon, while seven others exhibit incomplete termination codons. Furthermore, two protein-coding genes, namely NAD1 and NAD6, terminate with AGG and TAG stop codons, respectively. In the mitogenome of P. belcheri, the majority of transfer RNAs demonstrate the classical cloverleaf secondary structures, except for tRNA-serine, which lacks a DHU stem. Comparative analysis of conserved blocks within the control regions of two Psettodidae species unveiled that the CSB-II block extended to a length of 51 base pairs, surpassing the other blocks and encompassing highly variable sites. A comprehensive phylogenetic analysis using mitochondrial genomes (13 concatenated PCGs) categorized various Pleuronectiformes species, highlighting the basal position of the Psettodidae family and showed monophyletic clustering of Psettodes species. The approximate divergence time (35-10 MYA) between P. belcheri and P. erumei was estimated, providing insights into their separation and colonization during the early Miocene. The TimeTree analysis also estimated the divergence of two suborders, Psettodoidei and Pleuronectoidei, during the late Paleocene to early Eocene (56.87 MYA). The distribution patterns of Psettodes flatfishes were influenced by ocean currents and environmental conditions, contributing to their ecological speciation. In the face of climate change and anthropogenic activities, the conservation implications of Psettodes flatfishes are emphasized, underscoring the need for regulated harvesting and adaptive management strategies to ensure their survival in changing marine ecosystems. Overall, this study contributes to understanding the evolutionary history, genetic diversity, and conservation needs of Psettodes flatfishes globally. However, the multifaceted exploration of mitogenome and larger-scale genomic data of Psettodes flatfish will provide invaluable insights into their genetic characterization, evolutionary history, environmental adaptation, and conservation in the eastern Atlantic and Indo-West Pacific Oceans.

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.

A review of the pangenome: how it affects our understanding of genomic variation, selection and breeding in domestic animals?

As large-scale genomic studies have progressed, it has been revealed that a single reference genome pattern cannot represent genetic diversity at the species level. While domestic animals tend to have complex routes of origin and migration, suggesting a possible omission of some population-specific sequences in the current reference genome. Conversely, the pangenome is a collection of all DNA sequences of a species that contains sequences shared by all individuals (core genome) and is also able to display sequence information unique to each individual (variable genome). The progress of pangenome research in humans, plants and domestic animals has proved that the missing genetic components and the identification of large structural variants (SVs) can be explored through pangenomic studies. Many individual specific sequences have been shown to be related to biological adaptability, phenotype and important economic traits. The maturity of technologies and methods such as third-generation sequencing, Telomere-to-telomere genomes, graphic genomes, and reference-free assembly will further promote the development of pangenome. In the future, pangenome combined with long-read data and multi-omics will help to resolve large SVs and their relationship with the main economic traits of interest in domesticated animals, providing better insights into animal domestication, evolution and breeding. In this review, we mainly discuss how pangenome analysis reveals genetic variations in domestic animals (sheep, cattle, pigs, chickens) and their impacts on phenotypes and how this can contribute to the understanding of species diversity. Additionally, we also go through potential issues and the future perspectives of pangenome research in livestock and poultry.

Integration of Maps Enables a Cytogenomics Analysis of the Complete Karyotype in Solea senegalensis

The Pleuronectiformes order, which includes several commercially-important species, has undergone extensive chromosome evolution. One of these species is Solea senegalensis, a flatfish with 2n = 42 chromosomes. In this study, a cytogenomics approach and integration with previous maps was applied to characterize the karyotype of the species. Synteny analysis of S. senegalensis was carried out using two flatfish as a reference: Cynoglossus semilaevis and Scophthalmus maximus. Most S. senegalensis chromosomes (or chromosome arms for metacentrics and submetacentrics) showed a one-to-one macrosyntenic pattern with the other two species. In addition, we studied how repetitive sequences could have played a role in the evolution of S. senegalensis bi-armed (3, and 5–9) and acrocentric (11, 12 and 16) chromosomes, which showed the highest rearrangements compared with the reference species. A higher abundance of TEs (Transposable Elements) and other repeated elements was observed adjacent to telomeric regions on chromosomes 3, 7, 9 and 16. However, on chromosome 11, a greater abundance of DNA transposons was detected in interstitial BACs. This chromosome is syntenic with several chromosomes of the other two flatfish species, suggesting rearrangements during its evolution. A similar situation was also found on chromosome 16 (for microsatellites and low complexity sequences), but not for TEs (retroelements and DNA transposons). These differences in the distribution and abundance of repetitive elements in chromosomes that have undergone remodeling processes during the course of evolution also suggest a possible role for simple repeat sequences in rearranged regions.

Distinguishing Evolutionary Conservation from Derivedness

While the concept of “evolutionary conservation” has enabled biologists to explain many ancestral features and traits, it has also frequently been misused to evaluate the degree of changes from a common ancestor, or “derivedness”. We propose that the distinction of these two concepts allows us to properly understand phenotypic and organismal evolution. From a methodological aspect, “conservation” mainly considers genes or traits which species have in common, while “derivedness” additionally covers those that are not commonly shared, such as novel or lost traits and genes to evaluate changes from the time of divergence from a common ancestor. Due to these differences, while conservation-oriented methods are effective in identifying ancestral features, they may be prone to underestimating the overall changes accumulated during the evolution of certain lineages. Herein, we describe our recently developed method, “transcriptomic derivedness index”, for estimating the phenotypic derivedness of embryos with a molecular approach using the whole-embryonic transcriptome as a phenotype. Although echinoderms are often considered as highly derived species, our analyses with this method showed that their embryos, at least at the transcriptomic level, may not be much more derived than those of chordates. We anticipate that the future development of derivedness-oriented methods could provide quantitative indicators for finding highly/lowly evolvable traits.

Identification of Potential Blind-Side Hypermelanosis-Related lncRNA–miRNA–mRNA Regulatory Network in a Flatfish Species, Chinese Tongue Sole (Cynoglossus semilaevis)

Blind-side hypermelanosis has emerged as a major concern in commercial rearing environments of the flatfish aquaculture industry. To date, the underlying molecular mechanisms are not well understood. To fill this gap, in this study, whole transcriptomic sequencing and analyses were performed using normal skins and hypermelanic skins of the blind side of Chinese tongue sole (Cynoglossus semilaevis). Differentially expressed long non-coding RNAs (DElncRNAs), miRNAs (DEmiRNAs), and differentially expressed genes as well as their competing endogenous RNA (ceRNA) networks were identified. A total of 34 DElncRNAs, 226 DEmiRNAs, and 610 DEGs were identified. Finally, lncRNA–miRNA–mRNA regulatory networks (involving 29 DElncRNAs, 106 DEmiRNAs, and 162 DEGs) associated with blind-side hypermelanosis were constructed. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of 162 DEGs in ceRNA networks identified DEGs (e.g., oca2, mc1r, and ihhb) in pigmentation-related biological processes and DEGs (e.g., ca4, glul, and fut9) in nitrogen metabolism, glycosphingolipid biosynthesis, and folate biosynthesis pathways, as well as their corresponding DElncRNAs and DEmiRNAs to potentially play key regulatory roles in blind-side hypermelanosis. In conclusion, this is the first study on the ceRNA regulatory network associated with blind-side hypermelanosis in flatfish. These new findings expand the spectrum of non-coding regulatory mechanisms underpinning blind-side hypermelanosis, which facilitates the further exploration of molecular regulatory mechanisms of malpigmentation in flatfish.

Overview and New Insights Into the Diversity, Evolution, Role, and Regulation of Kisspeptins and Their Receptors in Teleost Fish

In the last two decades, kisspeptin (Kiss) has been identified as an important player in the regulation of reproduction and other physiological functions in vertebrates, including several fish species. To date, two ligands (Kiss1, Kiss2) and three kisspeptin receptors (Kissr1, Kissr2, Kissr3) have been identified in teleosts, likely due to whole-genome duplication and loss of genes that occurred early in teleost evolution. Recent results in zebrafish and medaka mutants have challenged the notion that the kisspeptin system is essential for reproduction in fish, in marked contrast to the situation in mammals. In this context, this review focuses on the role of kisspeptins at three levels of the reproductive, brain-pituitary-gonadal (BPG) axis in fish. In addition, this review compiled information on factors controlling the Kiss/Kissr system, such as photoperiod, temperature, nutritional status, sex steroids, neuropeptides, and others. In this article, we summarize the available information on the molecular diversity and evolution, tissue expression and neuroanatomical distribution, functional significance, signaling pathways, and gene regulation of Kiss and Kissr in teleost fishes. Of particular note are recent advances in understanding flatfish kisspeptin systems, which require further study to reveal their structural and functional diversity.

New Evidence for the Existence of Two Kiss/Kissr Systems in a Flatfish Species, the Turbot (Scophthalmus maximus), and Stimulatory Effects on Gonadotropin Gene Expression

Seasonal reproduction is generally controlled by the hypothalamus-pituitary-gonadal (HPG) axis in fish. Previous studies have demonstrated that the kisspeptin (Kiss)/kisspeptin receptor (Kissr) system, a positive regulator of the HPG axis, mediates the responses to environmental cues. Turbot (Scophthalmus maximus), a representative species of Pleuronectiformes, is one of the most commercially important fish species cultured in Europe and North China. However, the mechanisms by which the Kiss/Kissr system regulates the reproductive axis of turbot according to seasonal changes, especially photoperiod, have not been clearly characterized. In the current study, the cDNA sequences of kiss2/kissr2, along with kiss1/kissr3 which was thought to be lost in flatfish species, were cloned and functionally characterized. The kiss1, kiss2, and kissr3 transcripts were highly detected in the brain and gonad, while kissr2 mRNA was only abundantly expressed in the brain. Moreover, kiss/kissr mRNAs were further examined in various brain areas of both sexes. The kiss1, kissr2, kissr3 mRNAs were highly expressed in the mesencephalon, while a substantial degree of kiss2 transcripts were observed in the hypothalamus. During annual reproductive cycle, both kiss and kissr transcript levels declined significantly from the immature to mature stages and increased at the degeneration stage in the brains of both sexes, especially in the mesencephalon and hypothalamus. The ovarian kiss1, kiss2, and kissr2 mRNA levels were highest at the vitellogenic stage (mature stage), while expression of kissr3 was highest at the immature stage. The testicular kiss and kissr transcripts were highest in the immature and degeneration stages, and lowest at the mature stage. In addition, intraperitoneal injection of Kiss1-10 and Kiss2-10 significantly stimulated mRNA levels of pituitary lhβ, fhsβ, and gthα. In summary, two Kiss/Kissr systems were firstly proven in a flatfish species of turbot, and it has a positive involvement in controlling the reproduction of the Kiss/Kissr system in turbot. The results will provide preliminary information regarding how the Kiss/Kissr system controls seasonal reproduction in turbot broodstock.

Pattern of New Gene Origination in a Special Fish Lineage, the Flatfishes

Origination of new genes are of inherent interest of evolutionary geneticists for decades, but few studies have addressed the general pattern in a fish lineage. Using our recent released whole genome data of flatfishes, which evolved one of the most specialized body plans in vertebrates, we identified 1541 (6.9% of the starry flounder genes) flatfish-lineage-specific genes. The origination pattern of these flatfish new genes is largely similar to those observed in other vertebrates, as shown by the proportion of DNA-mediated duplication (1317; 85.5%), RNA-mediated duplication (retrogenes; 96; 6.2%), and de novo-origination (128; 8.3%). The emergence rate of species-specific genes is 32.1 per Mya and the whole average level rate for the flatfish-lineage-specific genes is 20.9 per Mya. A large proportion (31.4%) of these new genes have been subjected to selection, in contrast to the 4.0% in primates, while the old genes remain quite similar (66.4% vs. 65.0%). In addition, most of these new genes (70.8%) are found to be expressed, indicating their functionality. This study not only presents one example of systematic new gene identification in a teleost taxon based on comprehensive phylogenomic data, but also shows that new genes may play roles in body planning.

The complete mitochondrial genome of Toxotes chatareus (Toxotes; Toxotidae; Carangaria) assembled by the next-generation sequencing data and phylogenetic analysis of Carangaria

We present the complete mitochondrial genome of Toxotes chatareus yielded by the next-generation sequencing data in this study. The complete mitochondrial genome of T. chatareus has 16,543 bp and contained 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and a single control region (D-loop). The overall base composition was A 28.75%, C 29.80%, G 15.77%, T 25.68% and its gene arrangement was similar with other Carangaria mitochondrial genomes. Additionally, the phylogenetic relationships of 13 Carangaria species based on the complete mitochondrial genome was analyzed using the neighbor-joining method. The result showed T. chatareus was clustered with L. lactarius suggesting the close phylogenetic affinity they owned. Together, the complete mitochondrial genome of T. chatareus would be beneficial for the study of phylogenetic relationship, taxonomic classification and phylogeography of the Carangaria.

Developmental aspects of the hypothalamic-pituitary network related to reproduction in teleost fish

The hypothalamic-pituitary-gonadal axis is the main system that regulates reproduction in vertebrates through a complex network that involves different neuropeptides, neurotransmitters, and pituitary hormones. Considering that this axis is established early on life, the main goal of the present work is to gather information on its development and the actions of its components during early life stages. This review focuses on fish because their neuroanatomical characteristics make them excellent models to study neuroendocrine systems. The following points are discussed: i) developmental functions of the neuroendocrine components of this network, and ii) developmental disruptions that may impact adult reproduction. The importance of the components of this network and their susceptibility to external/internal signals that can alter their specific early functions and/or even the establishment of the reproductive axis, indicate that more studies are necessary to understand this complex and dynamic network.

Exon-capture data and locus screening provide new insights into the phylogeny of flatfishes (Pleuronectoidei)

There is an extensive collection of literature on the taxonomy and phylogenetics of flatfishes (Pleuronectiformes) that extends over two centuries, but consensus on many of their evolutionary relationships remains elusive. Phylogenetic uncertainty stems from highly divergent results derived from morphological and genetic characters, and between various molecular datasets. Deciphering relationships is complicated by rapid diversification early in the Pleuronectiformes tree and an abundance of studies that incompletely and inconsistently sample taxa and genetic markers. We present phylogenies based on a genome-wide dataset (4,434 nuclear markers via exon-capture) and wide taxon sampling (86 species spanning 12 of 16 families) of the largest flatfish suborder (Pleuronectoidei). Nine different subsets of the data and two tree construction approaches (eighteen phylogenies in total) are remarkably consistent with other recent molecular phylogenies, and show strong support for the monophyly of all families included except Pleuronectidae. Analyses resolved a novel phylogenetic hypothesis for the family Rhombosoleidae as being within the Pleuronectoidea rather than the Soleoidea, and failed to support the subfamily Hippoglossinae as a monophyletic group. Our results were corroborated with evidence from previous phylogenetic studies to outline regions of persistent phylogenetic uncertainty and identify groups in need of further phylogenetic inference.