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

Maternal developmental history alters transfer of circadian clock genes to offspring in Japanese quail (Coturnix japonica)

Maternal signals shape embryonic development, and in turn post-natal phenotypes. RNA deposition is one such method of maternal signalling and circadian rhythms are one trait thought to be maternally inherited, through this mechanism. These maternal circadian gene transcripts aid development of a functioning circadian system. There is increasing evidence that maternal signals can be modified, depending on prevailing environmental conditions to optimise offspring fitness. However, currently, it is unknown if maternal circadian gene transcripts, and consequently early embryonic gene transcription, are altered by maternal developmental conditions. Here, using avian mothers who experienced either pre-natal corticosterone exposure, and/or post-natal stress as juveniles we were able to determine the effects of the timing of stress on downstream circadian RNA deposition in offspring. We demonstrated that maternal developmental history does indeed affect transfer of offspring circadian genes, but the timing of stress was important. Avian mothers who experienced stress during the first 2 weeks of post-natal life increased maternally deposited transcript levels of two core circadian clock genes, BMAL1 and PER2. These differences in transcript levels were transient and disappeared at the point of embryonic genome transcription. Pre-natal maternal stress alone was found to elicit delayed changes in circadian gene expression. After activation of the embryonic genome, both BMAL1 and PER2 expression were significantly decreased. If both pre-natal and post-natal stress occurred, then initial maternal transcript levels of BMAL1 were significantly increased. Taken together, these results suggest that developmental stress differentially produces persistent transgenerational effects on offspring circadian genes.

Low incidence of cannibalism among brood parasitic cuckoo catfish embryos

Brood parasites have demanding needs of host resources. Brood parasitic offspring are highly competitive and frequently cause the failure of host broods and the survival of a single parasitic offspring. Accordingly, virulent brood parasites lay a single egg in the same host nest to avoid sibling competition. In the cuckoo catfish (Synodontis multipunctatus), which parasitize mouthbrooding cichlid fishes in Lake Tanganyika, the modes of host and parasite oviposition lead to frequent cases of multiple parasitism. We experimentally tested the prediction that multiple parasitism leads to frequent cannibalism among offspring. Cuckoo catfish embryos prey upon host offspring to obtain nourishment during their 3-week development in the host buccal cavity and may also consume conspecific embryos. The potential benefits of cannibalism in the system are, therefore, twofold; to decrease competition for limited resources (i.e., host brood with rich yolk sacs) and to directly obtain nourishment by consuming rivals. We found that cannibalism indeed provided measurable benefits in terms of increased growth of the cannibals, but cannibalism was rare and typically occurred once all host offspring had been consumed. This implies that cannibalism among cuckoo catfish embryos emerges to mitigate starvation rather than eliminate competition.

Maternal patterns of inheritance alter transcript expression in eggs

BACKGROUND

Modifications to early development can lead to evolutionary diversification. The early stages of development are under maternal control, as mothers produce eggs loaded with nutrients, proteins and mRNAs that direct early embryogenesis. Maternally provided mRNAs are the only expressed genes in initial stages of development and are tightly regulated. Differences in maternal mRNA provisioning could lead to phenotypic changes in embryogenesis and ultimately evolutionary changes in development. However, the extent that maternal mRNA expression in eggs can vary is unknown for most developmental models. Here, we use a species with dimorphic development- where females make eggs and larvae of different sizes and life-history modes-to investigate the extent of variation in maternal mRNA provisioning to the egg.

RESULTS

We find that there is significant variation in gene expression across eggs of different development modes, and that there are both qualitative and quantitative differences in mRNA expression. We separate parental effects from allelic effects, and find that both mechanisms contribute to mRNA expression differences. We also find that offspring of intraspecific crosses differentially provision their eggs based on the parental cross direction (a parental effect), which has not been previously demonstrated in reproductive traits like oogenesis.

CONCLUSION

We find that maternally controlled initiation of development is functionally distinct between eggs of different sizes and maternal genotypes. Both allele-specific effects and parent-of-origin effects contribute to gene expression differences in eggs. The latter indicates an intergenerational effect where a parent's genotype can affect gene expression in an egg made by the next generation.

Maternal body condition and season influence RNA deposition in the oocytes of alfalfa leafcutting bees (Megachile rotundata)

Maternal effects are an important source of phenotypic variance, whereby females influence offspring developmental trajectory beyond direct genetic contributions, often in response to changing environmental conditions. However, relatively little is known about the mechanisms by which maternal experience is translated into molecular signals that shape offspring development. One such signal may be maternal RNA transcripts (mRNAs and miRNAs) deposited into maturing oocytes. These regulate the earliest stages of development of all animals, but are understudied in most insects. Here we investigated the effects of female internal (body condition) and external (time of season) environmental conditions on maternal RNA in the maturing oocytes and 24-h-old eggs (24-h eggs) of alfalfa leafcutting bees. Using gene expression and WGCNA analysis, we found that females adjust the quantity of mRNAs related to protein phosphorylation, transcriptional regulation, and nuclease activity deposited into maturing oocytes in response to both poor body condition and shorter day lengths that accompany the late season. However, the magnitude of these changes was higher for time of season. Females also adjusted miRNA deposition in response to seasonal changes, but not body condition. We did not observe significant changes in maternal RNAs in response to either body condition or time of season in 24-h eggs, which were past the maternal-to-zygotic transition. Our results suggest that females adjust the RNA transcripts they provide for offspring to regulate development in response to both internal and external environmental cues. Variation in maternal RNAs may, therefore, be important for regulating offspring phenotype in response to environmental change.

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.

Developmental series of gene expression clarifies maternal mRNA provisioning and maternal-to-zygotic transition in a reef-building coral

BACKGROUND

Maternal mRNA provisioning of oocytes regulates early embryogenesis. Maternal transcripts are degraded as zygotic genome activation (ZGA) intensifies, a phenomenon known as the maternal-to-zygotic transition (MZT). Here, we examine gene expression over nine developmental stages in the Pacific rice coral, Montipora capitata, from eggs and embryos at 1, 4, 9, 14, 22, and 36 h-post-fertilization (hpf), as well as swimming larvae (9d), and adult colonies.

RESULTS

Weighted Gene Coexpression Network Analysis revealed four expression peaks, identifying the maternal complement, two waves of the MZT, and adult expression. Gene ontology enrichment revealed maternal mRNAs are dominated by cell division, methylation, biosynthesis, metabolism, and protein/RNA processing and transport functions. The first MZT wave occurs from ~4-14 hpf and is enriched in terms related to biosynthesis, methylation, cell division, and transcription. In contrast, functional enrichment in the second MZT wave, or ZGA, from 22 hpf-9dpf, includes ion/peptide transport and cell signaling. Finally, adult expression is enriched for functions related to signaling, metabolism, and ion/peptide transport. Our proposed MZT timing is further supported by expression of enzymes involved in zygotic transcriptional repression (Kaiso) and activation (Sox2), which peak at 14 hpf and 22 hpf, respectively. Further, DNA methylation writing (DNMT3a) and removing (TET1) enzymes peak and remain stable past ~4 hpf, suggesting that methylome programming occurs before 4 hpf.

CONCLUSIONS

Our high-resolution insight into the coral maternal mRNA and MZT provides essential baseline information to understand parental carryover effects and the sensitivity of developmental success under increasing environmental stress.

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.

Tails of reproduction: Regeneration leads to increased reproductive investment

Trade-offs between life-history traits are due to limited resources or constraints in the regulation of genetic and physiological networks. Tail autotomy, with subsequent regeneration, is a common anti-predation mechanism in lizards and is predicted to trade-off with life-history traits, such as reproduction. We utilize the brown anole lizard with its unusual reproductive pattern of single-egg clutches every 7-10 days to test for a trade-off in reproductive investment over 8 weeks of tail regeneration on a limited diet. In contrast to predictions, we found that investing in tissue regeneration had a positive effect on reproduction in terms of egg size (11.7% relative to controls) and hatchling size (11.5% relative to controls), and no effect on egg number or survival, with the increase in reproduction starting at peak regeneration. We discuss mechanistic hypotheses that the process of regeneration may cause increased energetic efficiency or utilized shared physiological pathways with reproductive investment.

Differences among families in craniofacial shape at early life-stages of Arctic charr (Salvelinus alpinus)

BACKGROUND

Organismal fitness can be determined at early life-stages, but phenotypic variation at early life-stages is rarely considered in studies on evolutionary diversification. The trophic apparatus has been shown to contribute to sympatric resource-mediated divergence in several taxa. However, processes underlying diversification in trophic traits are poorly understood. Using phenotypically variable Icelandic Arctic charr (Salvelinus alpinus), we reared offspring from multiple families under standardized laboratory conditions and tested to what extent family (i.e. direct genetic and maternal effects) contributes to offspring morphology at hatching (H) and first feeding (FF). To understand the underlying mechanisms behind early life-stage variation in morphology, we examined how craniofacial shape varied according to family, offspring size, egg size and candidate gene expression.

RESULTS

Craniofacial shape (i.e. the Meckel's cartilage and hyoid arch) was more variable between families than within families both across and within developmental stages. Differences in craniofacial morphology between developmental stages correlated with offspring size, whilst within developmental stages only shape at FF correlated with offspring size, as well as female mean egg size. Larger offspring and offspring from females with larger eggs consistently had a wider hyoid arch and contracted Meckel's cartilage in comparison to smaller offspring.

CONCLUSIONS

This study provides evidence for family-level variation in early life-stage trophic morphology, indicating the potential for parental effects to facilitate resource polymorphism.

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.

Comparative transcriptomics reveals candidate carotenoid color genes in an East African cichlid fish

BACKGROUND

Carotenoids contribute significantly to animal body coloration, including the spectacular color pattern diversity among fishes. Fish, as other animals, derive carotenoids from their diet. Following uptake, transport and metabolic conversion, carotenoids allocated to body coloration are deposited in the chromatophore cells of the integument. The genes involved in these processes are largely unknown. Using RNA-Sequencing, we tested for differential gene expression between carotenoid-colored and white skin regions of a cichlid fish, Tropheus duboisi "Maswa", to identify genes associated with carotenoid-based integumentary coloration. To control for positional gene expression differences that were independent of the presence/absence of carotenoid coloration, we conducted the same analyses in a closely related population, in which both body regions are white.

RESULTS

A larger number of genes (n = 50) showed higher expression in the yellow compared to the white skin tissue than vice versa (n = 9). Of particular interest was the elevated expression level of bco2a in the white skin samples, as the enzyme encoded by this gene catalyzes the cleavage of carotenoids into colorless derivatives. The set of genes with higher expression levels in the yellow region included genes involved in xanthophore formation (e.g., pax7 and sox10), intracellular pigment mobilization (e.g., tubb, vim, kif5b), as well as uptake (e.g., scarb1) and storage (e.g., plin6) of carotenoids, and metabolic conversion of lipids and retinoids (e.g., dgat2, pnpla2, akr1b1, dhrs). Triglyceride concentrations were similar in the yellow and white skin regions. Extracts of integumentary carotenoids contained zeaxanthin, lutein and beta-cryptoxanthin as well as unidentified carotenoid structures.

CONCLUSION

Our results suggest a role of carotenoid cleavage by Bco2 in fish integumentary coloration, analogous to previous findings in birds. The elevated expression of genes in carotenoid-rich skin regions with functions in retinol and lipid metabolism supports hypotheses concerning analogies and shared mechanisms between these metabolic pathways. Overlaps in the sets of differentially expressed genes (including dgat2, bscl2, faxdc2 and retsatl) between the present study and previous, comparable studies in other fish species provide useful hints to potential carotenoid color candidate genes.

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.

A way forward with eco evo devo: an extended theory of resource polymorphism with postglacial fishes as model systems

A major goal of evolutionary science is to understand how biological diversity is generated and altered. Despite considerable advances, we still have limited insight into how phenotypic variation arises and is sorted by natural selection. Here we argue that an integrated view, which merges ecology, evolution and developmental biology (eco evo devo) on an equal footing, is needed to understand the multifaceted role of the environment in simultaneously determining the development of the phenotype and the nature of the selective environment, and how organisms in turn affect the environment through eco evo and eco devo feedbacks. To illustrate the usefulness of an integrated eco evo devo perspective, we connect it with the theory of resource polymorphism (i.e. the phenotypic and genetic diversification that occurs in response to variation in available resources). In so doing, we highlight fishes from recently glaciated freshwater systems as exceptionally well-suited model systems for testing predictions of an eco evo devo framework in studies of diversification. Studies on these fishes show that intraspecific diversity can evolve rapidly, and that this process is jointly facilitated by (i) the availability of diverse environments promoting divergent natural selection; (ii) dynamic developmental processes sensitive to environmental and genetic signals; and (iii) eco evo and eco devo feedbacks influencing the selective and developmental environments of the phenotype. We highlight empirical examples and present a conceptual model for the generation of resource polymorphism - emphasizing eco evo devo, and identify current gaps in knowledge.

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.