- Open Access
Maternal mRNA input of growth and stress-response-related genes in cichlids in relation to egg size and trophic specialization
EvoDevovolume 9, Article number: 23 (2018)
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.
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).
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.
Parental investment is an important determinant of reproductive success that increases the fitness of the offspring at the cost of parental fitness . The evolutionary trade-off between the two results in parent–offspring evolutionary conflict, which is a cornerstone of kin selection theory . Parental investment is found across the broad taxonomic range of the animal kingdom and can be biparental or exclusively uniparental. It can be divided into two main categories: mating investment and rearing investment. In oviparous organisms, such as fish, maternal mating investment is significantly determined by egg size and clutch size trade-off, optimally modulated by environmental conditions . This investment can affect traits functionally, developmentally, physiologically related to offspring size and survival  and also influence later life stages . Thus, the maternal phenotype and the environment she experiences hold adaptive value as they can have important implications for the fitness of offspring.
Egg size is a life-history trait that affects both maternal and offspring fitness, simultaneously . Egg size is determined by a complex interplay of long-term local adaptation and short-term plasticity balancing egg size with brood size [4, 7, 8]. This makes it an evolutionarily important maternal effect [4, 9]. During the female reproductive cycle in oviparous fishes, multiple stages with similar basic patterns occur, such as follicular growth, oocyte maturation and ovulation as well as secondary growth due to incorporation of yolk material into the oocyte (vitellogenesis) . Changes in oocyte growth lead to differences in egg size, which is an important indicator of female energetic investment in reproduction and subsequently affect the development, growth, behaviour and fitness of the offspring [11,12,13,14,15,16,17,18,19]. A number of interconnected factors are implicated in variation of egg size including female age, size, fecundity, energy acquisition and batch sequence at the time of egg production [20,21,22,23,24] as well as environmental stressors like population density, water temperature, food availability, oxygen level, osmotic changes and predation risk [16, 25,26,27,28,29,30,31]. Life-history theories postulate that female investment in offspring quality is balanced with their risk of survival . The molecular mechanisms linking environmental adaptation, e.g. adaptive foraging strategies or stress response, and production of eggs variable in size are poorly understood. Importantly, it is also unclear whether eggs from different species with variable sizes carry different levels of maternally deposited mRNA transcripts triggering particular functions in development, growth and trophic adaptation.
Adaptive radiation is the rapid and simultaneous diversification of a lineage into an array of ecologically disparate species . The haplochromine cichlids from Lake Tanganyika (LT), Lake Malawi (LM) and Lake Victoria (LV) represent the most rapid and species-rich adaptive radiations known . The haplochromine cichlids not only exhibit novel and diverse morphological, ecological and behavioural adaptations, but several of these traits exhibit evolutionary parallelism . Thus, they provide an excellent framework for comparative studies on phenotypic evolution in triplicate radiations of varying age (LT 10–12 million years ago/MYA, LM 2.4–5 MYA, LV 15,000–100,000 years) . In addition, the riverine cichlid species in our study are known to be ancestral to the haplochromine lineages of the three lakes [36,37,38]. Remarkably, cichlids have highly complex breeding behaviours and all cichlids display a variety of parental investment in the form of egg and larvae care . Maternal mouthbrooding is one such key innovation that characterizes the haplochromines and strongly limits the number and size of eggs that can be incubated . This also means that maternal investment in reproduction is high. Interestingly, life-history traits, e.g. egg size and fecundity, in mouthbrooding cichlids were shown to evolve in parallel to different habitats across the 3 radiations . It was also recently reported that various life-history traits such as clutch size and egg mass in a riverine haplochromine, Astatotilapia burtoni, were diverging at the population level, suggestive of ‘ongoing’ adaptive radiation . Furthermore, haplochromines have evolved numerous trophic adaptations to acquire different food sources such as algae, snails, insects and shrimp . Since both the requirements of egg incubation and feeding are fitness-related traits that constrain size/shape of the buccal cavity, it is hypothesized that there is an evolutionary trade-off between the two at the species level . This raises the question of whether the observed traits, i.e. the divergence in trophic specialization and variations in the egg size, can be associated with differences in maternally supplied molecules, such as gene transcripts, with functions in environmental adaptation as well as development, growth and morphogenesis. Also it is interesting to know that maternal provisioning through the deposition of mRNA transcripts of genes with relevant functions can explain potential role for non-Mendelian inheritance in producing an evolutionary trade-off across species. Finally, since trophic specialization, egg size, mouth incubation time and size at independence might be correlated, it is to know if mothers are able to control this process by selective deposition of mRNA transcripts regulating embryonic growth.
The contribution of mRNA transcripts from the mother to the egg is a less studied aspect of maternal investment in offspring fitness. Maternal mRNA transcripts are among the main components deposited in eggs, which are not only critical for early embryonic development, but also affect developmental events after activation of zygotic gene transcription (or maternal to zygotic transition stage; MZT) and influence subsequent developmental life-history trajectories . At early stages of development, maternal RNAs control cellular programmes required for crucial events such as cleavage, blastula formation and gastrulation [45,46,47]. By the initiation of embryonic transcription during MZT up to half of the maternal RNAs are degraded [48, 49]; however, even at the peak zygotic transcription still majority of mRNA transcripts can have maternal origin in fish . Interestingly, a portion of maternal mRNAs can persist in embryos long after the MZT and function beyond this stage during embryonic morphogenesis [51, 52]. Notably, the products of certain maternal mRNAs might themselves control RNA degradation processes and by postponing the degradation time of specific maternal mRNAs exert their effects on later events of embryonic morphogenesis . It is worth emphasizing that the influence of maternal mRNAs on later developmental events after MZT might not necessarily depend on the presence of the mRNA molecules in later stages since their translated products can persist for much longer period. For instance, persisting activity of an enzyme (involved in retinoic acid production) translated from maternally deposited mRNA was found to be essential for pancreas development in zebrafish . In addition, the tuning of earlier developmental patterning by maternal mRNAs triggers a variety of molecular cascades and cellular processes that organize and guide later developmental events long after clearance of maternal mRNAs, as observed for highly interconnected maternal and zygotic control of dorsoventral patterning in zebrafish . Finally, maternally deposited mRNAs might affect the epigenetic programming during development which can later influence a variety of morphogenetic processes . The molecular mechanisms by which maternal mRNAs influence post-MZT, and particularly, late embryonic development and morphogenesis are poorly explored, although a growing number of studies have begun to focus the spotlight on these processes.
In annual killifish, for instance, different splice variants of maternally supplied gene transcripts have been shown to play a role in directing divergent developmental trajectories during somitogenesis and later affecting plastic responses of embryos to environmental stimuli . The level of maternal mRNA transcript for genes determining ventral embryonic fates was found to be important for the morphogenesis of retina and skeletal structures surrounding the eyes long after initiation of zygotic transcription in cavefish embryos . In round goby, the water temperature experienced by the mother before oviposition leads to selective deposition of maternal mRNAs for temperature-responsive genes which contribute to adaptation to temperature changes during and even beyond embryogenesis . Another recent study in rainbow trout has shown that thermal stress in mother regulates maternal mRNA deposition of genes contributing to neural development and the acquisition of neurocognitive function in embryos which can later shape inter-generational behaviour accordingly . The thermal-induced inter-generational adaptation driven by changes in maternal mRNA and egg size is also reported in stickleback . Maternally deposited transcripts of a nuclear receptor gene appeared to be able to regulate the epigenetic programming of development and indirectly influencing a range of morphogenetic processes in brain, heart, eye and skeletal system . Other examples for association of maternal mRNA deposition with later developmental phenotypes in fish include differences in the length of offspring that contain different levels of maternal mRNAs for specific growth-related genes in brook charr  and zebrafish .
The deposition of maternal mRNA transcripts occurs at different stages of oocyte formation and is thought to be correlated with certain properties of the egg such as size and its fertilization capacity . Inter-species comparison in amphibians has shown that evolutionary increase in egg size can result in changes in patterns of RNA localization in egg and early embryonic development [64, 65]. Moreover, various studies of teleost fishes have addressed the mechanisms linking egg size variation to changes in fecundity, egg number, follicular development and expression of growth-related genes [20, 66, 67]. Interestingly, several stages of the egg formation process are known to be highly conserved across vertebrates [68,69,70,71] and cross-species analysis of mature oocytes has identified the maternal deposition of conserved RNA transcripts across chordates [71, 72]. Growing evidence indicates that a stress-mediated molecular pathway, glucocorticoid (GC) signalling, can be a major determinant of egg size, as well as offspring growth and survival capacity in response to environmental changes [31, 73,74,75,76]. However, it is unclear if the GC-mediated changes in egg size are also accompanied with differential deposition of maternal mRNA transcripts for related genes which could then explain how the above-mentioned effects of GC pathway activation in mother during oogenesis are conveyed to the offspring. In addition, interconnected signals mediated by growth hormone and insulin-like growth factors, GH/IGF axis, are well known for their role in egg production under a regulatory mechanism of ghr expression (cognate GH receptor) . The axis can be self-limiting during oogenesis since GH-dependent activation of IGF can initiate a feedback mechanism controlling GH production [78, 79]. The level of maternal mRNA components of the GH/IGF axis has also been found to be correlated with egg properties and embryo survival [80, 81]. This indicates that distinct activity of a pathway affecting egg properties during oogenesis might coincide with differential deposition of maternal mRNAs for genes related to the pathway extending the effects to offspring.
In this study, we explored the level of maternal mRNA deposition of selected components of GC and GH/IGF pathways in the eggs of 15 haplochromine cichlid fish species with variable egg sizes and trophic specialization. The species included in this study belong to distinct trophic groups covering three independent radiations, Lake Tanganyika (LT), Lake Malawi (LM), Lake Victoria (LV), and riverine species (RV). The objectives of our study were, first, to investigate whether the eggs from closely related cichlid species contain different levels of maternal mRNAs for components of the two crucial pathways, second, to find out if the potential variations in the maternal mRNAs could be associated with differences in trophic niche specialization or/and egg size, and third, to provide technically accurate quantification method for further investigations in this topic using gene expression analysis. Our results provide first cross-species comparisons of maternal mRNA investment in relation to egg size, habitats and trophic niches in fish as well as a repository of validated reference genes with stable mRNA abundances across the eggs, for accurate normalization of qPCR data, for future studies.
Egg sampling and measurements
In this study we used 15 haplochromine cichlid species, endemic to different habitats in Eastern Africa, including 3 riverine species, 5 species from Lake Malawi, 3 species from Lake Victoria and 4 species from Lake Tanganyika (Fig. 1a). Within each lake, we had at least one herbivorous and one carnivorous species for trophic niche comparison, and a simplified representation of phylogenetic relationships between the species is depicted in Fig. 1a based on previous studies [35, 37]. The fish were raised in separate tanks of mixed-sex species-specific communities (10–14 individuals per tank, with higher ratio of females) with similar light and water conditions and receiving the same diet (Spirulina flakes for cichlids) until young adult stage when mating behaviour was first observed. From this time, we carefully monitored all individuals in each tank on an hourly basis every day, in order to identify any mating pair during the spawning period (up to 3 h). Immediately after the end of spawning, we removed the eggs from the mouth of the females using moderate manual pressure on their cheeks. A single female per species was used for sampling, and from each female up to 8 eggs with no deformities were selected and the females were measured for their standard length (Additional file 2). It should be emphasized that in this study we only look at cross-species comparisons and therefore a single female per species with at least 8 high-quality eggs was used from their first or second batch of laying eggs. Prior to weighing, eggs were checked under microscope for fertilization success and those in stages later than 2-cell cleavage were discarded [82, 83]. One by one, the eggs were quickly dried on a cotton pad in order to remove surface water then weighed to be immediately transferred to RNAlater solution (Ambion, USA) and kept at − 4 °C overnight. Each egg represented a biological replicate and therefore processed separately, giving a total of 8 replicates per breeding pair and species.
RNA isolation and cDNA synthesis
Each egg was transferred from RNAlater solution to a tube with 1.5 mL Qiazol lysis reagent (Qiagen, Hilden, Germany), and a 1.4-mm ceramic sphere was added in the tube to crush the egg. Then, we homogenized the eggs using FastPrep-24 Instrument (MP Biomedicals, Santa Ana, CA, USA) and we extracted RNA based on the protocol developed by the manufacturer suitable for tissues with high fat content. After three washes with 70% ethanol, the RNA was dissolved in 30 µl of nuclease-free and RNase-free water. The genomic DNA contamination was removed using DNase I enzyme (New England Biolabs), and the quantity and quality of each RNA samples were, respectively, checked by Nanophotometer (IMPLEN GmbH, Munich, Germany) and RNA ScreenTapes on an Agilent 2200 TapeStation (Agilent Technologies). First-strand cDNA was synthesized from 500 ng of RNA through High-Capacity cDNA Reverse Transcription kit (Applied Biosystems), following the protocol developed by the manufacturer. The final cDNA product for each sample was diluted 1:4 times in nuclease-free and RNase-free water and was used for qPCR.
Selection of candidate genes and primer design
We selected 8 candidate reference genes with highest maternal mRNA abundances known in zebrafish eggs, combined with least degradation from the fertilization stage onwards . In addition, we added five more classic reference genes, which are frequently used in studies of different tissues in African cichlids [85,86,87,88]. For target genes we selected 6 genes, including two paralogs of growth hormone receptor, ghr1 and ghr2; two insulin growth factors, igf1 and igf2; as well as two components of glucocorticoid pathway, gr (also known as nr3c1) and sgk1.
We designed the qPCR primers matching sections of the sequences that are conserved across East African cichlids, based on transcriptome data from 6 haplochromine species (Pundamilia nyererei, Simochromis diagramma, Tropheus duboisi, Gnathochromis pfefferi, Metriaclima zebra and Astatotilapia burtoni) and two species from more distant tribes (Oreochromis niloticus and Neolamprologus brichardi) [33, 89]. For this purpose, the coding sequences from all of the above-mentioned species were first aligned in CLC Genomic Workbench, version 7.5 (CLC Bio, Aarhus, Denmark), and then exon/exon junctions were identified based on the annotated genome of Oreochromis niloticus, which is available online in the Ensembl database (http://www.ensembl.org) . The primers were designed spanning over these junctions with an amplicon length inferior to 250 bp in order to remove potential effects of possible RNA degradation on qPCR quantification . Primer Express 3.0 (Applied Biosystems, CA, USA) and OligoAnalyzer 3.1 (Integrated DNA Technology) were used to design the primers and for the evaluation of their dimerization and secondary structures.
qPCR and data analysis
The protocol established by Maxima SYBR Green/ROX qPCR Master Mix (2X) (Thermo Fisher Scientific, Germany) was applied to prepare qPCR reactions, which were conducted in 96-well PCR plates and using ABI 7500 real-time PCR System (Applied Biosystems). Two technical replicates were allocated for each biological replicate following an experimental set-up described as sample maximization method to reach to optimal conditions for a qPCR run . The programme assigned for each qPCR run was initiated with steps of 2 min at 50 °C and a 10 min at 95 °C, continued by 40 cycles of amplification, consisting of 15 s at 95 °C and 1 min at 60 °C. In addition, a dissociation step (60–95 °C) was performed at the end of the amplification step. We also calculated primer efficiency for each primer pair value through LinRegPCR v11.0 program (http://LinRegPCR.nl) , and we redesigned primers for pairs showing efficiency values less than 0.9 (Additional file 1).
We used three software packages, with different algorithms, to rank the most stable reference genes: BestKeeper , NormFinder  and geNorm , as well as a ranking based on standard deviation (SD) using raw quantitation cycle (Cq) values. The higher r index, in BestKeeper ranking, indicates better stability, whereas in geNorm and NormFinder the lower values, respectively, identified by M (average expression stability values) and SV (stability values) suggest more stable reference genes. The average Cq values of the two most stable reference genes were used as normalization factor or Cqreference, and for each gene, ΔCq was calculated by subtracting Cq values of the target genes from the corresponding value of the reference genes (ΔCqtarget = Cq target − Cqreference). The data from all the egg samples were then normalized to the ΔCq value of a calibrator sample, to obtain a ΔΔCq value (ΔCqtarget − ΔCqcalibrator). The egg sample with the lowest mRNA abundance (highest ΔCq) was used as a calibrator sample for each gene across all of the species. Relative mRNA quantities (RQ) were calculated based on the expression level of the calibrator sample (E−ΔΔCq) , and then RQ values were transformed to logarithmic base 2 values (also known as fold differences; FD)  for statistical analysis. In order to identify species-specific, lake-specific and trophic niche-specific differences in mRNA abundances, we conducted ANOVA statistical tests, followed by Tukey’s HSD post hoc tests using FD values (Additional file 2). To assess correlation between egg weights and mRNA abundances, Pearson correlation coefficients (r) were calculated for each gene in respect of the egg weights across the species. Also, Pearson correlation coefficients (r) were used to calculate similarity in patterns of mRNA abundances between all gene pairs. All statistical analyses were implemented in R (http://www.r-project.org).
Weight differences of the haplochromine cichlid eggs
The egg measurements across all the species revealed a positive correlation between female length and egg weight; however, this was solely due to a single outlier, the LT species P. famula, and the removal of this species from the analysis obliterated the correlation between these variables (Fig. 1b). A comparison of corrected egg weights (using egg weight residuals against female standard lengths) across the lakes showed that LT species had the largest eggs among the selected species followed by LM species with the second largest eggs, while RV species appeared to have smallest eggs (Fig. 1c). In addition, LM and RV species displayed far less variations in egg weight than LT and LV species. Comparisons of egg weight in relation to trophic niche revealed larger eggs for herbivorous than carnivorous species across the lakes except for one of the carnivorous species in LV (Fig. 1d–f); notwithstanding that, more species per trophic niche from each lake are required to confirm this as a general trend in the eggs of haplochromine cichlids.
Suitable reference genes for quantification of maternal mRNA
In order to accurately measure mRNA abundance in eggs and be able to compare the quantities across the cichlid species, we needed to first identify reference genes with least variations in mRNA abundance among the eggs of different species . To select reference gene candidates, we retrieved transcriptome data available for zebrafish maternal mRNA and identified 8 genes displaying highest mRNA abundances with no significant variations from post-fertilization stage to the initiation of zygotic transcription (mid-blastula stage or MBT which happens at the beginning of MZT stage described above) . The study was conducted on a large number of zebrafish embryos harvested at five different developmental stages: 1-cell, 16-cell stage, 128-cell stage, MBT and post-MBT . In addition, we included 6 more candidates, actb1, ef1a, gapdh, rps11, rps18 and tbp, which are shown to have high expression in different tissues and were used as stable reference genes in a variety of studies on African cichlids [85,86,87,88]. Expectedly, the genes had a range of mRNA abundances, from highest levels of transcripts for mid1ip1, actb1 and ccnb2 to lowest levels for rps18, tbp and ef1a (Fig. 2a). When compared to zebrafish transcriptome data , we did not find significant correlation between the mRNA abundances of the genes in cichlid and zebrafish eggs (Fig. 2a). However, five genes, mid1ip1, mylipa, dvr1, tatdn2 and rps18 showed high to low mRNA abundances, respectively, in the cichlids and zebrafish, which are specified with triangles in Fig. 2a. It should be noted that the higher Cq values indicate lower mRNA abundances in the cichlid eggs, whereas higher RPKMs mean higher mRNA abundances in zebrafish; therefore, the opposite pattern between the two value types indicates similarity in mRNA abundance between the cichlids and zebrafish (Fig. 2a). In addition, ef1a with lowest mRNA abundance in the cichlid eggs had also below cut-off level of sequencing reads in zebrafish egg transcriptome , which might indicate partial conservation in the levels of maternal mRNA investment for certain genes across distant teleost taxa.
Among the candidate reference genes mid1ip1 had the lowest standard deviation (SD) in the eggs from different cichlid species followed by ccna1 and atf7ip as second and third most stable genes in this ranking (Table 1). After selecting only the top 8 genes with lowest SD to run on BestKeeper, atf7ip, ccnb2 and mid1ip1 were ranked as top three genes by this software. NormFinder ranked atf7ip, ccnb2 and actb1 as the top three most stable genes. On the other hand, geNorm suggested only atf7ip and mid1ip1 as suitable reference genes since their M value was below the recommended threshold (M < 1.5). Importantly, each of these two genes alone showed significant differences in at least one of the comparisons of the trophic niches and/or the lakes, but their average Cq values did not show differences in these comparisons (Additional file 3: Figure 1). Hence, we used average Cq value of atf7ip and mid1ip1 in each egg as normalization factor (NF) for further quantifications of mRNA abundance. It should be noted that NF was also ranked as most stable normalizing value when applied in all three ranking algorithms (Table 1).
Maternal mRNA abundance of candidate target genes
The mRNA abundance of 6 target genes, igf1, igf2, ghr1, ghr2, gr and sgk1, was examined across the eggs, and we found igf1 transcripts to be below the minimum detection level through qPCR in all of the species in this study. The lake-specific comparisons of mRNA abundances revealed lower mRNA investment of igf2 in LT with largest eggs (Fig. 2c). We also found lower mRNA abundance of ghr2, gr and sgk1 in LV compared to the other habitats. In contrast, ghr1 showed higher mRNA abundance in LV and LM compared to LT and RV. These results reveal clear variations in maternal mRNA investment for all of the candidate target genes, depending on the habitats of closely related haplochromine cichlid fishes. Moreover, the two paralogs of a growth hormone receptor, ghr1 and ghr2, displayed opposite expression patterns, reflecting contrasting maternal mRNA investments, which suggests lake-specific divergence of their role during early development. The lower transcripts of gr and sgk1 might imply distinct activation of GC pathway in LV species during deposition of maternal mRNA in the eggs.
The comparisons between herbivorous and carnivorous species revealed differences in mRNA abundance of the target genes in each lake (except for ghr1 in LV) (Fig. 3). A similar pattern of higher mRNA abundance was observed between a carnivorous and herbivorous species for ghr2 and gr in LM, and for gr and sgk1 in LT. Furthermore, a consistent pattern of higher mRNA abundance in carnivorous species was observed for igf2 in LV and LT, but this pattern was reversed in LM (with higher level of transcripts in herbivorous species). Among the target genes, only sgk1 showed a tendency towards higher mRNA abundance in carnivorous species across all lakes; however, the difference was not statistically significant in all of the species with contrasting trophic niches.
Correlation patterns of target gene maternal transcripts
We assessed the correlations between maternal mRNA abundance of the target genes and egg weights across all the haplochromine species in this study. Interestingly, the mRNA abundance of two genes, ghr2 and igf2, displayed significant positive and negative correlations with egg weights, respectively (Fig. 4a). This indicates inter-species differences in maternal mRNA deposition of the two growth-related genes, with respect to the egg mass in cichlids. Next, we conducted pairwise comparisons of mRNA abundance between the target genes in order to identify potential transcriptional regulatory connections linking the genes (Fig. 4b). Expectedly, we found a positive expression correlation between the mRNA abundances of gr and sgk1, which both belong to GC signalling pathway, confirming the accuracy of our qPCR analyses. Furthermore, both gr and sgk1 mRNA abundances appeared to have positive correlations with ghr2 mRNA abundance (but not with ghr1), suggesting a potential regulatory connection between ghr2 and activated GC pathway in the cichlid eggs. In contrast, ghr2 displayed negative correlation with igf2 and ghr1 mRNA abundances, while igf2 and ghr1 had a positive correlation across the eggs. These might be the result of transcriptional decoupling between ghr2 and GH-IGF regulatory axis due to a mechanism affecting only ghr2 transcription and maternal deposition.
Finally, we performed a hierarchical clustering of the target genes based on their mRNA abundances comparing the relationships between the species (Fig. 4c). Strikingly, the dendrogram showed the presence of two main clusters dividing LV and LM species from RV and LT species, implying similar patterns of maternal mRNA abundance in LT and RV species, while they are both distant from the other two lakes. In RV–LT cluster, the three LT herbivorous species (T. duboisi, P. famula and S. diagramma), with larger eggs, were sub-clustered together, whereas the LT carnivorous species (G. pfefferi), with smaller eggs, was sub-clustered with RV species. Two sub-clusters were also observed in LV–LM cluster, and in one sub-cluster, three of the LM species were grouped together, whereas in the other sub-cluster, the three LV species grouped with two other LM species. A closer look into the branching pattern of the dendrogram revealed that the two LM carnivorous species (A. hansbaenschi and S. fryeri) with an LM omnivorous species (C. afra) constitute one of the sub-clusters. Interestingly, in the other sub-cluster, the herbivorous species of both LM and LV (L. trewavasae and N. omnicaeruleus) are grouped together, whereas the two LV carnivorous species (P. sp. Red head and H. thereuterion) and a LM omnivorous species (M. zebra) constitute another group. These observations demonstrate that, based on maternal mRNA abundances of only the five target genes, a clear distinction is observed between LT and RV versus LM and LV. In addition, within each cluster a distinct grouping between herbivorous and carnivorous species is observed for each lake.
Differences in the level of maternal mRNA transcripts between eggs can be a result of distinct maternal provisioning for mRNA deposition during oogenesis. On the other hand, the amount of mRNA transcripts in eggs starts to change soon after fertilization and substantially declines from the initiation of zygotic transcription (mid-blastula transition) due to degradation process affecting most maternal gene transcripts . In zebrafish, a proportion of the maternal mRNA appeared to be degraded prior to mid-blastula transition (~ 3.5 h post-fertilization; hpf) . It should be noted, however, that the maternal mRNAs still account for a significant part of functional mRNA transcripts during zygotic transcription and some could persist in later stages of development as well (see introduction). The great diversity of cichlid fishes in trophic specialization, manifested in distinct morphology, growth, developmental patterning, egg size, provides an excellent model system to investigate potential links between maternal mRNA provisioning, life history and ecomorphology. Such studies can provide insight into the role of life-history traits and maternal effects in adaptive radiations. In cichlid fishes studied so far, the developmental stage comparable to mid-blastula transition occurs later in time than zebrafish, e.g. around 12 hpf in two African cichlids, Oreochromis niloticus and Astatotilapia burtoni [82, 100], and 4 hpf in two American cichlids, Symphysodon aequifasciatus and Amphilophus spp. [83, 101]. Therefore, we set out to investigate the maternal mRNA investment in cichlid eggs almost immediately after fertilization (up to 3 hpf) to minimize the effects of RNA degradation during early development. This early stage of cichlid embryogenesis also has implications for the developmental hourglass hypothesis where species are evolutionarily most divergent at the early and late stages of embryonic development, with a conserved stage in between .
Importance of reference genes for accurate quantification of maternal mRNA abundances in eggs
In order to accurately quantify maternal mRNA abundance of our target genes of interests we had to identify reference gene transcripts showing no signs of degradation in the early developmental stages. This is of paramount importance when studying relative gene expression in different tissues. Among the candidate reference genes, which were selected based on their resistance to RNA degradation and high level of maternal transcripts in zebrafish eggs , we found the cichlid homologous of two genes, atf7ip and mid1ip1, to have the most stable level of mRNA abundance across the eggs of haplochromine cichlids. The first gene, atf7ip (also known as mcaf1), encodes a transcriptional modulator participating in histone methylation process and chromatin organization  and contributes to the maintenance of X chromosomes inactivation during female cell development and division in mammals . In bovine oocyte, atf7ip has been implicated in the repression of oocyte-specific genes when the embryonic gene expression is beginning to take place . It is worth emphasizing that the maternal methylation imprinting can play an important role both during oogenesis and in later stages of embryonic development [106, 107]. In our study, the consistent mRNA abundance of atf7ip across all of the haplochromine species suggests its tight regulation during final stages of oogenesis when the maternal mRNA deposition takes place and its transcript stability in early development.
The second validated reference gene, mid1ip1 (mig12), encodes a protein required for efficient lipid biosynthesis, stabilization of microtubules in cytoskeleton and embryonic midline development [108, 109]. In zebrafish a mutation in this gene has demonstrated that the maternally provided mid1ip1 transcripts are crucial for cytoskeletal dynamics and membrane recycling in early cell divisions, and during the reorganization of the cytoskeleton at the egg-to-embryo transition . Therefore, the observed low variation in abundances of mid1ip1 transcripts across the eggs implies its tight regulation during late oogenesis, uniform maternal deposition across species and potential requirement in early development of the cichlid species. In addition to these candidate genes, we also tested the maternal abundance of 5 classic reference genes, frequently used in different studies of cichlid tissues, and we found that none of them appeared among the top validated reference genes with different ranking methods and therefore were not suitable for analysis of qPCR data obtained in this study. This once again underscores the necessity of careful validation of reference genes in each qPCR study, according to different experimental conditions, tissues and species . The reference genes validated here will be a useful resource for future studies on this topic.
Egg weight may be an indicator of adaptive divergence in cichlids
Given the importance of egg size as an indicator of maternal investment and life-history strategies, we were interested in investigating the relationship between mRNA abundance of a selected set of target genes, involved in early developmental growth and mediating the effects of environmental signals, and the egg size in haplochromine cichlids. Therefore, we first had to characterize the eggs based on their variations in weights across habitats and trophic niches. Indeed, the smallest eggs were found in the riverine haplochromines, while the largest eggs were observed in species of the oldest species flock in Lake Tanganyika. Also herbivorous species tended to have larger eggs than carnivorous species. However, the number of species with contrasting trophic niches per lake was not enough in our study to conclude with confidence that the observed differences are consistent throughout all haplochromine cichlids. A plausible explanation for the increase in egg size in all lake habitats, as well as for the increase with evolutionary age, might be adaptation to more predictable niches in lake habitats, as opposed to much more fluctuating environmental parameters in riverine environments.
Species of the youngest lake deposit the lowest amount of maternal mRNAs for stress-mediating genes
The glucocorticoid (GC) pathway can be considered as a prime candidate molecular signal affecting processes involved in ovulation, egg size determination, offspring growth and survival in response to environmental changes [31, 73,74,75,76, 112,113,114,115,116,117]. Concerning the target genes in our study, we examined the maternal mRNA abundance of two components of GC signalling pathway: a nuclear receptor, gr (or nr3c1), which also acts as an upstream transcriptional regulator of GC pathway, and a downstream effector of the pathway, sgk1, which is a kinase activating certain ion channels and mediating cellular stress response [118,119,120]. The activity of enzyme encoded by sgk1 is regulated by mTOR signalling pathway [121, 122], which is a modulator of oogenesis (primordial follicle activation and arrest) . In this study, we found a positive correlation between mRNA transcript abundances of gr and sgk1 genes across all the eggs confirming their regulatory connection during early development, prior to mid-blastula transition. Moreover, both genes had very similar patterns of transcript abundance across the lakes with lowest level of abundance in LV. This proposes a potential difference in response to environmental stressors, which can be reflected in maternal deposition of GC pathway gene transcripts, and possibly, their distinct effects on early development between the cichlid species in LV and those in the other habitats. Differences in sgk1 expression have been already suggested as an adaptive response to environmental stressors in vertebrates [124, 125], and for instance, it can be directly regulated by osmotic changes . This might be a highly conserved mechanism linking environmental cues to maternal provisioning since sgk1 has been found to be a regulator of fat storage (in eggs), embryonic growth and stress response in Caenorhabditis elegans as well [127, 128]. The maternal gr transcripts have been shown to be essential for developmental programming of different tissues, in particular for skeletal system development in zebrafish embryos . Thus, it would be interesting to know whether LV cichlid species display distinct skeletogenesis during embryonic development than other haplochromine species in the other lakes. Although the activation of GC pathway can affect egg size, we did not find any correlation between maternal abundance of gr and sgk1 with egg weight across the species, suggesting that such effect does not necessarily lead to changes in maternal deposition of the GC-related mRNAs, and thus studies of pre-ovulatory stages in the ovary are required to address this role in cichlids.
Larger eggs contain higher amount of maternal ghr2, but not ghr1 across haplochromine cichlids
There are two growth hormone (GH) receptors in teleost fish, ghr1 and ghr2, and even though both receptors show similarities in signal transduction and transcriptional regulation of certain downstream genes, there is evidence for their tissue-specific function and expression as well [129,130,131,132,133,134]. In particular, the overall expression level of ghr2 has been found to be higher than ghr1 in fish gonadal tissues [77, 131, 135], and during ovary reproduction cycle in cichlid, ghr2 shows higher expression at sexual recrudescent and regressed stages, whereas at sexual matured stage there is an increase in expression of both receptors without significant difference between them . This could indicate a stage-specific regulatory role of ghr2 during ovary reproduction cycle, as well as a regulatory role of both receptors in oocyte maturation. In our study, we found higher abundance of ghr1 maternal transcripts in LV and LM, whereas this pattern was reversed for ghr2 maternal transcripts (i.e. higher ghr2 in RV and LT eggs). This suggests a similarity between LT and RV haplochromine species in maternal deposition of ghr genes, and a possible evolutionary divergence of LV and LM from this pattern. These results support lineage-specific specialization of ghr genes in younger (LM and LV) versus older (LT) and seeding lineages (RV). Interestingly, ghr1 and ghr2 displayed negative correlation in maternal mRNA abundance across all the eggs, which can be a result of distinct regulatory effects of upstream factor(s) on their transcription during oocyte maturation and subsequent maternal mRNA deposition. It is already known that several factors with a role in oogenesis, such as cortisol, insulin and GH, also have distinct effects on ghr1 and ghr2 transcription in fish [131, 136]. For instance, insulin can repress ghr2 transcription, while inducing ghr1 transcription in cichlid liver cells . Furthermore, we found that larger eggs had higher abundance ghr2, but not ghr1, across the haplochromine cichlids. This is in agreement with the results of an intra-specific study of maternal mRNA investment in a haplochromine species with variable egg sizes in which ghr (not specified at the paralog level) had shown higher mRNA abundance in larger eggs . Interestingly, in the same study, individuals from smaller eggs with lower maternal ghr abundance later expressed higher ghr expression during larval development, a compensatory mechanism through which they could catch up in size with those individuals originating from larger eggs . Future inter-specific studies are required to know whether such a mechanism is ubiquitous in haplochromine cichlids. Our results, together with observations in other studies of GH receptors, suggest inter-specific and egg-size-dependent mechanism underlying maternal ghr2 transcript deposition which might contribute to adaptive developmental events in cichlids.
Larger eggs contain lower amount of maternal igf2 across haplochromine cichlids
Finally, we investigated the maternal investment for two insulin growth factors, igf1 and igf2, across the species. The proteins encoded by these paralog genes are functionally and structurally related to insulin hormone and play diverse roles in mediating growth and development in vertebrate including physiological processes controlling fish reproduction . Both factors, for instance, play a role in the formation of vitellogenic follicles and the promotion of oocyte maturation in fish ovaries [138,139,140,141]. The maternal mRNA abundances of igf1 and igf2 can differ in fish as greater abundance of igf2 mRNAs compared to igf1 was observed in follicles [44, 142] and post-fertilized eggs . Moreover, the oocyte investment of igf1 mRNA had been shown to be very little across distant teleost fish taxa [137, 142]. Similarly, we found maternal igf1 abundance to be below the minimum detection level by qPCR across the cichlid eggs, whereas maternal igf2 mRNA was detected in all of the eggs. We also found lower amounts of maternal igf2 transcripts in larger eggs which have never been reported in fish, proposing an unknown mechanism during maternal provisioning. To our knowledge, it has only been shown in porcine ovaries that the largest follicles contain lower igf2 expression, when compared to average or small size follicles ; however, the mechanism underlying this expression pattern and its potential impacts on embryonic development remains unknown. We observed negative correlation in mRNA abundance between igf2 and ghr2, which could implicate shared upstream regulator(s) influencing maternal mRNA deposition of both genes in an opposite manner. A candidate for such an upstream regulator could be the insulin signalling pathway, which is among the few signals that oocytes, at early stages of development, are responsive to, and while high levels of insulin can impair oocyte growth, low amounts can, on the contrary, promote oocyte growth . Interestingly, it is already shown that insulin can induce igf2 expression and repress ghr2 and igf1 expression in cichlids [136, 146]. It should be noted that insulin signalling is also involved in yolk formation of fish oocytes, through the regulation of vitellogenesis .
Adaptive trophic divergence may influence maternal mRNA deposition of growth- and stress-related genes in cichlids
Perhaps the most interesting finding of our study might be the clear clustering of species according to trophic niche, based on their maternal mRNA abundances of the target genes, when we combined the data together. For instance, the two carnivore species in each of LM and LV were clustered together, whereas the herbivore species of all 3 lakes were placed as separate clusters. Also, the three herbivore species in LT were clustered in one sub-branch, but the LT carnivore species was placed distantly, closer to the RV branch. Considering that distinct trophic niche clustering was observed using the maternal mRNA abundances of only 5 candidate genes, this could suggest potential predictability in differentiating herbivore and carnivore cichlid species based on their maternal mRNA investment. This points to disparate adaptive trajectories, which might influence subsequent developmental pathways in species specialized to a particular trophic niche. It should be noted that the haplochromine cichlids, with different trophic specializations, possess morphologically distinct skeletal structures features in dentition and jaw bones, which play an important role in their adaptive divergence [33, 148,149,150,151,152]. Our findings support the idea of the parallel evolution at the molecular level of trophic similar trophic niches and maternal mRNA investment across the parallel cichlid adaptive radiations . However, a more in-depth analysis of more genes and species is necessary to substantiate this claim. Knowing that the two signalling pathways, GC and IGF, are both involved in trophic skeletal formation during early development in fish , and particularly since recent studies in haplochromine cichlids and pupfish raise the possibility of IGF involvement in adaptive craniofacial divergence [153, 154], our study might implicate potential associations between maternal RNA deposition and differential early developmental patterning of diverse feeding apparatus in fishes. This is interesting because the pleiotropic involvement of both pathways in feeding and incubation of eggs could be key to controlling the trade-off between a morphological (jaws and feeding) and behavioural phenotype (mouthbrooding) in haplochromine cichlids .
In this study, we provide first insight into the molecular basis of an important life-history trait, namely by the identification of links between trophic specialization, habitats/lakes and maternal mRNA deposition of specific genes in East African cichlids. Our results show that maternal mRNA inputs vary substantially across cichlid adaptive radiations and with respect to trophic specializations, based on selected key genes involved in growth-related and environmental stress-mediating molecular pathways in association with egg size. We speculate that the evolutionary trade-off between mouthbrooding and feeding may be influenced by the same pleiotropic stress- and growth-related molecular pathways. However, further studies using large-scale transcriptional profiling with higher number of cichlid species with distinct feeding and breeding strategies are required to support this notion. Moreover, our findings shed light on the possible molecular interactions involved during this early and divergent stage of the hourglass hypothesis in cichlids and beyond, as the regulatory mechanisms causing the hourglass pattern still remain open to exploration [155, 156].
- ghr :
growth hormone receptor gene
- igf :
insulin-like growth factor gene
- gr :
glucocorticoid receptor gene
- sgk1 :
serum/glucocorticoid-regulated kinase 1 gene
insulin-like growth factor pathway
Trivers RL. Parental investment and sexual selection. In: Campbell, editor. Sexual selection and the descent of man. New York: Aldine de Gruyter; 1972.
Godfray HC. Evolutionary theory of parent-offspring conflict. Nature. 1995;376:133–8.
Hendry AP, Day T, Cooper AB. Optimal size and number of propagules: allowance for discrete stages and effects of maternal size on reproductive output and offspring fitness. Am Nat. 2001;157:387–407. https://doi.org/10.1086/319316.
Mousseau TA, Fox CW. The adaptive significance of maternal effects. Trends Ecol Evol. 1998;13:403–7.
Lindström J. Early development and fitness in birds and mammals. Trends Ecol Evol. 1999;14:343–8.
Sinervo B, Licht P. Proximate constraints on the evolution of egg size, number, and total clutch mass in lizards. Science. 1991;252:1300–2.
Roff D. The evolution of life histories. London: Chapman Hall; 1992.
Taborsky B. Mothers determine offspring size in response to own juvenile growth conditions. Biol Lett. 2006;2:225–8.
Bernardo J. Maternal effects in animal ecology. Am Zool. 1996;36:83–105. https://doi.org/10.1093/icb/36.2.83.
Tyler CR, Sumpter JP. Oocyte growth and development in teleosts. Rev Fish Biol Fish. 1996;6:287–318. https://doi.org/10.1007/BF00122584.
Benhaïm D, Skúlason S, Hansen BR. Behavioural variation in juvenile Arctic charr in relation to body size. J Fish Biol. 2003;62:1326–38.
Leblanc CA-L, Benhaïm D, Hansen BR, Kristjánsson BK, Skúlason S. The importance of egg size and social effects for behaviour of arctic charr juveniles: behaviour of charr juveniles. Ethology. 2011;117:664–74.
Leblanc CA-L, Kristjánsson BK, Skúlason S. The importance of egg size and egg energy density for early size patterns and performance of Arctic charr Salvelinus alpinus. Aquac Res. 2016;47:1100–11.
Bernardo J. The particular maternal effect of propagule size, especially egg size: patterns, models, quality of evidence and interpretations. Am Zool. 1996;36:216–36.
Mousseau TA, Fox CW, editors. Maternal effects as adaptations. New York: Oxford University Press; 1998.
Segers FHID, Taborsky B. Juvenile exposure to predator cues induces a larger egg size in fish. Proc Biol Sci. 2012;279:1241–8.
Einum S, Fleming IA. Highly fecund mothers sacrifice offspring survival to maximize fitness. Nature. 2000;405:565–7.
Räsänen K, Kruuk LEB. Maternal effects and evolution at ecological time-scales. Funct Ecol. 2007;21:408–21.
Kotrschal A, Heckel G, Bonfils D, Taborsky B. Life-stage specific environments in a cichlid fish: implications for inducible maternal effects. Evol Ecol. 2012;26:123–37. https://doi.org/10.1007/s10682-011-9495-5.
Forbes EL, Preston CD, Lokman PM. Zebrafish (Danio rerio) and the egg size versus egg number trade off: effects of ration size on fecundity are not mediated by orthologues of the Fec gene. Reprod Fertil Dev. 2010;22:1015–21.
Kamler E. Parent–egg–progeny relationships in teleost fishes: an energetics perspective. Rev Fish Biol Fish. 2005;15:399–421. https://doi.org/10.1007/s11160-006-0002-y.
McBride RS, Somarakis S, Fitzhugh GR, Albert A, Yaragina NA, Wuenschel MJ, et al. Energy acquisition and allocation to egg production in relation to fish reproductive strategies. Fish Fish. 2015;16:23–57. https://doi.org/10.1111/faf.12043.
Wang H-Y, Einhouse DW, Fielder DG, Rudstam LG, Vandergoot CS, VanDeValk AJ, et al. Maternal and stock effects on egg-size variation among walleye Sander vitreus stocks from the Great Lakes region. J Great Lakes Res. 2012;38:477–89.
Wootton RJ, Fletcher DA. Effect of spawning number and ration on reproductive performance of the batch-spawning three-spined stickleback Gasterosteus aculeatus. J Fish Biol. 2009;75:618–29.
Gregersen F, Haugen TO, Larsen ON. Egg size differentiation among sympatric demes of brown trout: possible effects of density-dependent interactions among fry. Ecol Freshw Fish. 2006;15:237–46. https://doi.org/10.1111/j.1600-0633.2006.00129.x.
Bownds C, Wilson R, Marshall DJ. Why do colder mothers produce larger eggs? An optimality approach. J Exp Biol. 2010;213:3796–801.
Gregersen F, Haugen TO, Vøllestad LA. Contemporary egg size divergence among sympatric grayling demes with common ancestors. Ecol Freshw Fish. 2008;17:110–8.
Finn RN, Fyhn HJ, Norberg B, Munholland J, Reith M. Oocyte hydration as a key feature in the adaptive evolution of teleost fishes to seawater. In: International symposium on reproductive physiology of fish; 2000.
Kucera CJ, Faulk CK, Holt GJ. The effect of spawning salinity on eggs of spotted seatrout (Cynoscion nebulosus, Cuvier) from two bays with historically different salinity regimes. J Exp Mar Biol Ecol. 2002;272:147–58.
Alderdice DF, Rao TR, Rosenthal H. Osmotic responses of eggs and larvae of the Pacific herring to salinity and cadmium. Helgoländer Wissenschaftliche Meeresuntersuchungen. 1979;32:508–38. https://doi.org/10.1007/BF02277992.
Giesing ER, Suski CD, Warner RE, Bell AM. Female sticklebacks transfer information via eggs: effects of maternal experience with predators on offspring. Proc Biol Sci. 2011;278:1753–9.
Schluter D. Ecological character displacement in adaptive radiation. Am Nat. 2000;156:S4–16. https://doi.org/10.1086/303412.
Brawand D, Wagner CE, Li YI, Malinsky M, Keller I, Fan S, et al. The genomic substrate for adaptive radiation in African cichlid fish. Nature. 2014;513:375–81.
Kocher TD. Adaptive evolution and explosive speciation: the cichlid fish model. Nat Rev Genet. 2004;5:288–98.
Irissari I, Singh P, Koblmüller S, Torres-Dowdall J, Henning F, Franchini P, et al. Anchored phylogenomics uncovers deep inter-tribal hybridizations in the Lake Tanganyika cichlid radiation and highlights adaptive loci shaping species’ ecology. Nat Commun. 2018;9:3159.
Mayer WE, Tichy H, Klein J. Phylogeny of African cichlid fishes as revealed by molecular markers. Heredity (Edinb). 1998;80:702–14. https://doi.org/10.1046/j.1365-2540.1998.00347.x.
Koblmüller S, Schliewen UK, Duftner N, Sefc KM, Katongo C, Sturmbauer C. Age and spread of the haplochromine cichlid fishes in Africa. Mol Phylogenet Evol. 2008;49:153–69.
Meier JI, Marques DA, Mwaiko S, Wagner CE, Excoffier L, Seehausen O. Ancient hybridization fuels rapid cichlid fish adaptive radiations. Nat Commun. 2017;8:14363. https://doi.org/10.1038/ncomms14363.
Goodwin NB, Balshine-Earn S, Reynolds JD. Evolutionary transitions in parental care in cichlid fish. Proc R Soc B Biol Sci. 1998;265:2265–72. https://doi.org/10.1098/rspb.1998.0569.
Duponchelle F, Paradis E, Ribbink AJ, Turner GF. Parallel life history evolution in mouthbrooding cichlids from the African Great Lakes. Proc Natl Acad Sci U S A. 2008;105:15475–80.
Parsons PJ, Bridle JR, Rüber L, Genner MJ. Evolutionary divergence in life history traits among populations of the Lake Malawi cichlid fish Astatotilapia calliptera. Ecol Evol. 2017;7:8488–506.
Fryer G, Iles T. The cichlid fishes of the Great Lakes of Africa: their biology and evolution. Edinburgh: Oliver Boyd; 1972. p. 656.
Tkint T, Verheyen E, De Kegel B, Helsen P, Adriaens D. Dealing with food and eggs in mouthbrooding cichlids: structural and functional trade-offs in fitness related traits. PLoS One. 2012;7:e31117.
Lubzens E, Bobe J, Young G, Sullivan CV. Maternal investment in fish oocytes and eggs: the molecular cargo and its contributions to fertility and early development. Aquaculture. 2017;472:107–43.
Baroux C, Autran D, Gillmor CS, Grimanelli D, Grossniklaus U. The maternal to zygotic transition in animals and plants. Cold Spring Harb Symp Quant Biol. 2008;73:89–100.
Harvey SA, Sealy I, Kettleborough R, Fenyes F, White R, Stemple D, et al. Identification of the zebrafish maternal and paternal transcriptomes. Development. 2013;140:2703–10.
Paranjpe SS, Jacobi UG, van Heeringen SJ, Veenstra GJC. A genome-wide survey of maternal and embryonic transcripts during Xenopus tropicalis development. BMC Genom. 2013;14:762.
Barckmann B, Simonelig M. Control of maternal mRNA stability in germ cells and early embryos. Biochim Biophys Acta Gene Regul Mech. 2013;1829:714–24.
Walser CB, Lipshitz HD. Transcript clearance during the maternal-to-zygotic transition. Curr Opin Genet Dev. 2011;21:431–43.
Lee MT, Bonneau AR, Takacs CM, Bazzini AA, DiVito KR, Fleming ES, et al. Nanog, Pou5f1 and SoxB1 activate zygotic gene expression during the maternal-to-zygotic transition. Nature. 2013;503:360–4.
Wagner DS, Dosch R, Mintzer KA, Wiemelt AP, Mullins MC. Maternal control of development at the midblastula transition and beyond. Dev Cell. 2004;6:781–90.
Mathavan S, Lee SGP, Mak A, Miller LD, Murthy KRK, Govindarajan KR, et al. Transcriptome analysis of zebrafish embryogenesis using microarrays. PLoS Genet. 2005;1:260–76.
Pikulkaew S, Benato F, Celeghin A, Zucal C, Skobo T, Colombo L, et al. The knockdown of maternal glucocorticoid receptor mRNA alters embryo development in zebrafish. Dev Dyn. 2011;240:874–89. https://doi.org/10.1002/dvdy.22586.
Alexa K, Choe S-K, Hirsch N, Etheridge L, Laver E, Sagerström CG. Maternal and zygotic aldh1a2 activity is required for pancreas development in zebrafish. PLoS One. 2009;4:e8261. https://doi.org/10.1371/journal.pone.0008261.
Langdon YG, Mullins MC. Maternal and zygotic control of zebrafish dorsoventral axial patterning. Annu Rev Genet. 2011;45:357–77.
Celeghin A, Benato F, Pikulkaew S, Rabbane MG, Colombo L, Dalla Valle L. The knockdown of the maternal estrogen receptor 2a (esr2a) mRNA affects embryo transcript contents and larval development in zebrafish. Gen Comp Endocrinol. 2011;172:120–9.
Romney AL, Podrabsky JE. Transcriptomic analysis of maternally provisioned cues for phenotypic plasticity in the annual killifish, Austrofundulus limnaeus. Evodevo. 2017;8:6.
Ma L, Strickler AG, Parkhurst A, Yoshizawa M, Shi J, Jeffery WR. Maternal genetic effects in Astyanax cavefish development. Dev Biol. 2018;441:209–20.
Adrian-Kalchhauser I, Walser J-C, Schwaiger M, Burkhardt-Holm P. RNA sequencing of early round goby embryos reveals that maternal experiences can shape the maternal RNA contribution in a wild vertebrate. BMC Evol Biol. 2018;18:34. https://doi.org/10.1186/s12862-018-1132-2.
Colson V, Cousture M, Zanerato-Damasceno D, Valotaire C, Nguyen T, Le Cam A, et al. Maternal temperature exposure triggers emotional and cognitive disorders and dysregulation of neurodevelopment genes in fish. PeerJ Prepr. 2018;6:e26910v1.
Shama LNS, Mark FC, Strobel A, Lokmer A, John U, Mathias Wegner K. Transgenerational effects persist down the maternal line in marine sticklebacks: gene expression matches physiology in a warming ocean. Evol Appl. 2016;9:1096–111. https://doi.org/10.1111/eva.12370.
Bougas B, Audet C, Bernatchez L. The influence of parental effects on transcriptomic landscape during early development in brook charr (Salvelinus fontinalis, Mitchill). Heredity (Edinb). 2013;110:484–91.
Amaral IPG, Johnston IA. Experimental selection for body size at age modifies early life-history traits and muscle gene expression in adult zebrafish. J Exp Biol. 2012;215:3895–904.
Elinson RP, Sabo MC, Fisher C, Yamaguchi T, Orii H, Nath K. Germ plasm in Eleutherodactylus coqui, a direct developing frog with large eggs. Evodevo. 2011;2:20.
Collazo A, Keller R. Early development of Ensatina eschscholtzii: an amphibian with a large, yolky egg. Evodevo. 2010;1:6.
Luckenbach JA, Kusakabe M, Swanson P, Young G. Unilateral ovariectomy increases egg size and reduces follicular atresia in the semelparous coho salmon, Oncorhynchus kisutch. J Exp Zool A Ecol Genet Physiol. 2008;309:468–76.
Segers FHID, Berishvili G, Taborsky B. Egg size-dependent expression of growth hormone receptor accompanies compensatory growth in fish. Proc Biol Sci. 2012;279:592–600.
Pepling ME, de Cuevas M, Spradling AC. Germline cysts: a conserved phase of germ cell development? Trends Cell Biol. 1999;9:257–62.
Kumano G. Polarizing animal cells via mRNA localization in oogenesis and early development. Dev Growth Differ. 2012;54:1–18.
Matova N, Cooley L. Comparative aspects of animal oogenesis. Dev Biol. 2001;231:291–320.
Song JL, Wessel GM. How to make an egg: transcriptional regulation in oocytes. Differentiation. 2005;73:1–17.
Evsikov AV, Graber JH, Brockman JM, Hampl A, Holbrook AE, Singh P, et al. Cracking the egg: molecular dynamics and evolutionary aspects of the transition from the fully grown oocyte to embryo. Genes Dev. 2006;20:2713–27.
Hayward LS, Wingfield JC. Maternal corticosterone is transferred to avian yolk and may alter offspring growth and adult phenotype. Gen Comp Endocrinol. 2004;135:365–71.
Leatherland JF, Li M, Barkataki S. Stressors, glucocorticoids and ovarian function in teleosts. J Fish Biol. 2010;76:86–111.
LaLone CA, Villeneuve DL, Olmstead AW, Medlock EK, Kahl MD, Jensen KM, et al. Effects of a glucocorticoid receptor agonist, dexamethasone, on fathead minnow reproduction, growth, and development. Environ Toxicol Chem. 2012;31:611–22.
Eriksen MS, Bakken M, Espmark A, Braastad BO, Salte R. Prespawning stress in farmed Atlantic salmon Salmo salar: maternal cortisol exposure and hyperthermia during embryonic development affect offspring survival, growth and incidence of malformations. J Fish Biol. 2006;69:114–29. https://doi.org/10.1111/j.1095-8649.2006.01071.x.
Ma X, Liu X, Zhang Y, Zhu P, Ye W, Lin H. Two growth hormone receptors in Nile tilapia (Oreochromis niloticus): molecular characterization, tissue distribution and expression profiles in the gonad during the reproductive cycle. Comp Biochem Physiol Part B Biochem Mol Biol. 2007;147:325–39.
Mommsen TP. Paradigms of growth in fish. Comp Biochem Physiol B Biochem Mol Biol. 2001;129:207–19.
Canosa LF, Chang JP, Peter RE. Neuroendocrine control of growth hormone in fish. Gen Comp Endocrinol. 2007;151:1–26.
Besseau L, Fuentès M, Sauzet S, Beauchaud M, Chatain B, Covès D, et al. Somatotropic axis genes are expressed before pituitary onset during zebrafish and sea bass development. Gen Comp Endocrinol. 2013;194:133–41.
Sullivan CV, Chapman RW, Reading BJ, Anderson PE. Transcriptomics of mRNA and egg quality in farmed fish: some recent developments and future directions. Gen Comp Endocrinol. 2015;221:23–30.
Fujimura K, Okada N. Development of the embryo, larva and early juvenile of Nile tilapia Oreochromis niloticus (Pisces: Cichlidae). Developmental staging system. Dev Growth Differ. 2007;49:301–24. https://doi.org/10.1111/j.1440-169X.2007.00926.x.
Kratochwil CF, Sefton MM, Meyer A. Embryonic and larval development in the Midas cichlid fish species flock (Amphilophus spp.): a new evo–devo model for the investigation of adaptive novelties and species differences. BMC Dev Biol. 2015;15:12.
Aanes H, Winata CL, Lin CH, Chen JP, Srinivasan KG, Lee SGP, et al. Zebrafish mRNA sequencing deciphers novelties in transcriptome dynamics during maternal to zygotic transition. Genome Res. 2011;21:1328–38.
Ahi EP, Sefc KM. A gene expression study of dorso-ventrally restricted pigment pattern in adult fins of Neolamprologus meeli, an African cichlid species. PeerJ. 2017;5:e2843.
Ahi EP, Richter F, Sefc KM. A gene expression study of ornamental fin shape in Neolamprologus brichardi, an African cichlid species. Sci Rep. 2017;7:17398.
Ahi EP, Sefc KM. Anterior-posterior gene expression differences in three Lake Malawi cichlid fishes with variation in body stripe orientation. PeerJ. 2017;5:e4080.
Yang CG, Wang XL, Tian J, Liu W, Wu F, Jiang M, et al. Evaluation of reference genes for quantitative real-time RT-PCR analysis of gene expression in Nile tilapia (Oreochromis niloticus). Gene. 2013;527:183–92.
Singh P, Börger C, More H, Sturmbauer C. The role of alternative splicing and differential gene expression in cichlid adaptive radiation. Genome Biol Evol. 2017;9:2764–81.
Zerbino DR, Achuthan P, Akanni W, Amode MR, Barrell D, Bhai J, et al. Ensembl 2018. Nucleic Acids Res. 2018;46:D754–61.
Fleige S, Pfaffl MW. RNA integrity and the effect on the real-time qRT-PCR performance. Mol Asp Med. 2006;27:126–39.
Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J. qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 2007;8:R19.
Ramakers C, Ruijter JM, Deprez RHL, Moorman AFM. Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neurosci Lett. 2003;339:62–6.
Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper–Excel-based tool using pair-wise correlations. Biotechnol Lett. 2004;26:509–15.
Andersen CL, Jensen JL, Ørntoft TF. Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004;64:5245–50.
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3:RESEARCH0034.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29:e45.
Bergkvist A, Rusnakova V, Sindelka R, Garda JMA, Sjögreen B, Lindh D, et al. Gene expression profiling—clusters of possibilities. Methods. 2010;50:323–35.
Kubista M, Andrade JM, Bengtsson M, Forootan A, Jonák J, Lind K, et al. The real-time polymerase chain reaction. Mol Aspects Med. 2006;27:95–125.
Woltering JM, Holzem M, Schneider RF, Nanos V, Meyer A. The skeletal ontogeny of Astatotilapia burtoni—a direct-developing model system for the evolution and development of the teleost body plan. BMC Dev Biol. 2018;18:8. https://doi.org/10.1186/s12861-018-0166-4.
da Mattos D, Cardoso LD, Fosse PJ, Radael MC, Filho JCF, de Manhães JVA, et al. Description of the ontogenic and larval period of discus fish (Symphysodon aequifasciatus). Zygote. 2015;23:460–6.
Duboule D. Temporal colinearity and the phylotypic progression: a basis for the stability of a vertebrate Bauplan and the evolution of morphologies through heterochrony. Dev Suppl. 1994;120:135–42.
Timms RT, Tchasovnikarova IA, Antrobus R, Dougan G, Lehner PJ. ATF7IP-mediated stabilization of the histone methyltransferase SETDB1 is essential for heterochromatin formation by the HUSH complex. Cell Rep. 2016;17:653–9.
Minkovsky A, Sahakyan A, Rankin-Gee E, Bonora G, Patel S, Plath K. The Mbd1-Atf7ip-Setdb1 pathway contributes to the maintenance of X chromosome inactivation. Epigenetics Chromatin. 2014;7:12. https://doi.org/10.1186/1756-8935-7-12.
McGraw S, Vigneault C, Sirard M-A. Temporal expression of factors involved in chromatin remodeling and in gene regulation during early bovine in vitro embryo development. Reproduction. 2007;133:597–608.
Marlow FL. Maternal control of development in vertebrates. Matern. Control Dev. Vertebr. My Mother Made Me Do It!. San Rafael: Morgan & Claypool Life Sciences; 2010.
Dosch R, Wagner DS, Mintzer KA, Runke G, Wiemelt AP, Mullins MC. Maternal control of vertebrate development before the midblastula transition: mutants from the zebrafish I. Dev Cell. 2004;6:771–80.
Berti C, Fontanella B, Ferrentino R, Meroni G. Mig12, a novel Opitz syndrome gene product partner, is expressed in the embryonic ventral midline and co-operates with Mid1 to bundle and stabilize microtubules. BMC Cell Biol. 2004;5:9. https://doi.org/10.1186/1471-2121-5-9.
Kim C-W, Moon Y-A, Park SW, Cheng D, Kwon HJ, Horton JD. Induced polymerization of mammalian acetyl-CoA carboxylase by MIG12 provides a tertiary level of regulation of fatty acid synthesis. Proc Natl Acad Sci U S A. 2010;107:9626–31.
Eno C, Solanki B, Pelegri F. Aura (mid1ip1l) regulates the cytoskeleton at the zebrafish egg-to-embryo transition. Development. 2016;143:1585–99.
Ahi EP, Guðbrandsson J, Kapralova KH, Franzdóttir SR, Snorrason SS, Maier VH, et al. Validation of reference genes for expression studies during craniofacial development in arctic charr. PLoS One. 2013;8:e66389.
Love OP, McGowan PO, Sheriff MJ. Maternal adversity and ecological stressors in natural populations: the role of stress axis programming in individuals, with implications for populations and communities. Funct Ecol. 2013;27:81–92. https://doi.org/10.1111/j.1365-2435.2012.02040.x.
Brann D, Mahesh V. Role of corticosteroids in female reproduction. FASEB J. 1991;5:2691–8.
Hillier S, Tetsuka M. An anti-inflammatory role for glucocorticoids in the ovaries? J Reprod Immunol. 1998;39:21–7.
Contreras-Sánchez WM, Schreck CB, Fitzpatrick MS, Pereira CB. Effects of stress on the reproductive performance of rainbow trout (Oncorhynchus mykiss)1. Biol Reprod. 1998;58:439–47. https://doi.org/10.1095/biolreprod58.2.439.
Love OP, Williams TD. The adaptive value of stress-induced phenotypes: effects of maternally derived corticosterone on sex-biased investment, cost of reproduction, and maternal fitness. Am Nat. 2008;172:E135–49. https://doi.org/10.1086/590959.
Monclús R, Tiulim J, Blumstein DT. Older mothers follow conservative strategies under predator pressure: the adaptive role of maternal glucocorticoids in yellow-bellied marmots. Horm Behav. 2011;60:660–5.
Firestone G, Giampaolo J, O’Keeffe B. Stimulus-dependent regulation of serum and glucocorticoid inducible protein kinase (SGK) transcription, subcellular localization and enzymatic activity. Cell Physiol Biochem. 2003;13:1–12.
Maiyar AC, Phu PT, Huang AJ, Firestone GL. Repression of glucocorticoid receptor transactivation and dna binding of a glucocorticoid response element within the serum/glucocorticoid-inducible protein kinase ( sgk ) gene promoter by the p53 tumor suppressor protein. Mol Endocrinol. 1997;11:312–29. https://doi.org/10.1210/mend.11.3.9893.
Galliher-Beckley AJ, Williams JG, Cidlowski JA. Ligand-independent phosphorylation of the glucocorticoid receptor integrates cellular stress pathways with nuclear receptor signaling. Mol Cell Biol. 2011;31:4663–75.
Lu M, Wang J, Jones KT, Ives HE, Feldman ME, Yao L, et al. mTOR complex-2 activates ENaC by phosphorylating SGK1. J Am Soc Nephrol. 2010;21:811–8.
Hong F, Larrea MD, Doughty C, Kwiatkowski DJ, Squillace R, Slingerland JM. mTOR-raptor binds and activates SGK1 to regulate p27 phosphorylation. Mol Cell. 2008;30:701–11.
Sánchez F, Smitz J. Molecular control of oogenesis. Biochim Biophys Acta. 2012;1822:1896–912.
Luca F, Kashyap S, Southard C, Zou M, Witonsky D, Di Rienzo A, et al. Adaptive variation regulates the expression of the human SGK1 gene in response to stress. PLoS Genet. 2009;5:e1000489. https://doi.org/10.1371/journal.pgen.1000489.
Anacker C, Cattaneo A, Musaelyan K, Zunszain PA, Horowitz M, Molteni R, et al. Role for the kinase SGK1 in stress, depression, and glucocorticoid effects on hippocampal neurogenesis. Proc Natl Acad Sci U S A. 2013;110:8708–13.
Evans TG. Co-ordination of osmotic stress responses through osmosensing and signal transduction events in fishes. J Fish Biol. 2010;76:1903–25.
Jones KT, Greer ER, Pearce D. Ashrafi K. Rictor/TORC2 regulates Caenorhabditis elegans fat storage, body size, and development through sgk-1. PLoS Biol. 2009;7:e60. https://doi.org/10.1371/journal.pbio.1000060.
Hertweck M, Göbel C, Baumeister RC. Elegans SGK-1 Is the critical component in the Akt/PKB kinase complex to control stress response and life span. Dev Cell. 2004;6:577–88.
Very NM, Kittilson JD, Norbeck LA, Sheridan MA. Isolation, characterization, and distribution of two cDNAs encoding for growth hormone receptor in rainbow trout (Oncorhynchus mykiss). Comp Biochem Physiol Part B Biochem Mol Biol. 2005;140:615–28.
Fukamachi S, Meyer A. Evolution of receptors for growth hormone and somatolactin in fish and land vertebrates: lessons from the lungfish and sturgeon orthologues. J Mol Evol. 2007;65:359–72. https://doi.org/10.1007/s00239-007-9035-7.
Jiao B, Huang X, Chan CB, Zhang L, Wang D, Cheng CHK. The co-existence of two growth hormone receptors in teleost fish and their differential signal transduction, tissue distribution and hormonal regulation of expression in seabream. J Mol Endocrinol. 2006;36:23–40.
Saera-Vila A, Calduch-Giner J-A, Pérez-Sánchez J. Duplication of growth hormone receptor (GHR) in fish genome: gene organization and transcriptional regulation of GHR type I and II in gilthead sea bream (Sparus aurata). Gen Comp Endocrinol. 2005;142:193–203.
Reindl KM, Kittilson JD, Sheridan MA. Differential ligand binding and agonist-induced regulation characteristics of the two rainbow trout GH receptors, Ghr1 and Ghr2, in transfected cells. J Endocrinol. 2009;202:463–71.
Kittilson JD, Jones E, Sheridan MA. ERK, Akt, and STAT5 are differentially activated by the two growth hormone receptor subtypes of a teleost fish (Oncorhynchus Mykiss). Front Endocrinol (Lausanne). 2011;2:30. https://doi.org/10.3389/fendo.2011.00030/abstract.
Gao F-Y, Lu M-X, Ye X, Huang Z-H, Wang H, Zhu H-P, et al. Identification and expression analysis of two growth hormone receptors in zanzibar tilapia (Oreochromis hornorum). Fish Physiol Biochem. 2011;37:553–65. https://doi.org/10.1007/s10695-010-9457-9.
Pierce AL, Breves JP, Moriyama S, Uchida K, Grau EG. Regulation of growth hormone (GH) receptor (GHR1 and GHR2) mRNA level by GH and metabolic hormones in primary cultured tilapia hepatocytes. Gen Comp Endocrinol. 2012;179:22–9.
Reinecke M. Insulin-like growth factors and fish reproduction. Biol Reprod. 2010;82:656–61.
Srivastava RK, Van Der Kraak G. Regulation of DNA synthesis in goldfish vitellogenic ovarian follicles by hormones and growth factors. J Exp Zool. 1994;270:263–72. https://doi.org/10.1002/jez.1402700305.
Kagawa H, Kobayashi M, Hasegawa Y, Aida K. Insulin and insulin-like growth factors I and II induce final maturation of oocytes of red seabream, pagrus major, in vitro. Gen Comp Endocrinol. 1994;95:293–300.
Negatu Z, Hsiao SM, Wallace RA. Effects of insulin-like growth factor-I on final oocyte maturation and steroid production in Fundulus heteroclitus. Fish Physiol Biochem. 1998;19:13–21. https://doi.org/10.1023/A:1007777926238.
Lokman PM, George KAN, Divers SL, Algie M, Young G. 11-Ketotestosterone and IGF-I increase the size of previtellogenic oocytes from shortfinned eel, Anguilla australis, in vitro. Reproduction. 2007;133:955–67.
Luckenbach JA, Dickey JT, Swanson P. Follicle-stimulating hormone regulation of ovarian transcripts for steroidogenesis-related proteins and cell survival, growth and differentiation factors in vitro during early secondary oocyte growth in coho salmon. Gen Comp Endocrinol. 2011;171:52–63.
Beck SV, Räsanen K, Ahi EP, Kristjánsson BK, Skúlason S, Jónsson ZO, et al. Gene expression in the phenotypically plastic Arctic charr (Salvelinus alpinus): a focus on growth and ossification at early stages of development. Evol Dev. 2018;1:1-15. https://doi.org/10.1111/ede.12275.
Liu J, Koenigsfeld AT, Cantley TC, Boyd CK, Kobayashi Y, Lucy MC. Growth and the initiation of steroidogenesis in porcine follicles are associated with unique patterns of gene expression for individual components of the ovarian insulin-like growth factor system1. Biol Reprod. 2000;63:942–52. https://doi.org/10.1095/biolreprod63.3.942.
Das D, Arur S. Conserved insulin signaling in the regulation of oocyte growth, development, and maturation. Mol Reprod Dev. 2017;84:444–59.
Pierce AL, Breves JP, Moriyama S, Hirano T, Grau EG. Differential regulation of Igf1 and Igf2 mRNA levels in tilapia hepatocytes: effects of insulin and cortisol on GH sensitivity. J Endocrinol. 2011;211:201–10.
Hiramatsu N, Todo T, Sullivan CV, Schilling J, Reading BJ, Matsubara T, et al. Ovarian yolk formation in fishes: molecular mechanisms underlying formation of lipid droplets and vitellogenin-derived yolk proteins. Gen Comp Endocrinol. 2015;221:9–15.
Albertson RC, Markert JA, Danley PD, Kocher TD. Phylogeny of a rapidly evolving clade: the cichlid fishes of Lake Malawi, East Africa. Proc Natl Acad Sci U S A. 1999;96:5107–10.
Muschick M, Indermaur A, Salzburger W. Convergent evolution within an adaptive radiation of cichlid fishes. Curr Biol. 2012;22:2362–8.
Albertson RC, Kocher TD. Genetic and developmental basis of cichlid trophic diversity. Heredity (Edinb). 2006;97:211–21.
Powder KE, Albertson RC. Cichlid fishes as a model to understand normal and clinical craniofacial variation. Dev Biol. 2016;415:338–46.
Ahi EP. Signalling pathways in trophic skeletal development and morphogenesis: insights from studies on teleost fish. Dev Biol. 2016;420:11–31.
Concannon MR, Albertson RC. The genetic and developmental basis of an exaggerated craniofacial trait in East African cichlids. J Exp Zool Part B Mol Dev Evol. 2015;324:662–70. https://doi.org/10.1002/jez.b.22641.
Lencer ES, Warren WC, Harrison R, McCune AR. The Cyprinodon variegatus genome reveals gene expression changes underlying differences in skull morphology among closely related species. BMC Genom. 2017;18:424. https://doi.org/10.1186/s12864-017-3810-7.
Irie N, Kuratani S. The developmental hourglass model: a predictor of the basic body plan? Development. 2014;141:4649–55.
Akhshabi S, Sarda S, Dovrolis C, Yi S. An explanatory evo–devo model for the developmental hourglass. 2014;F1000Research(3):156.
EPA, PS, LAL, WG and CS designed the study. EPA conducted the laboratory experiment, measurements and figure preparations. EPA, PS and LAL analysed the data, and EPA, PS, LAL and CS wrote the manuscript. WG performed fish breeding, egg sampling and photography of adult fishes shown in Fig. 1a. All authors reviewed the manuscript and approved its content. All authors read and approved the final manuscript.
The authors thank Holger Zimmermann and Stephan Koblmüller for sharing their precious knowledge on East African cichlid fishes. We would also like to extend special thanks to Isis Kowarik for her useful advices on statistical issues. The authors acknowledge Institute of Biology at University of Graz for providing fish breeding and laboratory facilities and the Austrian Science Fund for the financial support of our study.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analysed during this study are included in this published article.
Consent for publication
Ethics approval and consent to participate
Studies of fish eggs do not require ethics approval or consent to participate.
This study was funded by the Austrian Science Fund (Grant P22737).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.