Major events in the genome evolution of vertebrates: paranome age and size differ considerably between ray-finned fishes and land vertebrates

Exploration of Appetite Regulation in Yangtze Sturgeon (Acipenser dabryanus) During Weaning

Yangtze sturgeon is an endangered fish species. After weaning, some Yangtze sturgeon fry refuse to consume any food, which causes a low survival rate during the artificial breeding period. This study showed that the body length and body weight of failed weaning Yangtze sturgeons were significantly lower than those of successful weaning sturgeons. Since the brain is the center of appetite regulation, RNA-seq of the brain was employed to analyze the differentially expressed genes and their biological functions in successfully and unsuccessfully weaned fry. After that, 82,151 unigenes and 3222 DEGs were obtained. Based on the results of RNA-seq, appetite factors, including POMC, CART, NPY and AgRP, were cloned, and then a weaning experiment was designed to explore the changes in appetite after feeding a microcapsule diet (weaning group). The results showed that, during the weaning period, the expression of CART was increased on the 1st and 3rd days but decreased on the 5th, 6th, 8th and 10th days. The expression of AgRP was downregulated on the 1st and 3rd days but upregulated on the 5th, 6th, 8th and 10th days. These findings indicate that appetite was suppressed in the early and middle periods but enhanced in the latter period of weaning and that CART may play an important role in the appetite-suppressing effect.

Gone with the Species: From Gene Loss to Gene Extinction

BACKGROUND

Vertebrae protein-coding genes exhibit remarkable diversity and are organized into many gene families. These gene families have emerged through various gene duplication events, the most prominent being the two rounds of whole-genome duplication (WGD). The current research project analyzed a unique class of genes called "singletons". Notably, we introduce the concept of "super-singletons": genes that stand as the last representatives of their ancestral families and the sole representatives of their genetic makeup with no ortholog in any other species.

METHODS

We used the Ensembl/Biomart pipeline to identify duplicated and unduplicated protein-coding genes in different vertebrate species and found orthologs of human genes.

RESULTS

We showed the frequency of duplicated genes and singletons, demonstrating that singletons are more vulnerable to evolutionary loss than duplicated genes. Additionally, we found that contractions in vertebrate gene families are more prevalent than expansion.

CONCLUSION

Our study provides insight into the evolution of gene families and presents a novel scenario where the extinction of species would lead to the extinction of a gene, ultimately shifting the narrative from the impact of genetics on species extinction to the extinction of genes.

Zebrafish: unraveling genetic complexity through duplicated genes

The zebrafish is an invaluable model organism for genetic, developmental, and disease research. Although its high conservation with humans is often cited as justification for its use, the zebrafish harbors oft-ignored genetic characteristics that may provide unique insights into gene structure and function. Zebrafish, along with other teleost fish, underwent an additional round of whole genome duplication after their split from tetrapods-resulting in an abundance of duplicated genes when compared to other vertebrates. These duplicated genes have evolved in distinct ways over the ensuing 350 million years. Thus, each gene within a duplicated gene pair has nuanced differences that create a unique identity. By investigating both members of the gene pair together, we can elucidate the mechanisms that underly protein structure and function and drive the complex interplay within biological systems, such as signal transduction cascades, genetic regulatory networks, and evolution of tissue and organ function. It is crucial to leverage such studies to explore these molecular dynamics, which could have far-reaching implications for both basic science and therapeutic development. Here, we will review the role of gene duplications and the existing models for gene divergence and retention following these events. We will also highlight examples within each of these models where studies comparing duplicated genes in the zebrafish have yielded key insights into protein structure, function, and regulation.

Genome-wide identification of Gα family in grass carp (Ctenopharyngodon idella) and reproductive regulation functional characteristics of Cignaq

BACKGROUND

The Gα family plays a crucial role in the complex reproductive regulatory network of teleosts. However, the characterization and function of Gα family members, especially Gαq, remain poorly understood in teleosts. To analyze the characterization, expression, and function of grass carp (Ctenopharyngodon idella) Gαq, we identified the Gα family members in grass carp genome, and analyzed the expression, distribution, and signal transduction of Gαq/gnaq. We also explored the role of Gαq in the reproductive regulation of grass carp.

RESULTS

Our results showed that the grass carp genome contains 27 Gα genes with 46 isoforms, which are divided into four subfamilies: Gαs, Gαi/o, Gαq/11, and Gα12/13. The expression level of Cignaq in the testis was the highest and significantly higher than in other tissues, followed by the hypothalamus and brain. The luteinizing hormone receptor (LHR) was mainly localized to the nucleus in grass carp oocytes, with signals also present in follicular cells. In contrast, Gαq signal was mainly found in the cytoplasm of oocytes, with no signal in follicular cells. In the testis, Gαq and LHR were co-localized in the cytoplasm. Furthermore, the grass carp Gαq recombinant protein significantly promoted Cipgr expression.

CONCLUSIONS

These results provided preliminary evidence for understanding the role of Gαq in the reproductive regulation of teleosts.

Molecular Dating of the Teleost Whole Genome Duplication (3R) Is Compatible With the Expectations of Delayed Rediploidization

Vertebrate evolution has been punctuated by three whole genome duplication events that have been implicated causally in phenotypic evolution, from the origin of phenotypic novelties to explosive diversification. Arguably, the most dramatic of these is the 3R whole genome duplication event associated with the origin of teleost fishes which comprise more than half of all living vertebrate species. However, tests of a causal relationship between whole genome duplication and teleost diversification have proven difficult due to the challenge of establishing the timing of these phenomena. Here we show, based on molecular clock dating of concatenated gene alignments, that the 3R whole genome duplication event occurred in the early-middle Permian (286.18 to 267.20 million years ago; Ma), 52.02 to 12.84 million years (Myr) before the divergence of crown-teleosts in the latest Permian-earliest Late Triassic (254.36 to 234.16 Ma) and long before the major pulses of teleost diversification in Ostariophysi and Percomorpha (56.37 to 100.17 Myr and at least 139.24 to 183.29 Myr later, respectively). The extent of this temporal gap between putative cause and effect precludes 3R as a deterministic driver of teleost diversification. However, these age constraints remain compatible with the expectations of a prolonged rediploidization process following whole genome duplication which, through the effects of chromosome rearrangement and gene loss, remains a viable mechanism to explain the evolution of teleost novelties and diversification.

Whole-genome Sequencing Reveals Autooctoploidy in Chinese Sturgeon and Its Evolutionary Trajectories

The order Acipenseriformes, which includes sturgeons and paddlefishes, represents "living fossils" with complex genomes that are good models for understanding whole-genome duplication (WGD) and ploidy evolution in fishes. Here, we sequenced and assembled the first high-quality chromosome-level genome for the complex octoploid Acipenser sinensis (Chinese sturgeon), a critically endangered species that also represents a poorly understood ploidy group in Acipenseriformes. Our results show that A. sinensis is a complex autooctoploid species containing four kinds of octovalents (8n), a hexavalent (6n), two tetravalents (4n), and a divalent (2n). An analysis taking into account delayed rediploidization reveals that the octoploid genome composition of Chinese sturgeon results from two rounds of homologous WGDs, and further provides insights into the timing of its ploidy evolution. This study provides the first octoploid genome resource of Acipenseriformes for understanding ploidy compositions and evolutionary trajectories of polyploid fishes.

Conservation of Affinity Rather Than Sequence Underlies a Dynamic Evolution of the Motif-Mediated p53/MDM2 Interaction in Ray-Finned Fishes

The transcription factor and cell cycle regulator p53 is marked for degradation by the ubiquitin ligase MDM2. The interaction between these 2 proteins is mediated by a conserved binding motif in the disordered p53 transactivation domain (p53TAD) and the folded SWIB domain in MDM2. The conserved motif in p53TAD from zebrafish displays a 20-fold weaker interaction with MDM2, compared to the interaction in human and chicken. To investigate this apparent difference, we tracked the molecular evolution of the p53TAD/MDM2 interaction among ray-finned fishes (Actinopterygii), the largest vertebrate clade. Intriguingly, phylogenetic analyses, ancestral sequence reconstructions, and binding experiments showed that different loss-of-affinity changes in the canonical binding motif within p53TAD have occurred repeatedly and convergently in different fish lineages, resulting in relatively low extant affinities (KD = 0.5 to 5 μM). However, for 11 different fish p53TAD/MDM2 interactions, nonconserved regions flanking the canonical motif increased the affinity 4- to 73-fold to be on par with the human interaction. Our findings suggest that compensating changes at conserved and nonconserved positions within the motif, as well as in flanking regions of low conservation, underlie a stabilizing selection of "functional affinity" in the p53TAD/MDM2 interaction. Such interplay complicates bioinformatic prediction of binding and calls for experimental validation. Motif-mediated protein-protein interactions involving short binding motifs and folded interaction domains are very common across multicellular life. It is likely that the evolution of affinity in motif-mediated interactions often involves an interplay between specific interactions made by conserved motif residues and nonspecific interactions by nonconserved disordered regions.

Re-identification and characterization of grass carp Ctenopharyngodon idella TLR20

Toll-like receptors (TLRs) play a crucial role in the recognition of microbial-associated molecular patterns in the innate immune system. Fish TLRs have undergone significant gene expansion to adapt to complex aquatic environments. Among them, TLR20 from the TLR11 family actively responds to viral and bacterial invasions. Previous studies have reported two TLR20s in grass carp (Ctenopharyngodon idella), and in this study, we revised this conclusion. Based on the latest grass carp genome, we identified a new TLR20 member. These three TLR20s are arranged in tandem on chromosome 9, indicating that they are generated by gene duplication events. They were renamed CiTLR20.1 to CiTLR20.3 based on their chromosomal positions. The CiTLR20s in C. idella exhibit higher similarities with those in Danio rerio, Cyprinus carpio, and Megalobrama amblycephala, and lower similarities with those in other distantly related fish species. Selective pressure analysis revealed low conservation and negative evolution of TLR20s during evolution. The 3D structures of the three TLR20s showed significant differences, reflecting functional variations and different downstream adaptor molecule recruitment. Transcriptome data revealed tissue distribution differences of TLR20s, with TLR20.1 showing relatively low expression levels in all the tissues, while TLR20.2 and TLR20.3 showed higher expression in the head kidney, spleen, and gill. Additionally, TLR20.2 and TLR20.3 actively responded to GCRV-II infection, with higher upregulation of TLR20.2 in response to Aeromonas hydrophila challenge. In conclusion, this study corrected the number of grass carp TLR20 members and analyzed TLR20 from an evolutionary and structural perspective, exploring its role in antiviral and antibacterial defense. This study provides reference for future research on fish TLR20.

Neuropeptide B (NPB) and NPB receptor 2b (NPBWR2b) in the ricefield eel Monopterus albus: expression and potential involvement in the regulation of gonadotropins

The neuropeptide B/W signaling system is composed of neuropeptide B (NPB), neuropeptide W (NPW), and two cognate receptors, NPBWR1 and NPBWR2, which are involved in diverse physiological processes, including the central regulation of neuroendocrine axes in vertebrates. The components of this signaling system are not well conserved during vertebrate evolution, implicating its functional diversity. The present study characterized the ricefield eel neuropeptide B/W system, generated a specific antiserum against the neuropeptide B/W receptor, and examined the potential roles of the system in the regulation of adenohypophysial functions. The ricefield eel genome contains npba, npbb, and npbwr2b but lacks the npw, npbwr1, and npbwr2a genes. The loss of npw and npbwr1 probably occurred at the base of ray-finned fish radiation and that of npbwr2a species specifically in ray-finned fish. Npba and npbb genes are produced through whole-genome duplication (WGD) in ray-finned fish. The ricefield eel npba was expressed in the brain and some peripheral tissues, while npbb was predominantly expressed in the brain. The ricefield eel npbwr2b was also expressed in the brain and in some peripheral tissues, such as the pituitary, gonad, heart, and eye. Immunoreactive Npbwr2b was shown to be localized to Lh and Fsh cells but not to Gh or Prl cells in the pituitary of ricefield eels. Npba upregulated the expression of fshb and cga but not lhb mRNA in pituitary fragments of ricefield eels cultured in vitro. The results of the present study suggest that the NPB system of ricefield eels may be involved in the neuroendocrine regulation of reproduction.

Interlocus Gene Conversion, Natural Selection, and Paralog Homogenization

Following a duplication, the resulting paralogs tend to diverge. While mutation and natural selection can accelerate this process, they can also slow it. Here, we quantify the paralog homogenization that is caused by point mutations and interlocus gene conversion (IGC). Among 164 duplicated teleost genes, the median percentage of postduplication codon substitutions that arise from IGC rather than point mutation is estimated to be between 7% and 8%. By differentiating between the nonsynonymous codon substitutions that homogenize the protein sequences of paralogs and the nonhomogenizing nonsynonymous substitutions, we estimate the homogenizing nonsynonymous rates to be higher for 163 of the 164 teleost data sets as well as for all 14 data sets of duplicated yeast ribosomal protein-coding genes that we consider. For all 14 yeast data sets, the estimated homogenizing nonsynonymous rates exceed the synonymous rates.

Revealing the potential regulatory relationship between HSP70, HSP90 and HSF genes under temperature stress

Heat Shock Protein (HSPs) gene family members play fundamental roles in different environmental stress tolerances, protect the structure and function of cells, and perform a significant task in cellular homeostasis. In this study, we conducted a genome-wide identification, evolutionary relationship analysis and gene expression analysis of the HSP70, HSP90, and HSF gene families in Ruditapes philippinarum. We identified 83 RpHSP70, 6 RpHSP90, and 3 RpHSF genes in R. philippinarum. The structural characteristics, chromosomal localization, and the gene structure map were constructed to reveal the characteristics of protein structures. Furthermore, the expression profiling of transcriptome data showed the expression pattern of HSP70, HSP90 and HSF genes in Manila clam from different populations, and under high and low temperature stress. In addition, we performed protein-protein interaction network analysis between HSP70, HSP90, and HSF gene family which enabled us to recognize the regulatory relationship between the two HSP gene families and the HSF gene family. Furthermore, the predicted sub-cellular location revealed a diversified subcellular distribution of HSP70, HSP90, and HSF proteins, which may be directly or indirectly associated with functional diversification under heat stress condition.

The evolution and immunomodulatory role of Zc3h12 proteins in zebrafish (Danio rerio)

Zc3h12 family is an important RNA-binding protein family regulating mRNA of inflammatory cytokines in mammals. However, there are few studies on their post-transcriptional level regulation of inflammatory cytokines in fish. Here, we investigated the evolution of zebrafish Zc3h12 family and explored their immunomodulatory role. Phylogenetic and syntenic analysis indicated the number of zc3h12 family members had increased ranging from a single member in invertebrates to a single copy of four members in mammals. As the most evolutionarily diverse group of vertebrates, the number of zc3h12 family members was more complex and diverse in the teleost, each member experienced different fates and followed different rules in multiple rounds of whole-genome duplication events. Thereinto, zebrafish contained three zc3h12 genes, among which zc3h12aa and zc3h12ab were duplicated from the same gene. Zebrafish Zc3h12 family could recognize the 3'-UTR regions of inflammatory cytokines through binding to the specific RNA secondary structure and negatively regulate their expression. Deletion of either Zc3h12 domains or mutation of the key amino acid in RNAase domain attenuated their modulatory effect, suggesting both domain and RNAase activity are important to the immunomodulatory role. These results elucidated the evolution of Zc3h12 family and uncovered Zc3h12-mediated post-transcriptional regulation of cytokines in zebrafish.

Navigating the CoGe Online Software Suite for Polyploidy Research

The CoGe software suite at genomevolution.org hosts a number of tools that facilitate genomic research on plant and animal whole-genome multiplication-polyploidy. SynMap permits analysis and visualization of two-way syntenic dotplot alignments of genomes, includes many options and data/graphics download possibilities, and even permits three-genome synteny maps and interactive views. FractBias is a tool that operates within SynMap that permits calculation and graphic display of genome fragments (such as chromosomes) of one species mapped to another, displaying both blockwise homology depths and the extent of syntenic gene (syntelog) loss following polyploidy events. SynMap macrosynteny results can segue into the microsynteny tool GEvo, which provides genome-browser-like views of homologous genome blocks. CoGe FeatView allows call-up of given gene features already stored in the CoGe resource, and CoGeBlast permits searches for additional features that can be analyzed or downloaded further. Links from these tools can be fed into SynFind, which can find syntenic blocks surrounding a feature across multiple specified genomes while also simultaneously providing overall genome-wide syntenic depth calculations that can be interpreted to reflect polyploidy levels. Here, we describe basic use of these tools on the CoGe software suite.

Genome-wide identification and characterization of toll-like receptor 5 (TLR5) in fishes

Toll-like receptors 5 (TLR5), a member of the toll-like receptors (TLRs) family, is a class of pattern recognition receptors (PRRs) that recognize pathogen-associated molecular patterns (PAMPs). It responds to vertebrate recognition of bacterial flagellin and participates in innate immune responses. However, genome-wide identification and characterization of TLR5 in fishes have not been investigated. Here, three TLR5M isotypes (TLR5Ma, TLR5Mb1, and TLR5Mb2) and a TLR5S are all extracted from fish genomes on the basis of phylogenetic and synteny analyses. We confirmed that the non-teleost fishes have one TLR5M gene, as well as additional TLR5 genes (TLR5M and TLR5S) in teleost fishes. In addition, some special teleost fishes possess two to three TLR5 genes, which have undergone the fourth whole-genome duplication (WGD). According to our results, we inferred that the diversity of TLR5 genes in fishes seems to be the result of combinations of WGD and gene loss. Furthermore, TLR5 isoforms displayed differences at the flagellin interaction sites and viral binding sites, and showed lineage-specific, which indicated that TLR5 duplicates may generate functional divergence. Bacterial experiments also supported the idea that CiTLR5Ma and CiTLR5Mb are subfunctionalized to sense bacterial flagellin. In summary, our present comparative genomic survey will benefit for further functional investigations of TLR5 genes in fish.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

Two duplicated gsdf homeologs cooperatively regulate male differentiation by inhibiting cyp19a1a transcription in a hexaploid fish

Although evolutionary fates and expression patterns of duplicated genes have been extensively investigated, how duplicated genes co-regulate a biological process in polyploids remains largely unknown. Here, we identified two gsdf (gonadal somatic cell-derived factor) homeologous genes (gsdf-A and gsdf-B) in hexaploid gibel carp (Carassius gibelio), wherein each homeolog contained three highly conserved alleles. Interestingly, gsdf-A and gsdf-B transcription were mainly activated by dmrt1-A (dsx- and mab-3-related transcription factor 1) and dmrt1-B, respectively. Loss of either gsdf-A or gsdf-B alone resulted in partial male-to-female sex reversal and loss of both caused complete sex reversal, which could be rescued by a nonsteroidal aromatase inhibitor. Compensatory expression of gsdf-A and gsdf-B was observed in gsdf-B and gsdf-A mutants, respectively. Subsequently, we determined that in tissue culture cells, Gsdf-A and Gsdf-B both interacted with Ncoa5 (nuclear receptor coactivator 5) and blocked Ncoa5 interaction with Rora (retinoic acid-related orphan receptor-alpha) to repress Rora/Ncoa5-induced activation of cyp19a1a (cytochrome P450, family 19, subfamily A, polypeptide 1a). These findings illustrate that Gsdf-A and Gsdf-B can regulate male differentiation by inhibiting cyp19a1a transcription in hexaploid gibel carp and also reveal that Gsdf-A and Gsdf-B can interact with Ncoa5 to suppress cyp19a1a transcription in vitro. This study provides a typical case of cooperative mechanism of duplicated genes in polyploids and also sheds light on the conserved evolution of sex differentiation.

Panx1 channels promote both anti- and pro-seizure-like activities in the zebrafish via p2rx7 receptors and ATP signaling

The molecular mechanisms of excitation/inhibition imbalances promoting seizure generation in epilepsy patients are not fully understood. Evidence suggests that Pannexin1 (Panx1), an ATP release channel, modulates the excitability of the brain. In this report, we performed electrophysiological, behavioral, and molecular phenotyping experiments on zebrafish larvae bearing genetic or pharmacological knockouts of Panx1a and Panx1b channels, each homologous to human PANX1. When Panx1a function is lost, or both channels are under pharmacological blockade, seizures with ictal-like events and seizure-like locomotion are reduced in the presence of pentylenetetrazol. Transcriptome profiling by RNA-seq demonstrates a spectrum of distinct metabolic and cell signaling states which correlate with the loss of Panx1a. Furthermore, the pro- and anticonvulsant activities of both Panx1 channels affect ATP release and involve the purinergic receptor P2rx7. Our findings suggest a subfunctionalization of Panx1 enabling dual roles in seizures, providing a unique and comprehensive perspective to understanding seizure mechanisms in the context of this channel.

Diversification and Functional Evolution of HOX Proteins

Gene duplication and divergence is a major contributor to the generation of morphological diversity and the emergence of novel features in vertebrates during evolution. The availability of sequenced genomes has facilitated our understanding of the evolution of genes and regulatory elements. However, progress in understanding conservation and divergence in the function of proteins has been slow and mainly assessed by comparing protein sequences in combination with in vitro analyses. These approaches help to classify proteins into different families and sub-families, such as distinct types of transcription factors, but how protein function varies within a gene family is less well understood. Some studies have explored the functional evolution of closely related proteins and important insights have begun to emerge. In this review, we will provide a general overview of gene duplication and functional divergence and then focus on the functional evolution of HOX proteins to illustrate evolutionary changes underlying diversification and their role in animal evolution.

The first genome for the Cape Primrose Streptocarpus rexii (Gesneriaceae), a model plant for studying meristem-driven shoot diversity

Cape Primroses (Streptocarpus, Gesneriaceae) are an ideal study system for investigating the genetics underlying species diversity in angiosperms. Streptocarpus rexii has served as a model species for plant developmental research for over five decades due to its unusual extended meristem activity present in the leaves. In this study, we sequenced and assembled the complete nuclear, chloroplast, and mitochondrial genomes of S. rexii using Oxford Nanopore Technologies long read sequencing. Two flow cells of PromethION sequencing resulted in 32 billion reads and were sufficient to generate a draft assembly including the chloroplast, mitochondrial and nuclear genomes, spanning 776 Mbp. The final nuclear genome assembly contained 5,855 contigs, spanning 766 Mbp of the 929-Mbp haploid genome with an N50 of 3.7 Mbp and an L50 of 57 contigs. Over 70% of the draft genome was identified as repeats. A genome repeat library of Gesneriaceae was generated and used for genome annotation, with a total of 45,045 genes annotated in the S. rexii genome. Ks plots of the paranomes suggested a recent whole genome duplication event, shared between S. rexii and Primulina huaijiensis. A new chloroplast and mitochondrial genome assembly method, based on contig coverage and identification, was developed, and successfully used to assemble both organellar genomes of S. rexii. This method was developed into a pipeline and proved widely applicable. The nuclear genome of S. rexii and other datasets generated and reported here will be invaluable resources for further research to aid in the identification of genes involved in morphological variation underpinning plant diversification.

Divergent Evolution of Progesterone and Mineralocorticoid Receptors in Terrestrial Vertebrates and Fish Influences Endocrine Disruption

There is much concern about disruption of endocrine physiology regulated by steroid hormones in humans, other terrestrial vertebrates and fish by industrial chemicals, such as bisphenol A, and pesticides, such as DDT. These endocrine-disrupting chemicals influence steroid-mediated physiology in humans and other vertebrates by competing with steroids for receptor binding sites, disrupting diverse responses involved in reproduction, development and differentiation. Here I discuss that due to evolution of the progesterone receptor (PR) and mineralocorticoid receptor (MR) after ray-finned fish and terrestrial vertebrates diverged from a common ancestor, each receptor evolved to respond to different steroids in ray-finned fish and terrestrial vertebrates. In elephant shark, a cartilaginous fish that diverged before the separation between ray-finned fish and terrestrial vertebrates, both progesterone and 17,20β-dihydroxy-progesterone activate the PR. During the evolution of ray-finned fish and terrestrial vertebrates, the PR in terrestrial vertebrates continued responding to progesterone and evolved to weakly respond to 17,20β-dihydroxy-progesterone. In contrast, the physiological progestin for the PR in zebrafish and other ray-finned fish is 17,20β-dihydroxy-progesterone, and ray-finned fish PR responds weakly to progesterone. The MR in fish and terrestrial vertebrates also diverged to have different responses to progesterone. Progesterone is a potent agonist for elephant shark MR, zebrafish MR and other fish MRs, in contrast to progesterone's opposite activity as an antagonist for aldosterone, the physiological mineralocorticoid for human MR. These different physiological ligands for fish and terrestrial vertebrate PR and MR need to be considered in applying data for their disruption by chemicals in fish and terrestrial vertebrates to each other.

Genome-wide characterization of gap junction (connexins and pannexins) genes in turbot (Scophthalmus maximus L.): evolution and immune response following Vibrio anguillarum infection

Gap junction (GJ), a special intercellular junction between different cell types, directly connects the cytoplasm of adjacent cells, allows various molecules, ions and electrical impulses to pass through the intercellular regulatory gate, and plays vital roles in response to bacterial infection. Up to date, the information about the GJ in turbot (Scophthalmus maximus L.) is still limited. In current study, 43 gap junction genes were identified in turbot, phylogeny analysis suggested that gap junctions from turbot and other species were clustered into six groups, GJA, GJB, GJC, GJD, GJE and PANX, and turbot GJs together with respective GJs from Japanese flounder, half-smooth tongue sole and large yellow croaker, sharing same ancestors. In addition, these 43 GJ genes distributed in different chromosomes unevenly. According to gene structure and domain analysis, these genes (in GJA-GJE group) were highly conserved in that most of them contain the transmembrane area, connexin domain (CNX) and cysteine-rich domain (connexin CCC), while PANXs contain Pfam Innexin. Although only one tandem duplication was identified in turbot gap junction gene, 235 pairs of segmental duplications were identified in the turbot genome. To further investigate their evolutionary relationships, Ka/Ks was calculated, and results showed that most ratios were lower than 1, indicating they had undergone negative selection. Finally, expression analysis showed that gap junction genes were widely distributed in turbot tissues and significantly regulated after Vibrio anguillarum infection. Taken together, our research could provide valuable information for further exploration of the function of gap junction genes in teleost.

Genome-wide identification and transcriptome-based expression profiling of the Sox gene family in the spinyhead croaker (Collichthys lucidus)

Sox genes encode transcription factors with a high-mobility group (HMG) box, playing critical roles in the initiation and maintenance of a variety of developmental processes, such as sex determination and differentiation. In the present study, we identified 26 Sox genes in the genome of spinyhead croaker Collichthys lucidus (Richardson, 1844) with homology-based analysis of the HMG box. The transcriptome-based expression profiles revealed that the expression of the Sox gene in gonads began to differ between sexes when the body length was 2.74 ± 0.24 cm. At that time, three Sox genes (Sox11b, Sox8a and Sox19) were significantly upregulated, accompanied by the downregulation of 12 Sox genes in the ovary, and six Sox genes were temporarily significantly upregulated in the testis. Afterwards, the expression profile of Sox genes changed only with a small amplitude in both the ovary and testis. For adult tissues, huge differences were observed in the expression profiles of Sox genes between ovaries and testes, as well as small differences in somatic tissues between sexes. These results provide clues to further decipher the role of Sox genes in the processes of sex determination and differentiation in spinyhead croaker and other teleosts.

Fish Paralog Proteins RNASEK-a and -b Enhance Type I Interferon Secretion and Promote Apoptosis

Type I interferon and apoptosis elicit multifaceted effects on host defense and various diseases, such as viral infections and cancers. However, the gene/protein network regulating type I interferon and apoptosis has not been elucidated completely. In this study, we selected grass carp (Ctenopharyngodon idella) as an experimental model to investigate the modulation of RNASEK on the secretion of type I interferon and apoptosis. We first cloned two paralogs RNASEK-a and -b in grass carp, defined three exons in each gene, and found the length of both coding regions is 306 bp with 73.27% of protein homology. The protein sequences of the two paralogs are highly conserved across species. Two proteins were mainly localized in early and late endosomes and endoplasmic reticulum. Further, quantitative real-time PCR demonstrated that dsRNA poly I:C and grass carp reovirus upregulated RNASEK-a and -b in grass carp cells and tissues. Overexpression of RNASEK-a and -b individually induced type I interferon expression and the phosphorylation of IRF3/IRF7 shown by Western blot and immunofluorescent staining, increased Bax/Bcl-2 mRNA ratio, DNA fragmentations, TUNEL-positive cells, and the proportion of Annexin V-positive signals in flow cytometry, and activated eIF2α, opposite to that observed when RNASEK-a and -b were knocked down in multiple cell types. Taken together, we claim for the first time that fish paralog proteins RNASEK-a and -b enhance type I interferon secretion and promote apoptosis, which may be involved in the phosphorylation of IRF3/IRF7 and eIF2α, respectively. Our study reveals a previously unrecognized role of RNASEK as a new positive regulator of type I interferon and apoptosis.

A supernumerary "B-sex" chromosome drives male sex determination in the Pachón cavefish, Astyanax mexicanus

Sex chromosomes are generally derived from a pair of classical type-A chromosomes, and relatively few alternative models have been proposed up to now.^1,2 B chromosomes (Bs) are supernumerary and dispensable chromosomes with non-Mendelian inheritance found in many plant and animal species^3,4 that have often been considered as selfish genetic elements that behave as genome parasites.^5,6 The observation that in some species Bs can be either restricted or predominant in one sex^7-14 raised the interesting hypothesis that Bs could play a role in sex determination.¹⁵ The characterization of putative B master sex-determining (MSD) genes, however, has not yet been provided to support this hypothesis. Here, in Astyanax mexicanus cavefish originating from Pachón cave, we show that Bs are strongly male predominant. Based on a high-quality genome assembly of a B-carrying male, we characterized the Pachón cavefish B sequence and found that it contains two duplicated loci of the putative MSD gene growth differentiation factor 6b (gdf6b). Supporting its role as an MSD gene, we found that the Pachón cavefish gdf6b gene is expressed specifically in differentiating male gonads, and that its knockout induces male-to-female sex reversal in B-carrying males. This demonstrates that gdf6b is necessary for triggering male sex determination in Pachón cavefish. Altogether these results bring multiple and independent lines of evidence supporting the conclusion that the Pachón cavefish B is a "B-sex" chromosome that contains duplicated copies of the gdf6b gene, which can promote male sex determination in this species.

Chromosome-Level Genome Assembly of Chinese Sucker (Myxocyprinus asiaticus) Reveals Strongly Conserved Synteny Following a Catostomid-Specific Whole-Genome Duplication

Fishes of the family Catostomidae (“suckers”; Teleostei: Cypriniformes) are hypothesized to have undergone an allopolyploidy event approximately 60 Ma. However, genomic evidence has previously been unavailable to assess this hypothesis. We sequenced and assembled the first chromosome-level catostomid genome, Chinese sucker (Myxocyprinus asiaticus), and present clear evidence of a catostomid-specific whole-genome duplication (WGD) event (“Cat-4R”). Our results reveal remarkably strong, conserved synteny since this duplication event, as well as between Myxocyprinus and an unduplicated outgroup, zebrafish (Danio rerio). Gene content and repetitive elements are also approximately evenly distributed across homeologous chromosomes, suggesting that both subgenomes retain some function, with no obvious bias in gene fractionation or subgenome dominance. The Cat-4R duplication provides another independent example of genome evolution following WGD in animals, in this case at the extreme end of conserved genome architecture over at least 25.2 Myr since the duplication. The M. asiaticus genome is a useful resource for researchers interested in understanding genome evolution following WGD in animals.

Identification and evolution of transcription factors RHR gene family (NFAT and RBPJ) involving lamprey (Lethenteron reissneri) innate immunity

Nuclear factor of activated T cells (NFAT) and recombination signal binding protein (RBP) belong to the family of Rel homology region (RHR) transcription factors which regulate the expression of genes involved in different aspects of the immune response. To gain insights into the evolution and characterisation of RHR genes in lampreys, a jawless vertebrate, four RHR genes, including nuclear factor of activated T cells (NFAT) and recombination signal binding protein for immunoglobulin kappa J region (RBPJ), have been identified and cloned from the lamprey (Lethenteron reissneri) database. Evolutionary relationships of NFAT and RBPJ genes among different species were determined through molecular phylogenetic analysis. Motif, genetic structure, and tertiary structure analyses showed that NFATs and RBPJ are conserved and contain RHD and IPT domains. Moreover, synteny analysis showed that the neighbourhood genes of Lr-NFATs and Lr-RBPJ have undergone significant changes compared to jawed vertebrates. Real-time quantitative results demonstrated that the RHR gene family plays a significant role in immune defence. This study provides a new understanding of the origin and evolution of the RHR gene family in different species.

Molecular Evolution of clock Genes in Vertebrates

Circadian rhythms not only influence the overall daily routine of organisms but also directly affect life activities to varying degrees. Circadian locomotor output cycle kaput (Clock), the most critical gene in the circadian rhythm feedback system, plays an important role in the regulation of biological rhythms. Here, we aimed to elucidate the evolutionary history of the clock gene family in a taxonomically diverse set of vertebrates, providing novel insights into the evolution of the clock gene family based on 102 vertebrate genomes. Using genome-wide analysis, we extracted 264 clock sequences. In lobe-finned fishes and some basal non-teleost ray-finned fishes, only two clock isotypes were found (clock1 and clock2). However, the majority of teleosts possess three clock genes (two clock1 genes and one clock2 gene) owing to extra whole-genome duplication. The following syntenic analysis confirmed that clock1a, clock1b, and clock2 are conserved in teleost species. Interestingly, we discovered that osteoglossomorph fishes possess two clock2 genes. Moreover, protein sequence comparisons indicate that CLOCK protein changes among vertebrates were concentrated at the N-terminal and poly Q regions. We also performed a dN/dS analysis, and the results suggest that clock1 and clock2 may show distinct fates for duplicated genes between the lobe-finned and ray-finned fish clades. Collectively, these results provide a genome-wide insight into clock gene evolution in vertebrates.

Comprehensive phylogeny of Konosirus punctatus (Clupeiformes: Clupeidae) based on transcriptomic data

Konosirus punctatus is an important species for the structure of marine ecosystems. Meanwhile, it is a native species in the northwest Pacific Ocean and supports important commercial fishery. In the present study, we generated the whole transcriptome of K. punctatus from combined tissues (muscle, liver, gill, heart, kidney, swim bladder and sexual gland) using Illumina RNA-seq technology and a total of 46087110 clean reads were obtained, corresponding to 6531521430 nucleotides. Meanwhile, 10000 clean reads were randomly selected and compared with NT database to examine the possible contamination. Results showed that 6754 clean reads were distributed among some species closely related with K. punctatus, indicating no-pollution. De novo assembly was performed and all clean reads were assembled to produce 71610 longest unigenes with an N50 of 906 bp. Among all the unigenes, 43974 unigenes were annotated in at least one database and 3172 unigenes were annotated in all databases. All unigenes were further analyzed to predict the gene structure and we have obtained a total of 54864 coding sequences (CDSs) and 17326 simple sequence repeats (SSRs). Saturation analyses were applied to evaluate the accuracy of gene expression and we hypothesized that the detection of gene expression might be effective. Finally, single-copy orthologous genes were applied to construct the phylogenetic relationship of K. punctatus. Results showed that K. punctatus diverged from the common ancestor with Alosa alosa, Alosa pseudoharengus and Sardina pilchardus at approx. 61.16–92.52 MYA. The present study will provide a foundational molecular information for the biological research of K. punctatus.

Genetic Drift Dominates Genome-Wide Regulatory Evolution Following an Ancient Whole-Genome Duplication in Atlantic Salmon

Whole-genome duplications (WGD) have been considered as springboards that potentiate lineage diversification through increasing functional redundancy. Divergence in gene regulatory elements is a central mechanism for evolutionary diversification, yet the patterns and processes governing regulatory divergence following events that lead to massive functional redundancy, such as WGD, remain largely unknown. We studied the patterns of divergence and strength of natural selection on regulatory elements in the Atlantic salmon (Salmo salar) genome, which has undergone WGD 100-80 Ma. Using ChIPmentation, we first show that H3K27ac, a histone modification typical to enhancers and promoters, is associated with genic regions, tissue-specific transcription factor binding motifs, and with gene transcription levels in immature testes. Divergence in transcription between duplicated genes from WGD (ohnologs) correlated with difference in the number of proximal regulatory elements, but not with promoter elements, suggesting that functional divergence between ohnologs after WGD is mainly driven by enhancers. By comparing H3K27ac regions between duplicated genome blocks, we further show that a longer polyploid state post-WGD has constrained regulatory divergence. Patterns of genetic diversity across natural populations inferred from resequencing indicate that recent evolutionary pressures on H3K27ac regions are dominated by largely neutral evolution. In sum, our results suggest that post-WGD functional redundancy in regulatory elements continues to have an impact on the evolution of the salmon genome, promoting largely neutral evolution of regulatory elements despite their association with transcription levels. These results highlight a case where genome-wide regulatory evolution following an ancient WGD is dominated by genetic drift.

Growth genes are implicated in the evolutionary divergence of sympatric piscivorous and insectivorous rainbow trout (Oncorhynchus mykiss)

Background

Identifying ecologically significant phenotypic traits and the genomic mechanisms that underly them are crucial steps in understanding traits associated with population divergence. We used genome-wide data to identify genomic regions associated with key traits that distinguish two ecomorphs of rainbow trout (Oncorhynchus mykiss)—insectivores and piscivores—that coexist for the non-breeding portion of the year in Kootenay Lake, southeastern British Columbia. “Gerrards” are large-bodied, rapidly growing piscivores with high metabolic rates that spawn north of Kootenay Lake in the Lardeau River, in contrast to the insectivorous populations that are on average smaller in body size, with lower growth and metabolic rates, mainly forage on aquatic insects, and spawn in tributaries immediately surrounding Kootenay Lake. We used pool-seq data representing ~ 60% of the genome and 80 fish per population to assess the level of genomic divergence between ecomorphs and to identify and interrogate loci that may play functional or selective roles in their divergence.

Results

Genomic divergence was high between sympatric insectivores and piscivores (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{\text{ST}}$$\end{document}FST = 0.188), and in fact higher than between insectivorous populations from Kootenay Lake and the Blackwater River (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{\text{ST}}$$\end{document}FST = 0.159) that are > 500 km apart. A window-based \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{\text{ST}}$$\end{document}FST analysis did not reveal “islands” of genomic differentiation; however, the window with highest \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{\text{ST}}$$\end{document}FST estimate did include a gene associated with insulin secretion. Although we explored the use of the “Local score” approach to identify genomic outlier regions, this method was ultimately not used because simulations revealed a high false discovery rate (~ 20%). Gene ontology (GO) analysis identified several growth processes as enriched in genes occurring in the ~ 200 most divergent genomic windows, indicating many loci of small effect involved in growth and growth-related metabolic processes are associated with the divergence of these ecomorphs.

Conclusion

Our results reveal a high degree of genomic differentiation between piscivorous and insectivorous populations and indicate that the large body piscivorous phenotype is likely not due to one or a few loci of large effect. Rather, the piscivore phenotype may be controlled by several loci of small effect, thus highlighting the power of whole-genome resequencing in identifying genomic regions underlying population-level phenotypic divergences.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01795-9.

The American Paddlefish Genome Provides Novel Insights into Chromosomal Evolution and Bone Mineralization in Early Vertebrates

Sturgeons and paddlefishes (Acipenseriformes) occupy the basal position of ray-finned fishes, although they have cartilaginous skeletons as in Chondrichthyes. This evolutionary status and their morphological specializations make them a research focus, but their complex genomes (polyploidy and the presence of microchromosomes) bring obstacles and challenges to molecular studies. Here, we generated the first high-quality genome assembly of the American paddlefish (Polyodon spathula) at a chromosome level. Comparative genomic analyses revealed a recent species-specific whole-genome duplication event, and extensive chromosomal changes, including head-to-head fusions of pairs of intact, large ancestral chromosomes within the paddlefish. We also provide an overview of the paddlefish SCPP (secretory calcium-binding phosphoprotein) repertoire that is responsible for tissue mineralization, demonstrating that the earliest flourishing of SCPP members occurred at least before the split between Acipenseriformes and teleosts. In summary, this genome assembly provides a genetic resource for understanding chromosomal evolution in polyploid nonteleost fishes and bone mineralization in early vertebrates.

Functional and Evolutionary Diversification of Otx2 and Crx in Vertebrate Retinal Photoreceptor and Bipolar Cell Development

Otx family homeoproteins Otx2 and Crx are expressed in photoreceptor precursor cells and bind to the common DNA-binding consensus sequence, but these two proteins have distinct functions in retinal development. To examine the functional substitutability of Otx2 and Crx, we generate knockin mouse lines in which Crx is replaced by Otx2 and vice versa. We find that Otx2 and Crx cannot be substituted in photoreceptor development. Subsequently, we investigate the function of Otx2 in photoreceptor and bipolar cell development. High Otx2 levels induce photoreceptor cell fate but not bipolar cell fate, whereas reduced Otx2 expression impairs bipolar cell maturation and survival. Furthermore, we identify Otx2 and Crx in the lamprey genome by using synteny analysis, suggesting that the last common ancestor of vertebrates possesses both Otx2 and Crx. We find that the retinal Otx2 expression pattern is different between lampreys and mice, suggesting that neofunctionalization of Otx2 occurred in the jawed vertebrate lineage.

Biochemistry and physiology of zebrafish photoreceptors

All vertebrates share a canonical retina with light-sensitive photoreceptors in the outer retina. These photoreceptors are of two kinds: rods and cones, adapted to low and bright light conditions, respectively. They both show a peculiar morphology, with long outer segments, comprised of ordered stacks of disc-shaped membranes. These discs host numerous proteins, many of which contribute to the visual transduction cascade. This pathway converts the light stimulus into a biological signal, ultimately modulating synaptic transmission. Recently, the zebrafish (Danio rerio) has gained popularity for studying the function of vertebrate photoreceptors. In this review, we introduce this model system and its contribution to our understanding of photoreception with a focus on the cone visual transduction cascade.

Connexins during 500 Million Years-From Cyclostomes to Mammals

It was previously shown that the connexin gene family had relatively similar subfamily structures in several vertebrate groups. Still, many details were left unclear. There are essentially no data between tunicates, which have connexins that cannot be divided into the classic subfamilies, and teleosts, where the subfamilies are easily recognized. There are also relatively few data for the groups that diverged between the teleosts and mammals. As many of the previously analyzed genomes have been improved, and many more genomes are available, we reanalyzed the connexin gene family and included species from all major vertebrate groups. The major results can be summarized as follows: (i) The same connexin subfamily structures are found in all Gnathostomata (jawed vertebrates), with some variations due to genome duplications, gene duplications and gene losses. (ii) In contrast to previous findings, birds do not have a lower number of connexins than other tetrapods. (iii) The cyclostomes (lampreys and hagfishes) possess genes in the alpha, beta, gamma and delta subfamilies, but only some of the genes show a phylogenetic affinity to specific genes in jawed vertebrates. Thus, two major evolutionary transformations have occurred in this gene family, from tunicates to cyclostomes and from cyclostomes to jawed vertebrates.

Transposable Elements and Teleost Migratory Behaviour

Transposable elements (TEs) represent a considerable fraction of eukaryotic genomes, thereby contributing to genome size, chromosomal rearrangements, and to the generation of new coding genes or regulatory elements. An increasing number of works have reported a link between the genomic abundance of TEs and the adaptation to specific environmental conditions. Diadromy represents a fascinating feature of fish, protagonists of migratory routes between marine and freshwater for reproduction. In this work, we investigated the genomes of 24 fish species, including 15 teleosts with a migratory behaviour. The expected higher relative abundance of DNA transposons in ray-finned fish compared with the other fish groups was not confirmed by the analysis of the dataset considered. The relative contribution of different TE types in migratory ray-finned species did not show clear differences between oceanodromous and potamodromous fish. On the contrary, a remarkable relationship between migratory behaviour and the quantitative difference reported for short interspersed nuclear (retro)elements (SINEs) emerged from the comparison between anadromous and catadromous species, independently from their phylogenetic position. This aspect is likely due to the substantial environmental changes faced by diadromous species during their migratory routes.

The Evolution of Oxytocin and Vasotocin Receptor Genes in Jawed Vertebrates: A Clear Case for Gene Duplications Through Ancestral Whole-Genome Duplications

The neuronal and neuroendocrine peptides oxytocin (OT) and vasotocin (VT), including vasopressins, have six cognate receptors encoded by six receptor subtype genes in jawed vertebrates. The peptides elicit a broad range of responses that are specifically mediated by the receptor subtypes including neuronal functions regulating behavior and hormonal actions on reproduction and water/electrolyte balance. Previously, we have demonstrated that these six receptor subtype genes, which we designated VTR1A, VTR1B, OTR, VTR2A, VTR2B and VTR2C, arose from a syntenic ancestral gene pair, one VTR1/OTR ancestor and one VTR2 ancestor, through the early vertebrate whole-genome duplications (WGD) called 1R and 2R. This was supported by both phylogenetic and chromosomal conserved synteny data. More recently, other studies have focused on confounding factors, such as the OTR/VTR orthologs in cyclostomes, to question this scenario for the origin of the OTR/VTR gene family; proposing instead less parsimonious interpretations involving only one WGD followed by complex series of chromosomal or segmental duplications. Here, we have updated the phylogeny of the OTR/VTR gene family, including a larger number of vertebrate species, and revisited seven representative neighboring gene families from our previous conserved synteny analyses, adding chromosomal information from newer high-coverage genome assemblies from species that occupy key phylogenetic positions: the polypteriform fish reedfish (Erpetoichthys calabaricus), the cartilaginous fish thorny skate (Amblyraja radiata) and a more recent high-quality assembly of the Western clawed frog (Xenopus tropicalis) genome. Our analyses once again add strong support for four-fold symmetry, i.e., chromosome quadruplication in the same time window as the WGD events early in vertebrate evolution, prior to the jawed vertebrate radiation. Thus, the evolution of the OTR/VTR gene family can be most parsimoniously explained by two WGD events giving rise to the six ancestral genes, followed by differential gene losses of VTR2 genes in different lineages. We also argue for more coherence and clarity in the nomenclature of OT/VT receptors, based on the most parsimonious scenario.

The stress - Reproductive axis in fish: The involvement of functional neuroanatomical systems in the brain

The neuroendocrine-stress axis of nonmammalian species is evolutionarily conserved, which makes them useful to serve as important model systems for elucidating the function of the vertebrate stress response. The involvement of hypothalamo-pituitary-adrenal (HPA) axis hormones in regulation of stress and reproduction is well described in different vertebrates. However, the stress response is a complex process, which appears to be controlled by a number of neurochemicals in association with hypothalamo-pituitary-interrenal (HPI) axis or independent of HPI axis in fish. In recent years, the participation of neurohormones other than HPI axis in regulation of stress and reproduction is gaining more attention. This review mainly focuses on the involvement of functional neuroanatomical systems such as the catecholaminergic neurotransmitter dopamine (DA) and opioid peptides in regulation of the stress-reproductive axis in fish. Occurrences of DA and opioid peptides like β-endorphin, enkephalins, dynorphin, and endomorphins have been demonstrated in fish brain, and diverse roles such as pain modulation, social behaviour and reproduction are implicated for these hormones. Neuroanatomical studies using retrograde tracing, immunohistochemical staining and lesion methods have demonstrated that the neurons originating in the preoptic region and the nucleus lateralis tuberis directly innervate the pituitary gland and, therefore, the hypophysiotrophic role of these hormones. In addition, heightened synthetic and secretory activity of the opioidergic and the dopaminergic neurons in hypothalamic areas of the brain during stress exposure suggest potentially intricate relationship with the stress-reproductive axis in fish. Current evidence in early vertebrates like fish provides a novel insight into the underlying neuroendocrine mechanisms as additional pathways along the stress-reproductive axis that seem to be conserved during the course of evolution.

Uncovering patterns of the evolution of genomic sequence entropy and complexity

The lack of consensus concerning the biological meaning of entropy and complexity of genomes and the different ways to assess these data hamper conclusions concerning what are the causes of genomic entropy variation among species. This study aims to evaluate the entropy and complexity of genomic sequences of several species without using homologies to assess relationships among these variables and non-molecular data (e.g., the number of individuals) to seek a trigger of interspecific genomic entropy variation. The results indicate a relationship among genomic entropy, genome size, genomic complexity, and the number of individuals: species with a small number of individuals harbors large genome, and hence, low entropy but a higher complexity. We defined that the complexity of a genome relies on the entropy of each DNA segment within genome. Then, the entropy and complexity of a genome reflects its organization solely. Exons of vertebrates harbor smaller entropies than non-exon regions (likely by the repeats that accumulated from duplications), whereas other taxonomic groups do not present this pattern. Our findings suggest that small initial population might have defined current genomic entropy and complexity: actual genomes are less complex than ancestral ones. Besides, our data disagree with the relationship between phenotype and genomic entropies previously established. Finally, by establishing the relationship between genomic entropy/complexity with the number of individuals and genome size, under an evolutive perspective, ideas concerning the genomic variability may emerge.

Exon probe sets and bioinformatics pipelines for all levels of fish phylogenomics

Exon markers have a long history of use in phylogenetics of ray-finned fishes, the most diverse clade of vertebrates with more than 35,000 species. As the number of published genomes increases, it has become easier to test exons and other genetic markers for signals of ancient duplication events and filter out paralogues that can mislead phylogenetic analysis. We present seven new probe sets for current target-capture phylogenomic protocols that capture 1,104 exons explicitly filtered for paralogues using gene trees. These seven probe sets span the diversity of teleost fishes, including four sets that target five hyperdiverse percomorph clades which together comprise ca. 17,000 species (Carangaria, Ovalentaria, Eupercaria, and Syngnatharia + Pelagiaria combined). We additionally included probes to capture legacy nuclear exons and mitochondrial markers that have been commonly used in fish phylogenetics (despite some exons being flagged for paralogues) to facilitate integration of old and new molecular phylogenetic matrices. We tested these probes experimentally for 56 fish species (eight species per probe set) and merged new exon-capture sequence data into an existing data matrix of 1,104 exons and 300 ray-finned fish species. We provide an optimized bioinformatics pipeline to assemble exon capture data from raw reads to alignments for downstream analysis. We show that legacy loci with known paralogues are at risk of assembling duplicated sequences with target-capture, but we also assembled many useful orthologous sequences that can be integrated with many PCR-generated matrices. These probe sets are a valuable resource for advancing fish phylogenomics because targeted exons can easily be extracted from increasingly available whole genome and transcriptome data sets, and also may be integrated with existing PCR-based exon and mitochondrial data.

Mobile Elements in Ray-Finned Fish Genomes

Ray-finned fishes (Actinopterygii) are a very diverse group of vertebrates, encompassing species adapted to live in freshwater and marine environments, from the deep sea to high mountain streams. Genome sequencing offers a genetic resource for investigating the molecular bases of this phenotypic diversity and these adaptations to various habitats. The wide range of genome sizes observed in fishes is due to the role of transposable elements (TEs), which are powerful drivers of species diversity. Analyses performed to date provide evidence that class II DNA transposons are the most abundant component in most fish genomes and that compared to other vertebrate genomes, many TE superfamilies are present in actinopterygians. Moreover, specific TEs have been reported in ray-finned fishes as a possible result of an intricate relationship between TE evolution and the environment. The data summarized here underline the biological interest in Actinopterygii as a model group to investigate the mechanisms responsible for the high biodiversity observed in this taxon.

Collagen, stiffness, and adhesion: the evolutionary basis of vertebrate mechanobiology

The emergence of collagen I in vertebrates resulted in a dramatic increase in the stiffness of the extracellular environment, supporting long-range force propagation and the development of low-compliant tissues necessary for the development of vertebrate traits including pressurized circulation and renal filtration. Vertebrates have also evolved integrins that can bind to collagens, resulting in the generation of higher tension and more efficient force transmission in the extracellular matrix. The stiffer environment provides an opportunity for the vertebrates to create new structures such as the stress fibers, new cell types such as endothelial cells, new developmental processes such as neural crest delamination, and new tissue organizations such as the blood-brain barrier. Molecular players found only in vertebrates allow the modification of conserved mechanisms as well as the design of novel strategies that can better serve the physiological needs of the vertebrates. These innovations collectively contribute to novel morphogenetic behaviors and unprecedented increases in the complexities of tissue mechanics and functions.

Cloning, molecular characterization, and nutritional regulation of fatty acid-binding protein family genes in gold pompanos (Trachinotus ovatus)

Fatty acid-binding protein (Fabp) is an important protein family involved in fatty acid uptake and deposition. Elucidating the function and regulation of fabps could contribute to the efficient production of biologically relevant fatty acids, such as highly unsaturated fatty acids (HUFAs), from fish. Herein, five genes from Trachinotus ovatus named fabp4, fabp6a, fabp6b, fabp7a, and fabp7b coding 133, 127, 118, 132, and 132 amino acid residues were cloned and sequenced. The effect of dietary HUFA on the expression of these genes was also investigated. Multiple protein sequence alignment showed that these Fabps shared high identity to their orthologs from other fish and mammals. Two conserved domains, lipocalin and lipocalin 7, were predicted in the deduced protein sequence of fabp4 and fabp7 paralogs, whereas fabp6 paralogs did not present the lipocalin domain. The adipose tissue, spleen, gill, and intestine showed the highest levels of fabp6b expression. In the brain, fabp6b was weakly expressed, whereas the expression of fabp7a was at its highest. Conversely, fabp7a showed a lower mRNA level than the other fabps in the liver and heart. In the dorsal muscle and kidney, fabp6a was the most abundantly expressed gene. Increasing dietary HUFA from 1.0% to 2.1% increased the gene expression of hepatic fabp4 and fabp6a gene expression but decreased gene expression in the dorsal muscle. Similarly, the expression of fabp7a in the dorsal muscle also declined in the 2.1% HUFA group. This study lays the groundwork for further studies focused on the physiological function and regulation of fish fabps.

Regain of sex determination system and sexual reproduction ability in a synthetic octoploid male fish

Polyploids in vertebrates are generally associated with unisexual reproduction, but the direct consequences of polyploidy on sex determination system and reproduction mode remain unknown. Here, we synthesized a group of artificial octoploids between unisexual gynogenetic hexaploid Carassius gibelio and sexual tetraploid Carassius auratus. The synthetic octoploids were revealed to have more than 200 chromosomes, in which 50 chromosomes including the X/Y sex determination system were identified to transfer from sexual tetraploid C. auratus into the unisexual gynogenetic hexaploid C. gibelio. Significantly, a few synthetic octoploid males were found to be fertile, and one octoploid male was confirmed to regain sexual reproduction ability, which exhibits characteristics that are the same to sexual reproduction tetraploid males, such as 1:1 sex ratio occurrence, meiosis completion and euploid sperm formation in spermatogenesis, as well as normal embryo development and gene expression pattern during embryogenesis. Therefore, the current finding provides a unique case to explore the effect of sex determination system incorporation on reproduction mode transition from unisexual gynogenesis to sexual reproduction along with genome synthesis of recurrent polyploidy in vertebrates.

The lasting after-effects of an ancient polyploidy on the genomes of teleosts

The ancestor of most teleost fishes underwent a whole-genome duplication event three hundred million years ago. Despite its antiquity, the effects of this event are evident both in the structure of teleost genomes and in how the surviving duplicated genes still operate to drive form and function. I inferred a set of shared syntenic regions that survive from the teleost genome duplication (TGD) using eight teleost genomes and the outgroup gar genome (which lacks the TGD). I then phylogenetically modeled the TGD's resolution via shared and independent gene losses and applied a new simulation-based statistical test for the presence of bias toward the preservation of genes from one parental subgenome. On the basis of that test, I argue that the TGD was likely an allopolyploidy. I find that duplicate genes surviving from this duplication in zebrafish are less likely to function in early embryo development than are genes that have returned to single copy at some point in this species' history. The tissues these ohnologs are expressed in, as well as their biological functions, lend support to recent suggestions that the TGD was the source of a morphological innovation in the structure of the teleost retina. Surviving duplicates also appear less likely to be essential than singletons, despite the fact that their single-copy orthologs in mouse are no less essential than other genes.

Comprehensive Transcriptome Analysis Reveals Insights into Phylogeny and Positively Selected Genes of Sillago Species

Sillago species lives in the demersal environments and face multiple stressors, such as localized oxygen depletion, sulfide accumulation, and high turbidity. In this study, we performed transcriptome analyses of seven Sillago species to provide insights into the phylogeny and positively selected genes of this species. After de novo assembly, 82,024, 58,102, 63,807, 85,990, 102,185, 69,748, and 102,903 unigenes were generated from S. japonica, S. aeolus, S. sp.1, S. sihama, S. sp.2, S. parvisquamis, and S. sinica, respectively. Furthermore, 140 shared orthologous exon markers were identified and then applied to reconstruct the phylogenetic relationships of the seven Sillago species. The reconstructed phylogenetic structure was significantly congruent with the prevailing morphological and molecular biological view of Sillago species relationships. In addition, a total of 44 genes were identified to be positively selected, and these genes were potential participants in the stress response, material (carbohydrate, amino acid and lipid) and energy metabolism, growth and differentiation, embryogenesis, visual sense, and other biological processes. We suspected that these genes possibly allowed Sillago species to increase their ecological adaptation to multiple environmental stressors.

Origin and evolution of the Rax homeobox gene by comprehensive evolutionary analysis

Rax is one of the key transcription factors crucial for vertebrate eye development. In this study, we conducted comprehensive evolutionary analysis of Rax. We found that Bilateria and Cnidaria possess Rax, but Placozoa, Porifera, and Ctenophora do not, implying that the origin of the Rax gene dates back to the common ancestor of Cnidaria and Bilateria. The results of molecular phylogenetic and synteny analyses on Rax loci between jawed and jawless vertebrates indicate that segmental duplication of the Rax locus occurred in an early common ancestor of jawed vertebrates, resulting in two Rax paralogs in jawed vertebrates, Rax and Rax2. By analyzing 86 mammalian genomes from all four major groups of mammals, we found that at least five independent Rax2 gene loss events occurred in mammals. This study may provide novel insights into the evolution of the eye.

Phylogeny of teleost connexins reveals highly inconsistent intra- and interspecies use of nomenclature and misassemblies in recent teleost chromosome assemblies

BACKGROUND

Based on an initial collecting of database sequences from the gap junction protein gene family (also called connexin genes) in a few teleosts, the naming of these sequences appeared variable. The reasons could be (i) that the structure in this family is variable across teleosts, or (ii) unfortunate naming. Rather clear rules for the naming of genes in fish and mammals have been outlined by nomenclature committees, including the naming of orthologous and ohnologous genes. We therefore analyzed the connexin gene family in teleosts in more detail. We covered the range of divergence times in teleosts (eel, Atlantic herring, zebrafish, Atlantic cod, three-spined stickleback, Japanese pufferfish and spotted pufferfish; listed from early divergence to late divergence).

RESULTS

The gene family pattern of connexin genes is similar across the analyzed teleosts. However, (i) several nomenclature systems are used, (ii) specific orthologous groups contain genes that are named differently in different species, (iii) several distinct genes have the same name in a species, and (iv) some genes have incorrect names. The latter includes a human connexin pseudogene, claimed as GJA4P, but which in reality is Cx39.2P (a delta subfamily gene often called GJD2like). We point out the ohnologous pairs of genes in teleosts, and we suggest a more consistent nomenclature following the outlined rules from the nomenclature committees. We further show that connexin sequences can indicate some errors in two high-quality chromosome assemblies that became available very recently.

CONCLUSIONS

Minimal consistency exists in the present practice of naming teleost connexin genes. A consistent and unified nomenclature would be an advantage for future automatic annotations and would make various types of subsequent genetic analyses easier. Additionally, roughly 5% of the connexin sequences point out misassemblies in the new high-quality chromosome assemblies from herring and cod.

Building a Body Shape Morphospace of Teleostean Fishes

We present a dataset that quantifies body shape in three dimensions across the teleost phylogeny. Built by a team of researchers measuring easy-to-identify, functionally relevant traits on specimens at the Smithsonian National Museum of Natural History it contains data on 16,609 specimens from 6144 species across 394 families. Using phylogenetic comparative methods to analyze the dataset we describe the teleostean body shape morphospace and identify families with extraordinary rates of morphological evolution. Using log shape ratios, our preferred method of body-size correction, revealed that fish width is the primary axis of morphological evolution across teleosts, describing a continuum from narrow-bodied laterally compressed flatfishes to wide-bodied dorsoventrally flattened anglerfishes. Elongation is the secondary axis of morphological variation and occurs within the more narrow-bodied forms. This result highlights the importance of collecting shape on three dimensions when working across teleosts. Our analyses also uncovered the fastest rates of shape evolution within a clade formed by notothenioids and scorpaeniforms, which primarily thrive in cold waters and/or have benthic habits, along with freshwater elephantfishes, which as their name suggests, have a novel head and body shape. This unprecedented dataset of teleostean body shapes will enable the investigation of the factors that regulate shape diversification. Biomechanical principles, which relate body shape to performance and ecology, are one promising avenue for future research.

Genome-Wide Identification of hsp90 Gene in Rainbow Trout (Oncorhynchus mykiss) and Their Regulated Expression in Response to Heat Stress

In this study, we analyzed the gene structure, chemical characterizations, chromosome locations, evolutionary relationship, and expression profile of hsp90 genes with online database. In addition, the expression levels of hsp90s were also investigated under heat stress by quantitative real-time (qRT)-PCR. A total of eight hsp90 genes were identified from the rainbow trout genome. They were all distributed on chromosomes 2, 4, 8, and 13. The molecular weight ranged from 78.93 to 91.39 kDa, and the isoelectric point ranged from 4.84 to 4.96. The eight hsp90 genes were clustered into six subfamilies (A, B, C, D, E, and F). Genetic structure and conserved domain analysis revealed that all eight hsp90 genes had only one exon, and motif 1-motif 10 was shared by most genes. According to RNA-seq analysis of rainbow trout liver and head kidney, a total of seven out of eight genes were significantly upregulated under heat stress, and qRT-PCR was carried out on these seven genes; the expression levels of these genes were significantly upregulated under heat stress. The significantly regulated expressions of hsp90 genes under heat stress indicated that hsp90 genes are involved in heat stress response in rainbow trout. This study provides a theoretical basis for further study on the role of hsp90 in the heat stress tolerance of rainbow trout.