Towards a gene regulatory network shaping the fins of the Princess cichlid

Transcriptome analyses of betta fish (Betta splendens) provide novel insights into fin regeneration and color-related genes

The betta fish (Betta splendens), an important ornamental fish, haswell-developed and colorful fins.After fin amputation, betta fish can easily regenerate finssimilar to the originalsin terms of structureand color. The powerful fin regeneration ability and a variety of colors in the betta fish are fascinating. However, the underlying molecular mechanisms are still not fully understood. In this study, tail fin amputation and regeneration experiments were performed on two kinds of betta fish: red and white color betta fish. Then, transcriptome analyseswere conducted to screen out fin regeneration and color-relatedgenes in betta fish. Through enrichment analyses of differentially expressed genes (DEGs), we founda series of enrichment pathways and genes related to finregeneration, including cell cycle (i.e. plcg2), TGF-beta signaling pathway (i.e. bmp6), PI3K-Akt signaling pathway (i.e. loxl2aand loxl2b), Wnt signaling pathway(i.e. lef1), gap junctions (i.e. cx43), angiogenesis (i.e. foxp1), and interferon regulatory factor (i.e. irf8). Meanwhile, some fin color-related pathways and genes were identified in betta fish, especially melanogenesis (i.e. tyr, tyrp1a, tyrp1b, and mc1r) and carotenoid color genes (i.e. pax3, pax7, sox10, and ednrba). In conclusion, this studycan not only enrich the research onfish tissue regeneration, but also has a potential significance for the aquaculture and breeding of the betta fish.

A pituitary gene network linking vgll3 to regulators of sexual maturation in male Atlantic salmon

Age at maturity is a key life history trait and a significant contributor to life history strategy variation. The maturation process is influenced by genetic and environmental factors, but specific causes of variation in maturation timing remain elusive. In many species, the increase in the regulatory gonadotropin-releasing hormone 1 (GnRH1) marks the onset of puberty. Atlantic salmon, however, lacks the gnrh1 gene, suggesting gnrh3 and/or other regulatory factors are involved in the maturation process. Earlier research in Atlantic salmon has found a strong association between alternative alleles of vgll3 and maturation timing. Recently we reported strong induction of gonadotropin genes (fshb and lhb) in the pituitary of Atlantic salmon homozygous for the early maturation allele (E) of vgll3. The induction of gonadotropins was accompanied by increased expression of their direct upstream regulators, c-jun and sf1 (nr5a1b) but the regulatory connection between vgll3 and these regulators has never been investigated in any organism. In this study, we investigated the potential regulatory connection between vgll3 genotypes and these regulators through a stepwise approach of identifying a gene regulatory network (GRN) containing c-jun and sf1, and transcription factor motif enrichment analysis. We found a GRN containing c-jun with predicted upstream regulators, e2f1, egr1, foxj1 and klf4, to be differentially expressed in the pituitary. Finally, we suggest a vgll3 and Hippo pathway -dependent model for transcriptional regulation of c-jun and sf1 in the pituitary, which may have broader implications across vertebrates.

Gene expression patterns associated with caudal fin shape in the cichlid Lamprologus tigripictilis

Variation in fin shape is one of the most prominent features of morphological diversity among fish. Regulation of fin growth has mainly been studied in zebrafish, and it is not clear whether the molecular mechanisms underlying shape variation are equally diverse or rather conserved across species. In the present study, expression levels of 37 candidate genes were tested for association with fin shape in the cichlid fish Lamprologus tigripictilis. The tested genes included members of a fin shape-associated gene regulatory network identified in a previous study and novel candidates selected within this study. Using both intact and regenerating fin tissue, we tested for expression differences between the elongated and the short regions of the spade-shaped caudal fin and identified 20 genes and transcription factors (including angptl5, cd63, csrp1a, cx43, esco2, gbf1, and rbpj), whose expression patterns were consistent with a role in fin growth. Collated with available gene expression data of two other cichlid species, our study not only highlights several genes that were correlated with fin growth in all three species (e.g., angptl5, cd63, cx43, and mmp9), but also reveals species-specific gene expression and correlation patterns, which indicate considerable divergence in the regulatory mechanisms of fin growth across cichlids.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10750-022-05068-4.

Expression variations in ectodysplasin-A gene (eda) may contribute to morphological divergence of scales in haplochromine cichlids

BACKGROUND

Elasmoid scales are one of the most common dermal appendages and can be found in almost all species of bony fish differing greatly in their shape. Whilst the genetic underpinnings behind elasmoid scale development have been investigated, not much is known about the mechanisms involved in moulding of scales. To investigate the links between gene expression differences and morphological divergence, we inferred shape variation of scales from two different areas of the body (anterior and posterior) stemming from ten haplochromine cichlid species from different origins (Lake Tanganyika, Lake Malawi, Lake Victoria and riverine). Additionally, we investigated transcriptional differences of a set of genes known to be involved in scale development and morphogenesis in fish.

RESULTS

We found that scales from the anterior and posterior part of the body strongly differ in their overall shape, and a separate look on scales from each body part revealed similar trajectories of shape differences considering the lake origin of single investigated species. Above all, nine as well as 11 out of 16 target genes showed expression differences between the lakes for the anterior and posterior dataset, respectively. Whereas in posterior scales four genes (dlx5, eda, rankl and shh) revealed significant correlations between expression and morphological differentiation, in anterior scales only one gene (eda) showed such a correlation. Furthermore, eda displayed the most significant expression difference between species of Lake Tanganyika and species of the other two younger lakes. Finally, we found genetic differences in downstream regions of eda gene (e.g., in the eda-tnfsf13b inter-genic region) that are associated with observed expression differences. This is reminiscent of a genetic difference in the eda-tnfsf13b inter-genic region which leads to gain or loss of armour plates in stickleback.

CONCLUSION

These findings provide evidence for cross-species transcriptional differences of an important morphogenetic factor, eda, which is involved in formation of ectodermal appendages. These expression differences appeared to be associated with morphological differences observed in the scales of haplochromine cichlids indicating potential role of eda mediated signal in divergent scale morphogenesis in fish.

Transcriptional study reveals a potential leptin-dependent gene regulatory network in zebrafish brain

The signal mediated by leptin hormone and its receptor is a major regulator of body weight, food intake and metabolism. In mammals and many teleost fish species, leptin has an anorexigenic role and inhibits food intake by influencing the appetite centres in the hypothalamus. However, the regulatory connections between leptin and downstream genes mediating its appetite-regulating effects are still not fully explored in teleost fish. In this study, we used a loss of function leptin receptor zebrafish mutant and real-time quantitative PCR to assess brain expression patterns of several previously identified anorexigenic genes downstream of leptin signal under different feeding conditions (normal feeding, 7-day fasting, 2 and 6-h refeeding). These downstream factors include members of cart genes, crhb and gnrh2, as well as selected genes co-expressed with them based on a zebrafish co-expression database. Here, we found a potential gene expression network (GRN) comprising the abovementioned genes by a stepwise approach of identifying co-expression modules and predicting their upstream regulators. Among the transcription factors (TFs) predicted as potential upstream regulators of this GRN, we found expression pattern of sp3a to be correlated with transcriptional changes of the downstream gene network. Interestingly, the expression and transcriptional activity of Sp3 orthologous gene in mammals have already been implicated to be under the influence of leptin signal. These findings suggest a potentially conserved regulatory connection between leptin and sp3a, which is predicted to act as a transcriptional driver of a downstream gene network in the zebrafish brain.

Gene co-expression network analysis reveals mechanisms underlying ozone-induced carbamazepine toxicity in zebrafish (Danio rerio) embryos

Sewage effluent ozonation can reduce concentrations of chemical pollutants including pharmaceutical residues. However, the formation of potentially toxic ozonation byproducts (OBPs) is a matter of concern. This study sought to elucidate toxicity mechanisms of ozonated carbamazepine (CBZ), an anti-epileptic drug frequently detected in sewage effluents and surface water, in zebrafish embryos (Danio rerio). Embryos were exposed to ozonated and non-ozonated CBZ from 3 h post-fertilization (hpf) until 144 hpf. Embryotoxicity endpoints (proportion of dead and malformed embryos) were assessed at 24, 48, and 144 hpf. Heart rate was recorded at 48 hpf. Exposure to ozonated CBZ gave rise to cardiovascular-related malformations and reduced heart rate. Moreover, embryo-larvae exposed to ozonated CBZ displayed a lack of swim bladder inflation. Hence, the expression patterns of CBZ target genes involved in cardiovascular and embryonal development were investigated through a stepwise gene co-expression analysis approach. Two co-expression networks and their upstream transcription regulators were identified, offering mechanistic explanations for the observed toxicity phenotypes. The study presents a novel application of gene co-expression analysis elucidating potential toxicity mechanisms of an ozonated pharmaceutical with environmental relevance. The resulting data was used to establish a putative adverse outcome pathway (AOP).

Transcriptomics unravels molecular players shaping dorsal lip hypertrophy in the vacuum cleaner cichlid, Gnathochromis permaxillaris

BACKGROUND

Teleosts display a spectacular diversity of craniofacial adaptations that often mediates ecological specializations. A considerable amount of research has revealed molecular players underlying skeletal craniofacial morphologies, but less is known about soft craniofacial phenotypes. Here we focus on an example of lip hypertrophy in the benthivorous Lake Tangnayika cichlid, Gnathochromis permaxillaris, considered to be a morphological adaptation to extract invertebrates out of the uppermost layer of mud bottom. We investigate the molecular and regulatory basis of lip hypertrophy in G. permaxillaris using a comparative transcriptomic approach.

RESULTS

We identified a gene regulatory network involved in tissue overgrowth and cellular hypertrophy, potentially associated with the formation of a locally restricted hypertrophic lip in a teleost fish species. Of particular interest were the increased expression level of apoda and fhl2, as well as reduced expression of cyp1a, gimap8, lama5 and rasal3, in the hypertrophic lip region which have been implicated in lip formation in other vertebrates. Among the predicted upstream transcription factors, we found reduced expression of foxp1 in the hypertrophic lip region, which is known to act as repressor of cell growth and proliferation, and its function has been associated with hypertrophy of upper lip in human.

CONCLUSION

Our results provide a genetic foundation for future studies of molecular players shaping soft and exaggerated, but locally restricted, craniofacial morphological changes in fish and perhaps across vertebrates. In the future, we advocate integrating gene regulatory networks of various craniofacial phenotypes to understand how they collectively govern trophic and behavioural adaptations.

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.

Molecular mechanisms underlying nuchal hump formation in dolphin cichlid, Cyrtocara moorii

East African cichlid fishes represent a model to tackle adaptive changes and their connection to rapid speciation and ecological distinction. In comparison to bony craniofacial tissues, adaptive morphogenesis of soft tissues has been rarely addressed, particularly at the molecular level. The nuchal hump in cichlids fishes is one such soft-tissue and exaggerated trait that is hypothesized to play an innovative role in the adaptive radiation of cichlids fishes. It has also evolved in parallel across lakes in East Africa and Central America. Using gene expression profiling, we identified and validated a set of genes involved in nuchal hump formation in the Lake Malawi dolphin cichlid, Cyrtocara moorii. In particular, we found genes differentially expressed in the nuchal hump, which are involved in controlling cell proliferation (btg3, fosl1a and pdgfrb), cell growth (dlk1), craniofacial morphogenesis (dlx5a, mycn and tcf12), as well as regulators of growth-related signals (dpt, pappa and socs2). This is the first study to identify the set of genes associated with nuchal hump formation in cichlids. Given that the hump is a trait that evolved repeatedly in several African and American cichlid lineages, it would be interesting to see if the molecular pathways and genes triggering hump formation follow a common genetic track or if the trait evolved in parallel, with distinct mechanisms, in other cichlid adaptive radiations and even in other teleost fishes.

Transcriptional study of appetite regulating genes in the brain of zebrafish (Danio rerio) with impaired leptin signalling

The hormone leptin is a key regulator of body weight, food intake and metabolism. In mammals, leptin acts as an anorexigen and inhibits food intake centrally by affecting the appetite centres in the hypothalamus. In teleost fish, the regulatory connections between leptin and other appetite-regulating genes are largely unknown. In the present study, we used a zebrafish mutant with a loss of function leptin receptor to investigate brain expression patterns of 12 orexigenic and 24 anorexigenic genes under different feeding conditions (normal feeding, 7-day fasting, 2 and 6-hours refeeding). Expression patterns were compared to wild-type zebrafish, in order to identify leptin-dependent differentially expressed genes under different feeding conditions. We provide evidence that the transcription of certain orexigenic and anorexigenic genes is influenced by leptin signalling in the zebrafish brain. We found that the expression of orexigenic genes was not affected by impaired leptin signalling under normal feeding conditions; however, several orexigenic genes showed increased transcription during fasting and refeeding, including agrp, apln, galr1a and cnr1. This suggests an inhibitory effect of leptin signal on the transcription of these orexigenic genes during short-term fasting and refeeding in functional zebrafish. Most pronounced effects were observed in the group of anorexigenic genes, where the impairment of leptin signalling resulted in reduced gene expression in several genes, including cart family, crhb, gnrh2, mc4r, pomc and spx, in the control group. This suggests a stimulatory effect of leptin signal on the transcription of these anorexigenic genes under normal feeding condition. In addition, we found multiple gain and loss in expression correlations between the appetite-regulating genes, in zebrafish with impaired leptin signal, suggesting the presence of gene regulatory networks downstream of leptin signal in zebrafish brain. The results provide the first evidence for the effects of leptin signal on the transcription of various appetite-regulating genes in zebrafish brain, under different feeding conditions. Altogether, these transcriptional changes suggest an anorexigenic role for leptin signal, which is likely to be mediated through distinct set of appetite-regulating genes under different feeding conditions.

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.

Gene expression profiling suggests differences in molecular mechanisms of fin elongation between cichlid species

Comparative analyses of gene regulation inform about the molecular basis of phenotypic trait evolution. Here, we address a fin shape phenotype that evolved multiple times independently across teleost fish, including several species within the family Cichlidae. In a previous study, we proposed a gene regulatory network (GRN) involved in the formation and regeneration of conspicuous filamentous elongations adorning the unpaired fins of the Neolamprologus brichardi. Here, we tested the members of this network in the blockhead cichlid, Steatocranus casuarius, which displays conspicuously elongated dorsal and moderately elongated anal fins. Our study provided evidence for differences in the anatomy of fin elongation and suggested gene regulatory divergence between the two cichlid species. Only a subset of the 20 genes tested in S. casuarius showed the qPCR expression patterns predicted from the GRN identified in N. brichardi, and several of the gene-by-gene expression correlations differed between the two cichlid species. In comparison to N. brichardi, gene expression patterns in S. casuarius were in better (but not full) agreement with gene regulatory interactions inferred in zebrafish. Within S. casuarius, the dorsoventral asymmetry in ornament expression was accompanied by differences in gene expression patterns, including potential regulatory differentiation, between the anal and dorsal fin.