From 2R to 3R: evidence for a fish-specific genome duplication (FSGD)

The cellular expression patterns of gdnfa and gdnfb in the gonads of Nile tilapia and their differential response to retinoic acid

In mammals, glial cell derived neurotrophic factor (GDNF) plays a critical role in the self-renewal and maintenance of spermatogonial stem cells (SSCs) in testis and oogenesis in ovary, whilst retinoic acid (RA), the key factor of meiosis initiation, can downregulate its expression. Unlike mammals, two Gdnf replication genes are widely present in teleost fishes, however, our understanding of them is still poor. In the present study, two paralogous gdnf from Nile tilapia (Oreochromis niloticus), namely as Ongdnfa and Ongdnfb, were characterized, and then their cellular expression profiles in testis and ovary and responsiveness to RA treatment at the tissue and cellular levels were investigated. In phylogenetic tree, the Gdnfa and Gdnfb from teleost fishes were clustered into two different subclasses, respectively, and then clustered with the homologs from cartilaginous fish and tetrapods, suggesting that OnGdnfa and OnGdnfb are orthologous to GDNF and paralogous to each other. Ongdnfa is expressed in Sertoli cells and Leydig cells in testis and oocytes in ovary. The expression pattern of Ongdnfb is similar to Ongdnfa. In the ex vivo testicular organ culture, RA down-regulated the expression of Ongdnfa, whereas up-regulated the expression of Ongdnfb (P < 0.05), suggesting that they have differential responsiveness to RA signaling. RA treatment of the cultured cells derived from adult Nile tilapia testis which have the expression of RA receptors (RAR), Ongdnfa and Ongdnfb further confirmed the above result. Collectively, our study suggests that Ongdnfa and Ongdnfb have non-germline expression patterns in testis and germline expression patterns in ovary; furthermore, they have differential responsiveness to RA signaling, implying that they might have differential biological functions. This study broadens and enriches our understanding of fish GDNF homologs and lays foundation for the study of their biological functions in the future.

Characterization of the somatostatin system in tilapia: implications for growth and reproduction

Somatostatin (SST) plays diverse physiological roles in vertebrates, particularly in regulating growth hormone secretion from the pituitary. While the function of SST as a neuromodulator has been studied extensively, its role in fish and mammalian reproduction remains poorly understood. To address this gap, we investigated the involvement of the somatostatin system in the regulation of growth and reproductive hormones in tilapia. RNA sequencing of mature tilapia brain tissue revealed the presence of three SST peptides: SST6, SST3, and low levels of SST1. Four different isoforms of the somatostatin receptor (SSTR) subfamily were also identified in the tilapia genome. Phylogenetic and synteny analysis identified tiSSTR2-like as the root of the tree, forming two mega clades, with SSTR1 and SSTR4 in one and SSTR2a, SSTR3a, and SSTR5b in the other. Interestingly, the tiSSTR-5 isoforms 5x1, 5x2, and 5x3 were encoded in the sstr3b gene and were an artifact of misperception in the nomenclature in the database. RNA-seq of separated pituitary cell populations showed that SSTRs were expressed in gonadotrophs, with sstr3a enriched in luteinizing hormone (LH) cells and sstr3b significantly enriched in follicle-stimulating hormone (FSH) cells. Notably, cyclosomatostatin, an SSTR antagonist, induced cAMP activity in all SSTRs, with SSTR3a displaying the highest response, whereas octreotide, an SSTR agonist, showed a binding profile like that observed in human receptors. Binding site analysis of tiSSTRs from tilapia pituitary cells revealed the presence of canonical binding sites characteristic of peptide-binding class A G-protein-coupled receptors. Based on these findings, we explored the effect of somatostatin on gonadotropin release from the pituitary in vivo. Whereas cyclosomatostatin increased LH and FSH plasma levels at 2 h post-injection, octreotide decreased FSH levels after 2 h, but the LH levels remained unaffected. Overall, our findings provide important insights into the somatostatin system and its mechanisms of action, indicating a potential role in regulating growth and reproductive hormones. Further studies of the complex interplay between SST, its receptors, and reproductive hormones may advance reproductive control and management in cultured populations.

Genome-wide identification of Fgfr genes and function analysis of Fgfr4 in myoblasts differentiation of Lateolabrax maculatus

Fibroblast growth factor receptors (Fgfrs) are involved in cell proliferation, differentiation, and migration via complex signaling pathways in different tissues. Our previous studies showed that fibroblast growth factor receptor 4 (fgfr4) was detected in the most significant quantitative trait loci (QTL) for growth traits. However, studies focusing on the function of fgfr4 on the growth of bony fish are still limited. In this study, we identified seven fgfr genes in spotted sea bass (Lateolabrax maculatus) genome, namely fgfr1a, fgfr1b, fgfr2, fgfr3, fgfr4, fgfr5a, and fgfr5b. Phylogenetic analysis, syntenic analysis and gene structure analysis were conducted to further support the accuracy of our annotation and classification results. Additionally, fgfr4 showed the highest expression levels among fgfrs during the proliferation and differentiation stages of spotted sea bass myoblasts. To further study the function of fgfr4 in myogenesis, dual-fluorescence in situ hybridization (ISH) assay was conducted, and the results showed co-localization of fgfr4 with marker gene of skeletal muscle satellite cells. By treating differentiating myoblasts cultured in vitro with BLU-554, the mRNA expressions of myogenin (myog) and the numbers of myotubes formed by myoblasts increased significantly compared to negative control group. These results indicated that Fgfr4 inhibits the differentiation of myoblasts in spotted sea bass. Our findings contributed to filling a research gap on fgfr4 in bony fish myogenesis and the theoretical understanding of growth trait regulation of spotted sea bass.

Individual differences in the boldness of female zebrafish are associated with alterations in serotonin function

One of the most prevalent axes of behavioral variation in both humans and animals is risk taking, where individuals that are more willing to take risk are characterized as bold while those that are more reserved are regarded as shy. Brain monoamines (i.e. serotonin, dopamine and noradrenaline) have been found to play a role in a variety of behaviors related to risk taking. Using zebrafish, we investigated whether there was a relationship between monoamine function and boldness behavior during exploration of a novel tank. We found a correlation between serotonin metabolism (5-HIAA:5-HT ratio) and boldness during the initial exposure to the tank in female animals. The DOPAC:DA ratio correlated with boldness behavior on the third day in male fish. There was no relationship between boldness and noradrenaline. To probe differences in serotonergic function in bold and shy fish, we administered a selective serotonin reuptake inhibitor, escitalopram, and assessed exploratory behavior. We found that escitalopram had opposing effects on thigmotaxis in bold and shy female animals: the drug caused bold fish to spend more time near the center of the tank and shy fish spent more time near the periphery. Taken together, our findings indicate that variation in serotonergic function has sex-specific contributions to individual differences in risk-taking behavior.

Expansion and Functional Diversification of Long-Wavelength-Sensitive Opsin in Anabantoid Fishes

Gene duplication is one of the most important sources of novel genotypic diversity and the subsequent evolution of phenotypic diversity. Determining the evolutionary history and functional changes of duplicated genes is crucial for a comprehensive understanding of adaptive evolution. The evolutionary history of visual opsin genes is very dynamic, with repeated duplication events followed by sub- or neofunctionalization. While duplication of the green-sensitive opsins rh2 is common in teleost fish, fewer cases of multiple duplication events of the red-sensitive opsin lws are known. In this study, we investigate the visual opsin gene repertoire of the anabantoid fishes, focusing on the five lws opsin genes found in the genus Betta. We determine the evolutionary history of the lws opsin gene by taking advantage of whole-genome sequences of nine anabantoid species, including the newly assembled genome of Betta imbellis. Our results show that at least two independent duplications of lws occurred in the Betta lineage. The analysis of amino acid sequences of the lws paralogs of Betta revealed high levels of diversification in four of the seven transmembrane regions of the lws protein. Amino acid substitutions at two key-tuning sites are predicted to lead to differentiation of absorption maxima (λmax) between the paralogs within Betta. Finally, eye transcriptomics of B. splendens at different developmental stages revealed expression shifts between paralogs for all cone opsin classes. The lws genes are expressed according to their relative position in the lws opsin cluster throughout ontogeny. We conclude that temporal collinearity of lws expression might have facilitated subfunctionalization of lws in Betta and teleost opsins in general.

DNA methylation reprogramming in teleosts

Early embryonic development is crucially important but also remarkably diverse among animal taxa. Axis formation and cell lineage specification occur due to both spatial and temporal control of gene expression. This complex system involves various signaling pathways and developmental genes such as transcription factors as well as other molecular interactants that maintain cellular states, including several types of epigenetic marks. 5mC DNA methylation, the chemical modification of cytosines in eukaryotes, represents one such mark. By influencing the compaction of chromatin (a high-order DNA structure), DNA methylation can either repress or induce transcriptional activity. Mammals exhibit a reprogramming of DNA methylation from the parental genomes in the zygote following fertilization, and later in primordial germ cells (PGCs). Whether these periods of methylation reprogramming are evolutionarily conserved, or an innovation in mammals, is an emerging question. Looking into these processes in other vertebrate lineages is thus important, and teleost fish, with their extensive species richness, phenotypic diversity, and multiple rounds of whole genome duplication, provide the perfect research playground for answering such a question. This review aims to present a concise state of the art of DNA methylation reprogramming in early development in fish by summarizing findings from different research groups investigating methylation reprogramming patterns in teleosts, while keeping in mind the ramifications of the methodology used, then comparing those patterns to reprogramming patterns in mammals.

Transcriptional regulation of Znt family members znt4, znt5 and znt10 and their function in zinc transport in yellow catfish (Pelteobagrus fulvidraco)

The study characterized the transcriptionally regulatory mechanism and functions of three zinc (Zn) transporters (znt4, znt5 and znt10) in Zn2+ metabolism in yellow catfish (Pelteobagrus fulvidraco), commonly freshwater fish in China and other countries. We cloned the sequences of znt4 promoter, spanning from -1217 bp to +80 bp relative to TSS (1297 bp); znt5, spanning from -1783 bp to +49 bp relative to TSS (1832 bp) and znt10, spanning from -1923 bp to +190 bp relative to TSS (2113 bp). In addition, after conducting the experiments of sequential deletion of promoter region and mutation of potential binding site, we found that the Nrf2 binding site (-607/-621 bp) and Klf4 binding site (-5/-14 bp) were required on znt4 promoter, the Mtf-1 binding site (-1674/-1687 bp) and Atf4 binding site (-444/-456 bp) were required on znt5 promoter and the Atf4 binding site (-905/-918 bp) was required on znt10 promoter. Then, according to EMSA and ChIP, we found that Zn2+ incubation increased DNA affinity of Atf4 to znt5 or znt10 promoter, but decreased DNA affinity of Nrf2 to znt4 promoter, Klf4 to znt4 promoter and Mtf-1 to znt5 promoter. Using fluorescent microscopy, it was revealed that Znt4 and Znt10 were located in the lysosome and Golgi, and Znt5 was located in the Golgi. Finally, we found that znt4 knockdown reduced the zinc content of lysosome and Golgi in the control and zinc-treated group; znt5 knockdown reduced the zinc content of Golgi in the control and zinc-treated group and znt10 knockdown reduced the zinc content of Golgi in the zinc-treated group. High dietary zinc supplement up-regulated Znt4 and Znt5 protein expression. Above all, for the first time, we revealed that Klf4 and Nrf2 transcriptionally regulated the activities of znt4 promoter; Mtf-1 and Atf4 transcriptionally regulated the activities of znt5 promoter and Atf4 transcriptionally regulated the activities of znt10 promoter, which provided innovative regulatory mechanism of zinc transporting in yellow catfish. Our study also elucidated their subcellular location, and regulatory role of zinc homeostasis in yellow catfish.

Characterization of potential spermatogonia biomarker genes in the European eel (Anguilla anguilla)

Identification of specific molecular markers for spermatogonial stem cells in teleost is crucial for enhancing the efficacy of reproductive biotechnologies in aquaculture, such as transplantation and surrogate production in fishes. Since it is not yet possible to distinguish spermatogonial stem cells of European eel (Anguilla anguilla) using specific molecular markers, we isolated spermatogonial cells from immature European eels to find these potential markers. We attempted this by studying three candidate genes: vasa, nanos2, and dnd1. Two vasa (vasa1 and vasa2) genes, nanos2, and dnd1 were identified, characterized, and studied in the muscle, testis, and isolated spermatogonia. Our results showed that vasa1 and vasa2 had the highest levels of expression when measured by qPCR. In situ hybridization and immunochemistry assays showed that the four genes were localized explicitly in type A spermatogonia. However, vasa1 and vasa2 exhibited stronger signals in the immature testicular tissue than the other two potential markers. According to this, vasa1 and vasa2 were found to be the most effective markers for spermatogonial cells in the European eel.

Convergent gene losses and pseudogenizations in multiple lineages of stomachless fishes

The regressive evolution of independent lineages often results in convergent phenotypes. Several teleost groups display secondary loss of the stomach, and four gastric genes, atp4a, atp4b, pgc, and pga2 have been co-deleted in agastric (stomachless) fish. Analyses of genotypic convergence among agastric fishes showed that four genes, slc26a9, kcne2, cldn18a, and vsig1, were co-deleted or pseudogenized in most agastric fishes of the four major groups. kcne2 and vsig1 were also deleted or pseudogenized in the agastric monotreme echidna and platypus, respectively. In the stomachs of sticklebacks, these genes are expressed in gastric gland cells or surface epithelial cells. An ohnolog of cldn18 was retained in some agastric teleosts but exhibited an increased non-synonymous substitution when compared with gastric species. These results revealed novel convergent gene losses at multiple loci among the four major groups of agastric fish, as well as a single gene loss in the echidna and platypus.

Functional identification of two HMGB1 paralogues provides insights into autophagic machinery in teleost

High mobility group box 1 (HMGB1) is a multifunctional regulator that plays different roles in various physiological and pathological processes including cell development, autophagy, inflammation, tumor metastasis, and cell death based on its cellular localization. Unlike mammalian HMGB1, two HMGB1 paralogues (HMGB1a and HMGB1b) have been found in fathead minnow and other fish species and its function as an inflammatory cytokine has been well investigated. However, the role of fish HMGB1 in autophagy regulation has not been well clarified. In the present study, we generated HMGB1 paralogues single (HMGB1a-/- and HMGB1b-/-) and double knockout (DKO) epithelioma papulosum cyprini (EPC) cells from fathead minnow by CRISPR/Cas9 system, and the knockout efficiency of these genes was verified at both gene and protein levels. In this context, the effects of HMGB1 gene knockout on the protein expression of microtubule-associated protein 1 light chain 3 II (LC3-II), an autophagy marker, were determined, showing that single knockout of two HMGB1 paralogues significantly decreased the expression of LC3-II, and these inhibitory effects were further amplified in HMGB1 DKO cells under both basal and rapamycin treatment conditions, indicating the role of two HMGB1 paralogues in fish autophagy. In agreement with this notion, overexpression of HMGB1a or HMGB1b with Flag-tag markedly upregulated LC3-II protein expression. Interestingly, overexpressing two paralogues distributed in both cytoplasm and nucleus. Finally, the role of HMGB1-mediated autophagy was further explored, finding that HMGB1 could interact with Beclin1, a key initiation factor of autophagy. Taken together, these findings highlighted the role of HMGB1 paralogues as the autophagy regulator and increased our understanding of autophagic machinery in teleost.

Moderate replacement of fish oil with palmitic acid-stimulated mitochondrial fusion promotes β-oxidation by Mfn2 interacting with Cpt1α via its GTPase-domain

The mitochondrial matrix serves as the principal locale for the process of fatty acids (FAs) β-oxidation. Preserving the integrity and homeostasis of mitochondria, which is accomplished through ongoing fusion and fission events, is of paramount importance for the effective execution of FAs β-oxidation. There has been no investigation to date into whether and how mitochondrial fusion directly enhances FAs β-oxidation. The underlying mechanism of a balanced FAs ratio favoring hepatic lipid homeostasis remains largely unclear. To address such gaps, the present study was conducted to investigate the mechanism through which a balanced dietary FAs ratio enhances hepatic FAs β-oxidation. The investigation specifically focused on the involvement of Mfn2-mediated mitochondrial fusion in the regulation of Cpt1α in this process. In the present study, the yellow catfish (Pelteobagrus fulvidraco), recognized as a model organism for lipid metabolism, were subjected to eight weeks of in vivo feeding with six distinct diets featuring varying FAs ratios. Additionally, in vitro experiments were conducted to inhibit Mfn2-mediated mitochondrial fusion in isolated hepatocytes, achieved through the transfection of hepatocytes with si-mfn2. Further, deletion mutants for both Mfn2 and Cpt1α were constructed to elucidate the critical regions responsible for the interactions between these two proteins within the system. The key findings were: (1) Substituting palmitic acid (PA) for fish oil (FO) proved to be enhanced in reducing hepatic lipid accumulation. This beneficial effect was primarily attributed to the activation of mitochondrial FAs β-oxidation; (2) The balanced replacement of PA stimulated Mfn2-mediated mitochondrial fusion by diminishing Mfn2 ubiquitination, thereby enhancing its protein retention within the mitochondria; (3) Mfn2-mediated mitochondrial fusion promoted FAs β-oxidation through direct interaction between Mfn2 and Cpt1α via its GTPase-domains, which is essential for the maintenance of Cpt1 activity. Notably, the present research results unveil a previously undisclosed mechanism wherein Mfn2-mediated mitochondrial fusion promotes FAs β-oxidation by directly augmenting the capacity for FA transport into mitochondria (MT), in addition to expanding the mitochondrial matrix. This underscores the pivotal role of mitochondrial fusion in preserving hepatic lipid homeostasis. The present results further confirm that these mechanisms are evolutionarily conserved, extending their relevance from fish to mammals.

Uncoordinated 51-like kinase 1a/b and 2 in fish Megalobrama amblycephala: Molecular cloning, functional characterization, and their potential roles in glucose metabolism

Uncoordinated (Unc) 51-like kinase (ulk1) and ulk2 are closely involved in autophagy activation, but little is known about their roles in regulating glucose homeostasis. In this study, the genes of ulk1a, ulk1b and ulk2 were cloned and characterized in fish Megalobrama amblycephala. All the three genes shared the approximate N-terminal kinase domain and the C-terminal Atg1-like_tMIT domain structure, while the amino acid sequence identity of them are different between M. amblycephala and other vertebrates. Their transcripts were widely observed in various tissues (brain, muscle, gill, heart, spleen, eye, liver, intestine, abdominal adipose and kidney), but differed in tissue expression patterns. During the glucose tolerance test and the insulin tolerance test, the up-regulated transcriptions of ulk1a, ulk1b and ulk2 were all found despite some differences in the temporal patterns. At the same time, the activities of glycolytic enzymes like hexokinase and phosphofructokinase both showed parallel increases. Furthermore, the feeding of a high-carbohydrate diet decreased the transcriptions of ulk1a, ulk1b and ulk2. Collectively, this study demonstrated that ulk1a, ulk1b and ulk2 in M. amblycephala had similar molecular characterizations, but with different conservation and tissue expression patterns. In addition, ulk1/2 might play important roles in maintaining the glucose homeostasis in fish through regulating the glycolytic pathway.

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.

Caspase-1 activates gasdermin A in non-mammals

Gasdermins oligomerize to form pores in the cell membrane, causing regulated lytic cell death called pyroptosis. Mammals encode five gasdermins that can trigger pyroptosis: GSDMA, B, C, D, and E. Caspase and granzyme proteases cleave the linker regions of and activate GSDMB, C, D, and E, but no endogenous activation pathways are yet known for GSDMA. Here, we perform a comprehensive evolutionary analysis of the gasdermin family. A gene duplication of GSDMA in the common ancestor of caecilian amphibians, reptiles, and birds gave rise to GSDMA-D in mammals. Uniquely in our tree, amphibian, reptile, and bird GSDMA group in a separate clade than mammal GSDMA. Remarkably, GSDMA in numerous bird species contain caspase-1 cleavage sites like YVAD or FASD in the linker. We show that GSDMA from birds, amphibians, and reptiles are all cleaved by caspase-1. Thus, GSDMA was originally cleaved by the host-encoded protease caspase-1. In mammals the caspase-1 cleavage site in GSDMA is disrupted; instead, a new protein, GSDMD, is the target of caspase-1. Mammal caspase-1 uses exosite interactions with the GSDMD C-terminal domain to confer the specificity of this interaction, whereas we show that bird caspase-1 uses a stereotypical tetrapeptide sequence to confer specificity for bird GSDMA. Our results reveal an evolutionarily stable association between caspase-1 and the gasdermin family, albeit a shifting one. Caspase-1 repeatedly changes its target gasdermin over evolutionary time at speciation junctures, initially cleaving GSDME in fish, then GSDMA in amphibians/reptiles/birds, and finally GSDMD in mammals.

Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.

Molecular evolution of the Thrombospondin superfamily

Thrombospondins (TSPs) are multidomain, calcium-binding glycoproteins that have wide-ranging roles in vertebrates in cell interactions, extracellular matrix (ECM) organisation, angiogenesis, tissue remodelling, synaptogenesis, and also in musculoskeletal and cardiovascular functions. Land animals encode five TSPs, which assembly co-translationally either as trimers (subgroup A) or pentamers (subgroup B). The vast majority of research has focused on this canonical TSP family, which evolved through the whole-genome duplications that took place early in the vertebrate lineage. With benefit of the growth in genome- and transcriptome-predicted proteomes of a much wider range of animal species, examination of TSPs throughout metazoan phyla has revealed extensive conservation of subgroup B-type TSPs in invertebrates. In addition, these searches established that canonical TSPs are, in fact, one branch within a TSP superfamily that includes other clades designated mega-TSPs, sushi-TSPs and poriferan-TSPs. Despite the apparent simplicity of poriferans and cnidarians as organisms, these phyla encode a greater diversity of TSP superfamily members than vertebrates. We discuss here the molecular characteristics of the TSP superfamily members, current knowledge of their expression profiles and functions in invertebrates, and models for the evolution of this complex ECM superfamily.

Individual differences in the boldness of female zebrafish are associated with alterations in serotonin function

One of the most prevalent axes of behavioral variation in both humans and animals is risk taking, where individuals that are more willing to take risk are characterized as bold while those that are more reserved as shy. Brain monoamines (i.e., serotonin, dopamine, and norepinephrine) have been found to play a role in a variety of behaviors related to risk taking. Genetic variation related to monoamine function have also been linked to personality in both humans and animals. Using zebrafish, we investigated the relationship between monoamine function and boldness behavior during exploration of a novel tank. We found a sex-specific correlation between serotonin metabolism (5-HIAA:5-HT ratio) and boldness that was limited to female animals; there were no relationships between boldness and dopamine or norepinephrine. To probe differences in serotonergic function, we administered a serotonin reuptake inhibitor, escitalopram, to bold and shy fish, and assessed their exploratory behavior. We found that escitalopram had opposing effects on thigmotaxis in female animals with bold fish spending more time near the center of the tank and shy fish spent more time near the periphery. Taken together, our findings suggest that variation in serotonergic function makes sex-specific contributions to individual differences in risk taking behavior.

Evolution of two-pore domain potassium channels and their gene expression in zebrafish embryos

BACKGROUND

The two-pore domain potassium (K2P) channels are a major type of potassium channels that maintain the cell membrane potential by conducting passive potassium leak currents independent of voltage change. They play prominent roles in multiple physiological processes, including neuromodulation, perception of pain, breathing and mood control, and response to volatile anesthetics. Mutations in K2P channels have been linked to many human diseases, such as neuronal and cardiovascular disorders and cancers. Significant progress has been made to understand their protein structures, physiological functions, and pharmacological modifiers. However, their expression and function during embryonic development remain largely unknown.

RESULTS

We employed the zebrafish model and identified 23 k2p genes using BLAST search and gene cloning. We first analyzed vertebrate K2P channel evolution by phylogenetic and syntenic analyses. Our data revealed that the six subtypes of the K2P genes have already evolved in invertebrates long before the emergence of vertebrates. Moreover, the vertebrate K2P gene number increased, most likely due to two whole-genome duplications. Furthermore, we examined zebrafish k2p gene expression during early embryogenesis by in situ hybridization. Each subgroup's genes showed similar but distinct gene expression domains with some exceptions. Most of them were expressed in neural tissues consistent with their known function of neural excitability regulation. However, a few k2p genes were expressed temporarily in specific tissues or organs, suggesting that these K2P channels may be needed for embryonic development.

CONCLUSIONS

Our phylogenetic and developmental analyses of K2P channels shed light on their evolutionary history and potential roles during embryogenesis related to their physiological functions and human channelopathies.

Caspase-1 activates gasdermin A in non-mammals

Gasdermins oligomerize to form pores in the cell membrane, causing regulated lytic cell death called pyroptosis. Mammals encode five gasdermins that can trigger pyroptosis: GSDMA, B, C, D, and E. Caspase and granzyme proteases cleave the linker regions of and activate GSDMB, C, D, and E, but no endogenous activation pathways are yet known for GSDMA. Here, we perform a comprehensive evolutionary analysis of the gasdermin family. A gene duplication of GSDMA in the common ancestor of caecilian amphibians, reptiles and birds gave rise to GSDMA-D in mammals. Uniquely in our tree, amphibian, reptile and bird GSDMA group in a separate clade than mammal GSDMA. Remarkably, GSDMA in numerous bird species contain caspase-1 cleavage sites like YVAD or FASD in the linker. We show that GSDMA from birds, amphibians, and reptiles are all cleaved by caspase-1. Thus, GSDMA was originally cleaved by the host-encoded protease caspase-1. In mammals the caspase-1 cleavage site in GSDMA is disrupted; instead, a new protein, GSDMD, is the target of caspase-1. Mammal caspase-1 uses exosite interactions with the GSDMD C-terminal domain to confer the specificity of this interaction, whereas we show that bird caspase-1 uses a stereotypical tetrapeptide sequence to confer specificity for bird GSDMA. Our results reveal an evolutionarily stable association between caspase-1 and the gasdermin family, albeit a shifting one. Caspase-1 repeatedly changes its target gasdermin over evolutionary time at speciation junctures, initially cleaving GSDME in fish, then GSDMA in amphibians/reptiles/birds, and finally GSDMD in mammals.

The Importance of Monitoring Non-clonal Chromosome Aberrations (NCCAs) in Cancer Research

Cytogenetic analysis has traditionally focused on the clonal chromosome aberrations, or CCAs, and considered the large number of diverse non-clonal chromosome aberrations, or NCCAs, as insignificant noise. Our decade-long karyotype evolutionary studies have unexpectedly demonstrated otherwise. Not only the baseline of NCCAs is associated with fuzzy inheritance, but the frequencies of NCCAs can also be used to reliably measure genome or chromosome instability (CIN). According to the Genome Architecture Theory, CIN is the common driver of cancer evolution that can unify diverse molecular mechanisms, and genome chaos, including chromothripsis, chromoanagenesis, and polypoidal giant nuclear and micronuclear clusters, and various sizes of chromosome fragmentations, including extrachromosomal DNA, represent some extreme forms of NCCAs that play a key role in the macroevolutionary transition. In this chapter, the rationale, definition, brief history, and current status of NCCA research in cancer are discussed in the context of two-phased cancer evolution and karyotype-coded system information. Finally, after briefly describing various types of NCCAs, we call for more research on NCCAs in future cytogenetics.

Manipulation of Proteostasis Networks in Transgenic ZAAT Zebrafish via CRISPR-Cas9 Gene Editing

The CRISPR-Cas9 genome editing system is used to induce mutations in genes of interest resulting in the loss of functional protein. A transgenic zebrafish α1-antitrypsin deficiency (AATD) model displays an unusual phenotype, in that it lacks the hepatic accumulation of the misfolding Z α1-antitrypsin (ZAAT) evident in human and mouse models. Here we describe the application of the CRISPR-Cas9 system to generate mutant zebrafish with defects in key proteostasis networks likely to be involved in the hepatic processing of ZAAT in this model. We describe the targeting of the atf6a and man1b1 genes as examples.

Gene Characterization of Nocturnin Paralogues in Goldfish: Full Coding Sequences, Structure, Phylogeny and Tissue Expression

The aim of this work is the full characterization of all the nocturnin (noc) paralogues expressed in a teleost, the goldfish. An in silico analysis of the evolutive origin of noc in Osteichthyes is performed, including the splicing variants and new paralogues appearing after teleostean 3R genomic duplication and the cyprinine 4Rc. After sequencing the full-length mRNA of goldfish, we obtained two isoforms for noc-a (noc-aa and noc-ab) with two splice variants (I and II), and only one for noc-b (noc-bb) with two transcripts (II and III). Using the splicing variant II, the prediction of the secondary and tertiary structures renders a well-conserved 3D distribution of four α-helices and nine β-sheets in the three noc isoforms. A synteny analysis based on the localization of noc genes in the patrilineal or matrilineal subgenomes and a phylogenetic tree of protein sequences were accomplished to stablish a classification and a long-lasting nomenclature of noc in goldfish, and valid to be extrapolated to allotetraploid Cyprininae. Finally, both goldfish and zebrafish showed a broad tissue expression of all the noc paralogues. Moreover, the enriched expression of specific paralogues in some tissues argues in favour of neo- or subfunctionalization.

The evolution of NLRC3 subfamily genes in Sebastidae teleost fishes

BACKGROUND

With more than 36,000 valid fish species, teleost fishes constitute the most species-rich vertebrate clade and exhibit extensive genetic and phenotypic variation, including diverse immune defense strategies. NLRC3 subfamily genes, which are specific to fishes, play vital roles in the immune system of teleosts. The evolution of teleosts has been impacted by several whole-genome duplication (WGD) events, which might be a key reason for the expansions of the NLRC3 subfamily, but detailed knowledge of NLRC3 subfamily evolution in the family Sebastidae is still limited.

RESULTS

Phylogenetic inference of NLRC3 subfamily protein sequences were conducted to evaluate the orthology of NLRC3 subfamily genes in black rockfish (Sebastes schlegilii), 13 other fish species from the families Sebastidae, Serranidae, Gasterosteidae and Cyclopteridae, and three species of high vertebrates (bird, reptile and amphibian). WGD analyses were used to estimate expansions and contractions of the NLRC3 subfamily, and patterns of expression of NLRC3 subfamily genes in black rockfish following bacterial infections were used to investigate the functional roles of these genes in the traditional and mucosal immune system of the Sebastidae. Different patterns of gene expansions and contractions were observed in 17 fish and other species examined, and one and two whole-genome duplication events were observed in two members of family Sebastidae (black rockfish and honeycomb rockfish, Sebastes umbrosus), respectively. Subsequently, 179 copy numbers of NLRC3 genes were found in black rockfish and 166 in honeycomb rockfish. Phylogenetic analyses corroborated the conservation and evolution of NLRC3 orthologues between Sebastidae and other fish species. Finally, differential expression analyses provided evidence of the immune roles of NLRC3 genes in black rockfish during bacterial infections and gene ontology analysis also indicated other functional roles.

CONCLUSIONS

We hypothesize that NLRC3 genes have evolved a variety of different functions, in addition to their role in the immune response, as a result of whole genome duplication events during teleost diversification. Importantly, this study had underscored the importance of sampling across taxonomic groups, to better understand the evolutionary patterns of the innate immunity system on which complex immunological novelties arose. Moreover, the results in this study could extend current knowledge of the plasticity of the immune system.

Chromosome-level genome of butterflyfish unveils genomic features of unique colour patterns and morphological traits

Chaetodontidae, known as butterflyfishes, are typical fish in coral ecosystems, exhibiting remarkable interspecific differences including body colour patterns and feeding ecology. In this study, we report genomes of three butterflyfish species (Chelmon rostratus, Chaetodon trifasciatus and Chaetodon auriga) and a closely related species from the Pomacanthidae family, Centropyge bicolour, with an average genome size of 65,611 Mb. Chelmon rostratus, comprising 24 chromosomes assembled to the chromosome level, could be served as a reference genome for butterflyfish. By conducting a collinearity analysis between butterflyfishes and several fishes, we elucidated the specific and conserved genomic features of butterflyfish, with particular emphasis on novel genes arising from tandem duplications and their potential functions. In addition to the two melanocyte-specific tyr genes commonly found in fish, we found the gene tyrp3, a new tyrosinase-related proteins gene in the reef fish, including butterflyfish and clownfish, implicating their involvement in the pigmentation diversity of fish. Additionally, we observed a tandem duplication expansion of three copies of nell1 gene in C. rostratus genome, which likely contribute to its unique jaw development and distinctive morphology of its sharp mouth. These results provided valuable genomic resources for further investigations into the genetic diversity and evolutionary adaptations of reef fish.

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.

Identification, characterization, and transcription of serotonin receptors in rainbow trout (Oncorhynchus mykiss) in response to bacterial infection and salinity changes

Serotonergic system is involved in the regulation of physiological functions and behavioral traits including cognition, memory, aggression, stress coping, appetite and immunomodulation. Serotonin exerts its functions via binding distinct serotonin receptors which are classified into 7 groups. Salmonid exhibits expanded functional gene copies due to salmonid-specific whole genome duplication. However, serotonin receptor (htr) repertoire is not fully identified in rainbow trout (Oncorhynchus mykiss). In this study, we identified 39 htr genes, including 14 htr1, 4 htr2, 4 htr2 like, 3 htr3, 4 htr4, 2 htr5, 2 htr6, and 6 htr7 subtypes. We investigated physiological functions of serotonin receptors in response to bacterial pathogens exposure and salinity changes. We showed htr1, htr2, htr4 and htr7 subtypes were associated with immunomodulation in response to Vibrio anguillarum or Aeromonas salmonicida infection. Saltwater (salinity of 15) transfer significantly altered htr1, htr2, htr4, and htr7 subtypes, suggesting trout Htr was associated with osmoregulation. We further showed residues interacted with inverse agonist (methiothepin) and serotonin analogue (5-Carboxamidotryptamine) were conserved between trout and human, suggesting exogenous ligands targeting human HTRs might have a role in aquaculture. This study showed duplicated trout Htrs might be physiologically neofunctionalized and potentially exhibit pleiotropic effects in regulating immunomodulation and osmoregulation.

Two CcCCL19bs orchestrate an antibacterial immune response in Yellow River carp (Cyprinus carpio haematopterus)

Chemokines are a group of chemotactic cytokines with an essential role in homeostasis as well as immunity via specific G protein-coupled receptors and atypical receptors. In our study, two Yellow River carp (Cyprinus carpio haematopterus) CCL19b genes (CcCCL19bs), tentatively named CcCCL19b_a and CcCCL19b_b, were cloned. The open reading frames (ORFs) of CcCCL19b_a and CcCCL19b_b were both 333 bp that encoded a 12 kDa protein with 110 amino acid residues. CcCCL19bs contained a signal peptide and a SCY domain with four typical conserved cysteine residues. The two CcCCL19b proteins shared high similarities with each other in both secondary and three-dimensional structure. Phylogenetic analysis showed that CcCCL19bs and other CCL19bs from tetraploid cyprinid fish were clustered into one clade. CcCCL19bs were highly expressed in gill and intestine in healthy fish, and a significant up-regulation of gene expression after Aeromonas hydrophila infection and poly(I:C) stimulation was observed in gill, liver, and head kidney. Furthermore, chemotaxis and antibacterial activity of CcCCL19bs were studied. The results indicated that recombinant CcCCL19b_a and CcCCL19b_b protein (rCcCCL19b_a and rCcCCL19b_b) exhibited significant attraction to primary head kidney leukocytes (HKLs). Meanwhile, both of rCcCCL19bs could promote the proliferation of HKLs, and significantly up-regulate the expressions of IL-1β, CCR7, and IL-6, and down-regulate the expression of IL-10 in primary HKLs. In vitro, rCcCCL19bs could bind and aggregate A. hydrophila and Staphylococcus aureus. The rCcCCL19bs exhibited significant antibacterial activity against A. hydrophila, but not S. aureus. Moreover, they inhibited the growth of A. hydrophila and S. aureus. In vivo, overexpression of CcCCL19bs contributed to the bacterial clearance. These studies suggested that CcCCL19bs orchestrate an antibacterial immune response.

Interleukin-12 induces IFN-γ secretion and STAT signaling implying its potential regulation of Th1 cell response in Nile tilapia

As a pleiotropic cytokine consisting of IL-12p35 and IL-12p40, Interleukin-12 (IL-12) features in inflammation regulation and anti-bacterial immunity. While IL-12 homologs have been identified in non-mammalian species, the precise mechanisms by which IL-12 contributes to early adaptive immune responses in vertebrates remain incompletely understood. Herein, an evolutionary conserved Oreochromis niloticus IL-12 (defined as OnIL-12) was identified by synteny characterization, structural comparisons and phylogenetic pattern of IL-12p35b and IL-12p40a. IL-12p35b and IL-12p40a exhibited widespread expression in lymphoid-related tissues of tilapia, while their mRNA expression in head-kidney demonstrated a significant increase after Edwardsiella piscicida infection. Compared with other lymphocytes, recombinant OnIL-12 (rOnIL-12) displayed stronger affinity binding to T cells. Although stimulation of lymphocytes with the p35b or p40a subunit resulted in a significant induction of IFN-γ expression, rOnIL-12 showed stronger potential to promote IFN-γ expression than these subunits. rOnIL-12 not only elevated the mRNA expression level Th1 cell-associated transcription factor T-bet in lymphocytes, but also increased the proportion of CD4-1+IFN-γ+ lymphocytes. Moreover, the mRNA and phosphorylation levels of STAT1, STAT3, STAT4 and STAT5 were enhanced by rOnIL-12. These findings will offer previous evidence for further exploration into the regulatory mechanisms of Th1 cellular immunity in early vertebrates.

ncOrtho: efficient and reliable identification of miRNA orthologs

MicroRNAs (miRNAs) are post-transcriptional regulators that finetune gene expression via translational repression or degradation of their target mRNAs. Despite their functional relevance, frameworks for the scalable and accurate detection of miRNA orthologs are missing. Consequently, there is still no comprehensive picture of how miRNAs and their associated regulatory networks have evolved. Here we present ncOrtho, a synteny informed pipeline for the targeted search of miRNA orthologs in unannotated genome sequences. ncOrtho matches miRNA annotations from multi-tissue transcriptomes in precision, while scaling to the analysis of hundreds of custom-selected species. The presence-absence pattern of orthologs to 266 human miRNA families across 402 vertebrate species reveals four bursts of miRNA acquisition, of which the most recent event occurred in the last common ancestor of higher primates. miRNA families are rarely modified or lost, but notable exceptions for both events exist. miRNA co-ortholog numbers faithfully indicate lineage-specific whole genome duplications, and miRNAs are powerful markers for phylogenomic analyses. Their exceptionally low genetic diversity makes them suitable to resolve clades where the phylogenetic signal is blurred by incomplete lineage sorting of ancestral alleles. In summary, ncOrtho allows to routinely consider miRNAs in evolutionary analyses that were thus far reserved to protein-coding genes.

New approaches concerning the testis of Astyanax lacustris (Characidae): immunohistochemical studies

Astyanax lacustris, locally known as lambari-do-rabo-amarelo, is a study model for Neotropical fish. Testis of A. lacustris shows deep morphophysiological changes throughout the annual reproductive cycle. This work analyzed the distribution of claudin-1, actin, and cytokeratin as elements of the cytoskeleton in germinal epithelium and interstitium; the distribution of type I collagen, fibronectin, and laminin as extracellular matrix compounds; and the localization of androgen receptor in the testis of this species. Claudin-1, cytokeratin, and actin were present in the Sertoli cells and modified Sertoli cells, and actin was also detected in peritubular myoid cells. Type I collagen were in the interstitial tissue, laminin in the basement membrane of germinal epithelium and endothelium, but fibronectin was additionally detected in the germinal epithelium compartment. The labeling of androgen receptor was higher in peritubular myoid cells and undifferentiated spermatogonia, and weaker labeling was detected in type B spermatogonia. Therefore, the present work highlights new aspects of the biology of the testis of A. lacustris, and contribute to amplify the understanding of this organ.

Molecular evolution of the hemoglobin gene family across vertebrates

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

The rete mirabile: a possible control site for swimbladder function

In a recent study, a large number of transport proteins was detected in the transcriptome and proteome of saline perfused rete mirabile tissue of the European eel. In this study, the data set was reanalyzed for the presence of receptor proteins and proteins involved in intracellular signaling pathways. A large number of expressed receptor proteins and proteins involved in intracellular signal transduction was detected. Several G-protein-coupled receptor signal pathways were significantly enriched in their expression level, in particular receptors and signaling pathways involved in the control of blood flow. The enriched signaling pathways also include pathways involved in trafficking of crucial transport proteins like, monocarboxylate transporters, V-ATPase, and aquaporin. The data, therefore, suggest that the rete mirabile has the capacity to control swimbladder function by regulating blood flow and by modifying countercurrent multiplication.

Roles of the 5-HT2C receptor on zebrafish sociality

Serotonin (5-HT) receptors have been implicated in social behavior in vertebrates. Zebrafish (Danio rerio) have been increasingly being used behavioral neuroscience to study the neurobiological correlates of behavior, including sociality. Nonetheless, the role of 5-HT_2C receptors in different social functions were not yet studied in this species. Zebrafish were treated with the agonist MK-212 (2 mg/kg) or the antagonist RS-102221 (2 mg/kg) and tested in the social interaction and social novelty tests, conditional approach test, or mirror-induced aggressive displays. MK-212 increased preference for an unknown conspecific in the social investigation test, but also increased preference for the known conspecific in the social novelty test; RS-102221, on the other hand, decreased preference in the social investigation test but increased preference for the novel conspecific in the social novelty test. MK-212 also decreased predator inspection in the conditional approach test. While RS-102221 decreased time in the display zone in the mirror-induced aggressive display test, it increased display duration. Overall, these results demonstrate the complex role of 5-HT_2C receptors in different social contexts in zebrafish, revealing a participation in social plasticity in vertebrates.

Fatty acid-binding protein genes in gilthead seabream: molecular cloning and nutritional regulation under low water temperatures

The molecular characteristics and tissue disruption of 10 fatty acid-binding protein (fabp) genes in gilthead seabream (Sparus aurata) were investigated, and their expression levels were found in the fish fed diets with different vegetable oil (VO) sources, which may explore the potential function of fabp genes in S. aurata. For this purpose, the open reading frames of fabp genes involved in the transport and ß-oxidation of fatty acids (FA) were molecularly cloned and characterized. S. aurata was then exposed to a two-staged feeding trial (the grow-out period following a wash-out period) at low water temperatures. In the grow-out period, the fish were fed diets containing 50% and 100% ratios of various VOs for 60 days, and in the wash-out period, the fish were fed a diet containing 100% fish oil (FO) for 30 days. It has been determined that (a) S. aurata and vertebrate fabp/FABP genes are orthologues; (b) spatio-temporal differences in tissue-specific patterns of fabp genes differ importantly; for instance, the difference between the highest and lowest values reaches 13 × 10⁵ -fold in the fabp10a; and (c) VO-based diets upregulated fabp transcript levels in the liver and muscle with some exceptions, such as liver fabp11a and muscle fabp7a. Gene expressions of only the hepatic fabp7b and fabp10a genes were diminished at the end of the wash-out period. In this study, the authors provide further evidence that dietary FAs affect fabp mRNA expressions in S. aurata. This might be useful in the nutritional control of fabp genes to maintain lipid homeostasis in marine fish fed VO-based diets at low water temperatures.

Long-term waterborne Cu2+ exposure affects collagen metabolism in fish

Copper pollution might have a negative effect on collagen metabolism in fish. To test this hypothesis, we exposed an important economical fish, silver pomfret (Pampus argenteus), to three concentrations of Cu²⁺ for up to 21 days to simulate natural exposure to copper. With increasing copper exposure concentration and time, hematoxylin and eosin staining and picrosirius red staining revealed extensive vacuolization, cell necrosis, and tissue structure destruction, and a change of type and abnormal accumulation of collagen in the liver, intestine, and muscle tissues. To further study the mechanism of collagen metabolism disorder caused by copper exposure, we cloned and analyzed a key collagen metabolism regulation gene, timp, of silver pomfret. The full-length timp2b cDNA was 1035 bp with an open reading frame of 663 bp, encoding a protein of 220 amino acids. Copper treatment significantly increased the expression of akts, erks, and fgfs genes and decreased the mRNA and protein expression of Timp2b and MMPs. Finally, we constructed a silver pomfret muscle cell line (PaM) for the first time and used PaM Cu²⁺ exposure models (450 μM Cu²⁺ exposure for 9 h) to examine regulation function of the timp2b-mmps system. We knocked down or overexpressed timp2b in the model, and found that downregulation of mmps expression and upregulation of akt/erk/fgf were further aggravated in the timp2b^- group (subjected to RNA interference), whereas some recovery was achieved in the timp2b⁺ group (overexpression). These results indicated that long-term excessive copper exposure can lead to tissue damage and abnormal collagen metabolism in fish, which might be caused by the alteration of akt/erk/fgf expression, which disrupts the effects of the timp2b-mmps system on extracellular matrix balance. The present study assessed the impact of copper on the collagen of fish and clarified its regulatory mechanism, providing a basis for toxicity of copper pollution study.

N-Glycosylation of LRP6 by B3GnT2 Promotes Wnt/β-Catenin Signalling

Reception of Wnt signals by cells is predominantly mediated by Frizzled receptors in conjunction with a co-receptor, the latter being LRP6 or LRP5 for the Wnt/β-catenin signalling pathway. It is important that cells maintain precise control of receptor activation events in order to properly regulate Wnt/β-catenin signalling as aberrant signalling can result in disease in humans. Phosphorylation of the intracellular domain (ICD) of LRP6 is well known to regulate Wntβ-catenin signalling; however, less is known for regulatory post-translational modification events within the extracellular domain (ECD). Using a cell culture-based expression screen for functional regulators of LRP6, we identified a glycosyltransferase, B3GnT2-like, from a teleost fish (medaka) cDNA library, that modifies LRP6 and regulates Wnt/β-catenin signalling. We provide both gain-of-function and loss-of-function evidence that the single human homolog, B3GnT2, promotes extension of polylactosamine chains at multiple N-glycans on LRP6, thereby enhancing trafficking of LRP6 to the plasma membrane and promoting Wnt/β-catenin signalling. Our findings further highlight the importance of LRP6 as a regulatory hub in Wnt signalling and provide one of the few examples of how a specific glycosyltransferase appears to selectively target a signalling pathway component to alter cellular signalling events.

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.

Efficiently Editing Multiple Duplicated Homeologs and Alleles for Recurrent Polyploids

Research on the evolutionary fate of duplicated genes in recurrent polyploids is scarce due to the difficulties in disentangling the different homeologs and alleles of duplicated genes. This chapter describes the detailed procedures to identify different homeologs and alleles of duplicated genes, to analyze their molecular characteristics, and to reveal their functional divergence by gene editing with CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/CRISPR-associated system 9). Using the gene editing approach, we efficiently constructed multiple knockout mutant lines with single or simultaneously disrupted different homeologs or alleles in a recurrent polyploid fish, demonstrating its usability for targeting and mutating multiple divergent homeologs and alleles in recurrent duplicated genomes.

Improving Orthologous Signal and Model Fit in Datasets Addressing the Root of the Animal Phylogeny

There is conflicting evidence as to whether Porifera (sponges) or Ctenophora (comb jellies) comprise the root of the animal phylogeny. Support for either a Porifera-sister or Ctenophore-sister tree has been extensively examined in the context of model selection, taxon sampling, and outgroup selection. The influence of dataset construction is comparatively understudied. We re-examine five animal phylogeny datasets that have supported either root hypothesis using an approach designed to enrich orthologous signal in phylogenomic datasets. We find that many component orthogroups in animal datasets fail to recover major lineages as monophyletic with the exception of Ctenophora, regardless of the supported root. Enriching these datasets to retain orthogroups recovering ≥3 major lineages reduces dataset size by up to 50% while retaining underlying phylogenetic information and taxon sampling. Site-heterogeneous phylogenomic analysis of these enriched datasets recovers both Porifera-sister and Ctenophora-sister positions, even with additional constraints on outgroup sampling. Two datasets which previously supported Ctenophora-sister support Porifera-sister upon enrichment. All enriched datasets display improved model fitness under posterior predictive analysis. While not conclusively rooting animals at either Porifera or Ctenophora, we do see an increase in signal for Porifera-sister and a decrease in signal for Ctenophore-sister when data are filtered for orthologous signal. Our results indicate that dataset size and construction as well as model fit influence animal root inference.

Parvalbumin: A Major Fish Allergen and a Forensically Relevant Marker

Parvalbumins (PVALBs) are low molecular weight calcium-binding proteins. In addition to their role in many biological processes, PVALBs play an important role in regulating Ca2+ switching in muscles with fast-twitch fibres in addition to their role in many biological processes. The PVALB gene family is divided into two gene types, alpha (α) and beta (β), with the β gene further divided into two gene types, beta1 (β1) and beta2 (β2), carrying traces of whole genome duplication. A large variety of commonly consumed fish species contain PVALB proteins which are known to cause fish allergies. More than 95% of all fish-induced food allergies are caused by PVALB proteins. The authentication of fish species has become increasingly important as the seafood industry continues to grow and the growth brings with it many cases of food fraud. Since the PVALB gene plays an important role in the initiation of allergic reactions, it has been used for decades to develop alternate assays for fish identification. A brief review of the significance of the fish PVALB genes is presented in this article, which covers evolutionary diversity, allergic properties, and potential use as a forensic marker.

Klf4-Sirt3/Pparα-Lcad pathway contributes to high phosphate-induced lipid degradation

BACKGROUND

Phosphorus commonly reduces lipid deposition in the vertebrates. However, the underlying mechanisms involved in the process remain unclear.

METHODS

Yellow catfish were given three experimental diets with dietary phosphate levels of 3.22, 6.47 and 7.99 g Pi kg^- 1, respectively, for 8 weeks. The contents of triglyceride, non-esterified free fatty acids, adenosine triphosphate, nicotinamide adenine dinucleotide, nicotinamide adenine dinucleotide, enzymatic activities, mRNA and protein expression were determined in the intestinal tissues. Hematoxylin and eosin, Oil Red O staining, and transmission electron microscope were performed for intestinal tissues. Primary intestinal epithelial cells were isolated from yellow catfish intestine. Western blot analysis, Immunoprecipitation assays, Immunofluorescence staining, and RNA extraction and quantitative real-time PCR were decided. Luciferase reporter assays and electrophoretic mobility shift assay were used to evaluate the function of Sirt3, PPARα and Lcad promoters.

RESULTS

High dietary phosphate intake activated intestinal phosphate absorption and excretion, and reduced lipid deposition through increasing lipolysis in the intestine. Moreover, phosphate incubation increased the mRNA and protein expression of krüppel like factor 4 (klf4), silent mating-type information regulation 2 homolog 3 (sirt3), peroxisome proliferator activated receptor alpha (pparα) and long chain acyl-CoA dehydrogenase (lcad) in the intestinal epithelial cells (IECs), and klf4 knockdown attenuated the phosphate-induced increase of protein levels of Sirt3, Pparα and Lcad. Further investigation found that Klf4 overexpression increased the activity of sirt3 and pparα promoters, which in turn reduced the acetylation and protein level of Lcad.

CONCLUSION

Dietary Pi excess induced lipid degradation by the activation of the Klf4-Sirt3/Pparα-Lcad pathway in the intestine and primary IECs. Video Abstract.

Integrative omics-analysis of lipid metabolism regulation by peroxisome proliferator-activated receptor a and b agonists in male Atlantic cod

Lipid metabolism is essential in maintaining energy homeostasis in multicellular organisms. In vertebrates, the peroxisome proliferator-activated receptors (PPARs, NR1C) regulate the expression of many genes involved in these processes. Atlantic cod (Gadus morhua) is an important fish species in the North Atlantic ecosystem and in human nutrition, with a highly fatty liver. Here we study the involvement of Atlantic cod Ppar a and b subtypes in systemic regulation of lipid metabolism using two model agonists after in vivo exposure. WY-14,643, a specific PPARA ligand in mammals, activated cod Ppara1 and Ppara2 in vitro. In vivo, WY-14,643 caused a shift in lipid transport both at transcriptional and translational level in cod. However, WY-14,643 induced fewer genes in the fatty acid beta-oxidation pathway compared to that observed in rodents. Although GW501516 serves as a specific PPARB/D ligand in mammals, this compound activated cod Ppara1 and Ppara2 as well as Pparb in vitro. In vivo, it further induced transcription of Ppar target genes and caused changes in lipid composition of liver and plasma. The integrative approach provide a foundation for understanding how Ppars are engaged in regulating lipid metabolism in Atlantic cod physiology. We have shown that WY-14,643 and GW501516 activate Atlantic cod Ppara and Pparb, affect genes in lipid metabolism pathways, and induce changes in the lipid composition in plasma and liver microsomal membranes. Particularly, the combined transcriptomic, proteomics and lipidomics analyses revealed that effects of WY-14,643 on lipid metabolism are similar to what is known in mammalian studies, suggesting conservation of Ppara functions in mediating lipid metabolic processes in fish. The alterations in the lipid profiles observed after Ppar agonist exposure suggest that other chemicals with similar Ppar receptor affinities may cause disturbances in the lipid regulation of fish. Model organism: Atlantic cod (Gadus morhua). LSID: urn:lsid:zoobank.org:act:389BE401-2718-4CF2-BBAE-2E13A97A5E7B. COL Identifier: 6K72F.

The evolutionary history and spectral tuning of vertebrate visual opsins

Many animals depend on the sense of vision for survival. In eumetazoans, vision requires specialized, light-sensitive cells called photoreceptors. Light reaches the photoreceptors and triggers the excitation of light-detecting proteins called opsins. Here, we describe the story of visual opsin evolution from the ancestral bilaterian to the extant vertebrate lineages. We explain the mechanisms determining color vision of extant vertebrates, focusing on opsin gene losses, duplications, and the expression regulation of vertebrate opsins. We describe the sequence variation both within and between species that has tweaked the sensitivities of opsin proteins towards different wavelengths of light. We provide an extensive resource of wavelength sensitivities and mutations that have diverged light sensitivity in many vertebrate species and predict how these mutations were accumulated in each lineage based on parsimony. We suggest possible natural and sexual selection mechanisms underlying these spectral differences. Understanding how molecular changes allow for functional adaptation of animals to different environments is a major goal in the field, and therefore identifying mutations affecting vision and their relationship to photic selection pressures is imperative. The goal of this review is to provide a comprehensive overview of our current understanding of opsin evolution in vertebrates.

Sex Lethal Gene Manipulates Gonadal Development of Medaka, Oryzias latipes, through Estrogenic Interventions

Germ cells are pivotal for gonadal sexuality maintenance and reproduction. Sex lethal (sxl), the somatic sex determining gene of Drosophila, is the known regulator and initiator of germ cell femininity in invertebrates. However, the role of the Sxl homologue has rarely been investigated in vertebrates. So, we used medaka to clarify the role of sxl in vertebrate gonadogenesis and sexuality and identified two Sxl homologues, i.e., Sxl1a and Sxl1b. We found that sxl1a specifically expresses in the primordial germ cells (PGC), ovary, (early gonia and oocytes), while sxl1b distributions are ubiquitous. An mRNA overexpression of sxl1a accelerated germ cell numbers in 10 DAH XY fish, and sxl1a knockdown (KD), on the other hand, induced PGC mis-migration, aberrant PGC structuring and ultimately caused significant germ cell reduction in XX fish. Using an in vitro promoter analysis and in vivo steroid treatment, we found a strong link between sxl1a and estrogenic germ cell-population maintenance. Further, using sxl1a-KD and erβ2-knockout fish, we determined that sxl1 acts through erβ2 and controls PGC sexuality. Cumulatively, our study highlights the novel role of sxl1a in germ cell maintenance and sexual identity assignment and thus might become a steppingstone to understanding the commonalities of animal sexual development.

Identification and functional characteristics of two TLR5 subtypes in S. grahami

TLR5, as a member of Toll-like receptors (TLRs) family in mammals, is responsible for recognizing bacterial flagellin and initiating innate immunity, but its function is still unclear in fish species. In this study, two family members of TLR5 were cloned and identified from Sinocyclocheilus grahami (S. grahami), named sgTLR5a and sgTLR5b. The length of coding sequence of sgTLR5a and sgTLR5b is 2,622 bp and 2,658 bp, encoding 873 and 885 amino acids, respectively. Molecular phylogenetic analysis indicates that sgTLR5a and sgTLR5b have the closest genetic relationship with TLR5M (membrane-type) of Cyprinus carpio and Schizothorax prenanti, respectively. sgTLR5a and sgTLR5b were widely expressed in various tested tissues, of which the expression levels were the highest in skin tissue. After stimulations of Aeromonas hydrophila (A. hydrophila) and flagellin, the expression levels of sgTLR5a and sgTLR5b in liver, spleen and head kidney tissues were strongly up-regulated, but LPS stimulation only increased the expression of sgTLR5b in these tissues. The luciferase reporter assay displayed that sgTLR5a and sgTLR5b could specifically recognize bacterial flagellin and A. hydrophila and activate the downstream NF-κB signaling pathway in HEK293T cells. Moreover, the overexpression of sgTLR5a and sgTLR5b in EPC cells up-regulated the expression levels of IL-8 and TNF. sgTLR5a and sgTLR5b were observed to locate in the intracellular region by confocal microscope. Interestingly, we found that the NF-κB signaling pathway was positively regulated by co-transfecting sgTLR5a or sgTLR5b with TLR trafficking chaperone sgUNC93B1. In conclusion, our results reveal sgTLR5a and sgTLR5b may play an important role in antibacterial response by activating the NF-κB signaling pathway.

Three IRF4 paralogs act as negative regulators of type Ⅰ IFN responses in yellow catfish (Pelteobagrus fulvidraco)

IRF4 is a master member of the interferon regulatory factor (IRF) family playing vital regulatory roles in immune system development and function. Tetrapods have a single-copy IRF4 gene, while teleosts harbor duplicated IRF4 genes. This work describes three IRF4 paralogs from yellow catfish (Pelteobagrus fulvidraco), designated PfIRF4A, PfIRF4B and PfIRF4B-like. These genes all contain a typical IRF structural architecture. Phylogenic and synteny analyses indicate that they should arise from the teleost-specific whole-genome duplication. PfIRF4 genes are abundantly expressed in the immune-related tissues and upregulated by PolyI:C, LPS, and Edwardsiella ictaluri. Ectopic expression of these genes inhibits the activation of fish type Ⅰ IFN promoters and downregulates the transcription levels of IFN-responsive genes, thus allowing the efficient replication of a fish rhabdovirus, spring viremia of carp virus (SVCV). PfIRF4s possess a repressive effect on MyD88-mediated activation of IFN and NF-κB. Some differences are observed between each individual paralog. PfIRF4B is the main form expressed across the tissues and the most up-regulated one after pathogen induction. It exerts a stronger inhibitory effect on IFN antiviral response than the other two paralogs. PfIRF4A and PfIRF4B-like are primarily present in the nucleus, while PfIRF4B displays colocalization and direct associations with MyD88 in the cytoplasm. Overall, the data demonstrate that three PfIRF4 paralogs show shared and individual functional properties in the negative regulation of type Ⅰ IFN response.

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.

WGDTree: a phylogenetic software tool to examine conditional probabilities of retention following whole genome duplication events

BACKGROUND

Multiple processes impact the probability of retention of individual genes following whole genome duplication (WGD) events. In analyzing two consecutive whole genome duplication events that occurred in the lineage leading to Atlantic salmon, a new phylogenetic statistical analysis was developed to examine the contingency of retention in one event based upon retention in a previous event. This analysis is intended to evaluate mechanisms of duplicate gene retention and to provide software to generate the test statistic for any genome with pairs of WGDs in its history.

RESULTS

Here a software package written in Python, 'WGDTree' for the analysis of duplicate gene retention following whole genome duplication events is presented. Using gene tree-species tree reconciliation to label gene duplicate nodes and differentiate between WGD and SSD duplicates, the tool calculates a statistic based upon the conditional probability of a gene duplicate being retained after a second whole genome duplication dependent upon the retention status after the first event. The package also contains methods for the simulation of gene trees with WGD events. After running simulations, the accuracy of the placement of events has been determined to be high. The conditional probability statistic has been calculated for Phalaenopsis equestris on a monocot species tree with a pair of consecutive WGD events on its lineage, showing the applicability of the method.

CONCLUSIONS

A new software tool has been created for the analysis of duplicate genes in examination of retention mechanisms. The software tool has been made available on the Python package index and the source code can be found on GitHub here: https://github.com/cnickh/wgdtree .

Genome-wide identification, evolutionary analysis, and antimicrobial activity prediction of CC chemokines in allotetraploid common carp, Cyprinus carpio

Chemokines are a group of secreted small molecules which are essential for cell migration in physiological and pathological conditions by binding to specific chemokine receptors. They are structurally classified into five groups, namely CXC, CC, CX3C, XC and CX. CC chemokine group is the largest one among them. In this study, we identified and characterized 61 CC chemokines from allotetraploid common carp (Cyprinus carpio). The sequence analyses showed that the majority of CC chemokines had an N-terminal signal peptide, and an SCY domain, and all CC chemokines were located in the extracellular region. Phylogenetic, evolutionary and syntenic analyses confirmed that CC chemokines were annotated as 11 different types (CCL19, CCL20, CCL25, CCL27, CCL32, CCL33, CCL34, CCL35, CCL36, CCL39, and CCL44), which exhibited unique gene arrangement pattern and chromosomal location respectively. Furthermore, genome synteny analyses between common carp and four representative teleost species indicated expansion of common carp CC chemokines resulted from the whole genome duplication (WGD) event. Additionally, the continuous evolution of gene CCL25s in teleost afforded a novel viewpoint to explain the WGD event in teleost. Then, we predicted the three-dimensional structures and probable function regions of common carp CC chemokines. All the CC chemokines core structures were constituted of an N-loop, a three-stranded β-sheet, and a C-terminal helix. Finally, 43 CC chemokines were predicted to have probable general antimicrobial activity. Their tertiary structures, cationic and amphiphilic physicochemical property supported the viewpoint. To verify the prediction, six recombinant CCL19s proteins were prepared and the antibacterial activity against Escherichia coli and Aeromonas hydrophila were verified. The results supported our prediction that rCCL19a.1s (rCCL19a.1_a, rCCL19a.1_b) and rCCL19bs (rCCL19b_a, rCCL19b_b), especially rCCL19bs, exhibited extremely significant inhibition to the growth of both E. coli and A. hydrophila. On the contrary, two rCCL19a.2s had no significant inhibitory effect. These studies suggested that CC chemokines were essential in immune system evolution and not monofunctional during pathogen infection.