Development of the embryo, larva and early juvenile of Nile tilapia Oreochromis niloticus (Pisces: Cichlidae). Developmental staging system

From hatching to juvenile: Larval development of Vieja fenestrata (Teleostei: Cichlidae)

This study delves into the early development of Vieja fenestrata (Cichlidae), with a specific focus on the description of external morphological and morphometric changes, and growth patterns from hatching to the loss of larval characters under controlled laboratory conditions at a temperature of 28°C. Asynchronous hatching was observed between 58 and 60 h postfertilization, with the posterior body emerging first. Over 14 days, significant morphological, physiological, and behavioral changes were observed, revealing a complex developmental trajectory. The initial developmental phases were characterized by rapid vascularization, fin differentiation, and heightened activity, and the subsequent days witnessed the flexion of the notochord, emergence of swim bladder functionality, and transition to exogenous feeding. Maturation progressed with the absorption of the yolk sac, regression of cement glands, and fin ray development, culminating in metamorphosis by 14 days post-hatching. Throughout this period, evolving pigmentation patterns and structural adaptations highlight the species' adaptive strategies. During the larval period of V. fenestrata, substantial changes in morphological proportions were observed. Before the inflection, tail length, trunk length, and body depth had negative allometric growth, and head length, eye diameter, and snout length had positive allometric growth. After the inflection, body depth and snout length showed positive allometric growth; head length and trunk length exhibited isometric growth, whereas tail length and eye diameter demonstrated negative allometric growth. These findings contribute insights into the intricate developmental dynamics of V. fenestrata. Moreover, further research may explore these developmental dynamics' ecological and evolutionary implications.

Embryonic temperature has long-term effects on muscle circRNA expression and somatic growth in Nile tilapia

Embryonic temperature has a lasting impact on muscle phenotype in vertebrates, involving complex molecular mechanisms that encompass both protein-coding and non-coding genes. Circular RNAs (circRNAs) are a class of regulatory RNAs that play important roles in various biological processes, but the effect of variable thermal conditions on the circRNA transcriptome and its long-term impact on muscle growth plasticity remains largely unexplored. To fill this knowledge gap, we performed a transcriptomic analysis of circRNAs in fast muscle of Nile tilapia (Oreochromis niloticus) subjected to different embryonic temperatures (24°C, 28°C and 32°C) and then reared at a common temperature (28°C) for 4 months. Nile tilapia embryos exhibited faster development and subsequently higher long-term growth at 32°C compared to those reared at 28°C and 24°C. Next-generation sequencing data revealed a total of 5,141 unique circRNAs across all temperature groups, of which 1,604, 1,531, and 1,169 circRNAs were exclusively found in the 24°C, 28°C and 32°C groups, respectively. Among them, circNexn exhibited a 1.7-fold (log2) upregulation in the 24°C group and a 1.3-fold (log2) upregulation in the 32°C group when compared to the 28°C group. Conversely, circTTN and circTTN_b were downregulated in the 24°C groups compared to their 28°C and 32°C counterparts. Furthermore, these differentially expressed circRNAs were found to have multiple interactions with myomiRs, highlighting their potential as promising candidates for further investigation in the context of muscle growth plasticity. Taken together, our findings provide new insights into the molecular mechanisms that may underlie muscle growth plasticity in response to thermal variation in fish, with important implications in the context of climate change, fisheries and aquaculture.

Ontogeny and social context regulate the circadian activity patterns of Lake Malawi cichlids

Activity patterns tend to be highly stereotyped and critical for executing many different behaviors including foraging, social interactions, and predator avoidance. Differences in the circadian timing of locomotor activity and rest periods can facilitate habitat partitioning and the exploitation of novel niches. As a consequence, closely related species often display highly divergent activity patterns, suggesting that shifts from diurnal to nocturnal behavior, or vice versa, are critical for survival. In Africa's Lake Malawi alone, there are over 500 species of cichlids, which inhabit diverse environments and exhibit extensive phenotypic variation. We have previously identified a substantial range in activity patterns across adult Lake Malawi cichlid species, from strongly diurnal to strongly nocturnal. In many species, including fishes, ecological pressures differ dramatically across life-history stages, raising the possibility that activity patterns may change over ontogeny. To determine if rest-activity patterns change across life stages, we compared the locomotor patterns of six Lake Malawi cichlid species. While total rest and activity did not change between early juvenile and adult stages, rest-activity patterns did, with juveniles displaying distinct activity rhythms that are more robust than adults. One distinct difference between juveniles and adults is the emergence of complex social behavior. To determine whether social context is required for activity rhythms, we next measured locomotor behavior in group-housed adult fish. We found that when normal social interactions were allowed, locomotor activity patterns were restored, supporting the notion that social interactions promote circadian regulation of activity in adult fish. These findings reveal a previously unidentified link between developmental stage and social interactions in the circadian timing of cichlid activity.

Water properties and quality of the largest rice production region in Japan and their influence on the reproduction of zebrafish

The health of freshwater aquarium fish and their breeding success depend critically on the quality of tap water. In general, tap water in Japan is potable, although the properties of tap water vary among regions in Japan. The city of Niigata is located in the largest rice production region of Japan. We have faced challenges concerning the reproduction of freshwater aquarium fish in Niigata. To determine whether water properties and quality affect the reproduction of aquarium fish in Niigata, we investigated the chemical properties of water and raised zebrafish in water from three different sources, namely tap water of Niigata in May, artificial freshwater (i.e., prepared via reverse osmosis), and natural spring water of Gosen, to document any effects on their sexual maturation and reproduction. We found that the tap water of Niigata was not stable throughout a year (median electrical conductivity = 147.1 μS/cm; SD = 25.6), with springtime lower than the first quartile. We also found that low concentrations of four pesticides in the tap water have been detected in May (max. concentration in 2020, bromobutide 2,000 ng/L, butachlor 600 ng/L, pyraclonil 200 ng/L, ipfencarbazone 20 ng/L). Moreover, rearing zebrafish in tap water negatively influenced both fish growth and reproduction: The sex ratio of adults was male biased (proportion of F0 male 70.8%); the average total length (30.5 mm) and weight (182 mg) of F0 males was decreased; the GSI of F0 females (9.7%) was decreased; the fecundity (the mating success 58.7%; the number of F1 eggs 63.1) of adults was reduced. Rearing in artificial freshwater could improve these outcomes (the sex ratio 55.7%; the total length of F0 males 31.8 mm; the weight of F0 males 211 mg; the GSI of F0 females 11.7%; the mating success 72.6%; the number of F1 eggs 99.0), whereas rearing in natural spring water from Gosen could improve the sex ratio (56.3%) and the weight of F0 males (200 mg), but not the others. Therefore, artificial freshwater made via reverse osmosis should be used for breeding freshwater aquarium fish in rice production region like Niigata. Finally, our results demonstrate that the reproduction of freshwater aquarium fish can serve as a bioindicator of low levels of organic pollutants in tap water and thus provide a basis for evaluating the safety of tap water for human consumption.

Comparative proteomic profiling represents an inhibition of protein synthesis to regulate osmotic stress in Nile tilapia (Oreochromis niloticus) embryos

Seawater (SW)-acclimated Nile tilapia, Oreochromis niloticus, can tolerate up to 30 g.L-1 SW but rarely produce offspring. The embryos of SW-acclimated O. niloticus survived equally well from 0- to 10-g.L-1 environment but not under 20-g. L-1. However, when the embryos were incubated under 10 g.L-1 during days 0-3, and then the salinity was suddenly shifted to and maintained at 20 g.L-1 during days 4-6, their survival rate was comparable to those incubated under 0 and 10 g.L-1. To elucidate a molecular adaptation of the embryos that survived different salinity environments, the proteomic profiles of the newly hatched embryos, or early larvae, hatched under 0 g.L-1, 10 g.L-1, and those being incubated at 10 g.L-1 during days 0-3 followed at 20 g.L-1 during days 4-6 were compared. Total proteins extracted from the samples were identified with a gel-free shot-gun proteomics approach using the Nile tilapia protein database. The early larvae from the three groups expressed 2295 proteins, and 279 proteins showed statistically different expressions among groups. Downregulation of the 182 proteins in the larvae hatched under 10 and 20 g.L-1 was found to include 22 proteins that are responsible for cellular responses to osmotic stress. This adaptation may be a crucial factor in reducing cellular metabolism and ion transport between the intra- and extra-cellular environment to stabilize cellular osmolality. In addition, some of these proteins suppress cellular damage from oxygen free radicals generated from the osmotic stress. Eighty-seven proteins significantly changed in the larvae hatched under 20 g.L-1 were clustered. Nineteen of the cellular stress response proteins, which were considered to be mortality induction, were described.

Epigenetics and the evolution of form: Experimental manipulation of a chromatin modification causes species-specific changes to the craniofacial skeleton

A central question in biology is the molecular origins of phenotypic diversity. While genetic changes are key to the genotype-phenotype relationship, alterations to chromatin structure and the physical packaging of histone proteins may also be important drivers of vertebrate divergence. We investigate the impact of such an epigenetic mechanism, histone acetylation, within a textbook example of an adaptive radiation. Cichlids of Lake Malawi have adapted diverse craniofacial structures, and here we investigate how histone acetylation influences morphological variation in these fishes. Specifically, we assessed the effect of inhibiting histone deacetylation using the drug trichostatin A (TSA) on developing facial structures. We examined this during three critical developmental windows in two cichlid species with alternate adult morphologies. Exposure to TSA during neural crest cell (NCC) migration and as postmigratory NCCs proliferate in the pharyngeal arches resulted in significant changes in lateral and ventral shape in Maylandia, but not in Tropheops. This included an overall shortening of the head, widening of the lower jaw, and steeper craniofacial profile, all of which are paedomorphic morphologies. In contrast, treatment with TSA during early chondrogenesis did not result in significant morphological changes in either species. Together, these data suggest a sensitivity to epigenetic alterations that are both time- and species-dependent. We find that morphologies are due to nonautonomous or potentially indirect effects on NCC development, including in part a global developmental delay. Our research bolsters the understanding that proper histone acetylation is essential for early craniofacial development and identifies a species-specific robustness to developmental change. Overall, this study demonstrates how epigenetic regulation may play an important role in both generating and buffering morphological variation.

First Investigation of the Optimal Timing of Vaccination of Nile Tilapia (Oreochromis niloticus) Larvae against Streptococcus agalactiae

To investigate early immune responses and explore the optimal vaccination periods, Nile tilapia at 1, 7, 14, 21, 28, 35, and 42 days after yolk sac collapse (DAYC) were immersed in formalin-killed Streptococcus agalactiae vaccine (FKV-SA). A specific IgM was first detected via ELISA in the 21 DAYC larvae (0.108 g) at 336 h after vaccination (hav), whereas in the 28-42 DAYC larvae (0.330-0.580 g), the specific IgM could be initially detected at 24 hav. qRT-PCR analysis of the TCRβ, CD4, MHCIIα, IgHM, IgHT, and IgHD genes in 21-42 DAYC larvae immunized with the FKV-SA immersion route for 24, 168, and 336 hav revealed that the levels of most immune-related genes were significantly higher in the vaccinated larvae at all DAYCs than in the control larvae (p < 0.05) at 336 hav. Immunohistochemistry demonstrated stronger IgM signals in the gills, head kidney, and intestine tissues at 21, 28, and 35 DAYC in all vaccinated larvae compared with the control. Interestingly, at all DAYCs, FKV-SA larvae exhibited significantly higher survival rates and an increased relative percent survival (RPS) than the control after challenge with viable S. agalactiae, particularly in larvae that were immunized with FKV-SA at 168 and 336 hav (p < 0.05).

From egg to slaughter: monitoring the welfare of Nile tilapia, Oreochromis niloticus, throughout their entire life cycle in aquaculture

The primary aim of this study was to comprehensively evaluate the welfare of Nile tilapia (Oreochromis niloticus) throughout their entire life cycle within aquaculture, spanning from reproduction to slaughter. The methodology was structured to identify welfare indicators closely aligned with the principles of animal freedoms defined by the Farm Animal Council, encompassing environmental, health, nutritional, behavioral, and psychological freedom. Notably, psychological freedom was inherently considered within the behavioral and physical analyses of the animals. To accomplish this, an integrative systematic literature review was conducted to define precise indicators and their corresponding reference values for each stage of tilapia cultivation. These reference values were subsequently categorized using a scoring system that assessed the deviation of each indicator from established ideal (score 1), tolerable (score 2), and critical (score 3) ranges for the welfare of the target species. Subsequently, a laboratory experiment was executed to validate the pre-selected health indicators, specifically tailored for the early life stages of tilapia. This test facilitated an assessment of the applicability of these indicators under operational conditions. Building on the insights gained from this experimentation, partial welfare indices (PWIs) were computed for each assessed freedom, culminating in the derivation of a general welfare index (GWI). Mathematical equations were employed to calculate these indices, offering a quantitative and standardized measure of welfare. This approach equips tilapia farmers and processors with the tools necessary for the continuous monitoring and enhancement of their production systems and stimulate the adoption of more sustainable and ethical practices within the tilapia farming.

The importance of thyroid hormone signaling during early development: Lessons from the zebrafish model

The zebrafish is an optimal experimental model to study thyroid hormone (TH) involvement in vertebrate development. The use of state-of-the-art zebrafish genetic tools available for the study of the effect of gene silencing, cell fate decisions and cell lineage differentiation have contributed to a more insightful comprehension of molecular, cellular, and tissue-specific TH actions. In contrast to intrauterine development, extrauterine embryogenesis observed in zebrafish has facilitated a more detailed study of the development of the hypothalamic-pituitary-thyroid axis. This model has also enabled a more insightful analysis of TH molecular actions upon the organization and function of the brain, the retina, the heart, and the immune system. Consequently, zebrafish has become a trendy model to address paradigms of TH-related functional and biomedical importance. We here compilate the available knowledge regarding zebrafish developmental events for which specific components of TH signaling are essential.

Morphological and temporal variation in early embryogenesis contributes to species divergence in Malawi cichlid fishes

The cichlid fishes comprise the largest extant vertebrate family and are the quintessential example of rapid "explosive" adaptive radiations and phenotypic diversification. Despite low genetic divergence, East African cichlids harbor a spectacular intra- and interspecific morphological diversity, including the hyper-variable, neural crest (NC)-derived traits such as coloration and craniofacial skeleton. Although the genetic and developmental basis of these phenotypes has been investigated, understanding of when, and specifically how early, in ontogeny species-specific differences emerge, remains limited. Since adult traits often originate during embryonic development, the processes of embryogenesis could serve as a potential source of species-specific variation. Consequently, we designed a staging system by which we compare the features of embryogenesis between three Malawi cichlid species-Astatotilapia calliptera, Tropheops sp. 'mauve' and Rhamphochromis sp. "chilingali"-representing a wide spectrum of variation in pigmentation and craniofacial morphologies. Our results showed fundamental differences in multiple aspects of embryogenesis that could underlie interspecific divergence in adult adaptive traits. First, we identified variation in the somite number and signatures of temporal variation, or heterochrony, in the rates of somite formation. The heterochrony was also evident within and between species throughout ontogeny, up to the juvenile stages. Finally, the identified interspecific differences in the development of pigmentation and craniofacial cartilages, present at the earliest stages of their overt formation, provide compelling evidence that the species-specific trajectories begin divergence during early embryogenesis, potentially during somitogenesis and NC development. Altogether, our results expand our understanding of fundamental cichlid biology and provide new insights into the developmental origins of vertebrate morphological diversity.

Embryonic and Larval Development of Stinging Catfish, Heteropneustes fossilis, in Relation to Climatic and Water Quality Parameters

In terms of hatchery-based seed production, one of the most important aquaculture species in Bangladesh is the stinging catfish (Heteropneustes fossilis). Scientific and evidence-based embryonic and larval development research on this fish species in the context of climate change is limited. This experimental study was conducted via induced breeding of stinging catfish using a conventional hatchery system, rearing the larvae in hapas placed in ponds. A series of microscopic observations using a trinocular digital microscope and an analysis of the relationship between larval growth and climate-driven water quality parameters such as temperature, pH, dissolved oxygen, total dissolved solids, alkalinity, and ammonia were performed. During embryonic development, the first cleavage was observed between 30 and 35 min of post-fertilization. Embryonic development (ranging from the 2-cell to the pre-hatching stage) took 21:00 h. Hatching occurred at 22:30 to 23:00 h after fertilization, with an average larvae length of 2.78 ± 0.04 mm. In the post-hatching stage, four pairs of tiny barbels appeared at 36:00 h, and the larvae started feeding exogenously after 72:00 h. These larvae fully absorbed their yolk sacs on the 6th day and attained an average length of 6.44 ± 0.06 mm. Aerial respiration of the larvae was investigated through naked-eye observation on the 10th day of hatching. The average length of the larvae was 32.00 ± 2.0 mm at the end of the 30-day post-hatching observation period. Bivariate correlation analysis showed significant correlations between key climatic variables and water quality parameters under hapa-based larval-rearing conditions. According to canonical correlation analysis, the first canonical function revealed the highest significant correlation between the two sets of variables (r1 = 0.791). The response variable weight of larvae (6.607) was linked to two explanatory variables: pH (0.321) and dissolved oxygen (0.265). For the second canonical correlation function, a positive correlation (0.431) was observed between the two sets of variables. Larval weight (-18.304) was observed to be linked to climatic variables, including air temperature (-0.316) and surface pressure (0.338). Results of this study reveal the subtle correlation between larval growth and water quality driven by climatic variables.

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.

Origin, form and function of extraembryonic structures in teleost fishes

Teleost eggs have evolved a highly derived early developmental pattern within vertebrates as a result of the meroblastic cleavage pattern, giving rise to a polar stratified architecture containing a large acellular yolk and a small cellular blastoderm on top. Besides the acellular yolk, the teleost-specific yolk syncytial layer (YSL) and the superficial epithelial enveloping layer are recognized as extraembryonic structures that play critical roles throughout embryonic development. They provide enriched microenvironments in which molecular feedback loops, cellular interactions and mechanical signals emerge to sculpt, among other things, embryonic patterning along the dorsoventral and left-right axes, mesendodermal specification and the execution of morphogenetic movements in the early embryo and during organogenesis. An emerging concept points to a critical role of extraembryonic structures in reinforcing early genetic and morphogenetic programmes in reciprocal coordination with the embryonic blastoderm, providing the necessary boundary conditions for development to proceed. In addition, the role of the enveloping cell layer in providing mechanical, osmotic and immunological protection during early stages of development, and the autonomous nutritional support provided by the yolk and YSL, have probably been key aspects that have enabled the massive radiation of teleosts to colonize every ecological niche on the Earth. This article is part of the theme issue 'Extraembryonic tissues: exploring concepts, definitions and functions across the animal kingdom'.

Testosterone nanoemulsion produced masculinized Nile tilapia (Oreochromis niloticus)

The objective of this work was to develop a food additive for the sex reversal of Nile tilapia (Oreochromis niloticus) based on a simple oil in water (O/W) nanoemulsion with testosterone propionate for incorporation into commercial feed. Oil screening and evaluation of the organoleptic and physicochemical characteristics were carried out to determine the best formulation. A palatability test was also performed. Sex reversal test was assayed using 5 experimental groups: negative control - macerated feed without hormone; free testosterone - macerated feed with 60 mg/kg of testosterone propionate diluted in ethanol; and macerated feed with testosterone propionate nanoemulsion at a concentration of 30, 60, and 90 mg/kg. Stable nanoemulsions (size 76-210 nm) with testosterone propionate were produced. All nanoemulsion-added feed was palatable to tilapia. We obtained sex reversal values of ≈65, 75, and 72% in the groups of 30, 60, and 90 mg/kg, respectively. We can conclude that the nanoemulsion showed promising results; it is capable of inducing sex reversal in tilapia, is suitable as a commercial product, and has the potential to promote safety for rural staff and reduce the environmental impact of hormones.

Multigenerational inspections of environmental thermal perturbations promote metabolic trade-offs in developmental stages of tropical fish

Global warming both reduces global temperature variance and increases the frequency of extreme weather events. In response to these ambient perturbations, animals may be subject to trans- or intra-generational phenotype modifications that help to maintain homeostasis and fitness. Here, we show how temperature-associated transgenerational plasticity in tilapia affects metabolic trade-offs during developmental stages under a global warming scenario. Tropical tilapia reared at a stable temperature of 27 °C for a decade were divided into two temperature-experience groups for four generations of breeding. Each generation of one group was exposed to a single 15 °C cold-shock experience during its lifetime (cold-experienced CE group), and the other group was kept stably at 27 °C throughout their lifetimes (cold-naïve CN group). The offspring at early life stages from the CE and CN tilapia were then assessed by metabolomics-based profiling, and the results implied that parental cold-experience might affect energy provision during reproduction. Furthermore, at early life stages, progeny may be endowed with metabolic traits that help the animals cope with ambient temperature perturbations. This study also applied the feature rescaling and Uniform Manifold Approximation and Projection (UMAP) to visualize metabolic dynamics, and the result could effectively decompose the complex omic-based datasets to represent the energy trade-off variability. For example, the carbohydrate to free amino acid conversion and enhanced compensatory features appeared to be hypothermic-responsive traits. These multigenerational metabolic effects suggest that the tropical ectothermic tilapia may exhibit transgenerational phenotype plasticity, which could optimize energy allocation under ambient temperature challenges. Knowledge about such metabolism-related transgenerational plasticity effects in ectothermic aquatic species may allow us to better predict how adaptive mechanisms will affect fish populations in a climate with narrow temperature variation and frequent extreme weather events.

Nile tilapia TLR3 recruits MyD88 and TRIF as adaptors and is involved in the NF-κB pathway in the immune response

TLR3 plays a crucial role in innate immunity. In the present study, OnTLR3 was identified in the Nile tilapia Oreochromis niloticus, with a conserved LRR domain and a C-terminal TIR domain. OnTLR3 was broadly expressed in all tissues tested, with the highest expression levels in the blood and the lowest in the kidney. TLR3 mRNA could be detected from pharyngula (2.5 dpf) to late larva (8.5 dpf) during embryonic and larval development. Moreover, the expression level of OnTLR3 was clearly altered in all five tissues after Streptococcus agalactiae infection in vivo and could be induced by LPS, poly(I:C), S. agalactiae WC1535 and △CPS in Nile tilapia macrophages. When OnTLR3 was overexpressed in 293 T cells, it was distributed in the cytoplasm and could significantly increase NF-κB activation. The pulldown assays showed that OnTLR3 interacted with both OnMyD88 and OnTRIF. The binding assays revealed the specificity of OnTLR3 for pathogen-associated molecular patterns (PAMPs) and bacteria that included S. agalactiae, Aeromonas hydrophila and poly(I:C), LPS and PGN. Taken together, these findings suggest that OnTLR3, as a pattern recognition receptor (PRR), might play an important role in the immune response to pathogen invasion.

A show of Hands: Novel and conserved expression patterns of teleost hand paralogs during craniofacial, heart, fin, peripheral nervous system and gut development

BACKGROUND

Hand genes are required for the development of the vertebrate jaw, heart, peripheral nervous system, limb, gut, placenta, and decidua. Two Hand paralogues, Hand1 and Hand2, are present in most vertebrates, where they mediate different functions yet overlap in expression. In ray-finned fishes, Hand gene expression and function is only known for the zebrafish, which represents the rare condition of having a single Hand gene, hand2. Here we describe the developmental expression of hand1 and hand2 in the cichlid Copadichromis azureus.

RESULTS

hand1 and hand2 are expressed in the cichlid heart, paired fins, pharyngeal arches, peripheral nervous system, gut, and lateral plate mesoderm with different degrees of overlap.

CONCLUSIONS

Hand gene expression in the gut, peripheral nervous system, and pharyngeal arches may have already been fixed in the lobe- and ray-finned fish common ancestor. In other embryonic regions, such as paired appendages, hand2 expression was fixed, while hand1 expression diverged in lobe- and ray-finned fish lineages. In the lateral plate mesoderm and arch associated catecholaminergic cells, hand1 and hand2 swapped expression between divergent lineages. Distinct expression of cichlid hand1 and hand2 in the epicardium and myocardium of the developing heart may represent the ancestral pattern for bony fishes.

Glucose injection into the yolk influences intermediary metabolism in adult Nile tilapia fed with high levels of carbohydrates

Nutritional programming is a concept proposed to be applied in the field of fish nutrition to improve the use of new diets in aquaculture. This study aimed to investigate for the first time the effects of a glucose injection into the yolk at the alevin stage on intermediary metabolism and growth in adult Nile tilapia (Oreochromis niloticus) at 32-37 weeks later in the life. The early stimulus was performed through direct microinjection of 2 M glucose into yolk sacs of Nile tilapia alevin. Subsequently, in adult tilapia, the long-term effects of glucose stimulus on growth performance, blood metabolites, chemical composition in the liver and muscle, expression of genes involved in glucose transport and metabolism (glycolysis and gluconeogenesis) and related pathways (amino acid catabolism and lipogenesis) were investigated. Our results showed that, even though early glucose injection had no effect on growth performance in adult fish, very few significant effects on glucose metabolism were observed. Furthermore, to evaluate the potential metabolic programming after a dietary challenge, a 2 × 2 factorial design with two early stimuli (0.85% NaCl or 2 M glucose) and two different dietary carbohydrate intakes (medium-carbohydrate diet, CHO-M; high-carbohydrate diet, CHO-H) was performed between weeks 33 and 37. As expected, compared with the CHO-M diet, the CHO-H diet led to decreased growth performance, higher glyceamia and triglyceridemia, higher glycogen and lipid levels in the liver as well as down-regulation of gluconeogenesis and amino acid catabolism gene expressions. More interestingly, although early glucose injection had no significant effect on growth performance, it enhanced the capacities for lipogenesis, glycolysis and gluconeogenesis, particularly in fish that were fed the CHO-H diet. Thus, the nutritional programming of tilapia linked to glucose injection into the yolk of alevins is always visible at the adult stage albeit less intense than what we previously observed in juvenile.

Genome Editing Using the CRISPR-Cas9 System to Generate a Solid-Red Germline of Nile Tilapia (Oreochromis niloticus)

In recent years, there has been increasing demand for red tilapia, which are commercial strains of hybrids of different tilapiine species or red variants of highly inbred Nile tilapia. However, red tilapia phenotypes are genetically unstable and affected by environmental factors, resulting in nonuniform coloration with black or dark-red color blotches that reduce their market value. Solute carrier family 45 member 2 (SLC45A2) is a membrane transporter that mediates melanin biosynthesis and is evolutionarily conserved from fish to humans. In the present study, we describe the generation of a stable and heritable red tilapia phenotype by inducing loss-of-function mutations in the slc45a2 gene. For this purpose, we identified the slc45a2 gene in Nile tilapia and designed highly specific guide RNAs (gRNA) for its genomic sequence. Multiplex microinjection of slc45a2-specific ribonucleoproteins to Nile tilapia zygotes induced up to 97-99% albinism, including loss of melanin in the eye. Next-generation sequencing of the injected zygotes demonstrated that all injected fish carried mutant alleles with variable mutagenesis efficiencies. Sanger sequencing of the genomic target region in the slc45a2 gene from fin clips, sperm, and F1 offspring of a highly mutant male identified various genomic indels and germline transmission of the sperm-identified indels. Overall, this work demonstrates the generation of somatic and germline slc45a2 mutant alleles, which leads to complete albinism in Nile tilapia.

Distant hybridization and gynogenesis between Nile tilapia Oreochromis niloticus and Jaguar cichlid Parachromis managuensis

To investigate the distant hybridization and gynogenesis between Nile tilapia Oreochromis niloticus and Jaguar cichlid Parachromis managuensis, reciprocal crossing was first performed between the two species. No offspring, however, were viable when there were these hybridizations. Gynogenesis was induced in O. niloticus and P. managuensis using ultraviolet (UV)-irradiated spermatozoa from P. managuensis and O. niloticus, respectively. The morphology during embryonic development indicated gynogenetic offspring of both O. niloticus and the P. managuensis were normal and deformed, and the results from flow cytometric analysis indicated normal fry were diploid and deformed fry were haploid. Gynogenetic O. niloticus and P. managuensis had the same DNA content and chromosome number as their species of origin, indicating that gynogenetic individuals were produced in both species. The presence of only females for both gynogenetic P. managuensis and O. niloticus was indicative of an XX genotype in the female P. managuensis and O. niloticus. Results from studies on genetic diversity indicated the average heterozygosity of the gynogenetic diploid population of O. niloticus were less than that of the cultured population, but the genetic homozygosity of the gynogenetic diploid population of O. niloticus was greater than that of the cultured population after one generation of gynogenesis, which achieved the goal of rapidly establishing genetic homozygosity.

Stages of embryonic development in the live-bearing fish, Gambusia holbrooki

BACKGROUND

Divergent morphology and placentation of Poeciliids make them suitable model for investigating how evolutionary selection has altered and conserved the developmental mechanisms. However, there is limited description of their embryonic staging, despite representing a key evolutionary node that shares developmental strategy with placental vertebrates. Here, we describe the embryonic developmental stages of Gambusia holbrooki from zygote to parturition using freshly harvested embryos.

RESULTS

We defined 40 embryonic stages using a numbered (stages 0-39; zygote to parturition, respectively) and named (grouped into seven periods, ie, zygote, cleavage, blastula, gastrula, segmentation, pharyngula, and parturition) staging system. Two sets of quantitative (ie, egg diameter, embryonic total length, otic vesicle closure index, heart rates, the number of caudal fin rays and elements) and qualitative (ie, three-dimensional analysis of images and key morphological criteria) data were acquired and used in combination to describe each stage. All 40 stages are separated by well-defined morphological traits, revealing developmental novelties that are influenced by narrow perivitelline space, placentation, internal gestation, and sex differentiation.

CONCLUSIONS

The principal diagnostic features described are quick, reliable, and easy to apply. This system will benefit researchers investigating molecular ontogeny, particularly sexual differentiation mechanisms in G. holbrooki.

Nile Tilapia: A Model for Studying Teleost Color Patterns

The diverse color patterns of cichlid fishes play an important role in mate choice and speciation. Here we develop the Nile tilapia (Oreochromis niloticus) as a model system for studying the developmental genetics of cichlid color patterns. We identified 4 types of pigment cells: melanophores, xanthophores, iridophores and erythrophores, and characterized their first appearance in wild-type fish. We mutated 25 genes involved in melanogenesis, pteridine metabolism, and the carotenoid absorption and cleavage pathways. Among the 25 mutated genes, 13 genes had a phenotype in both the F0 and F2 generations. None of F1 heterozygotes had phenotype. By comparing the color pattern of our mutants with that of red tilapia (Oreochromis spp), a natural mutant produced during hybridization of tilapia species, we found that the pigmentation of the body and eye is controlled by different genes. Previously studied genes like mitf, kita/kitlga, pmel, tyrb, hps4, gch2, csf1ra, pax7b, and bco2b were proved to be of great significance for color patterning in tilapia. Our results suggested that tilapia, a fish with 4 types of pigment cells and a vertically barred wild-type color pattern, together with various natural and artificially induced color gene mutants, can serve as an excellent model system for study color patterning in vertebrates.

Gene coexpression networks reveal molecular interactions underlying cichlid jaw modularity

Background

The oral and pharyngeal jaw of cichlid fishes are a classic example of evolutionary modularity as their functional decoupling boosted trophic diversification and contributed to the success of cichlid adaptive radiations. Most studies until now have focused on the functional, morphological, or genetic aspects of cichlid jaw modularity. Here we extend this concept to include transcriptional modularity by sequencing whole transcriptomes of the two jaws and comparing their gene coexpression networks.

Results

We show that transcriptional decoupling of gene expression underlies the functional decoupling of cichlid oral and pharyngeal jaw apparatus and the two units are evolving independently in recently diverged cichlid species from Lake Tanganyika. Oral and pharyngeal jaw coexpression networks reflect the common origin of the jaw regulatory program as there is high preservation of gene coexpression modules between the two sets of jaws. However, there is substantial rewiring of genetic architecture within those modules. We define a global jaw coexpression network and highlight jaw-specific and species-specific modules within it. Furthermore, we annotate a comprehensive in silico gene regulatory network linking the Wnt and AHR signalling pathways to jaw morphogenesis and response to environmental cues, respectively. Components of these pathways are significantly differentially expressed between the oral and pharyngeal jaw apparatus.

Conclusion

This study describes the concerted expression of many genes in cichlid oral and pharyngeal jaw apparatus at the onset of the independent life of cichlid fishes. Our findings suggest that – on the basis of an ancestral gill arch network—transcriptional rewiring may have driven the modular evolution of the oral and pharyngeal jaws, highlighting the evolutionary significance of gene network reuse. The gene coexpression and in silico regulatory networks presented here are intended as resource for future studies on the genetics of vertebrate jaw morphogenesis and trophic adaptation.

Supplementary Information

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

Influence of temperature on larval ontogenesis of Geophagus brasiliensis (Quoy & Gaimard, 1824) (Pisces: Cichlidae)

Studies on the larval development of fish are essential for conservation and improvements in cultivation techniques. Geophagus brasiliensis popularly known as Cará has potential as a fish of interest in ornamental aquaculture. Wild adults of G. brasiliensis were kept in an aquarium for spontaneous reproduction. Newly hatched larvae were transferred to 5-litre aquaria at 22, 26 and 30°C until total yolk sac resorption. Histological slides were used for biometric analysis and monitoring of larval ontogenesis at different temperatures. Histologically, from the first to the fourth days it was possible to identify myomeres, optic vesicle, yolk syncytial layer, brain, heart and differentiation of the eye layers. From the fourth to the seventh days, it was possible to identify mandibular and gill cartilages, swim bladder, liver, prismatic epithelium with striated border in intestine and renal epithelium. All biometric measurements increased over the days, except height and length of the yolk sac, which gradually decreased until the complete resorption of the yolk sac that occurred on the fifth day at a temperature of 30°C, the sixth day at 26°C and the seventh day at 22°C. Morphological events at 30°C such as the reabsorption of the yolk sac, the appearance of cartilage in the branchial arches and differentiation of the layers of the eyes occurred faster compared with the other temperatures tested. Opening of the mouth and digestive tract occurred at a similar time on the fourth day in all temperatures.

The involvement of Nile tilapia (Oreochromis niloticus) Neu4 sialidase in neural differentiation during early ontogenesis

Neurogenesis is an important process for the formation of the central nervous system during ontogenesis. Mammalian sialidases are involved in neurogenesis through desialylation of sialo-glycoconjugates. However, the significance of fish sialidases, unlike that of mammals, in neurogenesis has not been investigated. The present study focuses on Nile tilapia (Oreochromis niloticus) because of its unique profiles of sialidases related to enzymatic properties, subcellular localization, and tissue-specific gene expression. First, the fish were cultured under aphotic condition, which is known to cause the delayed development of the retina and brain in various fish. Next, we investigate the effect of aphotic condition on the levels of tilapia sialidases. Our results revealed that the tilapia showed a decrease in the number of ganglion cell in the retina. The expression level of neu4 mRNA is up-regulated in the eyes from tilapia reared in Dark accompanied by the increase of retinal differentiation markers. These results indicated that tilapia Neu4 is involved in retinal development in Nile tilapia. Furthermore, we tried to clarify the function of tilapia Neu4 in the neuronal cells using two neuroblast cell lines (SH-SY5Y and Neuro2a cell lines). Tilapia Neu4 decreased sialic acid level of both nuclear glycoproteins as well as glycolipids. Moreover, tilapia Neu4 accelerated neurite formation in both two neural cell lines and, increased the acetylcholinesterase activity, but it did not affect cell proliferation. Collectively, these results suggest that Neu4 accelerates neurite differentiation during ontogenesis in tilapia.

Exploring the Expression of Cardiac Regulators in a Vertebrate Extremophile: The Cichlid Fish Oreochromis (Alcolapia) alcalica

Although it is widely accepted that the cellular and molecular mechanisms of vertebrate cardiac development are evolutionarily conserved, this is on the basis of data from only a few model organisms suited to laboratory studies. Here, we investigate gene expression during cardiac development in the extremophile, non-model fish species, Oreochromis (Alcolapia) alcalica. We first characterise the early development of O. alcalica and observe extensive vascularisation across the yolk prior to hatching. We further investigate heart development by identifying and cloning O. alcalica orthologues of conserved cardiac transcription factors gata4, tbx5, and mef2c for analysis by in situ hybridisation. Expression of these three key cardiac developmental regulators also reveals other aspects of O. alcalica development, as these genes are expressed in developing blood, limb, eyes, and muscle, as well as the heart. Our data support the notion that O. alcalica is a direct-developing vertebrate that shares the highly conserved molecular regulation of the vertebrate body plan. However, the expression of gata4 in O. alcalica reveals interesting differences in the development of the circulatory system distinct from that of the well-studied zebrafish. Understanding the development of O. alcalica embryos is an important step towards providing a model for future research into the adaptation to extreme conditions; this is particularly relevant given that anthropogenic-driven climate change will likely result in more freshwater organisms being exposed to less favourable conditions.

Establishment and characterization of Neu1-knockout zebrafish and its abnormal clinical phenotypes

Mammalian sialidase Neu1 is involved in various physiological functions, including cell adhesion, differentiation, cancer metastasis, and diabetes through lysosomal catabolism and desialylation of glycoproteins at the plasma membrane. Various animal models have been established to further explore the functions of vertebrate Neu1. The present study focused on zebrafish (Danio rerio) belonging to Cypriniformes as an experimental animal model with neu1 gene deficiency. The results revealed that the zebrafish Neu1 desialyzed both α2-3 and α2-6 sialic acid linkages from oligosaccharides and glycoproteins at pH 4.5, and it is highly conserved with other fish species and mammalian Neu1. Furthermore, Neu1-knockout zebrafish (Neu1-KO) was established through CRISPR/Cas9 genome editing. Neu1-KO fish exhibited slight abnormal embryogenesis with the accumulation of pleural effusion; however, no embryonic lethality was observed. Although Neu1-KO fish were able to be maintained as homozygous, they showed smaller body length and weight than the wild-type (WT) fish, and muscle atrophy and curvature of the vertebra were observed in adult Neu1-KO fish (8 months). The expression patterns of myod and myog transcription factors regulating muscle differentiation varied between Neu1-KO and WT fish embryo. Expression of lysosomal-related genes, including ctsa, lamp1a, and tfeb were up-regulated in adult Neu1-KO muscle as compared with WT. Furthermore, the expression pattern of genes involved in bone remodeling (runx2a, runx2b, and mmp9) was decreased in Neu1-KO fish. These phenotypes were quite similar to those of Neu1-KO mice and human sialidosis patients, indicating the effectiveness of the established Neu1-KO zebrafish for the study of vertebrate Neu1 sialidase.

Study on embryonic and larval developmental stages of Sucker head Garra gotyla (Gray 1830; Teleostei; Cyprinidae)

Garra gotyla is an indigenous coldwater fish of the cyprinid family and has wide geographical distribution in India as well as in other countries of Asia and Africa. Induced breeding in G. gotyla was carried out successfully for the first time and an attempt has been made to document developmental stages chronologically from the first minute of fertilization, through all stages of embryonic development until the fifth day post hatching. This experiment was carried out at 22-24°C water temperature at the Directorate of Coldwater Fisheries Research, Bhimtal, India. During the breeding trial, the fertilization rate was observed as 70-75% and hatching rate was 85-90%. The mature fertilized ova were measured as 0.8-1.0 mm in diameter and the perivitelline membrane became thick soon after fertilization and formation of the germ pole. The periods taken for complete developmental stages were recorded; cleavage stage 111 min (min post fertilization (pf)), blastulation stage 580 min (pf), neurulation and segmentation 1250 min (pf) and hatching was completed after 1420 min. The sac fry was measured as 3 mm in length and took almost 3 days for complete absorption of the yolk content. The major structural and differential changes observed are in head, tail, fins, alimentary canal, rudiments of each organ and appearance of melanophore pigmentation in the whole body. The 5-day-old larvae were measured as 6 mm in length with almost every organ fully differentiated. The present study will be utilized for large-scale production of fingerlings for stock enhancement in rivers, lakes and possibilities of genetic improvement and manipulation at the embryonic stage.

Comparative Toxicity of Nitrate to Common and Imperiled Freshwater Mussel Glochidia and Larval Fishes

Growing human populations and increasingly intensive agriculture have resulted in widespread aquatic nitrate pollution. Freshwater mussel populations have been in decline for decades but often are underrepresented in data used for the development of ambient water quality criteria and acute toxicity of nitrate to mussel glochidia has not yet been established. Additionally, toxicity testing with aquatic species often is limited to a few model species; however, relatively little is known about how representative model species are of imperiled species. Therefore, to better define the acute toxicity of nitrate to common and rare aquatic species, we conducted 24-h nitrate acute toxicity tests with glochidia of four species of freshwater mussels, including a federally threatened species (Hamiota altilis) and 7-day tests with larval fish of three species: fathead minnow (Pimephales promelas), tricolor shiner (Cyprinella trichroistia), and tilapia (Oreochromis spp.), across a range of water hardness. Median effective concentrations (EC₅₀s) in freshwater mussel glochidia ranged from 524 to 904 mg/L NO₃-N in moderately hard water. In fish, median lethal concentrations (LC₅₀s) ranged from 228 to 1725 mg/L NO₃-N and varied with water hardness. Of the species tested, generally sensitivity of the common species was similar to the rare species, although relative sensitivity varied with water hardness. Based on these results, we can conclude that acute lethal effects are unlikely for the fish and mussel species considered here at current environmental levels, but the results of these standardized tests are useful for the development of ambient water quality criteria and other regulatory and management decisions regarding acute nitrate exposures.

Glucose Injection Into Yolk Positively Modulates Intermediary Metabolism and Growth Performance in Juvenile Nile Tilapia (Oreochromis niloticus)

The aim of this study was to explore for the first time in omnivorous fish the concept of nutritional programming. A nutritional stimulus was accomplished by microinjecting 2 M glucose into yolk reserves during the alevin stage in Nile tilapia (Oreochromis niloticus). At the molecular level in fry, at 1 week post-injection, glucose stimuli were associated with the up-regulation of genes involved in glycolysis (pklr, hk1, hk2, and pkma), glucose transport (glut4) pathways and down-regulation of genes related to gluconeogenesis (g6pca1, g6pca2, and pck1) and amino acid catabolism (asat, alat) (P < 0.05), demonstrating that the larvae well received the glucose stimulus at a molecular level. Moreover, 20 weeks after glucose injection, early glucose stimuli were always linked to permanent effects in juvenile fish, as reflected by a higher level of glycolytic enzymes [gck, hk1 and hk2 at both mRNA and enzymatic levels and pyruvate kinase (PK) activity]. Finally, the effects of the glucose stimulus history were also examined in fish fed with two different dietary carbohydrate/protein levels (medium-carbohydrate diet, CHO-M; high-carbohydrate diet, CHO-H) in juvenile fish (during weeks 20-24). As expected, the CHO-H diet induced the expression of glycolytic and lipogenic genes (gck, pklr, hk1, hk2, fpkma, fasn, and g6pd) and suppressed the expression of gluconeogenic and amino acid catabolism genes (g6pca1, pck1, pck2, asat, alat, and gdh). Nevertheless, the early glucose stimulus led to persistent up-regulation of glycolytic enzymes (gck, pklr, hk1, and hk2) at both the mRNA and enzyme activity levels and glucose transporter glut4 as well as lower gluconeogenic pck1 gene expression (P < 0.05). More interestingly, the early glucose stimulus was associated with a better growth performance of juvenile fish irrespective of the diets. These permanent changes were associated with DNA hypomethylation in the liver and muscles, suggesting the existence of epigenetic mechanisms at the origin of programming. In conclusion, for the first time in tilapia, early glucose stimuli were found to be clearly associated with a positive metabolic programming effect later in life, improving the growth performance of the fish.

Appetite regulating genes may contribute to herbivory versus carnivory trophic divergence in haplochromine cichlids

Feeding is a complex behaviour comprised of satiety control, foraging, ingestion and subsequent digestion. Cichlids from the East African Great Lakes are renowned for their diverse trophic specializations, largely predicated on highly variable jaw morphologies. Thus, most research has focused on dissecting the genetic, morphological and regulatory basis of jaw and teeth development in these species. Here for the first time we explore another aspect of feeding, the regulation of appetite related genes that are expressed in the brain and control satiety in cichlid fishes. Using qPCR analysis, we first validate stably expressed reference genes in the brain of six haplochromine cichlid species at the end of larval development prior to foraging. We next evaluate the expression of 16 appetite related genes in herbivorous and carnivorous species from the parallel radiations of Lake Tanganyika, Malawi and Victoria. Interestingly, we find increased expression of two appetite-regulating genes (anorexigenic genes), cart and npy2r, in the brain of carnivorous species in all the three lakes. This supports the notion that appetite gene regulation might play a part in determining trophic niche specialization in divergent cichlid species, already prior to exposure to different diets. Our study contributes to the limited body of knowledge on the neurological circuitry that controls feeding transitions and adaptations in cichlids and other teleosts.

The impacts of seasonal variation on the immune status of Nile tilapia larvae and their response to different immunostimulants feed additives

Few data are available on the thermal tolerance of Nile tilapia fish larvae in relation to their immune status and survival. The aims of this work were to evaluate the immune status of one day old Nile tilapia (Oreochromis niloticus) larval stage collected at the beginning (March), middle (August) and at the end (October) of hatching season through morphometric assessment of the larvae parameters including yolk sac diameter, body length and width as well as the expression of some immune-related genes (rag, sacs and tlr), inflammatory (il1b and il8) and stress related genes (hsp27, hsp70). Also, to compare the effect of three different immunostimulants (β-glucan, Vitamin C, and methionine/lysine amino acids mix) on the expression of the studied genes at two variant temperatures (23 ± 1 °C and 30 ± 1 °C) in experimental study for 21 days. The immune status of Nile tilapia is affected by thermal fluctuation throughout the hatching season reflected by altered yolk sac size, length, and expression of the immune and stress related genes of the larvae, the best performances was observed at the beginning of the hatching season (March). High temperature (30 °C) suppress immune and stress responses throughout downregulation of all the genes under study, mask any effects for the immunostimulants, increased mortality in fish larvae suggesting narrow thermal tolerance range for the larvae compared with the adult fish. We recommend the use of amino acid mix as immunostimulant for Nile tilapia larvae, it reduces the mortality percentage and improve cellular response. Also, the use of β-glucan should be prohibited during this developmental stage of larvae, it induced the highest mortality percentage.

Developmental plasticity of epithelial stem cells in tooth and taste bud renewal

Nearly one-third of adults over the age of 65 have lost all their teeth. We set out to understand tooth renewal in animals that have replacement and regeneration capabilities. Using cichlid fishes and mouse models, we discovered plasticity between tooth and taste bud progenitor cell derivatives, mediated by BMP. Our results suggest that oral organs have surprising regenerative capabilities and can be manipulated to express characteristics of different tissue types.

In Lake Malawi cichlids, each tooth is replaced in one-for-one fashion every ∼20 to 50 d, and taste buds (TBs) are continuously renewed as in mammals. These structures are colocalized in the fish mouth and throat, from the point of initiation through adulthood. Here, we found that replacement teeth (RT) share a continuous band of epithelium with adjacent TBs and that both organs coexpress stem cell factors in subsets of label-retaining cells. We used RNA-seq to characterize transcriptomes of RT germs and TB-bearing oral epithelium. Analysis revealed differential usage of developmental pathways in RT compared to TB oral epithelia, as well as a repertoire of genome paralogues expressed complimentarily in each organ. Notably, BMP ligands were expressed in RT but excluded from TBs. Morphant fishes bathed in a BMP chemical antagonist exhibited RT with abrogated shh expression in the inner dental epithelium (IDE) and ectopic expression of calb2 (a TB marker) in these very cells. In the mouse, teeth are located on the jaw margin while TBs and other oral papillae are located on the tongue. Previous study reported that tongue intermolar eminence (IE) oral papillae of Follistatin (a BMP antagonist) mouse mutants exhibited dysmorphic invagination. We used these mutants to demonstrate altered transcriptomes and ectopic expression of dental markers in tongue IE. Our results suggest that vertebrate oral epithelium retains inherent plasticity to form tooth and taste-like cell types, mediated by BMP specification of progenitor cells. These findings indicate underappreciated epithelial cell populations with promising potential in bioengineering and dental therapeutics.

Bar, stripe and spot development in sand-dwelling cichlids from Lake Malawi

Background

Melanic patterns such as horizontal stripes, vertical bars and spots are common among teleost fishes and often serve roles in camouflage or mimicry. Extensive research in the zebrafish model has shown that the development of horizontal stripes depends on complex cellular interactions between melanophores, xanthophores and iridophores. Little is known about the development of horizontal stripes in other teleosts, and even less is known about bar or spot development. Here, we compare chromatophore composition and development of stripes, bars and spots in two cichlid species of sand-dwellers from Lake Malawi—Copadichromis azureus and Dimidiochromis compressiceps.

Results

(1) In D. compressiceps, stripes are made of dense melanophores underlaid by xanthophores and overlaid by iridophores. Melanophores and xanthophores are either loose or absent in interstripes, and iridophores are dense. In C. azureus, spots and bars are composed of a chromatophore arrangement similar to that of stripes but are separated by interbars where density of melanophores and xanthophores is only slightly lower than in stripes and iridophore density appears slightly greater. (2) Stripe, bar and spot chromatophores appear in the skin at metamorphosis. Stripe melanophores directly differentiate along horizontal myosepta into the adult pattern. In contrast, bar number and position are dynamic throughout development. As body length increases, new bars appear between old ones or by splitting of old ones through new melanophore appearance, not migration. Xanthophore and iridophore distributions follow melanophore patterns. (3) Metamorphic pigmentation arises in cichlids in a fashion similar to that described in zebrafish: melanophore progenitors derived from the medial route of neural crest migration migrate from the vicinity of the neural tube to the skin during metamorphosis.

Conclusion

The three pigment cell types forming stripes, bars and spots arise in the skin at metamorphosis. Stripes develop by differentiation of melanophores along horizontal myosepta, while bars do not develop along patent anatomical boundaries and increase in number in relation with body size. We propose that metamorphic melanophore differentiation and migratory arrest upon arrival to the skin lead to stripe formation, while bar formation must be supported by extensive migration of undifferentiated melanophores in the skin.

Hatching success and survival of fish early life stages in a chronic exposure to nevirapine: a case study of the Mozambique tilapia

The anti-retroviral nevirapine has been detected in surface waters throughout South Africa and its effects on non-target aquatic animals are still unknown. The aim was to investigate the potential effects of nevirapine on the hatching success and survival of Oreochromis mossambicus early life stages through a chronic exposure. The exposer started with newly fertilized O. mossambicus eggs and concluded 30 days after hatching. Environmental relevant concentration of nevirapine (1.48 µg/l) was used in a static renewal system and a controlled environment (27 ± 1°C; 14:10 day/night cycle). The main endpoints assessed included hatching success and survival; a morphological assessment was also done on whole individual on day 1 and 30 post-hatching to identify any physical abnormality. Nevirapine had no noticeable effects on the hatching success and survival of O. mossambicus larvae; no statistically significant differences were observed between the control and the nevirapine exposed fish (p > 0.05).

Divergence in larval jaw gene expression reflects differential trophic adaptation in haplochromine cichlids prior to foraging

Background

Understanding how variation in gene expression contributes to morphological diversity is a major goal in evolutionary biology. Cichlid fishes from the East African Great lakes exhibit striking diversity in trophic adaptations predicated on the functional modularity of their two sets of jaws (oral and pharyngeal). However, the transcriptional basis of this modularity is not so well understood, as no studies thus far have directly compared the expression of genes in the oral and pharyngeal jaws. Nor is it well understood how gene expression may have contributed to the parallel evolution of trophic morphologies across the replicate cichlid adaptive radiations in Lake Tanganyika, Malawi and Victoria.

Results

We set out to investigate the role of gene expression divergence in cichlid fishes from these three lakes adapted to herbivorous and carnivorous trophic niches. We focused on the development stage prior to the onset of exogenous feeding that is critical for understanding patterns of gene expression after oral and pharyngeal jaw skeletogenesis, anticipating environmental cues. This framework permitted us for the first time to test for signatures of gene expression underlying jaw modularity in convergent eco-morphologies across three independent adaptive radiations. We validated a set of reference genes, with stable expression between the two jaw types and across species, which can be important for future studies of gene expression in cichlid jaws. Next we found evidence of modular and non-modular gene expression between the two jaws, across different trophic niches and lakes. For instance, prdm1a, a skeletogenic gene with modular anterior-posterior expression, displayed higher pharyngeal jaw expression and modular expression pattern only in carnivorous species. Furthermore, we found the expression of genes in cichlids jaws from the youngest Lake Victoria to exhibit low modularity compared to the older lakes.

Conclusion

Overall, our results provide cross-species transcriptional comparisons of modularly-regulated skeletogenic genes in the two jaw types, implicating expression differences which might contribute to the formation of divergent trophic morphologies at the stage of larval independence prior to foraging.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1483-3) contains supplementary material, which is available to authorized users.

Peptide Transporters in the Primary Gastrointestinal Tract of Pre-Feeding Mozambique Tilapia Larva

Fish larvae differ greatly from the adult form in their morphology and organ functionality. The functionality of the gastrointestinal tract depends on the expression of various pumps, transporters, and channels responsible for feed digestion and nutrients absorption. During the larval period, the gastrointestinal tract develops from a simple closed tube, into its complex form with differentiated segments, crypts and villi, as found in the adult. In this study, we characterized the expression of three peptide transporters (PepT1a, PepT1b, and PepT2) in the gastrointestinal tract of Mozambique tilapia (Oreochromis mossambicus) larvae along 12 days of development, from pre-hatching to the completion of yolk sac absorption. Gene expression analysis revealed differential and complimentary time-dependent expression of the PepT1 variants and PepT2 along the larval development period. Immunofluorescence analysis showed differential protein localization of the three peptide transporters (PepTs) along the gastrointestinal tract, in a similar pattern to the adult. In addition, PepT1a was localized in mucosal cells in the larvae esophagus, in much higher abundance than in the adults. The results of this study demonstrate specialization of intestinal sections and absorbance potential of the enterocytes prior to the onset of active exogenous feeding, thus pointing to an uncharacterized function and role of the gastrointestinal tract and its transporters during the larval period.

Maternal ancestry analyses of red tilapia strains based on D-loop sequences of seven tilapia populations

Background

Many tilapia species or varieties have been widely introduced and have become an economically important food fish in China. Information on the genetic backgrounds of these populations is deficient and requires more research, especially for red tilapia strains.

Methods

In the present study, displacement loop (D-loop) sequences were used to evaluate the genetic relationship and diversity of seven tilapia populations that are widely cultured in China; this was done specifically to speculate on the maternal ancestry of red tilapia strains. Three red tilapia varieties of Oreochromis ssp., Taiwan (TW), Israel (IL), and Malaysia (MY) strains and other populations, including O. aureus (AR), O. niloticus (NL), O. mossambicus (MS), and the GIFT strain of O. niloticus, were collected and analyzed in this study.

Results

A total of 146 polymorphic sites and 32 haplotypes of D-loop sequences were detected among 332 fish and four major haplotypes were shared among the populations. The TW and NL populations had a greater number of haplotypes (20 and 8, respectively). The haplotype diversity (Hd) and nucleotide diversity (π) of each population ranged from 0.234 to 0.826, and 0 to 0.060, respectively. The significant positive Tajima’s D value of neutral test were detected in the NL, IL, and MY populations (P < 0.05), which indicated these populations might have not experienced historical expansion. According to the pairwise F-statistics, highly significant genetic differentiations were detected among populations (P < 0.01), with the exception of the IL and MY populations (P > 0.05). The nearest K2P genetic distance (D = 0.014) was detected between the MS and TW populations, whereas, the farthest (D = 0.101) was found between the GIFT and AR populations. The results from the molecular variance analysis (AMOVA) showed that there was an extremely significant genetic variation observed among the populations (P < 0.01), which contained 63.57% of the total variation. In view of the genetic relationship of red tilapia strains with other populations, TW and IL were detected with more similar genetic structures related to MS, and MY was more genetically similar to GIFT (or NL), which could provide more genetic evidence for the red tilapia strains maternal ancestry.

Maternal ancestry analyses of red tilapia strains based on D-loop sequences of seven tilapia populations

Background

Many tilapia species or varieties have been widely introduced and have become an economically important food fish in China. Information on the genetic backgrounds of these populations is deficient and requires more research, especially for red tilapia strains.

Methods

In the present study, displacement loop (D-loop) sequences were used to evaluate the genetic relationship and diversity of seven tilapia populations that are widely cultured in China; this was done specifically to speculate on the maternal ancestry of red tilapia strains. Three red tilapia varieties of Oreochromis ssp., Taiwan (TW), Israel (IL), and Malaysia (MY) strains and other populations, including O. aureus (AR), O. niloticus (NL), O. mossambicus (MS), and the GIFT strain of O. niloticus, were collected and analyzed in this study.

Results

A total of 146 polymorphic sites and 32 haplotypes of D-loop sequences were detected among 332 fish and four major haplotypes were shared among the populations. The TW and NL populations had a greater number of haplotypes (20 and 8, respectively). The haplotype diversity (Hd) and nucleotide diversity (π) of each population ranged from 0.234 to 0.826, and 0 to 0.060, respectively. The significant positive Tajima's D value of neutral test were detected in the NL, IL, and MY populations (P < 0.05), which indicated these populations might have not experienced historical expansion. According to the pairwise F-statistics, highly significant genetic differentiations were detected among populations (P < 0.01), with the exception of the IL and MY populations (P > 0.05). The nearest K2P genetic distance (D = 0.014) was detected between the MS and TW populations, whereas, the farthest (D = 0.101) was found between the GIFT and AR populations. The results from the molecular variance analysis (AMOVA) showed that there was an extremely significant genetic variation observed among the populations (P < 0.01), which contained 63.57% of the total variation. In view of the genetic relationship of red tilapia strains with other populations, TW and IL were detected with more similar genetic structures related to MS, and MY was more genetically similar to GIFT (or NL), which could provide more genetic evidence for the red tilapia strains maternal ancestry.

Profiling expression changes of genes associated with temperature and sex during high temperature-induced masculinization in the Nile tilapia brain

Fish sex-determining mechanisms can be classified as genotypic (GSD), temperature (TSD), or genotypic plus temperature effects (GSD+TE). Previous studies have shown that culturing water temperature during thermosensitive periods (TSP) could affect the expression of many genes in the gonad in some fish. However, few studies have focused on gene expression changes in the brain after temperature treatment during TSP in fish species. In this study, three families were developed by crossing XX neomales with XX females and one of them was used for transcriptome analysis. The results showed that a total of 105, 3164 and 4666 DEGs were respectively obtained in FC (female control) vs. FT (high temperature-treated females at TSP), FC vs. MC (male control), and MC vs. FT comparison groups. By profiling analysis, we show that the mRNA expression levels of 16 differentially expressed genes (DEGs) exhibited significant downregulation or upregulation after high temperature treatment and reached a similar level as that in MC. Among the 16 DEGs, LOC100699848 (lysine specific demethylase 6A) and Jarid2 contained JmjC domain, showing the possible important role of JmjC domain in response to temperature treatment in Nile tilapia. Kdm6b (lysine demethylase 6B) and Jarid2 have been shown to play important roles in reptile TSD, showing the relative conservation of underlying regulation mechanisms between TSD in reptile and TSD or GSD+TE in fish species. Finally, the transcriptome profiling was validated by quantitative real-time PCR in nine selected genes. These results provide a direction for investigating the GSD+TE molecular mechanism in fish species.

Frequent nonrandom shifts in the temporal sequence of developmental landmark events during teleost evolutionary diversification

Morphological transformations can be generated by evolutionary changes in the sequence of developmental events. In this study, we examined the evolutionary dynamics of the developmental sequence on a macroevolutionary scale in teleosts. Using the information from previous reports describing the development of 31 species, we extracted the developmental sequences of 19 landmark events involving the formation of phylogenetically conserved body parts; we then inferred ancestral developmental sequences by two different parsimony-based methods-event-pairing and continuous analysis. The phylogenetic comparisons of these sequences revealed event-dependent heterogeneity in the frequency of sequence changes. Most of the sequence changes occurred as exchanges of temporally neighboring events. These heterochronic changes in developmental sequences accumulated along evolutionary time, but the precise distribution of the changes over the teleostean phylogeny remains unclear due to technical limitations.

Sex biased expression of anti-Mullerian hormone (amh) gene in a live bearing fish, Gambusia holbrooki: Evolutionary implications and potential role in sex differentiation

The amh, a member of transforming growth factor-β (TGF-β) family, is known to play a critical role in vertebrate male sex differentiation, with its paralogue/s evolving to determine sex in few heterogametic (XX/XY) teleosts. However, it remains relatively unexplored in the reproductively unique live bearing teleosts. Therefore, this study comparatively examined the structure and content of G. holbrooki amh as well as characterised its expression. A paralogous Y-specific amh (amhy) was not detected, suggesting an unlikely role in sex determination. Two transcripts (1.4 and 1.5 kb) were detected in adults: the larger (1.5 kb) retaining intron 5, coding for a truncated AMH-N and no TGF-β domain. The small (1.4 kb) transcript, had both domains intact and clustered with members of Poeciliidae. In contrast to other vertebrates, a higher conservation between the N- rather than the C- terminus of amh in Poeciliidae was observed, suggesting an adaptation that may be unique to live bearing teleosts. The amh expression was 6 times higher in brain of both sexes and testis compared with ovaries (p = .001). Intriguingly, female splenic tissues showed 10 times higher expression (p = .006) and such female bias splenic expression has not been reported in any teleosts. Ontogenic expression was 25 times higher in male embryos at gastrulation stage (p = .001), much earlier than those reported in egg-laying teleosts. Such heightened expression in male embryos suggests a repressive role associated with proliferation and migration of primordial germ cells (PGCs) that are known to occur earlier at blastulation in teleosts-potentially influencing gonadal fate.

Maternal mRNA input of growth and stress-response-related genes in cichlids in relation to egg size and trophic specialization

BACKGROUND

Egg size represents an important form of maternal effect determined by a complex interplay of long-term adaptation and short-term plasticity balancing egg size with brood size. Haplochromine cichlids are maternal mouthbrooders showing differential parental investment in different species, manifested in great variation in egg size, brood size and duration of maternal care. Little is known about maternally determined molecular characters of eggs in fishes and their relation to egg size and trophic specialization. Here we investigate maternal mRNA inputs of selected growth- and stress-related genes in eggs of mouthbrooding cichlid fishes adapted to different trophic niches from Lake Tanganyika, Lake Malawi, Lake Victoria and compare them to their riverine allies.

RESULTS

We first identified two reference genes, atf7ip and mid1ip1, to be suitable for cross-species quantification of mRNA abundance via qRT-PCR in the cichlid eggs. Using these reference genes, we found substantial variation in maternal mRNA input for a set of candidate genes related to growth and stress response across species and lakes. We observed negative correlation of mRNA abundance between two of growth hormone receptor paralogs (ghr1 and ghr2) across all haplochromine cichlid species which also differentiate the species in the two younger lakes, Malawi and Lake Victoria, from those in Lake Tanganyika and ancestral riverine species. Furthermore, we found correlations between egg size and maternal mRNA abundance of two growth-related genes igf2 and ghr2 across the haplochromine cichlids as well as distinct clustering of the species based on their trophic specialization using maternal mRNA abundance of five genes (ghr1, ghr2, igf2, gr and sgk1).

CONCLUSIONS

These findings indicate that variations in egg size in closely related cichlid species can be linked to differences in maternal RNA deposition of key growth-related genes. In addition, the cichlid species with contrasting trophic specialization deposit different levels of maternal mRNAs in their eggs for particular growth-related genes; however, it is unclear whether such differences contribute to differential morphogenesis at later stages of development. Our results provide first insights into this aspect of gene activation, as a basis for future studies targeting their role during ecomorphological specialization and adaptive radiation.

An embryonic staging series up to hatching for Cyprinodon variegatus: An emerging fish model for developmental, evolutionary, and ecological research

Using multiple taxa to research development is necessary for making general conclusions about developmental patterns and mechanisms. We present a staging series for Cyprinodon variegatus as a basis for further study of the developmental biology of fishes in the genus Cyprinodon and for comparative work on teleost fishes beyond the standard models. Cyprinodon are small, euryhaline fishes, widely distributed in fresh, brackish, and hypersaline waters of southern and eastern North America. Cyprinodontids are closely related to fundulids, providing a comparative reference point to the embryological model, Fundulus heteroclitus. Ecologists and evolutionary biologists commonly study Cyprinodon, and we have been using Cyprinodon to study skull variation and its genetic basis among closely related species. We divided embryonic development of C. variegatus into 34 morphologically identifiable stages. We reference our staging series to that already defined for a related model species, Oryzias latipes (medaka) that is studied by a large community of researchers. We provide a description of the early chondrogenesis and ossification of skull and caudal fin bones during the latter stages of embryonic development. We show that Cyprinodon are tractable for studying development. Eggs can be obtained easily from breeding pairs and our study provides a staging system to facilitate future developmental studies.

Development and Validation of 58K SNP-Array and High-Density Linkage Map in Nile Tilapia (O. niloticus)

Despite being the second most important aquaculture species in the world accounting for 7.4% of global production in 2015, tilapia aquaculture has lacked genomic tools like SNP-arrays and high-density linkage maps to improve selection accuracy and accelerate genetic progress. In this paper, we describe the development of a genotyping array containing more than 58,000 SNPs for Nile tilapia (Oreochromis niloticus). SNPs were identified from whole genome resequencing of 32 individuals from the commercial population of the Genomar strain, and were selected for the SNP-array based on polymorphic information content and physical distribution across the genome using the Orenil1.1 genome assembly as reference sequence. SNP-performance was evaluated by genotyping 4991 individuals, including 689 offspring belonging to 41 full-sib families, which revealed high-quality genotype data for 43,588 SNPs. A preliminary genetic linkage map was constructed using Lepmap2 which in turn was integrated with information from the O_niloticus_UMD1 genome assembly to produce an integrated physical and genetic linkage map comprising 40,186 SNPs distributed across 22 linkage groups (LGs). Around one-third of the LGs showed a different recombination rate between sexes, with the female being greater than the male map by a factor of 1.2 (1632.9 to 1359.6 cM, respectively), with most LGs displaying a sigmoid recombination profile. Finally, the sex-determining locus was mapped to position 40.53 cM on LG23, in the vicinity of the anti-Müllerian hormone (amh) gene. These new resources has the potential to greatly influence and improve the genetic gain when applying genomic selection and surpass the difficulties of efficient selection for invasively measured traits in Nile tilapia.

Experimental evidence of the role of heterochrony in evolution of the Mesoamerican cichlids pigment patterns

The Mesoamerican cichlids display a spectacular diversity of pigment patterns, which serve a variety of functions and serve as a strong selective trait for this lineage. The development and variation of coloration in the Mesoamerican cichlids have been detailed by several groups. In particular, Říčan, Musilová, Muška, and Novák () and Říčan, Piálek, Dragová, and Novák () determined homology of pattern and revealed four alternative types of coloration and their ontogeny. In this work, this group posed an "ontogenetic timing hypothesis" proposing heterochronic shifts underlying major transitions in the evolution of the Mesoamerican cichlids. Here, we experimentally test this hypothesis by experimentally altering timing of pigment pattern formation in the convict cichlid Amatitlania nigrofasciata, a member of the Mesoamerican cichlids, via manipulations of thyroid hormone (TH) function. The response of different pigment cell lineages to TH-perturbations revealed that the transition from larval to juvenile coloration in the convict cichlid is under the control of TH-signaling. Importantly, hormonally induced changes in the timing of pigment cell lineages' development resulted in shifts of coloration ontogeny type observed between lineages and led to the appearance of phenotypes mimicking those in phylogenetically close and distant species. Thus, our findings support the hypothesis that simple changes in ontogenetic timing underlies species specific patterns in pigmentation and provide new perspectives for studying the role of endocrine signaling in the evolution of cichlids.