Search

High-altitude environments are inhospitable, but Triplophysa, the largest taxon among the three major fish groups in the Qinghai-Tibetan Plateau (QTP), is an exception. However, the evolutionary profiling of the common ancestor and its contribution to the adaptation of existing QTP native species is unclear.We researched the comparative genomics of Triplophysa species and found that the genome-wide genes of Triplophysa and its ancestry have the characteristics of rapid evolution.Moreover, the rapid evolution of the ancestral genes was caused by relaxed selection. Natural selection analysis showed that more ancestral relaxed selection genes were under strongly purifying selection and showed higher expression in QTP endemic Triplophysa species.The change in natural selection might be associated with the adaptation to QTP. It should be noted that SPT5 homolog, DSIF elongation factor subunit (supt5h) experienced relaxed selection in common ancestral populations of Triplophysa but under purifying selection in extant species, which might be related to hypoxia adaptation of QTP. In summary, the extant species in different environments were used to infer the evolutionary profile of the common ancestor and to identify candidate genes based on changes in natural selection. Our work might provide new clues for understanding adaptation to extreme environments., Competing Interests: Declarations. Competing interests: The authors declare no competing interests., (© 2025. The Author(s).)

Background: Silver arowana (Osteoglossum bicirrhosum) is a basal fish species with sexual monomorphism, while its sex determination mechanism has been poorly understood, posing a significant challenge to its captive breeding efforts., Results: We constructed two high-quality chromosome-level genome assemblies for both female and male silver arowana, with scaffold N50 values over 10 Mb. Combining re-sequencing data of 109 individuals, we identified a female-specific region, which was localized in a non-coding region, i.e., around 26-kb upstream of foxl2 gene (encoding forkhead box L2). Its strong interaction with the neighboring foxl2 on the same chromosome suggests foxl2 as a candidate sex-related gene in silver arowana. We subsequently propose a complex gene network in the sex determination process of silver arowana, with foxl2 acting as the central contributor. Transcriptome sequencing of gonads support our hypothesis that the regulation of foxl2 can be influenced by the spatial proximity of the female-specific fragment, thereby promoting ovarian function or inhibiting testicular function to stimulate gonadal differentiation. Furthermore, we found the sex chromosomes to be homomorphic with a potentially recent origin, as a linkage disequilibrium analysis proved minor recombination suppression., Conclusions: These results taken together serve as a crucial foundation for conducting extensive investigations on the evolution and differentiation of sex-determining mechanisms, as well as the emergence and development of sex chromosomes in various fishes., Competing Interests: Declarations. Ethics approval and consent to participate: Not applicable. Consent for publication: Not applicable. Competing interests: The authors declare no competing interests., (© 2025. The Author(s).)

Background: Spotted scat, a marine aquaculture fish, has variable body color development stages during their ontogenesis. However, the regulatory mechanism of body color patterns formation was poorly understood. Thyroid hormones (TH) function as an important endocrine factor in regulating metamorphosis. In this study, exogenous thyroid hormones 3,5,3'-L-triiodothyronine (T3) and its inhibitor thiourea (TU) were used to treat spotted scat juveniles during the metamorphosis stage (from 60 to 90 dpf). The function and molecular mechanism of thyroid hormone signaling in regulating body color patterns formation was revealed, using the micro-observation of pigments cells distribution, colorimetric evaluation and carotenoids concentration measurement by spectrophotometry, and comparative transcriptome analysis., Results: Spotted scat body color patterns consisted of whole body black color, black bar, black and red spots, and its final pattern was formed through the metamorphosis. When spotted scat were treated with the inhibitor TU to disrupt thyroid hormone signaling, the levels of T3 and T4 were significantly decreased, the melanophores numbers were significantly increased, as well as the expression of genes involved in melanin synthesis and melanophore differentiation (tyr, tyrp1, dct, mitf, pmel, oca2, slc24a5, and erbb3) was significantly increased. Besides, the expression of genes associated with carotenoids and pteridine metabolism (apod, pnpla2, rdh12, stard10, xdh, abca1, retsat, scarb1, rgs2, and gch1) and carotenoids accumulation were stimulated, when thyroid hormone signaling was disrupted by TU. On the contrary, the levels of T3 and T4 were significantly elevated in spotted scat treated with T3, which could weaken the skin redness and reduce the number of black spots and melanophores, as well as the number and diameter of larval erythrophores. Notably, unlike melanophores and erythrophores, the differentiation of iridophore was promoted by thyroid hormones, gene related to iridophore differentiation (fhl2-l, fhl2, ltk, id2a, alx4) and guanine metabolism (gmps, hprt1, ppat, impdh1b) were up-regulated after T3 treatment, but they were down-regulated after TU treatment., Conclusions: Above results showed that thyroid hormone signaling might play critical roles in regulation pigments synthesis and deposition, thereby affecting pigment cells (melanophores, iridophores and erythrophores) formation and body color patterns. The mechanisms of hyperthyroid and hypothyroid on different pigment cells development were different. Excess thyroid hormone might impact the rearrangement of melanophore by regulating cell cycle, resulting in the abnormalities of black spots in spotted scat. Meanwhile, the excessed thyroid hormone could reduce the number and diameter of larval erythrophores, as well as weaken the skin redness of juvenile erythrophores, but they were enhanced by the disruption of thyroid hormone. However, the formation of iridophore differentiation and guanine synthesis genes expression were stimulated by thyroid hormones. These findings provide new insights for exploring the formation of body color patterns in fish, and help to elucidate the molecular mechanism of thyroid hormone in regulating pigment cell development and body coloration, and may also contribute to selective breeding of ornamental fish., Competing Interests: Declarations. Ethics approval and consent to participate: The experiment and animals were approved by the Institutional Animal Care and Use Committee (IA-CUS) of Guangdong Ocean University and complied with China’s laws and regulations on biological research. Competing interests: The authors declare no competing interests., (© 2025. The Author(s).)

Background: The golden pompano (Trachinotus ovatus) is an economically significant warm-water aquaculture species in China. The time required for sexual maturity of T. ovatus is relatively long. Consequently, it has prompted researchers to investigate gonadal development process of this fish. To gain further insight into the function of competing endogenous RNA (ceRNA) in the gonads of T. ovatus and the regulatory mechanism of the ceRNA network, whole transcriptome libraries were constructed from the testes and ovaries., Results: Overall, a total of 96 differentially expressed microRNAs (DE-miRNAs), 2,338 differentially expressed messenger RNAs (DE-mRNAs), 973 differentially expressed long non-coding RNAs (DE-lncRNAs), and 94 differentially expressed circular RNAs (DE-circRNAs) were identified. Additionally, a ceRNA network was constructed, and enrichment analysis confirmed the involvement of numerous pathways in reproduction and gonadal development, including the TGF-β signaling pathway and GnRH signaling pathway. The ceRNA network analysis revealed that the oni-let-7d-1-p3 and PC-3p-112794_13 may play significant roles in T. ovatus gonadal development. And we have observed a possible relationship related to gonadal development involving R-spondin-1 (Rspo1), oni-let-7d-1-p3, and MSTRG.14909.1 (lnc-TGFβR). Dual-luciferase gene reporter system and fluorescence in situ hybridization analyses preliminary verified the regulation relationship between Rspo1 and oni-let-7d-1-p3, as well as lnc-TGFβR and oni-let-7d-1-p3 in the cytoplasm of sertoli cells., Conclusion: It is hypothesized that the lnc-TGFβR functions as a sponge for oni-let-7d-1-p3, participating in regulating the process of testis development. These findings could enhance our understanding of ncRNAs in gonadal development. It also provides new insight into the function of ncRNAs and the regulatory relationship of ceRNA between males and females. These results might contribute to discussions on the regulation of ncRNA during gametogenesis., Competing Interests: Declarations. Ethics approval and consent to participate: All animal care and use procedures were approved by the Committee of the Ethics of Animal Experiments of the Ocean University of China, and were conducted in compliance with the Chinese Guidelines for the Care and Use of Laboratory Animals (GB/T 35892‐2018). We confirm that this study adheres to the ARRIVE guidelines. Consent for publication: Not applicable. Competing interests: The authors declare no competing interests., (© 2025. The Author(s).)

Alanine-rich, alpha-helical type I antifreeze proteins (AFPs) in fishes are thought to have arisen independently in the last 30 Ma on at least four occasions. This hypothesis has recently been proven for flounder and sculpin AFPs, which both originated by gene duplication and divergence followed by substantial gene copy number expansion. Here, we examined the origins of the cunner (wrasse) and snailfish (liparid) AFPs. The cunner AFP has arisen by a similar route from the duplication and divergence of a GIMAP gene. The coding region for this AFP stems from an alanine-rich region flanking the GTPase domain of GIMAPa. The AFP gene has remained in the GIMAP gene locus and has undergone amplification there along with the GIMAPa gene. The AFP gene originated after the cunner diverged from its common ancestor with the closely related spotty and ballan wrasses, which exhibit similar gene synteny but lack AFP genes. Snailfish AFPs have also recently evolved because they are confined to a single genus of this family. In these AFP-producing species, the AFP locus does not share any similarity to functional genes. Instead, it is replete with repetitive DNAs and transposons, several stretches of which could encode alanine tracts with a dominant codon (GCC) that matches the bias observed in the AFP genes. All four known instances of type I AFPs occurring in fishes are independent evolutionary events that occurred soon after the onset of Northern Hemisphere Cenozoic glaciation events. Collectively, these results provide a remarkable example of convergent evolution to one AFP type., Competing Interests: Declarations. Ethics approval and consent to participate: Not applicable. Consent for publication: Not applicable. Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

Background: Animal venom systems are considered as valuable model for investigating the molecular mechanisms underlying phenotypic evolution. Stonefish are the most venomous and dangerous fish because of severe human envenomation and occasionally fatalities, whereas the genomic background of their venom has not been fully explored compared with that in other venomous animals., Results: In this study, we followed modern venomic pipelines to decode the Synanceia verrucosa venom components. A catalog of 478 toxin genes was annotated based on our assembled chromosome-level genome. Integrative analysis of the high-quality genome, the transcriptome of the venom gland, and the proteome of crude venom revealed mechanisms underlying the venom complexity in S. verrucosa. Six tandem-duplicated neoVTX subunit genes were identified as the major source for the neoVTX protein production. Further isoform sequencing revealed massive alternative splicing events with a total of 411 isoforms demonstrated by the six genes, which further contributed to the venom diversity. We then characterized 12 dominantly expressed toxin genes in the venom gland, and 11 of which were evidenced to produce the venom protein components, with the neoVTX proteins as the most abundant. Other major venom proteins included a presumed CRVP, Kuntiz-type serine protease inhibitor, calglandulin protein, and hyaluronidase. Besides, a few of highly abundant non-toxin proteins were also characterized and they were hypothesized to function in housekeeping or hemostasis maintaining roles in the venom gland. Notably, gastrotropin like non-toxin proteins were the second highest abundant proteins in the venom, which have not been reported in other venomous animals and contribute to the unique venom properties of S. verrucosa., Conclusions: The results identified the major venom composition of S. verrucosa, and highlighted the contribution of neoVTX genes to the diversity of venom composition through tandem-duplication and alternative splicing. The diverse neoVTX proteins in the venom as lethal particles are important for understanding the adaptive evolution of S. verrucosa. Further functional studies are encouraged to exploit the venom components of S. verrucosa for pharmaceutical innovation., Competing Interests: Declarations. Ethics approval and consent to participate: All the experiments were conducted in accordance with guidelines approved by the Animal Ethical and Welfare Committee of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou). Consent for publication: Not applicable. Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

Background: Primary transcripts are largely comprised of intronic sequences that are excised and discarded shortly after synthesis. In vertebrates, the shape of the intron size distribution is largely constant; however, most teleost fish have a diverged log-bimodal 'teleost distribution' (TD) that is seen only in teleosts. How the TD evolved and to what extent this was affected by adaptative or non-adaptive mechanisms is unknown., Results: Here, we show that the TD has evolved independently at least six times and that its appearance is linked to the loss of the aKRAB domain from PRDM9. We determined intron size distributions and identified PRDM9 orthologues from annotated genomes in addition to scanning 1193 teleost assemblies for the aKRAB domain. We show that a diverged form of PRDM9 ( β ) is predominant in teleosts whereas the α version is absent from most species. Only a subset of PRDM9- α proteins contain aKRAB, and hence, it is present only in a small number of teleost lineages. Almost all lineages lacking aKRAB (but no species with) had TDs., Conclusions: In mammals, PRDM9 defines the sites of meiotic recombination through a mechanism that increases structural variance and depends on aKRAB. The loss of aKRAB is likely to have shifted the locations of both recombination and structural variance hotspots. Our observations suggest that the TD evolved as a side-effect of these changes and link recombination to the evolution of intron size illustrating how genome architectures can evolve in the absence of selection., Competing Interests: Declarations. Ethics approval and consent to participate: Not applicable. Consent for publication: Not applicable. Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

Background: The short-snout icefish, Neosalanx brevirostris, a member of the Salangidae family, is an economically important fishery species in China. Understanding the mechanisms underlying sex determination in this species has crucial implications for conservation, ecology and evolution. Meanwhile, there is a shortage of rapid and cost-effective genetic methods for sex identification, which poses challenges in identifying the sex of immature individuals in sex determination mechanism studies and aquaculture breeding applications., Results: Based on whole genome resequencing data, sex-specific loci and regions were found to be concentrated in a region on chromosome 2. All sex-specific loci exhibited excess heterozygosity in females and complete homozygosity in males. This sex determining region contains seven genes, including cytochrome P450 aromatase CYP19B, which is involved in steroidogenesis and is associated with 24 sex-specific loci and two W-deletions. A haploid female-specific sequence was identified as paralogous to a diploid sequence with a significant length difference, making it suitable for rapid and cost-effective genetic sex identification by traditional PCR and agarose gel electrophoresis, which were further validated in 24 females and 24 males with known phenotypic sexes., Conclusions: Our results confirm that N. brevirostris exhibits a female heterogametic sex determination system (ZZ/ZW), with chromosome 2 identified as the putative sex chromosome containing a relatively small sex determining region (~ 48 Kb). The gene CYP19B is proposed as a candidate sex determining gene. Moreover, the development of PCR based method enables genetic sex identification at any developmental stage, thereby facilitating further studies on sex determination mechanisms and advancing aquaculture breeding applications for this species., Competing Interests: Declarations. Ethics approval and consent to participate: The sampling location is not privately owned or protected, and no specific permission is required for sample collection. All animal experimental procedures in this research were approved by the Animal Care and Use Committee of Institute of Oceanology, Chinese Academy of Sciences. The methods applied in this study were carried out in strict accordance with the Laboratory Animal Management Principles of China. Consent for publication: Not applicable. Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

Background: The white-blooded Antarctic icefishes is a representative organism that survive under the stenothermal conditions of the Southern Ocean without the hemoglobin genes. To compensate for inefficient oxygen transport, distinct features such as increased heart size, greater blood volume, and reduced hematocrit density enhance the amount of dissolved oxygen and the velocity of blood flow., Results: Here, we investigated these unique characteristics by comparing high-quality genomic data between white-blooded and red-blooded fishes and identified the loss of FAAP20, which is implicated in anemia. Although the gene region containing FAAP20 is conserved in notothenioids as shown through collinear analysis, only remnants of FAAP20 persist in several icefish species. Additionally, we observed the loss of SOAT1, which plays a pivotal role in cholesterol metabolism, providing a clue for further investigations into the unique mitochondrial form of the icefish., Conclusions: The loss of FAAP20, which is known to reduce erythrocyte counts under stress conditions in mice and humans, may provide a clue to understanding the genomic characteristics related to oxygen supply, such as low hematocrit, in Antarctic icefishes., Competing Interests: Declarations Ethics approval and consent to participate Not applicable. Consent for publication Not applicable. Competing interests The authors declare no competing interests., (© 2024. The Author(s).)

Fish early-stages constitute useful indicators of the states of marine ecosystems, as well as important fishery resources. Given the spectacular phenotypic changes during ontogeny, and the paucity of diagnostic morphological characters at the species level, the identification of fish early-stages is a challenging task. DNA barcoding, the use of the mitochondrial gene of the cytochrome c oxidase subunit I (COI) as an internal species tag, opened new perspectives for the identifications of both larval fish and fish eggs. However, the accuracy of the identifications assisted by DNA barcoding are dependent of the completeness of the DNA barcode reference libraries used to assigned unknown sequences to known species. Here, we built a DNA barcode reference library for 113 species of larval fish and 85 species of fish eggs involving the production of 741 newly generated DNA barcodes from South China Sea (63 localities). Together with 514 DNA barcodes mined from Genbank for 116 species from the South China Sea regions, a reference library including 1255 DNA barcodes for 308 species (248 locations) was assembled. The present study emphasizes the importance of integrating DNA barcoding to large scale inventories of early stages, as DNA-based species delimitation analyses delimited 305 molecular operational taxonomic units (MOTUs) and multiple cases of discordance with morphological identifications were detected. Cryptic diversity is detected with 14 species displaying two MOTUs and a total of 23 species were lumped into 11 MOTUs due to low interspecific divergence and/or mixed lineages., (© 2024. The Author(s).)

Background: Hydrogen peroxide (H 2 O 2 ), a novel water treatment agent, can be used for disinfection, water quality adjustment, and disease prevention, while excessive H 2 O 2 can injure farm animals, even leading to death. Hydrogen peroxide is a recommended disinfectant and bactericide for treating gill diseases and vibriosis in the greenfin horse-faced filefish Thamnaconus septentrionalis. However, its cumulative effect, toxic molecular mechanism and relevant signal transduction/metabolic networks in marine fishes are largely unknown., Results: We employed a multi-omics approach to investigate the detrimental effects of 50 mg/L H 2 O 2 exposure (2 h/d) on filefish for 2 d, 4 d, and 6 d. Transcriptome sequencing showed that differentially expressed genes (DEGs) were mainly classified into functions such as signal transduction, nervous system, liver and bile acid metabolism, energy metabolism, cell adhesion and communication, inflammation and immune response. Metabolomic analysis found that the significantly changed metabolites (SCMs) were involved in phenylalanine metabolism, inflammatory mediator regulation, linoleic acid metabolism, and necroptosis. The main SCMs were cholic acid, carnitine C12:1, dimethylmalonic acid, glutamic acid, L-lactic acid, shikimic acid, 2-methylsuccinic acid, and others. Moreover, H 2 O 2 -induced oxidative stress also disturbs the balance of the gut microbiota, altering the microbial composition and affecting digestive processes., Conclusions: Integrated multiomics analysis revealed that H 2 O 2 -induced detrimental impacts include mucosal damage, inflammatory and immune responses, altered energy metabolism, and gut microbiota disorders. These findings offer novel insights into the harmful effects and signal transduction/metabolic pathways triggered by H 2 O 2 exposure in marine fishes., (© 2024. The Author(s).)

Objective: Environmental DNA (eDNA) methods are crucial for monitoring populations, particularly rare and cryptic species. For confident eDNA application, rigorous assay validation is required including specificity testing with genomic DNA (gDNA). However, this critical step is often difficult to achieve as obtaining fresh tissue samples from at-risk species can be difficult, highly limited, or impossible. Natural history museum collections could serve as a valuable and ethical voucher specimen resource for eDNA assay validation. The present study demonstrates the effectiveness of whole genome amplification (WGA) in providing enough gDNA to assemble high quality mitogenomes from which robust targeted eDNA assays can be designed., Results: Using fresh and historical museum tissue samples from six species spanning fish, birds, and mammals, we successfully developed a WGA method with an average yield of 380 to 1,268 ng gDNA per 20 µL reaction. This gDNA was used for whole genome shotgun sequencing and subsequent assembly of high quality mitogenomes using mtGrasp. These mitogenomes were then used to develop six new robust, targeted quantitative real time polymerase chain reaction-based eDNA assays and 200 ng WGA-enriched yielded satisfactory C q values and near 100% detection frequencies for all assays tested. This approach offers a cost-effective and non-invasive alternative, streamlining eDNA research processes and aiding in conservation efforts., (© 2024. The Author(s).)

Background: Blotched snakehead (Channa maculata) displays significant sexual dimorphism, with males exhibiting faster growth rates and larger body sizes compared to females. The cultivation of the all-male population of snakeheads holds substantial economic and ecological value. Nonetheless, the intricate processes governing the development of bipotential gonads into either testis or ovary in C. maculata remain inadequately elucidated. Therefore, it is necessary to determine the critical time window of sex differentiation in C. maculata, providing a theoretical basis for sex control in production practices., Methods: The body length and weight of male and female C. maculata were measured at different developmental stages to reveal when sexual dimorphism in growth initially appears. Histological observations and spatiotemporal comparative transcriptome analyses were performed on ovaries and testes across various developmental stages to determine the crucial time windows for sex differentiation in each sex and the sex-related genes. Additionally, qPCR and MG2C were utilized to validate and locate sex-related genes, and levels of E 2 and T were quantified to understand sex steroid synthesis., Results: Sexual dimorphism in growth became evident starting from 90 dpf. Histological observations revealed that morphological sex differentiation in females and males occurred between 20 and 25 dpf or earlier and 30-35 dpf or earlier, respectively, corresponding to the appearance of the ovarian cavity or efferent duct anlage. Transcriptome analyses revealed divergent gene expression patterns in testes and ovaries after 30 dpf. The periods of 40-60 dpf and 60-90 dpf marked the initiation of molecular sex differentiation in females and males, respectively. Male-biased genes (Sox11a, Dmrt1, Amh, Amhr2, Gsdf, Ar, Cyp17a2) likely play crucial roles in male sex differentiation and spermatogenesis, while female-biased genes (Foxl2, Cyp19a1a, Bmp15, Figla, Er) could be pivotal in ovarian differentiation and development. Numerous biological pathways linked to sex differentiation and gametogenesis were also identified. Additionally, E 2 and T exhibited sexual dimorphism during sex differentiation and gonadal development. Based on these results, it is hypothesized that in C. maculata, the potential male sex differentiation pathway, Sox11a-Dmrt1-Sox9b, activates downstream sex-related genes (Amh, Amhr2, Gsdf, Ar, Cyp17a2) for testicular development, while the antagonistic pathway, Foxl2/Cyp19a1a, activates downstream sex-related genes (Bmp15, Figla, Er) for ovarian development., Conclusions: This study provides a comprehensive overview of gonadal dynamic changes during sex differentiation and gametogenesis in C. maculata, establishing a scientific foundation for sex control in this species., (© 2024. The Author(s).)

Resveratrol has been reported to promote immunity and decrease oxidative stress, but which demonstrates biphasic effects relied on the use concentration. In this study, the effects of diet supplement with a relative high concentration of resveratrol (0.32 mg/kg) on metabolism, antioxidation and apoptosis of liver were investigated in Siberian sturgeon. The results showed that resveratrol significantly increased the lipid synthesis and the apoptosis, but did not either activate the antioxidant NRF2/KEAP1 pathway or enhance the antioxidant enzyme activity. Transcriptome analysis revealed significant changes in regulatory pathways related to glycolipid, including PPAR signaling pathway, Insulin signaling pathway, Fatty acid biosynthesis, and Glycolysis/Gluconeogenesis. In addition, resveratrol significantly increased the lipid synthesis genes (accα and fas), fatty acid transport gene (fatp 6) and gluconeogenesis gene (gck), but decreased the survival-promoting genes (gadd45β and igf 1). These findings highlight a significant effect of resveratrol on glycolipid metabolism in Siberian sturgeon. Moreover, this study also demonstrated that 0.32 mg/kg resveratrol has physiological toxicity to the liver of Siberian sturgeon, indicating that this dose is too high for Siberian sturgeon. Thus, our study provides a valuable insight for future research and application of resveratrol in fish., (© 2024. The Author(s).)

Background: Daphnia galeata is a suitable model organism for investigating predator-induced defense. Genes and pathways exhibiting differential expression between fish kairomone-treated and untreated groups in D. galeata have been identified. However, understanding of the significance of alternative splicing, a crucial process of the regulation of gene expression in eukaryotes, to this mechanism remains limited. This study measured life-history traits and conducted short-read RNA sequencing and long-read isoform sequencing of two Korean D. galeata genotypes (KB1 and KE2) to uncover the genetic mechanism underlying their phenotypic plasticity under predation stress., Results: KB1 exhibited strategies to enhance fertility and decrease body length when exposed to fish kairomones, while KE2 deployed an adaptive strategy to increase body length. Full-length transcriptomes from KB1 and KE2 yielded 65,736 and 57,437 transcripts, respectively, of which 32 differentially expressed transcripts (DETs) were shared under predation stress across both genotypes. Prominent DETs common to both genotypes were related to energy metabolism and the immune system. Additionally, differential alternative splicing (DAS) events were detected in both genotypes in response to fish kairomones. DAS genes shared between both genotypes may indicate their significant role in the post-transcriptional stress response to fish predation. Calpain-3, involved in digestion and nutrient absorption, was identified as a DAS gene in both genotypes when exposed to fish kairomones. In addition, the gene encoding thymosin beta, which is related to growth, was found to be a statistically significant DAS only in KB1, while that encoding ultraspiracle protein, also associated with growth, was only identified in KE2. Moreover, transcripts encoding proteins such as EGF-like domain-containing protein, vitellogenin fused with superoxide dismutase, and others were identified overlapping between DAS events and DETs and potentially elucidating their association with the observed phenotypic variation in each genotype., Conclusions: Our findings highlight the crucial role of alternative splicing in modulating transcriptome landscape under predation stress in D. galeata, emphasizing the requirement for integrating gene expression and splicing analyses in evolutionary adaptation studies., (© 2024. The Author(s).)

Background: In aquaculture, sturgeons are generally maintained in the confined spaces, which not only hinders sturgeon movement, but also threatens their flesh quality that seriously concerned by aquaculture industry. As a typical antioxidant, resveratrol can improve the flesh quality of livestock and poultry. However, the mechanism of resveratrol's effect on the muscle of Siberian sturgeon is still unclear., Results: In this study, the dietary resveratrol increased the myofiber diameter, the content of the amino acids, antioxidant capacity markers (CAT, LDH and SOD) levels and the expression levels of mTORC1 and MYH9 in muscle of Siberian sturgeon. Further transcriptome analysis displayed that ROS production-related pathways ("Oxidative phosphorylation" and "Chemical carcinogenes-reactive oxygen species") were enriched in KEGG analysis, and the expression levels of genes related to the production of ROS (COX4, COX6A, ATPeF1A, etc.) in mitochondria were significantly down-regulated, while the expression levels of genes related to scavenging ROS (SOD1) were up-regulated., Conclusions: In summary, this study reveals that resveratrol may promote the flesh quality of Siberian sturgeon probably by enhancing myofiber growth, nutritional value and the antioxidant capacity of muscle, which has certain reference significance for the development of a new type of feed for Siberian sturgeon., (© 2024. The Author(s).)

Background: Genome-wide comparisons of populations are widely used to explore the patterns of nucleotide diversity and sequence divergence to provide knowledge on how natural selection and genetic drift affect the genome. In this study we have compared whole-genome sequencing data from Atlantic and Pacific herring, two sister species that diverged about 2 million years ago, to explore the pattern of genetic differentiation between the two species., Results: The genome comparison of the two species revealed high genome-wide differentiation but with islands of remarkably low genetic differentiation, as measured by an F ST analysis. However, the low F ST observed in these islands is not caused by low interspecies sequence divergence (d xy ) but rather by exceptionally high estimated intraspecies nucleotide diversity (π). These regions of low differentiation and elevated nucleotide diversity, termed high-diversity regions in this study, are not enriched for repeats but are highly enriched for immune-related genes. This enrichment includes genes from both the adaptive immune system, such as immunoglobulin, T-cell receptor and major histocompatibility complex genes, as well as a substantial number of genes with a role in the innate immune system, e.g. novel immune-type receptor, tripartite motif and tumor necrosis factor receptor genes. Analysis of long-read based assemblies from two Atlantic herring individuals revealed extensive copy number variation in these genomic regions, indicating that the elevated intraspecies nucleotide diversities were partially due to the cross-mapping of short reads., Conclusions: This study demonstrates that copy number variation is a characteristic feature of immune trait loci in herring. Another important implication is that these loci are blind spots in classical genome-wide screens for genetic differentiation using short-read data, not only in herring, likely also in other species harboring qualitatively similar variation at immune trait loci. These loci stood out in this study because of the relatively high genome-wide baseline for F ST values between Atlantic and Pacific herring., (© 2024. The Author(s).)

Background: Phylogenetic gaps of public databases of reference sequences are a major obstacle for comparative genomics and management of marine resources, particularly in the Global South, where economically important fisheries and conservation flagship species often lack closely-related references. We applied target-enrichment to obtain complete mitochondrial genomes of marine ichthyofauna from the Brazilian coast selected based on economic significance, conservation status and lack of phylogenetically-close references. These included sardines (Dorosomatidae, Alosidae), mackerels (Scombridae) croakers (Sciaenidae), groupers (Epinephelidae) and snappers (Lutjanidae)., Results: Custom baits were designed to enrich mitochondrial DNA across a broad phylogenetic range of fishes. Sequencing generated approximately 100k reads per sample, which were assembled in a total of 70 complete mitochondrial genomes and include fifty-two new additions to GenBank, including five species with no previous mitochondrial data. Departures from the typical gene content and order occurred in only three taxa and mostly involved tRNA gene duplications. Start-codons for all genes, except Cytochrome C Oxidase subunit I (COI), were consistently ATG, whilst a wide range of stop-codons deviated from the prevailing TAA. Phylogenetic analysis confirmed assembly accuracy and revealed signs of cryptic diversification within the Mullus genus. Lineage delimitation methods using Sardinella aurita and S. brasiliensis mitochondrial genomes support a single Operational Taxonomic Unit., Conclusions: Target enrichment was highly efficient, providing complete novel mitochondrial genomes with little sequencing effort. These sequences are deposited in public databases to enable subsequent studies in population genetics and adaptation of Latin American fish species and serve as a vital resource for conservation and management programs that rely on molecular data for species and genus-level identification., (© 2024. The Author(s).)

Habitat niches of fish species can exert a strong influence on population structure, even on a small geographical scale. In this scope, Pelasgus thesproticus is a great model species to study connectivity in riverine environments owing to its naturally patchy habitat distribution. Furthermore, it is important to conduct such studies in near-natural systems to avoid the impact of human disturbances on the river, such as fragmentation, morphological changes and habitat degradation. In this sense, the Vjosa in Albania is an excellent study area. A total of 204 individuals were sampled from five locations in the lower Vjosa and two tributaries and genotyped with 33 newly designed microsatellites loci using high throughput sequencing. The application of microsatellite genotyping by sequencing revealed genetic structure and some differentiation, even at a small spatial scale (< 65 river km). A total of 500 alleles were found with an average of 0.93 private alleles among sites with rather low F ST values (< 0.04). The extent of admixture observed in some populations indicate that the genetic structure is mainly influenced by upstream populations, either from the main river itself or from tributaries. In addition, the connection between a tributary and the other sites is disrupted by the flow regime, which is reflected in a high degree of divergence from the other populations. Our results indicate that hydrological conditions of the flowing river present strong barriers to gene flow, particularly in the upstream direction, but at the same time act as dispersal corridors in the downstream direction and exhibit source-sink dynamics in which upstream populations contribute disproportionately to downstream populations for this habitat specialist along the river. It is suggested that processes of colonization and reinforcement may play an important role in shaping the genetic structure of patchily distributed fish species in natural river systems. Future studies should increase the knowledge of dispersal factors, habitat heterogeneity, consequence of source-sink dynamics, and gene flow within the system, which will help to understand and maintain important processes related to metapopulation theory and the potential evolutionary consequences of habitat loss and fragmentation., (© 2023. The Author(s).)

Background: Streptococcus iniae is an important fish pathogen that cause significant economic losses to the global aquaculture industry every year. Although there have some reports on the genotype of S.iniae and its relationship with virulence, no genome-scale comparative analysis has been performed so far. In our previous work, we characterized 17 isolates of S.iniae from Trachinotus ovatus and divided them into two genotypes using RAPD and rep-PCR methods. Among them, BH15-2 was classified as designated genotype A (in RAPD) and genotype 1 (in rep-PCR), while BH16-24 was classified as genotype B and genotype 2. Herein, we compared the differences in growth, drug resistance, virulence, and genome between BH15-2 and BH16-24., Results: The results showed that the growth ability of BH16-24 was significantly faster than that of BH15-2 at the exponential stage. Antimicrobial tests revealed that BH15-2 was susceptible to most of the tested antibiotics except neomycin and gentamycin. In contrast, BH16-24 was resistant to 7 antibiotics including penicillin, sulfasomizole, compound sulfamethoxazole tablets, polymyxin B, spectinomycin, rifampin and ceftazidime. Intraperitoneal challenge of T.ovatus, showed that the LD 50 value of BH15-2 was 4.0 × 10 2 CFU/g, while that of BH16-24 was 1.2 × 10 5 CFU/g. The genome of S.iniae BH15-2 was 2,175,659 bp with a GC content of 36.80%. Meanwhile, the genome of BH16-24 was 2,153,918 bp with a GC content of 36.83%. Comparative genome analysis indicated that compared with BH15-2, BH16-24 genome had a large-scale genomic inversion fragment, at the location from 502,513 bp to 1,788,813 bp, resulting in many of virulence and resistance genes differentially expression. In addition, there was a 46 kb length, intact phage sequence in BH15-2 genome, which was absent in BH16-24., Conclusion: Comparative genomic studies of BH15-2 and BH16-24 showed that the main difference is a 1.28 Mbp inversion fragment. The inversion fragment may lead to abnormal expression of drug resistant and virulence genes, which is believed to be the main reason for the multiple resistance and weakened virulence of BH16-24. Our study revealed the potential mechanisms in underlying the differences of multidrug resistance and virulence among different genotypes of S.iniae., (© 2023. The Author(s).)

Background: Biodiversity exists at different levels of organisation: e.g. genetic, individual, population, species, and community. These levels of organisation all exist within the same system, with diversity patterns emerging across organisational scales through several key processes. Despite this inherent interconnectivity, observational studies reveal that diversity patterns across levels are not consistent and the underlying mechanisms for variable continuity in diversity across levels remain elusive. To investigate these mechanisms, we apply a spatially explicit simulation model to simulate the global diversification of tropical reef fishes at both the population and species levels through emergent population-level processes., Results: We find significant relationships between the population and species levels of diversity which vary depending on both the measure of diversity and the spatial partitioning considered. In turn, these population-species relationships are driven by modelled biological trait parameters, especially the divergence threshold at which populations speciate., Conclusions: To explain variation in multi-level diversity patterns, we propose a simple, yet novel, population-to-species diversity partitioning mechanism through speciation which disrupts continuous diversity patterns across organisational levels. We expect that in real-world systems this mechanism is driven by the molecular dynamics that determine genetic incompatibility, and therefore reproductive isolation between individuals. We put forward a framework in which the mechanisms underlying patterns of diversity across organisational levels are universal, and through this show how variable patterns of diversity can emerge through organisational scale., (© 2023. The Author(s).)

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., (© 2023. The Author(s).)

Background: Phenotypic plasticity is a crucial adaptive mechanism that enables organisms to modify their traits in response to changes in their environment. Predator-induced defenses are an example of phenotypic plasticity observed across a wide range of organisms, from single-celled organisms to vertebrates. In addition to morphology and behavior, these responses also affect life-history traits. The crustacean Daphnia galeata is a suitable model organism for studying predator-induced defenses, as it exhibits life-history traits changes under predation risk. To get a better overview of their phenotypic plasticity under predation stress, we conducted RNA sequencing on the transcriptomes of two Korean Daphnia galeata genotypes, KE1, and KB11, collected in the same environment., Results: When exposed to fish kairomones, the two genotypes exhibited phenotypic variations related to reproduction and growth, with opposite patterns in growth-related phenotypic variation. From both genotypes, a total of 135,611 unigenes were analyzed, of which 194 differentially expressed transcripts (DETs) were shared among the two genotypes under predation stress, which showed consistent, or inconsistent expression patterns in both genotypes. Prominent DETs were related to digestion and reproduction and consistently up-regulated in both genotypes, thus associated with changes in life-history traits. Among the inconsistent DETs, transcripts encode vinculin (VINC) and protein obstructor-E (OBST-E), which are associated with growth; these may explain the differences in life-history traits between the two genotypes. In addition, genotype-specific DETs could explain the variation in growth-related life-history traits between genotypes, and could be associated with the increased body length of genotype KE1., Conclusions: The current study allows for a better understanding of the adaptation mechanisms related to reproduction and growth of two Korean D. galeata genotypes induced by predation stress. However, further research is necessary to better understand the specific mechanisms by which the uncovered DETs are related with the observed phenotypic variation in each genotype. In the future, we aim to unravel the precise adaptive mechanisms underlying predator-induced responses., (© 2023. BioMed Central Ltd., part of Springer Nature.)

Background: The sturgeon group has been economically significant worldwide due to caviar production. Sturgeons consist of 27 species in the world. Mitogenome data could be used to infer genetic diversity and investigate the evolutionary history of sturgeons. A limited number of complete mitogenomes in this family were sequenced. Here, we annotated the mitochondrial Huso huso genome, which revealed new aspects of this species., Results: In this species, the mitochondrial genome consisted of 13 genes encoding proteins, 22tRNA and 2rRNA, and two non-coding regions that followed other vertebrates. In addition, H. huso had a pseudo-tRNA-Glu between ND6 and Cytb and a 52-nucleotide tandem repeat with two replications in 12S rRNA. This duplication event is probably related to the slipped strand during replication, which could remain in the strand due to mispairing during replication. Furthermore, an 82 bp repeat sequence with three replications was observed in the D-loop control region, which is usually visible in different species. Regulatory elements were also seen in the control region of the mitochondrial genome, which included termination sequences and conserved regulatory blocks. Genomic compounds showed the highest conservation in rRNA and tRNA, while protein-encoded genes and nonencoded regions had the highest divergence. The mitochondrial genome was phylogenetically assayed using 12 protein-encoding genes., Conclusions: In H. huso sequencing, we identified a distinct genome organization relative to other species that have never been reported. In recent years, along with the advancement in sequencing identified more genome rearrangements. However, it is an essential aspect of researching the evolution of the mitochondrial genome that needs to be recognized., (© 2023. BioMed Central Ltd., part of Springer Nature.)

Background: Genetic diversity and heterogeneous genomic signatures in marine fish populations may result from selection pressures driven by the strong effects of environmental change. Nearshore fishes are often exposed to complex environments and human activities, especially those with small ranges. However, studies on genetic diversity and population selection signals in these species have mostly been based on a relatively small number of genetic markers. As a newly recorded species of Sillaginidae, the population genetics and genomic selection signals of Sillago sinica are fragmented or even absent., Results: To address this theoretical gap, we performed whole-genome resequencing of 43 S. sinica individuals from Dongying (DY), Qingdao (QD) and Wenzhou (WZ) populations and obtained 4,878,771 high-quality SNPs. Population genetic analysis showed that the genetic diversity of S. sinica populations was low, but the genetic diversity of the WZ population was higher than that of the other two populations. Interestingly, the three populations were not strictly clustered within the group defined by their sampling location but showed an obvious geographic structure signal from the warm temperate to the subtropics. With further analysis, warm-temperate populations exhibited strong selection signals in genomic regions related to nervous system development, sensory function and immune function. However, subtropical populations showed more selective signalling for environmental tolerance and stress signal transduction., Conclusions: Genome-wide SNPs provide high-quality data to support genetic studies and localization of selection signals in S. sinica populations. The reduction in genetic diversity may be related to the bottleneck effect. Considering that low genetic diversity leads to reduced environmental adaptability, conservation efforts and genetic diversity monitoring of this species should be increased in the future. Differences in genomic selection signals between warm temperate and subtropical populations may be related to human activities and changes in environmental complexity. This study deepened the understanding of population genetics and genomic selection signatures in nearshore fishes and provided a theoretical basis for exploring the potential mechanisms of genomic variation in marine fishes driven by environmental selection pressures., (© 2023. BioMed Central Ltd., part of Springer Nature.)

Background: Smell abilities differ greatly among vertebrate species due to distinct sensory needs, with exceptional variability reported in the number of olfactory genes and the size of the odour-processing regions of the brain. However, key environmental factors shaping genomic and phenotypic changes linked to the olfactory system remain difficult to identify at macroevolutionary scales. Here, we investigate the association between diverse ecological traits and the number of olfactory chemoreceptors in approximately two hundred ray-finned fishes., Results: We found independent expansions producing large gene repertoires in several lineages of nocturnal amphibious fishes, generally able to perform active terrestrial exploration. We reinforced this finding with on-purpose genomic and transcriptomic analysis of Channallabes apus, a catfish species from a clade with chemosensory-based aerial orientation. Furthermore, we also detected an augmented information-processing capacity in the olfactory bulb of nocturnal amphibious fishes by estimating the number of cells contained in this brain region in twenty-four actinopterygian species., Conclusions: Overall, we report a convergent genomic and phenotypic magnification of the olfactory system in nocturnal amphibious fishes. This finding suggests the possibility of an analogous evolutionary event in fish-like tetrapod ancestors during the first steps of the water-to-land transition, favouring terrestrial adaptation through enhanced aerial orientation., (© 2023. The Author(s).)

Background: Teleost fishes comprise more than half of the vertebrate species. Within teleosts, most phylogenies consider the split between Osteoglossomorpha and Euteleosteomorpha/Otomorpha as basal, preceded only by the derivation of the most primitive group of teleosts, the Elopomorpha. While Osteoglossomorpha are generally species poor, the taxon contains the African weakly electric fish (Mormyroidei), which have radiated into numerous species. Within the mormyrids, the genus Campylomormyrus is mostly endemic to the Congo Basin. Campylomormyrus serves as a model to understand mechanisms of adaptive radiation and ecological speciation, especially with regard to its highly diverse species-specific electric organ discharges (EOD). Currently, there are few well-annotated genomes available for electric fish in general and mormyrids in particular. Our study aims at producing a high-quality genome assembly and to use this to examine genome evolution in relation to other teleosts. This will facilitate further understanding of the evolution of the osteoglossomorpha fish in general and of electric fish in particular., Results: A high-quality weakly electric fish (C. compressirostris) genome was produced from a single individual with a genome size of 862 Mb, consisting of 1,497 contigs with an N50 of 1,399 kb and a GC-content of 43.69%. Gene predictions identified 34,492 protein-coding genes, which is a higher number than in the two other available Osteoglossomorpha genomes of Paramormyrops kingsleyae and Scleropages formosus. A Computational Analysis of gene Family Evolution (CAFE5) comparing 33 teleost fish genomes suggests an overall faster gene family turnover rate in Osteoglossomorpha than in Otomorpha and Euteleosteomorpha. Moreover, the ratios of expanded/contracted gene family numbers in Osteoglossomorpha are significantly higher than in the other two taxa, except for species that had undergone an additional genome duplication (Cyprinus carpio and Oncorhynchus mykiss). As potassium channel proteins are hypothesized to play a key role in EOD diversity among species, we put a special focus on them, and manually curated 16 Kv1 genes. We identified a tandem duplication in the KCNA7a gene in the genome of C. compressirostris., Conclusions: We present the fourth genome of an electric fish and the third well-annotated genome for Osteoglossomorpha, enabling us to compare gene family evolution among major teleost lineages. Osteoglossomorpha appear to exhibit rapid gene family evolution, with more gene family expansions than contractions. The curated Kv1 gene family showed seven gene clusters, which is more than in other analyzed fish genomes outside Osteoglossomorpha. The KCNA7a, encoding for a potassium channel central for EOD production and modulation, is tandemly duplicated which may related to the diverse EOD observed among Campylomormyrus species., (© 2023. The Author(s).)

Biodiversity assessment approaches based on molecular biology techniques such as metabarcoding, RAD-seq, or SnaPshot sequencing are increasingly applied in assessing marine and aquatic ecosystems. Here we present a new reference database for fish meta-barcoding based on mitochondrial genes. The Mare-MAGE database contains quality-checked sequences of the mitochondrial 12S ribosomal RNA and Cytochrome c Oxidase I gene. All sequences were obtained from the National Center for Biotechnology Information- GenBank (NBCI-GenBank), the European Nucleotide Archive (ENA), AquaGene Database and BOLD database, and have undergone intensive processing. They were checked for false annotations and non-target anomalies, according to the Integrated Taxonomic Information System (ITIS) and FishBase. The dataset is compiled in ARB-Home, FASTA and Qiime2 formats, and is publicly available from the Mare-MAGE database website ( http://mare-mage.weebly.com/ ). It includes altogether 231,333 COI and 12S rRNA gene sequences of fish, covering 19,506 species of 4,058 genera and 586 families., (© 2023. The Author(s).)

An important food-producing sector in Egypt is aquaculture and fisheries; however, several pathogenic microorganisms lead to high mortalities and significant economic losses. The occurrence of Psychrobacter glacincola infection among 180 wild marine fishes collected from the Red sea at Hurghada, Egypt were investigated in the present study. The disease prevalence rate was 6.7%. The recovered isolates were subjected to biochemical and molecular identification. The study also investigated pathogenicity and the antibiogram profile of the recovered isolates. The clinical examination of the infected fish revealed various signs that included lethargy and sluggish movement, hemorrhages and ulcers on the body and the operculum, scale loss, and fin congestion and rot, especially at the tail fin. Furthermore, during postmortem examination, congestion of the liver, spleen, and kidney was observed. Interestingly, 12 isolates were recovered and were homogenous bacteriologically and biochemically. The phylogenetic analysis based on 16S rRNA gene confirmed that MRB62 identified strain was closely related the genus Psychrobacter and identified as P. glacincola and was pathogenic to Rhabdosargus haffara fish, causing 23.3% mortality combined with reporting a series of clinical signs similar to that found in naturally infected fishes. The present study also showed that P. glacincola isolates were sensitive to all antibiotics used for sensitivity testing. Our findings add to the body of knowledge regarding the occurrence of pathogenic P. glacincola infection in Egyptian marine fishes and its potential effects on fish. Future large-scale surveys exploring this bacterium among other freshwater and marine fishes in Egypt would be helpful for the implementation of effective strategies for the prevention and control of this infection are warranted., (© 2022. The Author(s).)

Background: Antibody affinity maturation in vertebrates requires the enzyme activation-induced cytidine deaminase (AID) which initiates secondary antibody diversification by mutating the immunoglobulin loci. AID-driven antibody diversification is conserved across jawed vertebrates since bony and cartilaginous fish. Two exceptions have recently been reported, the Pipefish and Anglerfish, in which the AID-encoding aicda gene has been lost. Both cases are associated with unusual reproductive behavior, including male pregnancy and sexual parasitism. Several cold water fish in the Atlantic cod (Gadinae) family carry an aicda gene that encodes for a full-length enzyme but lack affinity-matured antibodies and rely on antibodies of broad antigenic specificity. Hence, we examined the functionality of their AID., Results: By combining genomics, transcriptomics, immune responsiveness, and functional enzymology of AID from 36 extant species, we demonstrate that AID of that Atlantic cod and related fish have extremely lethargic or no catalytic activity. Through ancestral reconstruction and functional enzymology of 71 AID enzymes, we show that this enzymatic inactivation likely took place relatively recently at the emergence of the true cod family (Gadidae) from their ancestral Gadiformes order. We show that this AID inactivation is not only concordant with the previously shown loss of key adaptive immune genes and expansion of innate and cell-based immune genes in the Gadiformes but is further reflected in the genomes of these fish in the form of loss of AID-favored sequence motifs in their immunoglobulin variable region genes., Conclusions: Recent demonstrations of the loss of the aicda gene in two fish species challenge the paradigm that AID-driven secondary antibody diversification is absolutely conserved in jawed vertebrates. These species have unusual reproductive behaviors forming an evolutionary pressure for a certain loss of immunity to avoid tissue rejection. We report here an instance of catalytic inactivation and functional loss of AID rather than gene loss in a conventionally reproducing vertebrate. Our data suggest that an expanded innate immunity, in addition to lower pathogenic pressures in a cold environment relieved the pressure to maintain robust secondary antibody diversification. We suggest that in this unique scenario, the AID-mediated collateral genome-wide damage would form an evolutionary pressure to lose AID function., (© 2022. The Author(s).)

Background: Antarctica harbors the bulk of the species diversity of the dominant teleost fish suborder-Notothenioidei. However, the forces that shape their evolution are still under debate., Results: We sequenced the genome of an icefish, Chionodraco hamatus, and used population genomics and demographic modelling of sequenced genomes of 52 C. hamatus individuals collected mainly from two East Antarctic regions to investigate the factors driving speciation. Results revealed four icefish populations with clear reproduction separation were established 15 to 50 kya (kilo years ago) during the last glacial maxima (LGM). Selection sweeps in genes involving immune responses, cardiovascular development, and photoperception occurred differentially among the populations and were correlated with population-specific microbial communities and acquisition of distinct morphological features in the icefish taxa. Population and species-specific antifreeze glycoprotein gene expansion and glacial cycle-paced duplication/degeneration of the zona pellucida protein gene families indicated fluctuating thermal environments and periodic influence of glacial cycles on notothenioid divergence., Conclusions: We revealed a series of genomic evidence indicating differential adaptation of C. hamatus populations and notothenioid species divergence in the extreme and unique marine environment. We conclude that geographic separation and adaptation to heterogeneous pathogen, oxygen, and light conditions of local habitats, periodically shaped by the glacial cycles, were the key drivers propelling species diversity in Antarctica., (© 2022. The Author(s).)

Background: Ray-finned fishes (Actinopterygii) perceive their environment through a range of sensory modalities, including olfaction. Anatomical diversity of the olfactory organ suggests that olfaction is differentially important among species. To explore this topic, we studied the evolutionary dynamics of the four main gene families (OR, TAAR, ORA/VR1 and OlfC/VR2) coding for olfactory receptors in 185 species of ray-finned fishes., Results: The large variation in the number of functional genes, between 28 in the ocean sunfish Mola mola and 1317 in the reedfish Erpetoichthys calabaricus, is the result of parallel expansions and contractions of the four main gene families. Several ancient and independent simplifications of the olfactory organ are associated with massive gene losses. In contrast, Polypteriformes, which have a unique and complex olfactory organ, have almost twice as many olfactory receptor genes as any other ray-finned fish., Conclusions: We document a functional link between morphology of the olfactory organ and richness of the olfactory receptor repertoire. Further, our results demonstrate that the genomic underpinning of olfaction in ray-finned fishes is heterogeneous and presents a dynamic pattern of evolutionary expansions, simplifications, and reacquisitions., (© 2022. The Author(s).)

Background: In aquatic environments, pH, salinity, and ammonia concentration are extremely important for aquatic animals. NHE is a two-way ion exchange carrier protein, which can transport Na + into cells and exchange out H + , and also plays key roles in regulating intracellular pH, osmotic pressure, and ammonia concentration., Results: In the present study, ten NHEs, the entire NHE gene family, were identified from Coilia nasus genome and systemically analyzed via phylogenetic, structural, and synteny analysis. Different expression patterns of C. nasus NHEs in multiple tissues indicated that expression profiles of NHE genes displayed tissue-specific. Expression patterns of C. nasus NHEs were related to ammonia excretion during multiple embryonic development stages. To explore the potential functions on salinity challenge and ammonia stress, expression levels of ten NHEs were detected in C. nasus gills under hypotonic stress, hypertonic stress, and ammonia stress. Expression levels of all NHEs were upregulated during hypotonic stress, while they were downregulated during hypertonic stress. NHE2 and NHE3 displayed higher expression levels in C. nasus larvae and juvenile gills under ammonia stress., Conclusions: Our study revealed that NHE genes played distinct roles in embryonic development, salinity stress, and ammonia exposure. Syntenic analysis showed significant difference between stenohaline fish and euryhaline fishes. Our findings will provide insight into effects of C. nasus NHE gene family on ion transport and ammonia tolerance and be beneficial for healthy aquaculture of C. nasus., (© 2022. The Author(s).)

Objectives: Anisakis pegreffii is a zoonotic parasite requiring marine organisms to complete its life-history. Human infection (anisakiasis) occurs when the third stage larvae (L3) are accidentally ingested with raw or undercooked infected fish or squids. A new de novo transcriptome of A. pegreffii was here generated aiming to provide a robust bulk of data to be used for a comprehensive "ready-to-use" resource for detecting functional studies on genes and gene products of A. pegreffii involved in the molecular mechanisms of parasite-host interaction., Data Description: A RNA-seq library of A. pegreffii L3 was here newly generated by using Illumina TruSeq platform. It was combined with other five RNA-seq datasets previously gathered from L3 of the same species stored in SRA of NCBI. The final dataset was analyzed by launching three assembler programs and two validation tools. The use of a robust pipeline produced a high-confidence protein-coding transcriptome of A. pegreffii. These data represent a more robust and complete transcriptome of this species with respect to the actually existing resources. This is of importance for understanding the involved adaptive and immunomodulatory genes implicated in the "cross talk" between the parasite and its hosts, including the accidental one (humans)., (© 2022. The Author(s).)

Background: Aggression is an evolutionarily conserved behavior critical for animal survival. In the fish Betta splendens, across different stages of fighting interactions, fighting opponents suffer from various stressors, especially from the great demand for oxygen. Using RNA sequencing, we profiled differential alternative splicing (DAS) events in the brains of fish collected before fighting, during fighting, and after fighting to study the involvement of alternative splicing (AS) in the response to stress during the fight., Results: We found that fighting interactions induced the greatest increase in AS in the 'during-fighting' fish, followed by that of the 'after-fighting' fish. Intron retention (IR) was the most enriched type among all the basic AS events. DAS genes were mainly associated with synapse assembly, ion transport, and regulation of protein secretion. We further observed that IR events significantly differentiated between winners and losers for 19 genes, which were associated with messenger RNA biogenesis, DNA repair, and transcription machinery. These genes share many common features, including shorter intron length and higher GC content., Conclusions: This study is the first comprehensive view of AS induced by fighting interactions in a fish species across different stages of those interactions, especially with respect to IR events in winners and losers. Together, these findings facilitate future investigations into transcriptome complexity and AS regulation in response to stress under the context of aggression in vertebrates., (© 2022. The Author(s).)

Objectives: American shad (Alosa sapidissima) is an important migratory fish under Alosinae and has long been valued for its economic, nutritional and cultural attributes. Overfishing and barriers across the passage made it vulnerable to sustain. To protect this valuable species, aquaculture action plans have been taken though there are no published genetic resources prevailing yet. Here, we reported the first de novo assembled and annotated transcriptome of A. sapidissima using blood and brain tissues., Data Description: We generated 160,481 and 129,040 non-redundant transcripts from brain and blood tissues. The entire work strategy involved RNA extraction, library preparation, sequencing, de novo assembly, filtering, annotation and validation. Both coding and non-coding transcripts were annotated against Swissprot and Pfam datasets. Nearly, 83% coding transcripts were functionally assigned. Protein clustering with clupeiform and non-clupeiform taxa revealed ~ 82% coding transcripts retained the orthologue relationship which improved confidence over annotation procedure. This study will serve as a useful resource in future for the research community to elucidate molecular mechanisms for several key traits like migration which is fascinating in clupeiform shads., (© 2022. The Author(s).)

Background: Seahorses, seadragons, pygmy pipehorses, and pipefishes (Syngnathidae, Syngnathiformes) are among the most recognizable groups of fishes because of their derived morphology, unusual life history, and worldwide distribution. Despite previous phylogenetic studies and recent new species descriptions of syngnathids, the evolutionary relationships among several major groups within this family remain unresolved., Results: Here, we provide a reconstruction of syngnathid phylogeny based on genome-wide sampling of 1314 ultraconserved elements (UCEs) and expanded taxon sampling to assess the current taxonomy and as a basis for macroevolutionary insights. We sequenced a total of 244 new specimens across 117 species and combined with published UCE data for a total of 183 species of Syngnathidae, about 62% of the described species diversity, to compile the most data-rich phylogeny to date. We estimated divergence times using 14 syngnathiform fossils, including nine fossils with newly proposed phylogenetic affinities, to better characterize current and historical biogeographical patterns, and to reconstruct diversification through time. We present a phylogenetic hypothesis that is well-supported and provides several notable insights into syngnathid evolution. We found nine non-monophyletic genera, evidence for seven cryptic species, five potentially invalid synonyms, and identified a novel sister group to the seahorses, the Indo-Pacific pipefishes Halicampus macrorhynchus and H. punctatus. In addition, the morphologically distinct southwest Pacific seahorse Hippocampus jugumus was recovered as the sister to all other non-pygmy seahorses. As found in many other groups, a high proportion of syngnathid lineages appear to have originated in the Central Indo-Pacific and subsequently dispersed to adjoining regions. Conversely, we also found an unusually high subsequent return of lineages from southern Australasia to the Central Indo-Pacific. Diversification rates rose abruptly during the Middle Miocene Climate Transition and peaked after the closure of the Tethys Sea., Conclusions: Our results reveal a previously underappreciated diversity of syngnathid lineages. The observed biogeographic patterns suggest a significant role of the southern Australasian region as a source and sink of lineages. Shifts in diversification rates imply possible links to declining global temperatures, the separation of the Atlantic and Pacific faunas, and the environmental changes associated with these events., (© 2022. The Author(s).)

Background: Environmental DNA (eDNA) is an effective tool for the detection and monitoring of presence or absence of rare and invasive species. These techniques have been extended to quantify biomass in vertebrates, particularly in fish species. However, the efficacy of eDNA techniques to quantify biomass in invertebrate species has rarely been examined. This study tested whether eDNA could be used to determine the biomass of the world-wide invasive green crab, Carcinus maenas. In a controlled laboratory study, the relationship between biomass and C. maenas eDNA concentration was examined in the context of different biotic (activity) and abiotic (temperature) parameters., Results: When incubating different numbers of crabs in sterile saltwater for up to 7 days, a relationship between eDNA concentration and biomass was observed at temperatures of 6.7 ℃ and 18.7 ℃, but not at 12.8 ℃. Additionally, motor activity, aggression level, time of sampling, and features of organismal decay had significant impact on the concentration of C. maenas eDNA collected., Conclusions: We show that eDNA concentration did not correlate with biomass, and that biomass, temperature, organismal characteristics, and potentially many more parameters affect shedding and degradation rates for eDNA in this species, thus, impacting the recoverable eDNA concentration. Therefore, eDNA techniques are not likely to provide a reliable signal of biomass in the invasive invertebrate species C. maenas., (© 2022. The Author(s).)

Background: Hatching is identified as one of the most important events in the reproduction of oviparous vertebrates. The genes for hatching enzymes, which are vital in the hatching process, are conserved among vertebrates. However, especially in teleost, it is difficult to trace their molecular evolution in detail due to the presence of other C6astacins, which are the subfamily to which the genes for hatching enzymes belong and are highly diverged. In particular, the hatching enzyme genes are diversified with frequent genome translocations due to retrocopy., Results: In this study, we took advantage of the rapid expansion of whole-genome data in recent years to examine the molecular evolutionary process of these genes in vertebrates. The phylogenetic analysis and the genomic synteny analysis revealed C6astacin genes other than the hatching enzyme genes, which was previously considered to be retained only in teleosts, was also retained in the genomes of basal ray-finned fishes, coelacanths, and cartilaginous fishes. These results suggest that the common ancestor of these genes can be traced back to at least the common ancestor of the Gnathostomata. Moreover, we also found that many of the C6astacin genes underwent multiple gene duplications during vertebrate evolution, and the results of gene expression analysis in frogs implied that genes derived from hatching enzyme genes underwent neo-functionalization., Conclusions: In this study, we describe in detail the molecular evolution of the C6astacin gene in vertebrates, which has not been summarized previously. The results revealed the presence of the previously unknown C6astacin gene in the basal-lineage of jawed vertebrates and large-scale gene duplication of hatching enzyme genes in amphibians. The comprehensive investigation reported in this study will be an important basis for studying the molecular evolution of the vertebrate C6astacin genes, hatching enzyme, and its paralogous genes and for identifying these genes without the need for gene expression and functional analysis., (© 2022. The Author(s).)

Background: Fishes are the one of the most diverse groups of animals with respect to their modes of sex determination, providing unique models for uncovering the evolutionary and molecular mechanisms underlying sex determination and reversal. Here, we have investigated how sex is determined in a species of both commercial and ecological importance, the Siamese fighting fish Betta splendens., Results: We conducted association mapping on four commercial and two wild populations of B. splendens. In three of the four commercial populations, the master sex determining (MSD) locus was found to be located in a region of ~ 80 kb on LG2 which harbours five protein coding genes, including dmrt1, a gene involved in male sex determination in different animal taxa. In these fish, dmrt1 shows a male-biased gonadal expression from undifferentiated stages to adult organs and the knockout of this gene resulted in ovarian development in XY genotypes. Genome sequencing of XX and YY genotypes identified a transposon, drbx1, inserted into the fourth intron of the X-linked dmrt1 allele. Methylation assays revealed that epigenetic changes induced by drbx1 spread out to the promoter region of dmrt1. In addition, drbx1 being inserted between two closely linked cis-regulatory elements reduced their enhancer activities. Thus, epigenetic changes, induced by drbx1, contribute to the reduced expression of the X-linked dmrt1 allele, leading to female development. This represents a previously undescribed solution in animals relying on dmrt1 function for sex determination. Differentiation between the X and Y chromosomes is limited to a small region of ~ 200 kb surrounding the MSD gene. Recombination suppression spread slightly out of the SD locus. However, this mechanism was not found in the fourth commercial stock we studied, or in the two wild populations analysed, suggesting that it originated recently during domestication., Conclusions: Taken together, our data provide novel insights into the role of epigenetic regulation of dmrt1 in sex determination and turnover of SD systems and suggest that fighting fish are a suitable model to study the initial stages of sex chromosome evolution., (© 2021. The Author(s).)

Background: To determine the presence and abundance of an aquatic species in large waterbodies, especially when populations are at low densities, is highly challenging for conservation biologists. Environmental DNA (eDNA) has the potential to offer a noninvasive and cost-effective method to complement traditional population monitoring, however, eDNA has not been extensively applied to study large migratory species. Chinese sturgeon (Acipenser sinensis), is the largest anadromous migratory fish in the Yangtze River, China, and in recent years its population has dramatically declined and spawning has failed, bringing this species to the brink of extinction. In this study, we aim to test the detectability of eDNA methods to determine the presence and relative abundance of reproductive stock of the species and whether eDNA can be used as a tool to reflect behavioral patterns. Chinese sturgeon eDNA was collected from four sites along the spawning ground across an eight month period, to investigate the temporal and spatial distribution using droplet digital PCR (ddPCR)., Results: We designed a pair of specific primers for Chinese sturgeon and demonstrated the high sensitivity of ddPCR to detect and quantify the Chinese sturgeon eDNA concentration with the limit of detection 0.17 copies/μl, with Chinese sturgeon eDNA been intermittently detected at all sampling sites. There was a consistent temporal pattern among four of the sampling sites that could reflect the movement characteristics of the Chinese sturgeon in the spawning ground, but without a spatial pattern. The eDNA concentration declined by approximately 2-3 × between December 2018 and December 2019., Conclusions: The results prove the efficacy of eDNA for monitoring reproductive stock of the Chinese sturgeon and the e decreased eDNA concentration reflect that Chinese sturgeon may survive with an extremely small number of reproductive stock in the Yangtze River. Accordingly, we suggest future conservation measures should focus on both habitat restoration and matured fish restocking to ensure successful spawning. Overall, this study provides encouraging support for the application of eDNA methods to monitor endangered aquatic species., (© 2021. The Author(s).)

Objective: The Caribbean is an important global biodiversity hotspot. Adaptive radiations there lead to many speciation events within a limited period and hence are particularly prominent biodiversity generators. A prime example are freshwater fish of the genus Limia, endemic to the Greater Antilles. Within Hispaniola, nine species have been described from a single isolated site, Lake Miragoâne, pointing towards extraordinary sympatric speciation. This study examines the evolutionary history of the Limia species in Lake Miragoâne, relative to their congeners throughout the Caribbean., Results: For 12 Limia species, we obtained almost complete sequences of the mitochondrial cytochrome b gene, a well-established marker for lower-level taxonomic relationships. We included sequences of six further Limia species from GenBank (total N  = 18 species). Our phylogenies are in concordance with other published phylogenies of Limia. There is strong support that the species found in Lake Miragoâne in Haiti are monophyletic, confirming a recent local radiation. Within Lake Miragoâne, speciation is likely extremely recent, leading to incomplete lineage sorting in the mtDNA. Future studies using multiple unlinked genetic markers are needed to disentangle the relationships within the Lake Miragoâne clade., (© 2021. The Author(s).)

Background: The genetic control of sex determination in teleost species is poorly understood. This is partly because of the diversity of mechanisms that determine sex in this large group of vertebrates, including constitutive genes linked to sex chromosomes, polygenic constitutive mechanisms, environmental factors, hermaphroditism, and unisexuality. Here we use a de novo genome assembly of New Zealand silver trevally (Pseudocaranx georgianus) together with sex-specific whole genome sequencing data to detect sexually divergent genomic regions, identify candidate genes and develop molecular makers., Results: The de novo assembly of an unsexed trevally (Trevally&#95;v1) resulted in a final assembly of 579.4 Mb in length, with a N50 of 25.2 Mb. Of the assembled scaffolds, 24 were of chromosome scale, ranging from 11 to 31 Mb in length. A total of 28,416 genes were annotated after 12.8 % of the assembly was masked with repetitive elements. Whole genome re-sequencing of 13 wild sexed trevally (seven males and six females) identified two sexually divergent regions located on two scaffolds, including a 6 kb region at the proximal end of chromosome 21. Blast analyses revealed similarity between one region and the aromatase genes cyp19 (a1a/b) (E-value < 1.00E-25, identity > 78.8 %). Males contained higher numbers of heterozygous variants in both regions, while females showed regions of very low read-depth, indicative of male-specificity of this genomic region. Molecular markers were developed and subsequently tested on 96 histologically-sexed fish (42 males and 54 females). Three markers amplified in absolute correspondence with sex (positive in males, negative in females)., Conclusions: The higher number of heterozygous variants in males combined with the absence of these regions in females support a XY sex-determination model, indicating that the trevally&#95;v1 genome assembly was developed from a male specimen. This sex system contrasts with the ZW sex-determination model documented in closely related carangid species. Our results indicate a sex-determining function of a cyp19a1a-like gene, suggesting the molecular pathway of sex determination is somewhat conserved in this family. The genomic resources developed here will facilitate future comparative work, and enable improved insights into the varied sex determination pathways in teleosts. The sex marker developed in this study will be a valuable resource for aquaculture selective breeding programmes, and for determining sex ratios in wild populations., (© 2021. The Author(s).)

Background: Microsatellite repeats are ubiquitous in organism genomes and play an important role in the chromatin organization, regulation of gene activity, recombination and DNA replication. Although microsatellite distribution patterns have been studied in most phylogenetic lineages, they are unclear in fish species., Results: Here, we present the first systematic examination of microsatellite distribution in coding and non-coding regions of 14 fish genomes. Our study showed that the number and type of microsatellites displayed nonrandom distribution for both intragenic and intergenic regions, suggesting that they have potential roles in transcriptional or translational regulation and DNA replication slippage theories alone were insufficient to explain the distribution patterns. Our results showed that microsatellites are dominant in non-coding regions. The total number of microsatellites ranged from 78,378 to 1,012,084, and the relative density varied from 4925.76 bp/Mb to 25,401.97 bp/Mb. Overall, (A + T)-rich repeats were dominant. The dependence of repeat abundance on the length of the repeated unit (1-6 nt) showed a great similarity decrease, whereas more tri-nucleotide repeats were found in exonic regions than tetra-nucleotide repeats of most species. Moreover, the incidence of different repeated types appeared species- and genomic-specific. These results highlight potential mechanisms for maintaining microsatellite distribution, such as selective forces and mismatch repair systems., Conclusions: Our data could be beneficial for the studies of genome evolution and microsatellite DNA evolutionary dynamics, and facilitate the exploration of microsatellites structural, function, composition mode and molecular markers development in these species.

Background: In commercial fish, dominance effects could be exploited by predicting production abilities of the offspring that would be generated by different mating pairs and choosing those pairs that maximise the average offspring phenotype. Consequently, matings would be performed to reduce inbreeding depression. This can be achieved by applying mate selection (MS) that combines selection and mating decisions in a single step. An alternative strategy to MS would be to apply minimum coancestry mating (MCM) after selection based on estimated breeding values. The objective of this study was to evaluate, by computer simulations, the potential benefits that can be obtained by implementing MS or MCM based on genomic data for exploiting dominance effects when creating commercial fish populations that are derived from a breeding nucleus., Methods: The selected trait was determined by a variable number of loci with additive and dominance effects. The population consisted of 50 full-sib families with 30 offspring each. Males and females with the highest estimated genomic breeding values were selected in the nucleus and paired using the MCM strategy. Both MCM and MS were used to create the commercial population., Results: For a moderate number of SNPs, equal or even higher mean phenotypic values are obtained by selecting on genomic breeding values and then applying MCM than by using MS when the trait exhibited substantial inbreeding depression. This could be because MCM leads to high levels of heterozygosity across the whole genome, even for loci affecting the trait that are in linkage equilibrium with the SNPs. In contrast, MS specifically promotes heterozygosity for SNPs for which a dominance effect has been detected., Conclusions: In most scenarios, for the management of aquaculture breeding programs it seems advisable to follow the MCM strategy when creating the commercial population, especially for traits with large inbreeding depression. Moreover, MCM has the appealing property of reducing inbreeding levels, with a corresponding reduction in inbreeding depression for traits beyond those included in the selection objective.

Background: The electrosensory ampullary organs (AOs) and mechanosensory neuromasts (NMs) found in sturgeon and some other non-neopterygian fish or amphibians are both originated from lateral line placodes. However, these two sensory organs have characteristic morphological and physiological differences. The molecular mechanisms for the specification of AOs and NMs are not clearly understood., Results: We sequenced the transcriptome for neomycin treated sturgeon AOs and NMs in the early regeneration stages, and de novo assembled a sturgeon transcriptome. By comparing the gene expression differences among untreated AOs, NMs and general epithelia (EPs), we located some specific genes for these two sensory organs. In sturgeon lateral line, the voltage-gated calcium channels and voltage-gated potassium channels were predominant calcium and potassium channel subtypes, respectively. And by correlating gene expression with the regeneration process, we predicated several candidate key transcriptional regulation related genes might be involved in AOs and NMs regeneration., Conclusions: Genes with specific expression in the two lateral line sensory organs suggests their important roles in mechanoreceptor and electroreceptor formation. The candidate transcriptional regulation related genes may be important for mechano- and electro- receptor specification, in a "dosage-related" manner. These results suggested the molecular basis for specification of these two sensory organs in sturgeon.

Background: Hundreds of genomes and transcriptomes of fish species have been sequenced in recent years. However, fish scholarship currently lacks a comprehensive, integrated, and up-to-date collection of fish genomic data., Results: Here we present FishDB, the first database for fish multi-level omics data, available online at http://fishdb.ihb.ac.cn . The database contains 233 fish genomes, 201 fish transcriptomes, 5841 fish mitochondrial genomes, 88 fish gene sets, 16,239 miRNAs of 65 fishes, 1,330,692 piRNAs and 4852 lncRNAs of Danio rerio, 59,040 Mb untranslated regions (UTR) of 230 fishes, and 31,918 Mb coding sequences (CDS) of 230 fishes. Among these, we newly generated a total of 11 fish genomes and 53 fish transcriptomes., Conclusions: This release contains over 410,721.67 Mb sequences and provides search functionality, a BLAST server, JBrowse, and PrimerServer modules.

Background: Studying functional divergences between paralogs that originated from genome duplication is a significant topic in investigating molecular evolution. Genes that exhibit basal level cyclic expression patterns including circadian and light responsive genes are important physiological regulators. Temporal shifts in basal gene expression patterns are important factors to be considered when studying genetic functions. However, adequate efforts have not been applied to studying basal gene expression variation on a global scale to establish transcriptional activity baselines for each organ. Furthermore, the investigation of cyclic expression pattern comparisons between genome duplication created paralogs, and potential functional divergence between them has been neglected. To address these questions, we utilized a teleost fish species, Xiphophorus maculatus, and profiled gene expression within 9 organs at 3-h intervals throughout a 24-h diurnal period., Results: Our results showed 1.3-21.9% of genes in different organs exhibited cyclic expression patterns, with eye showing the highest fraction of cycling genes while gonads yielded the lowest. A majority of the duplicated gene pairs exhibited divergences in their basal level expression patterns wherein only one paralog exhibited an oscillating expression pattern, or both paralogs exhibit oscillating expression patterns, but each gene duplicate showed a different peak expression time, and/or in different organs., Conclusions: These observations suggest cyclic genes experienced significant sub-, neo-, or non-functionalization following the teleost genome duplication event. In addition, we developed a customized, web-accessible, gene expression browser to facilitate data mining and data visualization for the scientific community.

Background: EP300 is a conserved protein in vertebrates, which serves as a key mediator of cellular homeostasis. Mutations and dysregulation of EP300 give rise to severe human developmental disorders and malignancy. Danio rerio is a promising model organism to study EP300 related diseases and drugs; however, the effect of EP300 duplicates derived from teleost-specific whole genome duplication should not just be neglected., Results: In this study, we obtained EP300 protein sequences of representative teleosts, mammals and sauropsids, with which we inferred a highly supported maximum likelihood tree. We observed that Ep300 duplicates (Ep300a and Ep300b) were widely retained in teleosts and universally expressed in a variety of tissues. Consensus sequences of Ep300a and Ep300b had exactly the same distribution of conserved domains, suggesting that their functions should still be largely overlapped. We analyzed the molecular evolution of Ep300 duplicates in teleosts, using branch-site models, clade models and site models. The results showed that both duplicates were subject to strong positive selection; however, for an extant species, generally at most one copy was under positive selection. At the clade level, there were evident positive correlations between evolutionary rates, the number of positively selected sites and gene expression levels. In Ostariophysi, Ep300a were under stronger positive selection than Ep300b; in Neoteleostei, another species-rich teleost clade, the contrary was the case. We also modeled 3D structures of zf-TAZ domain and its flanking regions of Ep300a and Ep300b of D. rerio and Oryzias latipes and found that in either species the faster evolving copy had more short helixes., Conclusions: Collectively, the two copies of Ep300 have undoubtedly experienced directional divergence in main teleost clades. The divergence of EP300 between teleosts and mammals should be greater than the divergence between different teleost clades. Further studies are needed to clarify to what extent the EP300 involved regulatory network has diverged between teleosts and mammals, which would also help explain the huge success of teleosts.

Background: Through its ability to open pores in cell membranes, perforin-1 plays a key role in the immune system. Consistent with this role, the gene encoding perforin shows hallmarks of complex evolutionary events, including amplification and pseudogenization, in multiple species. A large proportion of these events occurred in phyla for which scarce genomic data were available. However, recent large-scale genomics projects have added a wealth of information on those phyla. Using this input, we annotated perforin-1 homologs in more than eighty species including mammals, reptiles, birds, amphibians and fishes., Results: We have annotated more than 400 perforin genes in all groups studied. Most mammalian species only have one perforin locus, which may contain a related pseudogene. However, we found four independent small expansions in unrelated members of this class. We could reconstruct the full-length coding sequences of only a few avian perforin genes, although we found incomplete and truncated forms of these gene in other birds. In the rest of reptilia, perforin-like genes can be found in at least three different loci containing up to twelve copies. Notably, mammals, non-avian reptiles, amphibians, and possibly teleosts share at least one perforin-1 locus as assessed by flanking genes. Finally, fish genomes contain multiple perforin loci with varying copy numbers and diverse exon/intron patterns. We have also found evidence for shorter genes with high similarity to the C2 domain of perforin in several teleosts. A preliminary analysis suggests that these genes arose at least twice during evolution from perforin-1 homologs., Conclusions: The assisted annotation of new genomic assemblies shows complex patterns of birth-and-death events in the evolution of perforin. These events include duplication/pseudogenization in mammals, multiple amplifications and losses in reptiles and fishes and at least one case of partial duplication with a novel start codon in fishes.