Anterior-posterior gene expression differences in three Lake Malawi cichlid fishes with variation in body stripe orientation

Comparative Transcriptome Analysis Reveals Differential Cutaneous Gene Expression in the Color Variation of Two Ornamental Discus, Red Melon and Red Cover

Red Melon (RM) and Red Cover (RC) discus (Symphysodon spp.) are ornamental fish varieties that were selectively bred from the wild parental lineages of the brown discus S. aquafaciatus over many generations, resulting in distinct cutaneous patterns from juveniles to adults. To better understand the underlying mechanisms, skin samples were collected from juveniles aged 60 days and adults aged 1 year from RM and RC for investigations. Microscopic observation detected xanthophores and erythrophores in all samples, except RC juveniles with no erythrophores. Melanophores were presented only in RC. The comparative analysis revealed that genes involved in pteridine synthesis (gch1 and zgc:153031), one-carbon metabolism (aldh1l2 and zgc153031), and lipid metabolism (apoda and klf1) were differentially expressed in RM juveniles, which may be associated with the development of erythrophores and xanthophores. The temporal inhibition of melanophore differentiation and development was observed in RM juveniles, coupled with elevated expression of notum2 and sost, two antagonist genes in Wnt-signaling, suggesting their roles in melanophore development. Distinct pigment pattern between RM and RC since the juvenile stage may be driven by the differential expression of multiple axial developmental genes, including GATA, ankyrin, and mitotic spindle orientation proteins. This is the first report to describe the differential growth of cutaneous pigments and the molecular processes involved in red discus. The results provided valuable insights into pigment pattern differences in an interesting ornamental fish model.

African cichlid fishes: morphological data and taxonomic insights from a genus-level survey of supraneurals, pterygiophores, and vertebral counts (Ovalentaria, Blenniiformes, Cichlidae, Pseudocrenilabrinae)

Background

The iconic freshwater cichlid fishes (Cichlidae) comprise about 1750 validly named species and hundreds more that are known, but not yet described and named. Cichlids are an important source of protein for millions of people on several continents, are model organisms in studies of evolution, speciation, ecology, development, behaviour and physiology and are popular as aquarium fishes. Yet, comparative studies of cichlid internal anatomy are rare. Even their osteology has not been taxonomically surveyed. The cichlid postcranial skeleton has been especially neglected.

New information

Here, I provide the first survey in cichlids of the considerable variation in numbers of vertebrae, supraneurals and dorsal- and anal-fin supports (pterygiophores), as well as the patterns with which the pterygiophores insert between the neural or haemal spines. The study includes some 1700 specimens of nearly 400 cichlid species. Focusing on the largest subfamily, the African cichlids or Pseudocrenilabrinae, the survey furnishes data from species in all but one of its 166 genera. Limited data from species in the other cichlid subfamilies (Etroplinae, Ptychochrominae and Cichlinae) and from the related leaffishes, Polycentridae, are also presented. Key examples of pterygiophore insertion patterns from throughout the range of variation are illustrated and discussed. Detailed analytical tables and all raw data are provided in supplementary files.A bizarre specialisation in Cyprichromis is noted, evidently for the first time. Uniquely in this Lake Tanganyikan genus, five to seven anal pterygiophores are abdominal in position, located anterior to the anal fin and inserting toward or between successive pairs of pleural ribs.Taxonomic changes: The most speciose tribe of African cichlids, currently known as Haplochromini, is correctly called Pseudocrenilabrini. Based chiefly on the molecular phylogenetic findings of other workers, I propose four pseudocrenilabrine subtribes, one occurring in rivers and three endemic to Lake Malawi. I also re-assign the Lake Tanganyikan tribe Tropheini as another subtribe of Pseudocrenilabrini, in line with numerous molecular studies placing tropheines firmly within this tribe. The remaining genera of Pseudocrenilabrini remain incertae sedis in this tribe pending clarification of their phylogenetic relationships.The character complex here surveyed is a promising source of taxonomically and phylogenetically informative characteristics distinguishing or uniting cichlid taxa at multiple hierarchical levels, from species through subfamily. This reference set of novel character data can also provide information for palaeontological studies of African cichlids. These attributes are skeletal features potentially available for study in well preserved fossils and may help determine their correct taxonomic placement.

What Makes a Mimic? Orange, Red, and Black Color Production in the Mimic Poison Frog (Ranitomeya imitator)

Aposematic organisms rely on their conspicuous appearance to signal that they are defended and unpalatable. Such phenotypes are strongly tied to survival and reproduction. Aposematic colors and patterns are highly variable; however, the genetic, biochemical, and physiological mechanisms producing this conspicuous coloration remain largely unidentified. Here, we identify genes potentially affecting color variation in two color morphs of Ranitomeya imitator: the orange-banded Sauce and the redheaded Varadero morphs. We examine gene expression in black and orange skin patches from the Sauce morph and black and red skin patches from the Varadero morph. We identified genes differentially expressed between skin patches, including those that are involved in melanin synthesis (e.g. mlana, pmel, tyrp1), iridophore development (e.g. paics, ppat, ak1), pteridine synthesis (e.g. gch1, pax3-a, xdh), and carotenoid metabolism (e.g. dgat2, rbp1, scarb2). In addition, using weighted correlation network analysis, we identified the top 50 genes with high connectivity from the most significant network associated with gene expression differences between color morphs. Of these 50 genes, 13 were known to be related to color production (gch1, gmps, gpr143, impdh1, mc1r, pax3-a, pax7, ppat, rab27a, rlbp1, tfec, trpm1, xdh).

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.

Developmental and Cellular Basis of Vertical Bar Color Patterns in the East African Cichlid Fish Haplochromis latifasciatus

The East African adaptive radiations of cichlid fishes are renowned for their diversity in coloration. Yet, the developmental basis of pigment pattern formation remains largely unknown. One of the most common melanic patterns in cichlid fishes are vertical bar patterns. Here we describe the ontogeny of this conspicuous pattern in the Lake Kyoga species Haplochromis latifasciatus. Beginning with the larval stages we tracked the formation of this stereotypic color pattern and discovered that its macroscopic appearance is largely explained by an increase in melanophore density and accumulation of melanin during the first 3 weeks post-fertilization. The embryonal analysis is complemented with cytological quantifications of pigment cells in adult scales and the dermis beneath the scales. In adults, melanic bars are characterized by a two to threefold higher density of melanophores than in the intervening yellow interbars. We found no strong support for differences in other pigment cell types such as xanthophores. Quantitative PCRs for twelve known pigmentation genes showed that expression of melanin synthesis genes tyr and tyrp1a is increased five to sixfold in melanic bars, while xanthophore and iridophore marker genes are not differentially expressed. In summary, we provide novel insights on how vertical bars, one of the most widespread vertebrate color patterns, are formed through dynamic control of melanophore density, melanin synthesis and melanosome dispersal.

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.

Genome of the Malawi golden cichlid fish (Melanochromis auratus) reveals exon loss of oca2 in an amelanistic morph

The tropical freshwater fish family Cichlidae is famous for its record-breaking rates of speciation and diversity in colors and color patterns. Here, we sequenced the genome of the Lake Malawi cichlid Melanochromis auratus to study the genetic basis of an amelanistic morph of this species that lacks the typical melanic stripes and markings. Genome sequencing of the amelanistic and wild-type morph revealed the loss of the second exon of the known pigmentation gene oculocutaneous albinism II (oca2), also known as p(ink-eyed dilution) gene or melanocyte-specific transporter gene. Additional genotyping confirms the complete association with this recessive Mendelian phenotype. The deletion results in a shorter transcript, lacking an acidic di-leucine domain that is crucial for trafficking of the Oca2 protein to melanosomes. The fact that oca2 is involved in a wide range of amelanistic morphs across vertebrates demonstrates its highly conserved function.

Variation in pigmentation gene expression is associated with distinct aposematic color morphs in the poison frog Dendrobates auratus

Background

Color and pattern phenotypes have clear implications for survival and reproduction in many species. However, the mechanisms that produce this coloration are still poorly characterized, especially at the genomic level. Here we have taken a transcriptomics-based approach to elucidate the underlying genetic mechanisms affecting color and pattern in a highly polytypic poison frog. We sequenced RNA from the skin from four different color morphs during the final stage of metamorphosis and assembled a de novo transcriptome. We then investigated differential gene expression, with an emphasis on examining candidate color genes from other taxa.

Results

Overall, we found differential expression of a suite of genes that control melanogenesis, melanocyte differentiation, and melanocyte proliferation (e.g., tyrp1, lef1, leo1, and mitf) as well as several differentially expressed genes involved in purine synthesis and iridophore development (e.g., arfgap1, arfgap2, airc, and gart).

Conclusions

Our results provide evidence that several gene networks known to affect color and pattern in vertebrates play a role in color and pattern variation in this species of poison frog.

Electronic supplementary material

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

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.

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

Variation in fin shape and size contributes to the outstanding morphological diversity of teleost fishes, but the regulation of fin growth has not yet been studied extensively outside the zebrafish model. A previous gene expression study addressing the ornamental elongations of unpaired fins in the African cichlid fish Neolamprologus brichardi identified three genes (cx43, mmp9 and sema3d) with strong and consistent expression differences between short and elongated fin regions. Remarkably, the expression patterns of these genes were not consistent with inferences on their regulatory interactions in zebrafish. Here, we identify a gene expression network (GRN) comprising cx43, mmp9, and possibly also sema3d 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 11 co-expressed genes, six TFs (foxc1, foxp1, foxd3, myc, egr2, irf8) showed expression patterns consistent with their cooperative transcriptional regulation of the gene network. Some of these TFs have already been implicated in teleost fish fin regeneration and formation. We particularly discuss the potential function of foxd3 as driver of the network and its role in the unexpected gene expression correlations observed in N. brichardi.