A high quality assembly of the Nile Tilapia (Oreochromis niloticus) genome reveals the structure of two sex determination regions

Reconstruction of the X and Y haplotypes in the genetically improved Abbassa nile tilapia genome assembly

The success of the Nile tilapia (Oreochromis niloticus) as an aquaculture species is partly the result of continuous selective breeding leading to high performing strains. These elite strains have been derived from breeding populations of diverse origins and crosses with other Oreochromis species. Owing to the complex and unique evolutionary histories of each strain, existing reference genomes of wild populations are unsuitable to implement genomic selection for beneficial traits such as growth or environmental resilience in aquaculture programmes. Here we generated a high-quality genome assembly and annotation of the WorldFish Genetically Improved Abbassa Nile tilapia (GIANT) elite strain using a combination of PacBio HiFi, and Omni-C Illumina sequencing. As a male Abbassa Nile tilapia was used for the generation of the genome assembly, we reconstructed both X and Y haplotypes, identifying both amhY and amhΔy on LG23 indicating that Abbassa likely shares the same sex determination system as GIFT, and thereby differs from the existing reference genome, whose sex determination loci are located on LG1.

Heat inducible nuclear translocation of Kdm6bb drives temperature dependent sex reversal in Nile tilapia

Understanding the primary molecular events driving temperature-dependent sex reversal (TSR) has proven challenging, particularly in distinguishing these from secondary effects of sexual differentiation. The mechanisms translating temperature into a sex-determining signal in fish are still largely unknown. Through combined transcriptomic and genome-wide histone methylation analyses of gonads in Nile tilapia (Oreochromis niloticus) exposed to normal and elevated temperatures, we observed significant upregulation of male-promoting genes (amh, dmrt1, gsdf) and suppression of female-promoting genes (wt1a and foxl3) at high temperature. These changes were correlated with methylation changes in H3K27 and H3K4 in the promoter regions of these genes. Among the histone methylation enzymes induced by high temperature, we identified the H3K27 demethylase Kdm6bb to be a key factor. Gene deletion and biochemical studies confirmed that Kdm6bb significantly impacts the H3K27 methylation level, that influences sex determination. Crucially, we discovered that the TSR function of Kdm6bb is mediated by the alternative inclusion of a previously unrecognized intron, enabling nuclear translocation of the demethylase to perform its function. Our findings refute the previously proposed "translation deficiency" mechanism of kdm6bb, and highlight the critical role of mRNA alternative splicing and subcellular localization of the demethylase in temperature-induced sex reversal.

Genome assembly variation and its implications for gene discovery in nematodes

Genome assemblers are a critical component of genome science, but the choice of assembly software and protocols can be daunting. Here, we investigate genome assembly variation and its implications for gene discovery across three nematode species-Caenorhabditis bovis, Haemonchus contortus, and Heligmosomoides bakeri-highlighting the critical interplay between assembly choice and downstream genomic analysis. Selecting commonly used genome assemblers, we generated multiple assemblies for each species, analyzing their structure, completeness, and effect on gene family analysis. Our findings demonstrate that assembly variations can significantly affect gene family composition, with notable differences in gene families important in anthelmintic discovery and immunomodulation. Despite broadly similar performance using various assembly metrics, comparisons of assemblies with a single species revealed underlying structural rearrangements and inconsistencies in gene content, which would affect downstream analyses. This emphasizes the need for continuous refinement of genome assemblies and their annotations.

The origin, evolution, and translocation of sex chromosomes in Silurus catfish mediated by transposons

BACKGROUND

Sex chromosome (SC) evolution is a longstanding topic of focus in evolutionary biology. Teleosts often exhibit rapid turnover of SCs and sex-determining (SD) genes, alongside a diverse range of SC differentiation mechanisms.

RESULTS

On the basis of new chromosome-scale assemblies of three Silurus species (S. microdorsalis, S. glanis, and S. lanzhouensis) and two outgroup species (Pterocryptis cochinchinensis and Kryptopterus bicirrhis), along with our previous assemblies of S. meridionalis and S. asotus, we traced the evolution of SC in the Silurus genus (Siluriformes), following the fate of the known SD gene amhr2y. Phylogenetic analysis showed that amhr2y occurred at least before the divergence of Pterocryptis, Kryptopterus, and Silurus and lost in P. cochinchinensis and K. bicirrhis. Chr24 has become the SC in the ancestor of five Silurus species due to the duplication-and-translocation of amhr2 mediated by LTR transposon. Then, a proto Y was formed and maintained with a shared 60 kb male-specific region of the Y chromosome (MSY) by transposable elements (TEs) expansion and gene gathering. Due to the continuous TEs accumulation, genes other than amhr2y in MSYs have degenerated or been lost, while non-recombinant regions continue to expend, forming MSYs of different sizes in different Silurus species (from 320 to 550 kb). Two turnover events, one homologous (from the left arm to the right arm of Chr24) and one nonhomologous (from Chr24 to Chr5), occurring among five Silurus species were possibly mediated by hAT and Helitron transposons.

CONCLUSIONS

Our results on the dynamic evolutionary trajectory of SD gene amhr2y, MSYs, and SCs in Silurus catfish indicated the variability and diversity of fish SCs and confirmed that frequent turnover is an important way to maintain the homology and low differentiation of fish SCs.

Identification of a Novel QTL on LG16 Associated with Acute Salt Tolerance in Red Tilapia (Oreochromis spp.) Using GWAS

Culturing saline tilapia has become a new trend in the aquaculture due to the scarcity of freshwater resources. In this study, the genetic basis controlling for salt tolerance were investigated by using a ddRAD-seq-based GWAS in 288 individuals with extreme salt tolerant traits from half-sib families of red tilapia. 12 genome-wide significant SNPs and 6 chromosome-wide significant SNPs associated with acute salt tolerance were identified. Two QTLs on LG18:25,593,701-7009020 and on LG16:19,735,164-21,231,391 were defined. It is noteworthy that the QTL on LG16 is a novel QTL associated with acute salt stress. Near the significant SNP sites, we identified candidate genes sik1, ltb4r2b, pnp5b and kirrel1b with differential transcript expression under salt stress. Furthermore, significant physiological differences in serum osmolality and ion concentrations were confirmed between the tolerant group and sensitive group under 4.5 h of 22 ppt stress. The sensitive group had much higher serum osmolality (osmolality: 642.20 ± 6.30 mOsm/kg) and higher concentrations of sodium and chloride ions (sodium: 317.67 ± 5.03 mmol/L and chloride: 316.43 ± 8.28 mmol/L) than the tolerant group (547.60 ± 15.44 mOsm/kg, p osmolality = 0.0002; sodium: 280.53 ± 9.13 mmol/L, p sodium < 0.0242; chloride: 266.00 ± 12.00 mmol/L, p chloride < 0.0184). However, the lowest bicarbonate concentration was detected in the sensitive group at 22 ppt (2.53 ± 0.30 mmol/L), which was significantly different from both the sensitive group at 0 ppt (p = 0.0008) and the tolerant group at 22 ppt (p = 0.0164). Our research laid the foundation for exploring the genetic mechanisms of acute salt tolerance and osmoregulation in red tilapia and for developing strains of red tilapia adapted to saltwater.

Dietary purslane (Portulaca oleracea L.) leaf powder maintains growth and intestinal health in Oreochromis niloticus under chronic water-borne cadmium exposure by strengthening the gut barriers, modulating the intestinal nutrient transporters, and relieving oxidative stress

High cadmium (Cd) concentrations pose a threat to aquatic life globally. This study examined the efficiency of adding purslane (Portulaca oleracea L.) leaf powder (PLP) to Oreochromis niloticus diets on Cd's negative effects. PLP was analyzed using high-performance liquid chromatography, and its main constituents were gallic acid, chlorogenic acid, and pyrocatechol. Nile tilapia (180 fish, 34.5 ± 0.5 g) were divided into four groups in triplicate. A basal diet was given to the control group. The PLP group received a basal diet containing 10 g PLP/kg diet. The Cd group was exposed to 50 µg/L water. The Cd + PLP group was exposed to Cd and fed diets containing PLP. Results showed that PLP significantly rescued Cd-induced effects. PLP improved fish survival, feed conversion ratio, and growth retardation caused by Cd. PLP also restored decreased activities of lipase, trypsin, and amylase in the intestine. Furthermore, PLP corrected disturbances in leptin and growth hormone levels induced by Cd. PLP mitigated pathological alterations, replenished antioxidants (SOD, CAT, and GSH), and reduced lipid peroxidation in the intestinal tissues. PLP supplementation depleted significant Cd accumulation in the intestine and muscles. Additionally, PLP corrected altered expressions of tight junction proteins (zo-1, zo-2, and claudin-4) and nutrient transporters (glut-1, slc15a2, slc26a6, and slc4a4) in Cd-exposed fish. Conclusively, PLP shows promise as a dietary supplement to mitigate Cd's harmful impacts on fish growth. Its antioxidant activity and regulation of intestinal tight junction proteins and nutrient transporters contribute to its effectiveness. PLP supplementation holds the potential for reducing the detrimental effects of Cd in aquaculture.

Towards Reliable Detection of Introgression in the Presence of Among-Species Rate Variation

The role of interspecific hybridization has recently seen increasing attention, especially in the context of diversification dynamics. Genomic research has now made it abundantly clear that both hybridization and introgression-the exchange of genetic material through hybridization and backcrossing-are far more common than previously thought. Besides cases of ongoing or recent genetic exchange between taxa, an increasing number of studies report "ancient introgression"- referring to results of hybridization that took place in the distant past. However, it is not clear whether commonly used methods for the detection of introgression are applicable to such old systems, given that most of these methods were originally developed for analyses at the level of populations and recently diverged species, affected by recent or ongoing genetic exchange. In particular, the assumption of constant evolutionary rates, which is implicit in many commonly used approaches, is more likely to be violated as evolutionary divergence increases. To test the limitations of introgression detection methods when being applied to old systems, we simulated thousands of genomic datasets under a wide range of settings, with varying degrees of among-species rate variation and introgression. Using these simulated datasets, we showed that some commonly applied statistical methods, including the D-statistic and certain tests based on sets of local phylogenetic trees, can produce false-positive signals of introgression between divergent taxa that have different rates of evolution. These misleading signals are caused by the presence of homoplasies occurring at different rates in different lineages. To distinguish between the patterns caused by rate variation and genuine introgression, we developed a new test that is based on the expected clustering of introgressed sites along the genome and implemented this test in the program Dsuite.

A pheromone receptor in cichlid fish mediates attraction to females but inhibits male parental care

Reproductive behaviors differ across species, but the mechanisms that control variation in mating and parental care systems remain unclear. In many animal species, pheromones guide mating and parental care. However, it is not well understood how vertebrate pheromone signaling evolution can lead to new reproductive behavior strategies. In fishes, prostaglandin F2α (PGF2α) drives mating and reproductive pheromone signaling in fertile females, but this pheromonal activity appears restricted to specific lineages, and it remains unknown how a female fertility pheromone is sensed for most fish species. Here, we utilize single-cell transcriptomics and CRISPR gene editing in a cichlid fish model to identify and test the roles of key genes involved in olfactory sensing of reproductive cues. We find that a pheromone receptor, Or113a, detects fertile cichlid females and thereby promotes male attraction and mating behavior, sensing a ligand other than PGF2α. Furthermore, while cichlid fishes exhibit extensive parental care, for most species, care is provided solely by females. We find that males initiate mouthbrooding parental care if they have disrupted signaling in ciliated sensory neurons due to cnga2b mutation or if or113a is inactivated. Together, these results show that distinct mechanisms of pheromonal signaling drive reproductive behaviors across taxa. Additionally, these findings indicate that a single pheromone receptor has gained a novel role in behavior regulation, driving avoidance of paternal care among haplochromine cichlid fishes. Lastly, a sexually dimorphic, evolutionarily derived parental behavior is controlled by central circuits present in both sexes, while olfactory signals gate this behavior in a sex-specific manner.

Genome assemblies for Chromidotilapia guntheri (Teleostei: Cichlidae) identify a novel candidate gene for vertebrate sex determination, RIN3

Advances in genome sequencing have greatly accelerated the identification of sex chromosomes in a variety of species. Many of these species have experienced structural rearrangements that reduce recombination between the sex chromosomes, allowing the accumulation of sequence differences over many megabases. Identification of the genes that are responsible for sex determination within these sometimes large regions has proved difficult. Here, we identify an XY sex chromosome system on LG19 in the West African cichlid fish Chromidotilapia guntheri in which the region of differentiation extends over less than 400 kb. We develop high-quality male and female genome assemblies for this species, which confirm the absence of structural variants, and which facilitate the annotation of genes in the region. The peak of differentiation lies within rin3, which has experienced several debilitating mutations on the Y chromosome. We suggest two hypotheses about how these mutations might disrupt endocytosis, leading to Mendelian effects on sexual development.

Transcriptomics and gut microbiome analysis of the edible herb Bidens pilosa as a functional feed additive to promote growth and metabolism in tilapia (Oreochromis spp.)

To reduce the use of antibiotics and chemicals in aquaculture, an edible herb, Bidens pilosa, has been selected as a multifunctional feed additive. Although there has been considerable research into the effects of B. pilosa on poultry, the wider effects of B. pilosa, particularly on the growth and gut microbiota of fish, remain largely unexplored. We aimed to investigate the interactive effects between the host on growth and the gut microbiota using transcriptomics and the gut microbiota in B. pilosa-fed tilapia. In this study, we added 0.5% and 1% B. pilosa to the diet and observed that the growth performance of tilapia significantly increased over 8 weeks of feeding. Comparative transcriptome analysis was performed on RNA sequence profiles obtained from liver and muscle tissues. Functional enrichment analysis revealed that B. pilosa regulates several pathways and genes involved in amino acid metabolism, lipid metabolism, carbohydrate metabolism, endocrine system, signal transduction, and metabolism of other amino acids. The expression of the selected growth-associated genes was validated by qRT-PCR. The qRT-PCR results indicated that B. pilosa may enhance growth performance by activating the expression of the liver igf1 and muscle igf1rb genes and inhibiting the expression of the muscle negative regulator mstnb. Both the enhancement of liver endocrine IGF1/IGF1Rb signaling and the suppression of muscle autocrine/paracrine MSTN signaling induced the expression of myogenic regulatory factors (MRFs), myod1, myog and mrf4 in muscle to promote muscle growth in tilapia. The predicted function of the gut microbiota showed several significantly different pathways that overlapped with the KEGG enrichment results of differentially expressed genes in the liver transcriptomes. This finding suggested that the gut microbiota may influence liver metabolism through the gut-liver axis in B. pilosa-fed tilapia. In conclusion, dietary B. pilosa can regulate endocrine IGF1 signaling and autocrine/paracrine MSTN signaling to activate the expression of MRFs to promote muscle growth and alter the composition of gut bacteria, which can then affect liver amino acid metabolism, carbohydrate metabolism, endocrine system, lipid metabolism, metabolism of other amino acids, and signal transduction in the host, ultimately enhancing growth performance. Our results suggest that B. pilosa has the potential to be a functional additive that can be used as an alternative to reduce antibiotic use as a growth promoter in aquaculture.

The regulatory effects of water probiotic supplementation on the blood physiology, reproductive performance, and its related genes in Red Tilapia (Oreochromis niloticus X O. mossambicus)

Probiotics are becoming increasingly popular as eco-friendly alternatives in aquaculture. However, there is limited research on their impacts on the reproductive efficiency of Red Tilapia (Oreochromis niloticus x O. mossambicus) broodstock. Therefore, this experiment aimed to explore the combined effects of selective probiotics Bacillus subtilis and B. licheniformis (BSL; 1:1) added to water on blood hematology, serum metabolites, gonadal histology, reproductive performance, and reproductive associated genes in Red Tilapia broodstock. Tilapia broodfish weighing 140-160 g were stocked in four treatment groups: control (T0), and the other three groups were added different levels of BSL to the water as follows: T1 (0.01 g/m3), T2 (0.02 g/m3), and T3 (0.03 g/m3), respectively. Results indicate that BSL administration significantly improved RBCs, hemoglobin, hematocrit, MCH, and MCHC, with the highest improvement seen in the T3 group (P < 0.05). BSL added to the fish water significantly enhanced serum protein fractions (total protein, albumin, and globulins), while AST, ALT, ALP, creatinine, uric acid, and glucose were significantly diminished in a dose-dependent way (P < 0.05). Adding 0.02-0.03 g/ m3 of BSL resulted in higher antioxidant status (superoxide dismutase and catalase) compared to other groups (P < 0.05). Testosterone levels were higher in T3 than in other groups (P < 0.05). All female hormones (LH, FSH, estradiol, and progesterone) were substantially augmented by the addition of BSL. Additionally, the BSL groups exhibited higher GSI, HSI, VSI (male only), egg diameter (mm), mean number of fry/fish, and mean fry weight (g) compared to the control group (P < 0.05). Expression of reproductive-associated genes (vasa, nanos1a, nanos2, dnd1, pum1, AMH, and vtg) were significantly up-regulated in the gonads of fish in the 0.03 g/m3 treatment. The histological gonadal structure exhibited that BSL improved gonad maturation in both genders of Tilapia fish. Overall, adding a mixture of B. subtilis and B. licheniformis (0.03 g/m3 water) can accelerate reproductive performance in Red Tilapia through up-regulation of reproductive genes and enhance the health profile.

Coping with salinity extremes: Gill transcriptome profiling in the black-chinned tilapia (Sarotherodon melanotheron)

Steeper and sometimes extreme salinity gradients increasingly affect aquatic organisms because of climate change. Hypersalinity habitats demand powerful physiological adaptive strategies. Few teleost species have the capacity to spend their whole life cycle in salinities way over seawater levels. Focusing on the multifunctional gill, we unraveled the tilapia S. melanotheron key strategies to cope with different environmental conditions, ranging from freshwater up to hypersaline habitats. De novo transcriptome assembly based on RNAseq allowed for the analysis of 40,967 annotated transcripts among samples collected in three wild populations at 0, 40 and 80 ‰. A trend analysis of the expression patterns revealed responses across the salinity gradient with different gene pathways involved. Genes linked to ion transport, pH regulation and cell surface receptor signaling were mainly upregulated in the high salinity habitat. We identified tight junction proteins that were critical in high salinity habitats and that were different from the well-known tightening junctional proteins identified and expressed in fresh water. Expression profiles also suggest a change in the vascular tone that could be linked to an osmorespiratory compromise not only in fresh water, but also in high salinity environments. A striking downregulation of genes linked to the immune system and to the heat shock response was observed suggesting an energetic trade-off between immunity and acclimation/adaptation in the hypersaline habitat. The high expression of transcripts coding for immune and heat shock response in the freshwater habitat suggests the establishment of powerful mechanisms to protect gills from environmental threats and to maintain protein integrity. Non-directional expression trends were also detected with an upregulation of genes only in the hypersaline habitat (80 ‰) or only in the marine habitat (40 ‰). Unravel physiological strategies in S. melanotheron populations will help to better understand the molecular basis of fish euryhalinity in salinity-contrasted environments.

slc26a12-A novel member of the slc26 family, is located in tandem with slc26a2 in coelacanths, amphibians, reptiles, and birds

Solute carrier family 26 (Slc26) is a family of anion exchangers with 11 members in mammals (named Slc26a1-a11). Here, we identified a novel member of the slc26 family, slc26a12, located in tandem with slc26a2 in the genomes of several vertebrate lineages. BLAST and synteny analyses of various jawed vertebrate genome databases revealed that slc26a12 is present in coelacanths, amphibians, reptiles, and birds but not in cartilaginous fishes, lungfish, mammals, or ray-finned fishes. In some avian and reptilian lineages such as owls, penguins, egrets, and ducks, and most turtles examined, slc26a12 was lost or pseudogenized. Phylogenetic analysis showed that Slc26a12 formed an independent branch with the other Slc26 members and Slc26a12, Slc26a1 and Slc26a2 formed a single branch, suggesting that these three members formed a subfamily in Slc26. In jawless fish, hagfish have two genes homologous to slc26a2 and slc26a12, whereas lamprey has a single gene homologous to slc26a2. African clawed frogs express slc26a12 in larval gills, skin, and fins. These results show that slc26a12 was present at least before the separation of lobe-finned fish and tetrapods; the name slc26a12 is appropriate because the gene duplication occurred in the distant past.

Genome sequencing, comparative analysis, and gene expression responses of cytochrome P450 genes in Oryzias curvinotus provide insights into environmental adaptation

The mangrove fish (Oryzias curvinotus) serves as a model for researching environmental adaptation and sexual development. To further such research, we sequenced and assembled a high-quality 842 Mb reference genome for O. curvinotus. Comparative genomic analysis revealed 891 expanded gene families, including significantly expanded cytochrome P450 (CYP) detoxification genes known to be involved in xenobiotic defense. We identified 69 O. curvinotus CYPs (OcuCYPs) across 18 families and 10 clans using multiple methods. Extensive RNA-seq and qPCR analysis demonstrated diverse spatiotemporal expression patterns of OcuCYPs by developmental stage, tissue type, sex, and pollutant exposure (17β-estradiol (E2) and testosterone (MT)). Many OcuCYPs exhibited sexual dimorphism in gonads, suggesting reproductive roles in steroidogenesis, while their responsiveness to model toxicants indicates their importance in environmental adaptation through enhanced detoxification. Pathway analysis highlighted expanded CYP genes in arachidonic acid metabolism, drug metabolism, and steroid hormone biosynthesis. This chromosome-level genomic resource provides crucial biological insights to elucidate the functional roles of expanded CYPs in environmental adaptation, sexual development, early life history, and conservation in the anthropogenically impacted mangrove habitats of O. curvinotus. It also enables future ecotoxicology research leveraging O. curvinotus as a pollution sentinel species.

Genomic insights into the seawater adaptation in Cyprinidae

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Whole genome assembly and annotation of the King Angelfish (Holacanthus passer) gives insight into the evolution of marine fishes of the Tropical Eastern Pacific

Holacanthus angelfishes are some of the most iconic marine fishes of the Tropical Eastern Pacific (TEP). However, very limited genomic resources currently exist for the genus. In this study we: (i) assembled and annotated the nuclear genome of the King Angelfish (Holacanthus passer), and (ii) examined the demographic history of H. passer in the TEP. We generated 43.8 Gb of ONT and 97.3 Gb Illumina reads representing 75× and 167× coverage, respectively. The final genome assembly size was 583 Mb with a contig N50 of 5.7 Mb, which captured 97.5% of the complete Actinoterygii Benchmarking Universal Single-Copy Orthologs (BUSCOs). Repetitive elements accounted for 5.09% of the genome, and 33,889 protein-coding genes were predicted, of which 22,984 were functionally annotated. Our demographic analysis suggests that population expansions of H. passer occurred prior to the last glacial maximum (LGM) and were more likely shaped by events associated with the closure of the Isthmus of Panama. This result is surprising, given that most rapid population expansions in both freshwater and marine organisms have been reported to occur globally after the LGM. Overall, this annotated genome assembly provides a novel molecular resource to study the evolution of Holacanthus angelfishes, while facilitating research into local adaptation, speciation, and introgression in marine fishes.

De Novo Genome Assembly of the Whitespot Parrotfish (Scarus forsteni): A Valuable Scaridae Genomic Resource

Scarus forsteni, a whitespot parrotfish from the Scaridae family, is a herbivorous fish inhabiting coral reef ecosystems. The deterioration of coral reefs has highly affected the habitats of the parrotfish. The decline in genetic diversity of parrotfish emphasizes the critical importance of conserving their genetic variability to ensure the resilience and sustainability of marine ecosystems for future generations. In this study, a genome of S. forsteni was assembled de novo through using Illumina and Nanopore sequencing. The 1.71-Gb genome of S. forsteni, was assembled into 544 contigs (assembly level: contig). It exhibited an N50 length of 17.97 Mb and a GC content percentage of 39.32%. Our BUSCO analysis revealed that the complete protein of the S. forsteni genome had 98.10% integrity. Combined with structure annotation data, 34,140 (74.81%) genes were functionally annotated out of 45,638 predicted protein-coding genes. Upon comparing the genome size and TE content of teleost fishes, a roughly linear relationship was observed between these two parameters. However, TE content is not a decisive factor in determining the genome size of S. forsteni. Population history analysis results indicate that S. forsteni experienced two major population expansions, both of which occurred before the last interglacial period. In addition, through a comparative genomic analysis of the evolutionary relationship of other species, it was found that S. forsteni had the closest relationship with Cheilinus undulatus, another member of the Labridae family. Our expansion and contraction analysis of the gene family showed that the expansion genes were mainly associated with immune diseases, organismal systems, and cellular processes. At the same time, cell transcription and translation, sex hormone regulation, and other related pathways were also more prominent in the positive selection genes. The genomic sequence of S. forsteni offers valuable resources for future investigations on the conservation, evolution, and behavior of fish species.

Convergent genomic signatures associated with vertebrate viviparity

BACKGROUND

Viviparity-live birth-is a complex and innovative mode of reproduction that has evolved repeatedly across the vertebrate Tree of Life. Viviparous species exhibit remarkable levels of reproductive diversity, both in the amount of care provided by the parent during gestation, and the ways in which that care is delivered. The genetic basis of viviparity has garnered increasing interest over recent years; however, such studies are often undertaken on small evolutionary timelines, and thus are not able to address changes occurring on a broader scale. Using whole genome data, we investigated the molecular basis of this innovation across the diversity of vertebrates to answer a long held question in evolutionary biology: is the evolution of convergent traits driven by convergent genomic changes?

RESULTS

We reveal convergent changes in protein family sizes, protein-coding regions, introns, and untranslated regions (UTRs) in a number of distantly related viviparous lineages. Specifically, we identify 15 protein families showing evidence of contraction or expansion associated with viviparity. We additionally identify elevated substitution rates in both coding and noncoding sequences in several viviparous lineages. However, we did not find any convergent changes-be it at the nucleotide or protein level-common to all viviparous lineages.

CONCLUSIONS

Our results highlight the value of macroevolutionary comparative genomics in determining the genomic basis of complex evolutionary transitions. While we identify a number of convergent genomic changes that may be associated with the evolution of viviparity in vertebrates, there does not appear to be a convergent molecular signature shared by all viviparous vertebrates. Ultimately, our findings indicate that a complex trait such as viviparity likely evolves with changes occurring in multiple different pathways.

A phylogenetics-based nomenclature system for steroid receptors in teleost fishes

Teleost fishes have emerged as tractable models for studying the neuroendocrine regulation of social behavior via molecular genetic techniques, such as CRISPR/Cas9 gene editing. Moreover, teleosts provide an opportunity to investigate the evolution of steroid receptors and their functions, as species within this lineage possess novel steroid receptor paralogs that resulted from a teleost-specific whole genome duplication. Although teleost fishes have grown in popularity as models for behavioral neuroendocrinology, there is not a consistent nomenclature system for steroid receptors and their genes, which may impede a clear understanding of steroid receptor paralogs and their functions. Here, we used a phylogenetic approach to assess the relatedness of protein sequences encoding steroid receptor paralogs in 18 species from 12 different orders of the Infraclass Teleostei. While most similarly named sequences grouped based on the established phylogeny of the teleost lineage, our analysis revealed several inconsistencies in the nomenclature of steroid receptor paralogs, particularly for sequences encoding estrogen receptor beta (ERβ). Based on our results, we propose a nomenclature system for teleosts in which Greek symbols refer to proteins and numbers refer to genes encoding different subtypes of steroid receptors within the five major groups of this nuclear receptor subfamily. Collectively, our results bridge a critical gap by providing a cohesive naming system for steroid receptors in teleost fishes, which will serve to improve communication, promote collaboration, and enhance our understanding of the evolution and function of steroid receptors across vertebrates.

Effect of high temperatures on sex ratio and differential expression analysis (RNA-seq) of sex-determining genes in Oreochromis niloticus from different river basins in Benin

The high temperature sex reversal process leading to functional phenotypic masculinization during development has been widely described in Nile tilapia (Oreochromis n iloticus) under laboratory or aquaculture conditions and in the wild. In this study, we selected five wild populations of O. niloticus from different river basins in Benin and produced twenty full-sib families of mixed-sex (XY and XX) by natural reproduction. Progenies were exposed to room temperature or high (36.5°C) temperatures between 10 and 30 days post-fertilization (dpf). In control groups, we observed sex ratios from 40% to 60% males as expected, except for 3 families from the Gobé region which showed a bias towards males. High temperature treatment significantly increased male rates in each family up to 88%. Transcriptome analysis was performed by RNA-sequencing (RNA-seq) on brains and gonads from control and treated batches of six families at 15 dpf and 40 dpf. Analysis of differentially expressed genes, differentially spliced genes, and correlations with sex reversal was performed. In 40 dpf gonads, genes involved in sex determination such as dmrt1, cyp11c1, amh, cyp19a1b, ara, and dax1 were upregulated. In 15 dpf brains, a negative correlation was found between the expression of cyp19a1b and the reversal rate, while at 40 dpf a negative correlation was found between the expression of foxl2, cyp11c1, and sf1 and positive correlation was found between dmrt1 expression and reversal rate. Ontology analysis of the genes affected by high temperatures revealed that male sex differentiation processes, primary male sexual characteristics, autophagy, and cilium organization were affected. Based on these results, we conclude that sex reversal by high temperature treatment leads to similar modifications of the transcriptomes in the gonads and brains in offspring of different natural populations of Nile tilapia, which thus may activate a common cascade of reactions inducing sex reversal in progenies.

Tilapia, a good model for studying reproductive endocrinology

The Nile tilapia (Oreochromis niloticus), with a system of XX/XY sex determination, is a worldwide farmed fish with a shorter sexual maturation time than that of most cultured fish. Tilapia show a spawning cycle of approximately 14 days and can be artificially propagated in the laboratory all year round to obtain genetically all female (XX) and all male (XY) fry. Its genome sequence has been opened, and a perfect gene editing platform has been established. With a moderate body size, it is convenient for taking enough blood to measure hormone level. In recent years, using tilapia as animal model, we have confirmed that estrogen is crucial for female development because 1) mutation of star2, cyp17a1 or cyp19a1a (encoding aromatase, the key enzyme for estrogen synthesis) results in sex reversal (SR) due to estrogen deficiency in XX tilapia, while mutation of star1, cyp11a1, cyp17a2, cyp19a1b or cyp11c1 affects fertility due to abnormal androgen, cortisol and DHP levels in XY tilapia; 2) when the estrogen receptors (esr2a/esr2b) are mutated, the sex is reversed from female to male, while when the androgen receptors are mutated, the sex cannot be reversed; 3) the differentiated ovary can be transdifferentiated into functional testis by inhibition of estrogen synthesis, and the differentiated testis can be transdifferentiated into ovary by simultaneous addition of exogenous estrogen and androgen synthase inhibitor; 4) loss of male pathway genes amhy, dmrt1, gsdf causes SR with upregulation of cyp19a1a in XY tilapia. Disruption of estrogen synthesis rescues the male to female SR of amhy and gsdf but not dmrt1 mutants; 5) mutation of female pathway genes foxl2 and sf-1 causes SR with downregulation of cyp19a1a in XX tilapia; 6) the germ cell SR of foxl3 mutants fails to be rescued by estrogen treatment, indicating that estrogen determines female germ cell fate through foxl3. This review also summarized the effects of deficiency of other steroid hormones, such as androgen, DHP and cortisol, on fish reproduction. Overall, these studies demonstrate that tilapia is an excellent animal model for studying reproductive endocrinology of fish.

DNA hydroxymethylation differences underlie phenotypic divergence of somatic growth in Nile tilapia reared in common garden

Phenotypic plasticity of metabolism and growth are essential for adaptation to new environmental conditions, such as those experienced during domestication. Epigenetic regulation plays a key role in this process but the underlying mechanisms are poorly understood, especially in the case of hydroxymethylation. Using reduced representation 5-hydroxymethylcytosine profiling, we compared the liver hydroxymethylomes in full-sib Nile tilapia with distinct growth rates (3.8-fold difference) and demonstrated that DNA hydroxymethylation is strongly associated with phenotypic divergence of somatic growth during the early stages of domestication. The 2677 differentially hydroxymethylated cytosines between fast- and slow-growing fish were enriched within gene bodies (79%), indicating a pertinent role in transcriptional regulation. Moreover, they were found in genes involved in biological processes related to skeletal system and muscle structure development, and there was a positive association between somatic growth and 5hmC levels in genes coding for growth factors, kinases and receptors linked to myogenesis. Single nucleotide polymorphism analysis revealed no genetic differentiation between fast- and slow-growing fish. In addition to unveiling a new link between DNA hydroxymethylation and epigenetic regulation of growth in fish during the initial stages of domestication, this study suggests that epimarkers may be applied in selective breeding programmes for superior phenotypes.

Anemonefishes: A model system for evolutionary genomics

Anemonefishes are an iconic group of coral reef fish particularly known for their mutualistic relationship with sea anemones. This mutualism is especially intriguing as it likely prompted the rapid diversification of anemonefish. Understanding the genomic architecture underlying this process has indeed become one of the holy grails of evolutionary research in these fishes. Recently, anemonefishes have also been used as a model system to study the molecular basis of highly complex traits such as color patterning, social sex change, larval dispersal and life span. Extensive genomic resources including several high-quality reference genomes, a linkage map, and various genetic tools have indeed enabled the identification of genomic features controlling some of these fascinating attributes, but also provided insights into the molecular mechanisms underlying adaptive responses to changing environments. Here, we review the latest findings and new avenues of research that have led to this group of fish being regarded as a model for evolutionary genomics.

Nutritional programming in Nile tilapia (Oreochromis niloticus): Effect of low dietary protein on growth and the intestinal microbiome and transcriptome

Nutritional programming is the idea that early nutrient contributions can influence organismal structure or function and is documented in a variety of vertebrates, yet studies in fish are largely lacking. Tilapia are an important foodfish, with global production having increased rapidly since the 1990s. They exhibit high disease-resistance and grow well on formulated feeds which makes them an ideal aquaculture species, however incorporating high quality proteins into feeds can be costly. As feed constitutes 50-70% of total production costs in aquaculture, reducing protein content could curb these costs and increase revenue. Thus, we examined the effects of feeding Nile tilapia (O. niloticus) fry a restricted protein diet for the first 7-21 days on growth, gut microbial flora, and the intestinal transcriptome. Fish were fed either a 25% restricted or 48% control crude protein starter (ST) diet for up to 21 days and then switched to a 25% or 38% control crude protein growout (GO) diet. Fish fed a 25% ST diet for 14 days followed by a 38% GO diet had significantly higher lengths and weights and better feed efficiency than fish fed the control 48% ST and 38% GO diet after 56 days of culture. Growth of fry on the 25% ST, 7-day/38% GO and the 25% ST,7-day/25% GO diets did not differ from the those fed the control protein diets, while fish fed the 25% ST diet for 21 days had significantly lower growth and survival rates. We observed no significant differences in either alpha or beta diversity of the gut microbial flora between diets, however species richness (Shannon Index) was higher in fry fed the 25% protein ST diet regardless of the GO diet. Similarly, fish fed the 25% ST diet for 14 days followed by the 38% GO diet had minimal changes to the intestinal transcriptome relative to fish fed the control 48% ST and 38% GO diet. However, those fed 25% ST and GO diets for the entire 56 days exhibited substantial differences in the gut transcriptome from other groups showing gene expression profiles characteristic of detrimental changes to gut physiology, protein metabolism and immune function. Results suggest protein restriction for up to 14 days early in development leads to enhanced growth and feed efficiency with minimal effects on gut microbes or intestinal function. Protein restriction beyond this period appears detrimental to fish growth and health as underscored by expression of disease related genes and higher mortality rates.

Mapping of a major QTL for increased robustness and detection of genome assembly errors in Asian seabass (Lates calcarifer)

BACKGROUND

For Asian seabass (Lates calcarifer, Bloch 1790) cultured at sea cages various aquatic pathogens, complex environmental and stress factors are considered as leading causes of disease, causing tens of millions of dollars of annual economic losses. Over the years, we conducted farm-based challenges by exposing Asian seabass juveniles to complex natural environmental conditions. In one of these challenges, we collected a total of 1,250 fish classified as either 'sensitive' or 'robust' individuals during the 28-day observation period.

RESULTS

We constructed a high-resolution linkage map with 3,089 SNPs for Asian seabass using the double digest Restriction-site Associated DNA (ddRAD) technology and a performed a search for Quantitative Trait Loci (QTL) associated with robustness. The search detected a major genome-wide significant QTL for increased robustness in pathogen-infected marine environment on linkage group 11 (ASB_LG11; 88.9 cM to 93.6 cM) with phenotypic variation explained of 81.0%. The QTL was positioned within a > 800 kb genomic region located at the tip of chromosome ASB_LG11 with two Single Nucleotide Polymorphism markers, R1-38468 and R1-61252, located near to the two ends of the QTL. When the R1-61252 marker was validated experimentally in a different mass cross population, it showed a statistically significant association with increased robustness. The majority of thirty-six potential candidate genes located within the QTL have known functions related to innate immunity, stress response or disease. By utilizing this ddRAD-based map, we detected five mis-assemblies corresponding to four chromosomes, namely ASB_LG8, ASB_LG9, ASB_LG15 and ASB_LG20, in the current Asian seabass reference genome assembly.

CONCLUSION

According to our knowledge, the QTL associated with increased robustness is the first such finding from a tropical fish species. Depending on further validation in other stocks and populations, it might be potentially useful for selecting robust Asian seabass lines in selection programs.

Chromatin accessibility in gill tissue identifies candidate genes and loci associated with aquaculture relevant traits in tilapia

The Nile tilapia (Oreochromis niloticus) accounts for ∼9% of global freshwater finfish production however, extreme cold weather and decreasing freshwater resources has created the need to develop resilient strains. By determining the genetic bases of aquaculture relevant traits, we can genotype and breed desirable traits into farmed strains. We generated ATAC-seq and gene expression data from O. niloticus gill tissues, and through the integration of SNPs from 27 tilapia species, identified 1168 highly expressed genes (4% of all Nile tilapia genes) with highly accessible promoter regions with functional variation at transcription factor binding sites (TFBSs). Regulatory variation at these TFBSs is likely driving gene expression differences associated with tilapia gill adaptations, and differentially segregate in freshwater and euryhaline tilapia species. The generation of novel integrative data revealed candidate genes e.g., prolactin receptor 1 and claudin-h, genetic relationships, and loci associated with aquaculture relevant traits like salinity and osmotic stress acclimation.

Evolution of MHC IIB Diversity Across Cichlid Fish Radiations

The genes of the major histocompatibility complex (MHC) are among the most polymorphic genes in vertebrates and crucial for their adaptive immune response. These genes frequently show inconsistencies between allelic genealogies and species phylogenies. This phenomenon is thought to be the result of parasite-mediated balancing selection maintaining ancient alleles through speciation events (trans-species polymorphism [TSP]). However, allele similarities may also arise from postspeciation mechanisms, such as convergence or introgression. Here, we investigated the evolution of MHC class IIB diversity in the cichlid fish radiations across Africa and the Neotropics by a comprehensive review of available MHC IIB DNA sequence information. We explored what mechanism explains the MHC allele similarities found among cichlid radiations. Our results showed extensive allele similarity among cichlid fish across continents, likely due to TSP. Functionality at MHC was also shared among species of the different continents. The maintenance of MHC alleles for long evolutionary times and their shared functionality may imply that certain MHC variants are essential in immune adaptation, even in species that diverged millions of years ago and occupy different environments.

Simplified Genomic Data Revealing the Decline of Aleuritopteris grevilleoides Population Accompanied by the Uplift of Dry-Hot Valley in Yunnan, China

Understanding the evolutionary history of endangered species is crucial for identifying the main reasons for species endangerment in the past and predicting the changing trends and evolutionary directions of their future distribution. In order to study the impact of environmental changes caused by deep valley incision after the uplift of the Qinghai-Tibet Plateau on endangered species, we collected 23 samples belonging to four populations of Aleuritopteris grevilleoides, an endangered fern endemic to the dry-hot valleys (DHV) of Yunnan. Single-nucleotide variation sites (SNPs) were obtained by the genotyping-by-sequencing (GBS) method, and approximately 8085 SNP loci were identified. Through the reconstruction and analysis of genetic diversity, population structure, population dynamics, evolution time, and ancestral geographical distribution, combined with geological historical events such as the formation of dry-hot valleys, this study explores the formation history, current situation, reasons for endangerment and scientifically sound measures for the protection of A. grevilleoides. In our study, A. grevilleoides had low genetic diversity (Obs_Het = 0.16, Exp_Het = 0.32, Pi = 0.33) and a high inbreeding coefficient (Fis = 0.45). The differentiation events were 0.18 Mya, 0.16 Mya, and 0.11 Mya in the A. grevilleoides and may have been related to the formation of terraces within the dry-hot valleys. The history of population dynamics results shows that the diversion of the river resulted in a small amount of gene flow between the two clades, accompanied by a rapid increase in the population at 0.8 Mya. After that, the effective population sizes of A. grevilleoides began to contract continuously due to topographic changes resulting from the continuous expansion of dry-hot valleys. In conclusion, we found that the environmental changes caused by geological events might be the main reason for the changing population size of A. grevilleoides.

No evidence for sex chromosomes in natural populations of the cichlid fish Astatotilapia burtoni

Sex determination (SD) is not conserved among teleost fishes and can even differ between populations of the same species. Across the outstandingly species-rich fish family Cichlidae, more and more SD systems are being discovered. Still, the picture of SD evolution in this group is far from being complete. Lake Tanganyika and its affluent rivers are home to Astatotilapia burtoni, which belongs to the extremely successful East African cichlid lineage Haplochromini. Previously, in different families of an A. burtoni laboratory strain, an XYW system and an XY system have been described. The latter was also found in a second laboratory strain. In a laboratory-reared family descending from a population of the species' southern distribution, a second XY system was discovered. Yet, an analysis of sex chromosomes for the whole species distribution is missing. Here, we examined the genomes of 11 natural populations of A. burtoni, encompassing a wide range of its distribution, for sex-linked regions. We did not detect signs of differentiated sex chromosomes and also not the previously described sex chromosomal systems present in laboratory lines, suggesting different SD systems in the same species under natural and (long-term) artificial conditions. We suggest that SD in A. burtoni is more labile than previously assumed and consists of a combination of non-genetic, polygenic, or poorly differentiated sex chromosomes.

Genome-wide hydroxymethylation profiles in liver of female Nile tilapia with distinct growth performance

The mechanisms underlying the fast genome evolution that occurs during animal domestication are poorly understood. Here, we present a genome-wide epigenetic dataset that quantifies DNA hydroxymethylation at single nucleotide resolution among full-sib Nile tilapia (Oreochromis niloticus) with distinct growth performance. In total, we obtained 355 million, 75 bp reads from 5 large- and 5 small-sized fish on an Illumina NextSeq500 platform. We identified several growth-related genes to be differentially hydroxymethylated, especially within gene bodies and promoters. Previously, we proposed that DNA hydroxymethylation greatly affects the earliest responses to adaptation and potentially drives genome evolution through its targeted enrichment and elevated nucleotide transversion rates. This dataset can be analysed in various contexts (e.g., epigenetics, evolution and growth) and compared to other epigenomic datasets in the future, namely DNA methylation and histone modifications. With forthcoming advancements in genome research, this hydroxymethylation dataset will also contribute to better understand the epigenetic regulation of key genomic features, such as cis-regulatory and transposable elements.

East African cichlid fishes

Cichlid fishes are a very diverse and species-rich family of teleost fishes that inhabit lakes and rivers of India, Africa, and South and Central America. Research has largely focused on East African cichlids of the Rift Lakes Tanganyika, Malawi, and Victoria that constitute the biodiversity hotspots of cichlid fishes. Here, we give an overview of the study system, research questions, and methodologies. Research on cichlid fishes spans many disciplines including ecology, evolution, physiology, genetics, development, and behavioral biology. In this review, we focus on a range of organismal traits, including coloration phenotypes, trophic adaptations, appendages like fins and scales, sensory systems, sex, brains, and behaviors. Moreover, we discuss studies on cichlid phylogenies, plasticity, and general evolutionary patterns, ranging from convergence to speciation rates and the proximate and ultimate mechanisms underlying these processes. From a methodological viewpoint, the last decade has brought great advances in cichlid fish research, particularly through the advent of affordable deep sequencing and advances in genetic manipulations. The ability to integrate across traits and research disciplines, ranging from developmental biology to ecology and evolution, makes cichlid fishes a fascinating research system.

Genome-wide association study for Streptococcus iniae in Nile tilapia (Oreochromis niloticus) identifies a significant QTL for disease resistance

Streptococcus iniae is a problematic gram-positive bacterium negatively affecting Nile tilapia (Oreochromis niloticus), one of the main aquacultural species produced worldwide. The aim of this study was to identify the genetic architecture of survival to S. iniae and identify single nucleotide polymorphism (SNPs) linked to quantitative trait loci (QTL) related to survival to S. iniae challenge. With this purpose, Nile tilapia from the Spring Genetics breeding program were sent to a controlled S. iniae challenge test where phenotypes were scored as dead for fish that died during challenge test and survivors for the fish alive at the termination of the test. Additionally, fin-clip samples from all fish in the test were collected for DNA extraction. Out of 1904 fish in the challenge test, tissue samples of 321 fish were sent for genotyping using double digest restriction site associated DNA sequencing (ddRADseq). After quality control and filtering, 9,085 SNPs were used to perform a genome-wide association study (GWAS). A significant signal in LG8 was observed indicating association with survival to S. iniae challenge, with SNPs explaining from 12% to 26% of the genetic variance. To demonstrate the usefulness of marker assisted selection (MAS) to selectively breed fish for survival to S. iniae, offspring of breeding candidates classified as "resistant" and "susceptible" based on haplotypes of the four most significant markers were sent to a controlled S. iniae challenge test. At the end of the test, the differences in mortality between the two groups were strikingly different with a final cumulative percent mortality of less than 1% and 73% for offspring from "resistant" and "susceptible" parents, respectively. These results demonstrate that MAS for improved resistance to S. iniae is feasible.

Chromosome-level genome sequence of the Genetically Improved Farmed Tilapia (GIFT, Oreochromis niloticus) highlights regions of introgression with O. mossambicus

BACKGROUND

The Nile tilapia (Oreochromis niloticus) is the third most important freshwater fish for aquaculture. Its success is directly linked to continuous breeding efforts focusing on production traits such as growth rate and weight. Among those elite strains, the Genetically Improved Farmed Tilapia (GIFT) programme initiated by WorldFish is now distributed worldwide. To accelerate the development of the GIFT strain through genomic selection, a high-quality reference genome is necessary.

RESULTS

Using a combination of short (10X Genomics) and long read (PacBio HiFi, PacBio CLR) sequencing and a genetic map for the GIFT strain, we generated a chromosome level genome assembly for the GIFT. Using genomes of two closely related species (O. mossambicus, O. aureus), we characterised the extent of introgression between these species and O. niloticus that has occurred during the breeding process. Over 11 Mb of O. mossambicus genomic material could be identified within the GIFT genome, including genes associated with immunity but also with traits of interest such as growth rate.

CONCLUSION

Because of the breeding history of elite strains, current reference genomes might not be the most suitable to support further studies into the GIFT strain. We generated a chromosome level assembly of the GIFT strain, characterising its mixed origins, and the potential contributions of introgressed regions to selected traits.

Meiotic pairing and double-strand break formation along the heteromorphic threespine stickleback sex chromosomes

Double-strand break repair during meiosis is normally achieved using the homologous chromosome as a repair template. Heteromorphic sex chromosomes share little sequence homology, presenting unique challenges to the repair of double-strand breaks. Our understanding of how heteromorphic sex chromosomes behave during meiosis has been focused on ancient sex chromosomes, where the X and Y differ markedly in overall structure and gene content. It remains unclear how more recently evolved sex chromosomes that share considerably more sequence homology with one another pair and form double-strand breaks. One possibility is barriers to pairing evolve rapidly. Alternatively, recently evolved sex chromosomes may exhibit pairing and double-strand break repair that more closely resembles that of their autosomal ancestors. Here, we use the recently evolved X and Y chromosomes of the threespine stickleback fish (Gasterosteus aculeatus) to study patterns of pairing and double-stranded break formation using molecular cytogenetics. We found that the sex chromosomes of threespine stickleback fish did not pair exclusively in the pseudoautosomal region. Instead, the chromosomes fully paired in a non-homologous fashion. To achieve this, the X chromosome underwent synaptic adjustment during pachytene to match the axis length of the Y chromosome. Double-strand break formation and repair rate also matched that of the autosomes. Our results highlight that recently evolved sex chromosomes exhibit meiotic behavior that is reminiscent of autosomes and argues for further work to identify the homologous templates that are used to repair double-strand breaks on the X and Y chromosomes.

Genetic and epigenetic regulation of growth, reproduction, disease resistance and stress responses in aquaculture

Major progress has been made with genomic and genetic studies in aquaculture in the last decade. However, research on epigenetic regulation of aquaculture traits is still at an early stage. It is apparent that most, if not all, aquaculture traits are regulated at both genetic and epigenetic levels. This paper reviews recent progress in understanding of genetic and epigenetic regulation of important aquaculture traits such as growth, reproduction, disease resistance, and stress responses. Although it is challenging to make generalized statements, DNA methylation is mostly correlated with down-regulation of gene expression, especially when at promoters and enhancers. As such, methylation of growth factors and their receptors is negatively correlated with growth; hypomethylation of genes important for stress tolerance is correlated with increased stress tolerance; hypomethylation of genes important for male or female sex differentiation leads to sex differentiation into males or females, respectively. It is apparent that environmental regulation of aquaculture traits is mediated at the level of epigenetic regulation, and such environment-induced epigenetic changes appeared to be intergenerationally inherited, but evidences for transgenerational inheritance are still limited.

Population genomics of an icefish reveals mechanisms of glacier-driven adaptive radiation in Antarctic notothenioids

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.

Symmetric subgenomes and balanced homoeolog expression stabilize the establishment of allopolyploidy in cyprinid fish

BACKGROUND

Interspecific postzygotic reproduction isolation results from large genetic divergence between the subgenomes of established hybrids. Polyploidization immediately after hybridization may reset patterns of homologous chromosome pairing and ameliorate deleterious genomic incompatibility between the subgenomes of distinct parental species in plants and animals. However, the observation that polyploidy is less common in vertebrates raises the question of which factors restrict its emergence. Here, we perform analyses of the genome, epigenome, and gene expression in the nascent allotetraploid lineage (2.95 Gb) derived from the intergeneric hybridization of female goldfish (Carassius auratus, 1.49 Gb) and male common carp (Cyprinus carpio, 1.42 Gb), to shed light on the changes leading to the stabilization of hybrids.

RESULTS

We firstly identify the two subgenomes derived from the parental lineages of goldfish and common carp. We find variable unequal homoeologous recombination in somatic and germ cells of the intergeneric F1 and allotetraploid (F22 and F24) populations, reflecting high plasticity between the subgenomes, and rapidly varying copy numbers between the homoeolog genes. We also find dynamic changes in transposable elements accompanied by genome merger and duplication in the allotetraploid lineage. Finally, we observe the gradual decreases in cis-regulatory effects and increases in trans-regulatory effects along with the allotetraploidization, which contribute to increases in the symmetrical homoeologous expression in different tissues and developmental stages, especially in early embryogenesis.

CONCLUSIONS

Our results reveal a series of changes in transposable elements, unequal homoeologous recombination, cis- and trans-regulations (e.g. DNA methylation), and homoeologous expression, suggesting their potential roles in mediating adaptive stabilization of regulatory systems of the nascent allotetraploid lineage. The symmetrical subgenomes and homoeologous expression provide a novel way of balancing genetic incompatibilities, providing a new insight into the early stages of allopolyploidization in vertebrate evolution.

Chromosome-level genome assembly of largemouth bass (Micropterus salmoides) using PacBio and Hi-C technologies

The largemouth bass (Micropterus salmoides) has become a cosmopolitan species due to its widespread introduction as game or domesticated fish. Here a high-quality chromosome-level reference genome of M. salmoides was produced by combining Illumina paired-end sequencing, PacBio single molecule sequencing technique (SMRT) and High-through chromosome conformation capture (Hi-C) technologies. Ultimately, the genome was assembled into 844.88 Mb with a contig N50 of 15.68 Mb and scaffold N50 length of 35.77 Mb. About 99.9% assembly genome sequences (844.00 Mb) could be anchored to 23 chromosomes, and 98.03% assembly genome sequences could be ordered and directed. The genome contained 38.19% repeat sequences and 2693 noncoding RNAs. A total of 26,370 protein-coding genes from 3415 gene families were predicted, of which 97.69% were functionally annotated. The high-quality genome assembly will be a fundamental resource to study and understand how M. salmoides adapt to novel and changing environments around the world, and also be expected to contribute to the genetic breeding and other research.

CRP Involved in Nile Tilapia (Oreochromis niloticus) against Bacterial Infection

Simple Summary

C-reactive protein (CRP) is an acute-phase protein that can be used as an early diagnostic marker for inflammation. Few CRPs have been isolated from teleost, and the specific immunological functions and mechanism of fish CRP have not been well-studied. Therefore, in this research, a CRP gene from Nile tilapia was identified, and its roles during bacterial infection were investigated. The current results revealed that CRP participated in anti-bacterial immune response through agglutinating bacterial, regulating phagocytosis and inflammation. Hopefully, our data might be beneficial in further study to understand the protective mechanism of fish CRP against bacterial infection.

Abstract

C-reactive protein (CRP) is an acute-phase protein that can be used as an early diagnostic marker for inflammation, which is also an evolutionarily conserved protein and has been identified from arthropods to mammals. However, the roles of CRP during the immune response of Nile tilapia (Oreochromis niloticus) remain unclear. In this study, a CRP gene from Nile tilapia (On-CRP) was identified, and its roles in response to bacterial infection were investigated in vivo or in vitro. On-CRP was found to contain an open reading frame of 675 bp, encoding a polypeptide of 224 amino acids with the conservative pentraxin domain. On-CRP shares more than 50% of its identity with other fish species, and 30% of its identity with mammals. The transcriptional level of On-CRP was most abundant in the liver and its transcripts can be remarkably induced following Streptococcus agalactiae and Aeromonas hydrophila infection. Furthermore, in vitro analysis indicated that the recombinant protein of On-CRP improved phagocytic activity of monocytes/macrophages, and possessed a bacterial agglutination activity in a calcium-dependent manner. Both in vivo and in vitro experiments indicated that On-CRP could promote inflammation and activate the complement pathway. However, a direct relationship between CRP and several immune pathways could not be confirmed. The present data lays a theoretical foundation to further explore the mechanism of how CRP protects fish against bacterial infection.

Comparing Transcriptomes Reveals Key Metabolic Mechanisms in Superior Growth Performance Nile Tilapia (Oreochromis niloticus)

Metabolic capacity is intrinsic to growth performance. To investigate superior growth performance in Nile tilapia, three full-sib families were bred and compared at the biochemical and transcriptome levels to determine metabolic mechanisms involved in significant growth differences between individuals under the same culture environment and feeding regime. Biochemical analysis showed that individuals in the higher growth group had significantly higher total protein, total triglyceride, total cholesterol, and high- and low-density lipoproteins, but significantly lower glucose, as compared with individuals in the lower growth group. Comparative transcriptome analysis showed 536 differentially expressed genes (DEGs) were upregulated, and 622 DEGs were downregulated. These genes were significantly enriched in three key pathways: the tricarboxylic acid cycle (TCA cycle), fatty acid biosynthesis and metabolism, and cholesterol biosynthesis and metabolism. Conjoint analysis of these key pathways and the biochemical parameters suggests that Nile tilapia with superior growth performance have higher ability to consume energy substrates (e.g., glucose), as well as higher ability to biosynthesize fatty acids and cholesterol. Additionally, the fatty acids biosynthesized by the superior growth performance individuals were less active in the catabolic pathway overall, but were more active in the anabolic pathway, and might be used for triglyceride biosynthesis to store excess energy in the form of fat. Furthermore, the tilapia with superior growth performance had lower ability to convert cholesterol into bile acids, but higher ability to convert it into sterols. We discuss the molecular mechanisms of the three key metabolic pathways, map the pathways, and note key factors that may impact the growth of Nile tilapia. The results provide an important guide for the artificial selection and quality enhancement of superior growth performance in tilapia.

Absence of Figla-like Gene Is Concordant with Femaleness in Cichlids Harboring the LG1 Sex-Determination System

Oreochromis niloticus has been used as a reference genome for studies of tilapia sex determination (SD) revealing segregating genetic loci on linkage groups (LGs) 1, 3, and 23. The master key regulator genes (MKR) underlying the SD regions on LGs 3 and 23 have been already found. To identify the MKR in fish that segregate for the LG1 XX/XY SD-system, we applied short variant discovery within the sequence reads of the genomic libraries of the Amherst hybrid stock, Coptodon zillii and Sarotherodon galilaeus, which were aligned to a 3-Mbp-region of the O. aureus genome. We obtained 66,372 variants of which six were concordant with the XX/XY model of SD and were conserved across these species, disclosing the male specific figla-like gene. We further validated this observation in O. mossambicus and in the Chitralada hybrid stock. Genome alignment of the 1252-bp transcript showed that the figla-like gene’s size was 2664 bp, and that its three exons were capable of encoding 99 amino acids including a 45-amino-acid basic helix–loop–helix domain that is typical of the ovary development regulator—factor-in-the-germline-alpha (FIGLA). In Amherst gonads, the figla-like gene was exclusively expressed in testes. Thus, the figla-like genomic presence determines male fate by interrupting the female developmental program. This indicates that the figla-like gene is the long-sought SD MKR on LG1.

Integration of mRNA and miRNA Profiling Reveals Heterosis in Oreochromis niloticus × O. aureus Hybrid Tilapia

Heterosis is a widespread biological phenomenon in fishes, in which hybrids have superior traits to parents. However, the underlying molecular basis for heterosis remains uncertain. Heterosis in growth and survival rates is apparent in hybrid tilapia (Oreochromis niloticus ♀ × O. aureus ♂). Comparisons of growth and hematological biochemical characteristics and mRNA and miRNA transcriptional analyses were performed in hybrid and parents tilapia stocks to investigate the underlying molecular basis for heterosis. Growth characteristics and hematological glucose and cholesterol parameters were significantly improved in hybrids. Of 3097 differentially expressed genes (DEGs) and 120 differentially expressed miRNAs (DEMs) identified among three stocks (O. niloticus, O. aureus, and hybrids), 1598 DEGs and 62 DEMs were non-additively expressed in hybrids. Both expression level dominance and overdominance patterns occurred for DEGs and DEMs, indicating that dominance and overdominance models are widespread in the transcriptional and post-transcriptional regulation of genes involved in growth, metabolism, immunity, and antioxidant capacity in hybrid tilapia. Moreover, potential negative regulation networks between DEMs and predicted target DEGs revealed that most DEGs from miRNA-mRNA pairs are up-regulated. Dominance and overdominance models in levels of transcriptome and miRNAome facilitate the integration of advantageous parental alleles into hybrids, contributing to heterosis of growth and improved survival. The present study provides new insights into molecular heterosis in hybrid tilapia, advancing our understanding of the complex mechanisms involved in this phenomenon in aquatic animals.

Homozygous Mutation of gsdf Causes Infertility in Female Nile Tilapia (Oreochromis niloticus)

Gonadal somatic cell-derived factor (Gsdf) is a member of the TGF-β superfamily, which exists mainly in fishes. Homozygous gsdf mutations in Japanese medaka and zebrafish resulted in infertile females, and the reasons for their infertility remain unknown. This study presents functional studies of Gsdf in ovary development using CRISPR/Cas9 in Nile tilapia (Oreochromis niloticus). The XX wild type (WT) female fish regularly reproduced from 12 months after hatching (mah), while the XX gsdf-/- female fish never reproduced and were infertile. Histological observation showed that at 24 mah, number of phase IV oocyte in the XX gsdf-/- female fish was significantly lower than that of the WT fish, although their gonadosomatic index (GSI) was similar. However, the GSI of the XX gsdf-/- female at 6 mah was higher than that of the WT. The mutated ovaries were hyperplastic with more phase I oocytes. Transcriptome analysis identified 344 and 51 up- and down-regulated genes in mutants compared with the WT ovaries at 6 mah. Some TGF-β signaling genes that are critical for ovary development in fish were differentially expressed. Genes such as amh and amhr2 were up-regulated, while inhbb and acvr2a were down-regulated in mutant ovaries. The cyp19a1a, the key gene for estrogen synthesis, was not differentially expressed. Moreover, the serum 17β-estradiol (E2) concentrations between XX gsdf-/- and WT were similar at 6 and 24 mah. Results from real-time PCR and immunofluorescence experiments were similar and validated the transcriptome data. Furthermore, Yeast-two-hybrid assays showed that Gsdf interacts with TGF-β type II receptors (Amhr2 and Bmpr2a). Altogether, these results suggest that Gsdf functions together with TGF-β signaling pathway to control ovary development and fertility. This study contributes to knowledge on the function of Gsdf in fish oogenesis.

Identification and characterization of mkk genes and their expression profiles in rainbow trout (Oncorhynchus mykiss) symptomatically or asymptomatically infected with Vibrio anguillarum

Mitogen-activated protein kinase kinases (MKKs) are intermediate kinases of mitogen-activated protein kinases (MAPKs) signaling pathways. MKKs are activated by mitogen-activated protein kinase kinase kinase (MKKK) and then the activated MKKs trigger the activation of downstream MAPKs. MAPK signaling pathways play an important role in regulating immune functions including apoptosis and inflammation. However, studies on identification and characterization of mkk repertoire in rainbow trout (Oncorhynchus mykiss) are still limited. Trout experienced 4 rounds (4R) of whole genome duplication (WGD), thus exhibiting increased paralogs of mkks with potentially functional diversity. In this study, we identified 17 mkk genes in trout and the following bacterial challenge (Vibrio anguillarum) studies showed functional diversity of different mkk subtypes. Vibrio anguillarum infection resulted in significantly up-regulated mkk2 subtypes in spleen and liver, and mkk4b3 in spleen, suggesting immunomodulation was regulated by activation of ERK, p38 and JNK pathways. Compared to other mkk subtypes, mkk6s were down-regulated in symptomatic group, rather than asymptomatic group. The organisms present negative feedback on MAPK activation, thus reducing extra damage to cells. We observed down-regulated mkk6s with up-regulated genes (dusp1 & dusp2) involved in negative feedback of MAPK activation. Based on these results, we might propose the distinct expression patterns of genes associated with MAPK pathways resulted in different phenotypes and symptoms of trout in response to bacterial challenge.

A Chromosome-Level Genome Assembly of Mozambique Tilapia (Oreochromis mossambicus) Reveals the Structure of Sex Determining Regions

The Mozambique tilapia (Oreochromis mossambicus) is a fascinating taxon for evolutionary and ecological research. It is an important food fish and one of the most widely distributed tilapias. Because males grow faster than females, genetically male tilapia are preferred in aquaculture. However, studies of sex determination and sex control in O. mossambicus have been hindered by the limited characterization of the genome. To address this gap, we assembled a high-quality genome of O. mossambicus, using a combination of high coverage of Illumina and Nanopore reads, coupled with Hi-C and RNA-Seq data. Our genome assembly spans 1,007 Mb with a scaffold N50 of 11.38 Mb. We successfully anchored and oriented 98.6% of the genome on 22 linkage groups (LGs). Based on re-sequencing data for male and female fishes from three families, O. mossambicus segregates both an XY system on LG14 and a ZW system on LG3. The sex-patterned SNPs shared by two XY families narrowed the sex determining regions to ∼3 Mb on LG14. The shared sex-patterned SNPs included two deleterious missense mutations in ahnak and rhbdd1, indicating the possible roles of these two genes in sex determination. This annotated chromosome-level genome assembly and identification of sex determining regions represents a valuable resource to help understand the evolution of genetic sex determination in tilapias.

Effects of long-term cadmium exposure on growth, antioxidant defense and DNA methylation in juvenile Nile tilapia (Oreochromis niloticus)

Cadmium (Cd) is a ubiquitous environmental contaminant, posing serious threats to aquatic organisms. The aims of the present study were to investigate the effects of long-term Cd exposure on the growth, GH/IGF axis, antioxidant defense and DNA methylation in juvenile Nile tilapia (Oreochromis niloticus). To this end, juvenile Nile tilapia were exposed to 0, 10 and 50 µg∙L-1 Cd for 45 and 90 days. The obtained results revealed that exposure to high concentrations of Cd significantly decreased body mass and body length, and down-regulated mRNA levels of GHRs, IGF-I and IGF-II in the liver of Nile tilapia. Cd exposure induced oxidative stress including the reduction of antioxidant activities and increases of malondialdehyde (MDA) and 8-hydroxydeoxyguanosine (8-OHdG) contents. Beside, the global DNA methylation levels significantly decreased with increasing Cd concentration and exposure time, which might result from increased oxidative DNA damage, the down-regulated expression of DNMT3a and DNMT3b and up-regulated expression of TET1 and TET2. In conclusion, long-term Cd exposure could inhibit growth, reduce antioxidant capacity and lead to oxidative damages to lipid and DNA, and decrease global DNA methylation level in juvenile Nile tilapia.

Whole-Genome Sequencing and Genome-Wide Studies of Spiny Head Croaker (Collichthys lucidus) Reveals Potential Insights for Well-Developed Otoliths in the Family Sciaenidae

Spiny head croaker (Collichthys lucidus), belonging to the family Sciaenidae, is a small economic fish with a main distribution in the coastal waters of Northwestern Pacific. Here, we constructed a nonredundant chromosome-level genome assembly of spiny head croaker and also made genome-wide investigations on genome evolution and gene families related to otolith development. A primary genome assembly of 811.23 Mb, with a contig N50 of 74.92 kb, was generated by a combination of 49.12-Gb Illumina clean reads and 35.24 Gb of PacBio long reads. Contigs of this draft assembly were further anchored into chromosomes by integration with additional 185.33-Gb Hi-C data, resulting in a high-quality chromosome-level genome assembly of 817.24 Mb, with an improved scaffold N50 of 26.58 Mb. Based on our phylogenetic analysis, we observed that C. lucidus is much closer to Larimichthys crocea than Miichthys miiuy. We also predicted that many gene families were significantly expanded (p-value <0.05) in spiny head croaker; among them, some are associated with "calcium signaling pathway" and potential "inner ear functions." In addition, we identified some otolith-related genes (such as otol1a that encodes Otolin-1a) with critical deletions or mutations, suggesting possible molecular mechanisms for well-developed otoliths in the family Sciaenidae.

Correspondence of aCGH and long-read genome assembly for detection of copy number differences: A proof-of-concept with cichlid genomes

Copy number variation is an important source of genetic variation, yet data are often lacking due to technical limitations for detection given the current genome assemblies. Our goal is to demonstrate the extent to which an array-based platform (aCGH) can identify genomic loci that are collapsed in genome assemblies that were built with short-read technology. Taking advantage of two cichlid species for which genome assemblies based on Illumina and PacBio are available, we show that inter-species aCGH log2 hybridization ratios correlate more strongly with inferred copy number differences based on PacBio-built genome assemblies than based on Illumina-built genome assemblies. With regard to inter-species copy number differences of specific genes identified by each platform, the set identified by aCGH intersects to a greater extent with the set identified by PacBio than with the set identified by Illumina. Gene function, according to Gene Ontology analysis, did not substantially differ among platforms, and platforms converged on functions associated with adaptive phenotypes. The results of the current study further demonstrate that aCGH is an effective platform for identifying copy number variable sequences, particularly those collapsed in short read genome assemblies.

Breeding Strategy Shapes the Composition of Bacterial Communities in Female Nile Tilapia Reared in a Recirculating Aquaculture System

In industrial animal production, breeding strategies are essential to produce offspring of better quality and vitality. It is also known that host microbiome has a bearing on its health. Here, we report for the first time the influence of crossbreeding strategy, inbreeding or outbreeding, on the buccal and intestinal bacterial communities in female Nile tilapia (Oreochromis niloticus). Crossbreeding was performed within a family and between different fish families to obtain the inbred and outbred study groups, respectively. The genetic relationship and structure analysis revealed significant genetic differentiation between the inbred and outbred groups. We also employed a 16S rRNA gene sequencing technique to understand the significant differences between the diversities of the bacterial communities of the inbred and outbred groups. The core microbiota composition in the mouth and the intestine was not affected by the crossbreeding strategy but their abundance varied between the two groups. Furthermore, opportunistic bacteria were abundant in the buccal cavity and intestine of the outbred group, whereas beneficial bacteria were abundant in the intestine of the inbred group. The present study indicates that crossbreeding can influence the abundance of beneficial bacteria, core microbiome and the inter-individual variation in the microbiome.

The effect of benzo[a]pyrene on the gut microbiota of Nile tilapia (Oreochromis niloticus)

Benzo[a]pyrene (BaP) is a highly toxic and carcinogenic polycyclic aromatic hydrocarbon (PAH) whose toxicological effects in the gut microbiota of aquatic organisms have not yet been fully revealed. Therefore, in this study, we used high-throughput 16S rRNA gene sequencing to evaluate the effects of BaP in the gut microbiome of Oreochromis niloticus, including its possible participation in the process of detoxification and its ability to recover. The fish were injected with a single intraperitoneal dose of 20 mg kg^-1 of BaP, and the effects in the microbiome were evaluated at 24, 72, and 120 h post-injection. The results indicate a clear dysbiosis (in composition, relative abundance, diversity, and interaction networks) of the gut microbiota during 24 h post-injection, dominated by Fusobacteria and Bacteroidetes and a decrease in Proteobacteria and Spirochaetae. Interestingly, a slight recovery of the microbiome begins at 72 h and stabilises at 120 h post-injection. Pathway analysis revealed the participation of the gut microbiome in PAH degradation mainly at 24 h post-injection. This study provides new insights in the toxicology of BaP in O. niloticus and the first evidence of the ability of the gut microbiome to recovery after a chemical disturbance. KEY POINTS: • Benzo[a]pyrene caused a dysbiosis in the gut microbiota of Oreochromis niloticus. • We observed an enrichment of bacteria involved in the metabolism of xenobiotics. • The gut microbiota was recovered after exposure to benzo[a]pyrene.