Skip to main content

Revelation of candidate genes and molecular mechanism of reproductive seasonality in female rohu (Labeo rohita Ham.) by RNA sequencing

Abstract

Background

Carp fish, rohu (Labeo rohita Ham.) is important freshwater aquaculture species of South-East Asia having seasonal reproductive rhythm. There is no holistic study at transcriptome level revealing key candidate genes involved in such circannual rhythm regulated by biological clock genes (BCGs). Seasonality manifestation has two contrasting phases of reproduction, i.e., post-spawning resting and initiation of gonadal activity appropriate for revealing the associated candidate genes. It can be deciphered by RNA sequencing of tissues involved in BPGL (Brain-Pituitary-Gonad-Liver) axis controlling seasonality. How far such BCGs of this fish are evolutionarily conserved across different phyla is unknown. Such study can be of further use to enhance fish productivity as seasonality restricts seed production beyond monsoon season.

Result

A total of ~ 150 Gb of transcriptomic data of four tissues viz., BPGL were generated using Illumina TruSeq. De-novo assembled BPGL tissues revealed 75,554 differentially expressed transcripts, 115,534 SSRs, 65,584 SNPs, 514 pathways, 5379 transcription factors, 187 mature miRNA which regulates candidate genes represented by 1576 differentially expressed transcripts are available in the form of web-genomic resources. Findings were validated by qPCR. This is the first report in carp fish having 32 BCGs, found widely conserved in fish, amphibian, reptile, birds, prototheria, marsupials and placental mammals. This is due to universal mechanism of rhythmicity in response to environment and earth rotation having adaptive and reproductive significance.

Conclusion

This study elucidates evolutionary conserved mechanism of photo-periodism sensing, neuroendocrine secretion, metabolism and yolk synthesis in liver, gonadal maturation, muscular growth with sensory and auditory perception in this fish. Study reveals fish as a good model for research on biological clock besides its relevance in reproductive efficiency enhancement.

Peer Review reports

Background

Reproductive rhythm of fish varies from prolific breeders to seasonal breeders. Though prolific breeder like common carp and zebra fish breed independent of season but vast majority of the freshwater fishes of Indian subcontinent breed seasonally during the monsoon season (June–August) when rainfall is heaviest [1]. This rhythm reflects evolution of differential biological clock.

In seasonal reproductive behaviour, BPG/ BPGL (Brain-Pituitary-Gonad-Liver) axis and fish biological clock is well known having entrainment of its biological clock with environment, temperature and food availability [2, 3]. Most active liver metabolism has been found in fish having potential to acquire bigger size thus liver transcriptome should be included in study of gonadal maturation by BPG (Brain-Pituitary-Gonad) axis [4]. Biological rhythm of seasonal breeding in fish is regulated by neuro-endocrine system, controlled by BPGL axis [5]. In order to decipher candidate genes regulating seasonal fish breeding, RNA sequencing of associated tissues are imperative. Such data can be used to construct gene regulatory networks (GRNs) to depict process of cell differentiation, metabolism, cell cycle and signal transduction in the biological system. Understanding dynamics of these networks can shed light on mechanism of reproductive seasonality.

Since circannual biological rhythm of seasonal reproduction is controlled by biological clock, regulating BPGL axis, thus RNA seq approach can also be used to discover biological clock genes. Biological clock genes originated 2.5 billion years ago and they are structurally conserved across phyla i.e., from lower vertebrates to higher vertebrates and can be used for their in silico identification. Differential expression of such genes across different tissues can be used to study control of reproductive seasonality and adaptation during species divergence [6]. Major carp of South-East Asia like Labeo rohita (Hamilton), which is an important freshwater aquaculture species having economic, ecological and cultural attributes would be a best model to study how seasonal manifestation of reproductive behavior is controlled by circannual biological clock.

Earlier attempts to investigate gene expression associated with seasonal breeding in carp fish, rohu were confined to limited genes (< 200) by microarray approach, limited ESTs (< 5000) comparing preparatory and post-spawning stages along with limited genomic resources in form of SSR (< 250) but without SNP [7]. In another such investigation by RNA Seq approach [8] limited reproductive genes were identified. This study enlists reproductive genes which are based on pooled RNA samples (of twelve tissues related to digestive and reproductive system) without differentially expressed genes (DEGs), thus lacks in identification of key candidate genes associated with reproductive seasonality. Even in one of the largest global consortium model fish transcriptome projects, FishT1K (https://www.fisht1k.org/) which is covering 124 species representing 46 orders, having six cypriniformes species but there is no data on carp fish, Labeo rohita which can be used as genomic resource to discover seasonality controlling genes along with biological clock genes (BCGs). RNA seq data of BPGL axis involved tissues can be used for evolutionary studies. A tissue-wise RNA sequencing comparison between two physiologically contrasting phases i.e., post-spawning resting (PSR) and initiation of gonadal activity (IGA) can decipher the entire mechanism of reproductive seasonality. Such transcriptomic studies can reveal key candidate genes and their regulatory networks controlling neuroendocrine and physiological pathway operating in four different tissues, namely, BPGL.

The present work aims at revealing the candidate genes of seasonal biological rhythm of breeding and gene regulatory networks. It further aims to investigate differentially expressed biological clock genes (BCGs), along with their evolutionary conservation across jawless fish, cartilaginous fish (whale shark), bony fishes (zebrafish and European carp), amphibia, bird, reptile, monotreme (egg laying mammal), marsupials (young born) and placental mammals by comparative genomics. This study also aims to develop a genomic resource database of conserved BCGs of each phylum along with seasonality associated candidate genes in female carp fish, Labeo rohita (rohu).

Results and discussion

Reproductive phase identification by GSI estimation with ovarian morphology and histology

Ovarian morphology: The ovaries in December (Fig. 1A) appeared as two unequal elongated, cylindrical lobes or thread like structures lying on each side of the air bladder and joining together posteriorly with a common opening as a pore above the anus. The total weight of the two lobes varies between 2 to 10 g according to body weight and age. The color of the lobes brownish to light pink, opaque with very less vascularization in the tissue. By February end (Fig. 1B), the ovaries grown bigger in size with flattened lobes at the anterior end, yellowish brown in color, increased vascularization and little more transparent at the surfaces. The weight increases to 7 to 25 g. GSI ranges between 0.43 to 1.04 (average 0.78) in December as against 0.76 to 1.39 (average 1.14) in February during initiating phase (Table 1).

Fig. 1
figure1

Morphological and histological view of rohu ovary in December and February (1A. Rohu ovary in December; 1B. Rohu ovary in February; 1C. Histological view of rohu ovary in December at 10X; 1D. Histological view of rohu ovary in February at 10X)

Table 1 Average Gonadosomatic index, Oocyte diameter and % immature and mature oocytes observed in the ovarian histology of rohu during resting and initiating phases

Histology of ovaries

Resting phase (December): ovary is predominantly occupied (~ 98%) by immature oocytes of both type I and II and they constitute 80 and 20% respectively, of the total immature oocyte count. Both type oocytes are having larger nucleus. However, as the size of type II cells grow larger than type I, the nucleus become prominent with a greater number of nucleoli arranged along the inner margin of the nuclear membrane. The average diameter of type I oocytes is 75.77 μm. Some type-II oocytes are also seen which are larger in size (average diameter 140.59 μm) but less in number. In this phase neither any maturing oocytes (type-III) nor any mature oocytes are observed. Immature type I oocytes are seen towards the periphery whereas large type II cells are located towards the center (Fig. 1C).

Initiating phase (February): Histologically, ovaries show more type-II oocytes (60%) than type I immature oocytes (40%) and both of which constitute 98.58% of the total oocyte count. The type II cells are at perinucleolar stage. The average diameter is 89.19 μm and 151.77 μm of type I and type II oocytes, respectively. Rarely, type-III maturing oocytes may be visible (1.42%) with yolk vesicles appear along the periphery in the ooplasm but no matured follicles could be observed in this phase also (Fig. 1D).

Histological analysis of ovarian tissue of PSR and IGA phases revealed the changes in response to seasonality. Primary ovarian follicles (POF), smaller in size were present during PSR dormancy (lower GSI). During IGA phase having higher GSI maturation of secondary ovarian follicles (SOF) with larger size was seen in the histological photomicrographs, revealing conspicuous seasonal differences (Fig. 1A, B, C, D). It confirmed the accuracy of tissues collected in the investigation. Similar histological attributes of ovarian development stages have also been reported in other fishes like zebrafish [9].

Like other Indian teleosts, the reproductive cycle of L. rohita may be divided into four stages, viz., preparatory period (February–April), pre-spawning period (May–June), spawning period (July–August), and post-spawning period (September–January) and at each stage gonads show discrete changes [10]. Generally, in an adult rohu female, the ovary remains quiescent and almost inactive during last part of post spawning period also called post spawning resting (PSR). With gradual increase in day length and temperature, there is increase of gonadal activity (IGA) for preparatory period. This confirmed that these two reproductive contrasting phases, viz., PSR and IGA as most appropriate transcriptomic data points to study the changes that happens at the molecular level during onset of ovary maturation in this tropical species manifesting biological rhythm of seasonal breeding.

Pre-processing and assembly

A total of ~ 150 Gb of transcriptomic data of four tissues viz., BPGL were pre-processed before assembly by removing 1,098,340 poor quality reads. The remaining 305,461,765 high quality reads were subjected to de novo assembly to obtain DEGs. A total of 440,665 transcripts were generated with k-mer size 25. The minimum contig length was 201 bp and maximum of 19,331 bp with average length of 867 bp. N50 was found to be 1850 bp and total GC content as 42.94% (Table 2). Figure 2 shows the sequence length distribution of the data under study. A total of 45,019 unigenes were involved in the formation of isoforms. Numbers of isoforms were found in single unigenes ranging from 2 to 30.

Table 2 De novo assembled statistics of Rohu transcriptome (All IGA vs. all PSR)
Fig. 2
figure2

Sequence length distribution of the data

Identification of candidate genes

A total of 35,840, 21,767, 31,097 and 14,558 DEGs were identified for IGA [BR] vs. PSR [BR], IGA [PIT] vs. PSR [PIT], IGA [OVA] vs. PSR [OVA] and IGA [LIV] vs. PSR [LIV] respectively. The upregulated and downregulated DEGs from the datasets are presented in Table 3. Expression profiling and hierarchical clustering of DEGs of all tissue specific comparison were shown in the form of heatmaps. MA plots and volcano plots of differentially expressed genes were generated by edgeR tool which showed the red color dots are differentially expressed and black color dot represents the non-DEGs and all these plots were generated on the basis of log2 fold change values (Supplementary file 1).

Table 3 Number of upregulated and downregulated DEGs from IGA [BR] vs. PSR [BR], IGA [PIT] vs. PSR [PIT], IGA [OVA] vs. PSR [OVA] and IGA [LIV] vs. PSR [LIV]

Venn diagram was constructed to identify the unique and common unigenes among all tissue specific candidate genes (Fig. 3). A total of 705 unigenes were found common in all datasets of tissue specific DEGs. We found 20,845, 9388, 17,739 and 5758 unigenes as unique in IGA [BR] vs. PSR [BR], IGA [PIT] vs. PSR [PIT], IGA [OVA] vs. PSR [OVA] and IGA [LIV] vs. PSR [LIV], respectively.

Fig. 3
figure3

Venn diagram to identify the unique and common unigenes among all tissue specific differentially expressed genes

Annotations and gene ontology classification of DEGs

We found 79.27, 84.38, 82.87 and 86.44% of differentially expressed genes of IGA [BR] vs. PSR [BR], IGA [PIT] vs. PSR [PIT], IGA [OVA] vs. PSR [OVA] and IGA [LIV] vs. PSR [LIV] having blast match, respectively (Supplementary file 2: Uploaded at “Supplements” tab of http://webtom.cabgrid.res.in/lrsatdb/). Table 4 described the number of sequences which showed the blast hits, mapping, interproscan and annotations statistics for dataset under study.

Table 4 Sequence similarity to know sequence, annotation and mapping unigenes

DEGs from all tissue specific comparison were categorized into three sub division i.e. biological process, molecular function and cellular component (Supplementary file 3). Gene ontology and KEGG (Kyoto Encyclopaedia of Genes and Genomes) analysis were performed by using Blast2Go pro software v 4.0 [11]. Tissue specific DEGs were mapped to KEGG pathways and 131, 126, 130 and 124 pathways were found in IGA [BR] vs. PSR [BR], IGA [PIT] vs. PSR [PIT], IGA [OVA] vs. PSR [OVA] and IGA [LIV] vs. PSR [LIV], respectively (Supplementary file 4).

Prediction of transcription factors controlling seasonality genes

We found 2458, 1685, 2354 and 1225 transcriptions factor (TFs) in sets IGA [BR] vs. PSR [BR], IGA [PIT] vs. PSR [PIT], IGA [OVA] vs. PSR [OVA] and IGA [LIV] vs. PSR [LIV], respectively (Supplementary file 5). Out of these transcription factors, 53 were found common in all the datasets. We also found 1234, 674, 1169 and 458 transcription factors as unique in all the tissue specific sets i.e., IGA [BR] vs. PSR [BR], IGA [PIT] vs. PSR [PIT], IGA [OVA] vs. PSR [OVA] and IGA [LIV] vs. PSR [LIV], respectively (Fig. 4). The molecular pathway controlled by these differentially expressed TFs are predicted in this study for all the four sets of data of both the phases (Supplementary file 5). Many of our reported TFs are already reported to regulate biological clock genes in flies and humans by acting as time switch genes [12]. We found some of these are well known to control BCGs and are also conserved even in this carp fish. Example of such TFs are clock, foxl, bmal, kiss, esr, prdm, ofx, myc, srebf, bhlh and PAS. These TF and/or BCGs regulates gene regulatory network (GRN) of BPGL axis manifesting circannual rhythm of breeding [13].

Fig. 4
figure4

Venn diagram of transcription factors in all the tissue specific

Prediction of miRNA controlling seasonality genes

We identified miRNA target prediction of all tissue specific differential expressed genes with 1637 mature miRNAs of nine fish species with cut-off threshold of total score at 1000 and total energy > − 14 kcal/mol. We obtained 790, 641, 764 and 629 mature miRNAs that regulate 1444, 770, 1229 and 1044 differential expressed genes at tissue specific stages IGA [BR] vs. PSR [BR], IGA [LIV] vs. PSR [LIV], IGA [OVA] vs. PSR [OVA] and IGA [PIT] vs. PSR [PIT] (Table 5) (Supplementary Table 6). From all the tissue specific stages, we identified several miRNAs which were involved in various specific functions, like, gonadal development, oocyte and early embryo, nervous system already reported in nine other fish species (Supplementary file 7).

Table 5 List of mature miRNAs of nine fish species which regulate transcripts at tissue specific stages IGA [BR] vs. PSR [BR], IGA [PIT] vs. PSR [PIT], IGA [OVA] vs. PSR [OVA] and IGA [LIV] vs. PSR [LIV]

Our finding of predicted miRNA involved in reproductive seasonality fills the gap as there is no such report in rohu fish. This might be due to difficulty in isolation by cloning due to low expression, temporal expression and stability, tissue especially for smaller size RNA. Thus, computational identification of miRNAs binding site can be used for prediction and classification of the functional attributes. Instead of isolating small RNA afresh, the same larger RNA can be used to predict specific binding site to have information without any additional cost and time. Our enlisted predicted miRNAs can be used for exploration and validation of future biomarker research in this species. Such use has been reported in other fish species as marker for egg quality and embryonic development is a potential biomarker in rainbow trout [14]. These miRNAs are also reported to be involved in oocyte and early embryo [15,16,17] and gonadal development [18].

Discovery of putative markers of seasonality associated genes

A total of 115,534 putative SSR markers were mined from de novo assembly of transcriptome using MISA perl script. We found 89,133, 35,931, 13,645, 3312, 73, 51 mononucleotides, dinucleotides, trinucleotides, tetranucleotides, pentanucleotides and hexanucleotides, respectively (Table 6, Supplementary file 8: Uploaded at “Supplements” tab of http://webtom.cabgrid.res.in/lrsatdb/).

Table 6 SSR marker of all tissues specific DEGs

We obtained a total of 123,994 and 120,652 variants from IGA and PSR tissues samples, respectively against Danio rerio GRCz10 reference genome. A total of 99,981 variants were found common in IGA and PSR transcriptomic sets. In heterologous mapping with species Danio rerio, maximum numbers of variants were found on chromosome number 5. Maximum numbers of variants were 7617 and 7238 in IGA and PSR, respectively. In de novo-based mining, a total of 45,098 and 43,420 SNPs and Indels were obtained in IGA and PSR samples, respectively. Among this, 22,934 SNPs and Indels were common. In IGA samples, unigenes number c144686_g1_i1 (polyprotein) were found maximum number of times i.e., 116 and in PSR samples, c144749_g1_i1 unigene was found abundantly i.e., 77. Unigene “c144749_g1_i1” might be a novel transcript, as it did not show any specific matched with NCBI blastx (uncharacterized protein). Stringent filtering criteria was used for the detection of SNPs and Indels, i.e., minimum 15x [19] read depth and quality of 30 were used. Annotations of detected variants were performed against Danio rerio GRCz10 reference and identified the region where maximum number of variants were present (Supplementary file 9: Uploaded at “Supplements” tab of http://webtom.cabgrid.res.in/lrsatdb/).

Molecular mediation of reproductive seasonality and key candidate genes

We report a total of 137 key candidate genes controlling reproductive seasonality and its major pathways operating in four tissues, viz., brain, pituitary, gonad and liver. Since major molecular pathways of reproductive seasonality is orchestrated by interplay among these four tissues, thus, we computed co-expressional network to unveil this in the form of gene regulatory networks (GRN). A GRN depicts set of genes that interact with each other to control a specific cell function pertaining to a particular trait. Such studies are important to understand development, differentiation and response to environmental stimuli. By making such network of large number of DEGs can be further narrowed down [20]. SNPs of genes involved in GRN affects phenotype/trait [21] thus are valuable for eQTL discovery. Since prediction of the behaviour of GRN are rapid and economical than wet lab experiments, thus computational methods can be a valuable research tool [22]. GRN studies has been used to identify candidate genes associated with the traits [23]. GRN is further controlled by transcription factor and microRNA [24].

The GRN models of brain, pituitary, gonad and liver tissues were represented by 21, 35, 46 and 25 genes, respectively. Some of them were predicted as hub genes having similar protein-protein interaction. Magnitude of such genes expression can be correlated with trait of interest [25]. As not much information is available on seasonality aspect of these tissues, thus limited sample size transcriptomic data was used to develop “logical model” to have a basic picture of GRN reflecting qualitative attributes of genes involved in the network rather than precise quantitative relationship [26]. Similar BPGL/BPG axis has been reported with various candidate genes controlling reproduction and its synchronization with food availability and photoperiodism in other farm fishes [27].

Seasonality genes and their pathways in brain

Transcriptome analysis of brain revealed 21 hub genes which are further controlling network of > 80 genes to mediate signals of seasonality events in reproduction. GRN in brain controlling seasonality was constructed using top 100 co-expressed genes to depict protein-protein interactive network. Genes were selected at a cut off degree value of 167. A total of 20 hub genes were predicted having 8 up, 13 down regulated genes shown in Supplementary file 10. Since fish neuroendocrine system and brain is having very high similarity (82%) with human brain [28], thus available literature can be used to understand and explain the GRN operating in brain manifesting seasonal breeding. A set of 21 key genes which are differentially expressed in brain of rohu fish are involved in mediation of seasonality events. Major physiological and molecular events are summarised in Table 7 along with the respective key genes investigated in this tissue. Our co-expressional network reveals that brain is the key tissue right from sensing of daylength by visual stimuli, neuro modulation, neuro translation, neurotransmitter synthesis and release, integration of light and temperature regulated circadian gene expression, melanocytes and growth and learning regulation of feeding behaviour.

Table 7 Role of different hub genes identified in rohu fish brain

Above finding reveals that seasonality is synchronized by temperate latitudes, photoperiod, abundance of food, temperature and other environmental cues which are input to neuroendocrine output. Fish brain perceives the signal of photoperiod by melatonin pathway. Eyes are the only photoreceptive organs operating through pineal gland producing melatonin. Melatonin receptors are present in suprachiasmatic nuclei (SCN) which are the transducer of such signal acting as circadian rhythm pacemaker. Saccus vasculosus (SV) acts as a sensor of such seasonal changes in day length in fish as its removal affects the BPGL axis [29]. Photoperiod is sensed by SCN and transduced to give signal to thyroid for seasonal physiological events thus thyroid-stimulating hormone (TSH) and DIO gene also play role in seasonality [30]. There are two steps in cyclical phase of seasonality viz., PSR and IGA with slow and faster growth, respectively which is controlled by BPGL axis in rohu fish [31]. Such biological rhythm is even connected with fish auditory function for environmental cue required for reproduction and breeding especially in seasonal breeding fish [32].

Seasonality genes and their pathways in pituitary

Transcriptome analysis of pituitary revealed 35 hub genes which are further controlling network of > 65 genes involved in mediation of reproductive seasonality. To depict the integration of molecular events of seasonality operating in pituitary, we constructed co-expressional network. This protein-protein interactive network of pituitary GRN was constructed with cut off degree of 175. In pituitary of rohu fish, we found 35 hub genes with 15 and 20 up and down regulated genes, respectively (Supplementary file 10).

Our analysis further reveals that major molecular events operating in pituitary of rohu fish. Role of rohu specific hub genes is summarized in the Table 8. Our GRN depicts how MAPK signalling pathway is controlling environmental signal mediation, gonadogenesis, retenoic acid pathway controlling gametogenesis, ovarian growth and gametogenesis, increase in cell size and proliferation, oocyte maturation and ovulation, control of BPGL axis for liver metabolism, skeletal muscle development and immunity control.

Table 8 Role of different hub genes identified in rohu fish pituitary gland

Our above findings are duly supported by various other reports having role of pituitary in manifestation of seasonality of breeding in fish. The fish anterior pituitary secretes two gonadotrophin (GTH) hormones viz., GTH-I and GTH-II which are structurally similar to FSH and LH, respectively. GTH-II in goldfish is reported to stimulate gonadal growth, steroidogenesis, ovulation and sperm release with seasonality [33]. Pituitary regulates seasonality by variation in pituitary GnRH receptor activity and gonadal feedback mechanism in catfish and goldfish [34]. Teleost have two gonadotropins, viz., GtH beta I and II. Beta subunit synthesis is regulated by hypothalamic hormones and gonadal steroids. In pre-spawning and early maturing stage of fish beta I and in post spawning beta II predominates with seasonality. In tilapia, this activity is mediated by cAMP-PKA and PKC pathways [35] FSH and LH upregulates genes involved in steroid synthesis in fish [36]. This leads to growth of primary oocyte into secondary oocyte. FSH level remains increased in ovary along with FSH receptor gene expression in ovary from pre-spawning to post spawning transforming primary. During post spawning, i.e., vitellogenesis FSH continue to rise but during ovulation it declines with surge of LH [37]. Our finding of IGA and PSR stages, GSI and histology are in consonance of these reports.

Seasonality genes and their pathways in ovary

Transcriptome analysis of ovary revealed 46 hub genes which are further controlling network of > 54 genes involved in mediation of reproductive seasonality. In order to depict integrated molecular pathway operating in ovary, we constructed co-expressed genes network with a cut off degree of 175. Our study reveals that in rohu fish ovary, there are 46 hub genes having 23 up and 23 down regulated genes (Supplementary file 10). Our result shows rohu fish specific molecular events which are operated by neuroendocrine mechanism synchronized with physiological events. The role of various candidate genes/ hub genes identified in rohu fish ovary are shown in Table 9.

Table 9 Role of different hub genes identified in rohu fish gonad

Seasonal breeding fish, Labeo rohita has seasonal variation in ovarian development having seven distinct stages (phases) viz., Virgin/Immature, primary growth, perinucleolar, previtellogenic (yolk vesicle), vitellogenic (post-vitellogenic), germinal vesicle break down and spawning stage [38]. In ovary, there is coordinated activities for various regulators of ovulation like proteases, protease inhibitors, progestational steroids, eicosanoids, catecholamines and vasoactive peptides [39]. Our study supports the molecular events observed in two reproductive phases of rohu with respect to ovarian development and vitellogenesis with increase in body weight.

Our revealed ovarian hub genes are majorly involved in controlling mitogenic action of estrogen in grandulosa cell of ovary, growth of ovarian follicles, oocyte maturation, chondrocyte differentiation, prostaglandin synthesis, collagen synthesis and receptor mediated endocytosis. These events are reported in study of various other fishes. Neuroendocrine mechanism regulating reproduction through hypothalamo-pituitary-gonadal (HPG) axis is evolutionarily conserved in vertebrates [40]. Neuroendocrine control of BPGL axis and its associated gene regulatory network operating in brain, pituitary, gonad and liver are well coordinated to manifest seasonality in fish breeding by sensing photoperiod using biological clock genes along with temperature and food abundance [3]. We observed LH controlled DEG (poliovirus receptor and Bloodthirsty) involved in gonad development and maturation associated genes (ehmt2 and racgap1) which is also reported in rainbow trout. FSH modulates steroidogenic pathway as well as early germ cell proliferation and differentiation. We found four major pathways and their genes which are differentially expressed in BPGL axis viz., IGF pathway (insulin gene enhancer protein, insulin receptor, insulin-induced gene, precursor of insulin and its receptor), the TGF pathway (amh, inha and fstl3), the Wnt pathway (wisp1), and pleiotrophin (mdka) which are also reported in rainbow trout [3]. Insulin gene enhancer protein isl with its 3 isoforms families (Isl-1,2,3), observed in our dataset which is reported to play role in regulation of IGF pathway [41]. SV organ acts as a sensor of seasonality by sensing the day length/photoperiod by rhodopsin molecule. Thus, rhodopsin gene family must be one of the major hub in brain tissue acting as a transducer for neuroendocrine output in form of GnRH [29].

Seasonality genes and their pathways in liver

Transcriptome analysis of liver revealed 25 hub genes which are further controlling network of > 75 genes involved in metabolism associated with reproductive seasonality. The co-expressional network was constructed with a cut off degree 165. Among 25 hub genes, 9 were up and 16 were down regulated as shown in Supplementary file 10.

Our study reveals that in liver of rohu fish, JAK/STAT signaling is major pathway controlling entire liver metabolism to manifest reproductive seasonality with synchronized events in ovarian development. This study further reveals how rohu specific genes are involved in modulation of glycolysis, gluconeogenesis, oxidative phosphorhylation, apolipo protein production, MHC response, retenoic acid, betain and GABA, immune response, cholesterol and lipid transport while meanifesting reproductive seasonality. List of key candidate genes controlling these events in rohu fish liver with their specific functions are summarized in Table 10.

Table 10 Role of different hub genes identified in rohu fish liver

We also found upregulation of Wnt and Notch signaling genes in all the four tissues of BPGL axis having multifarious role including higher perceptivity of auditory system in breeding season. This is in response to estrogen induced seasonal changes [42]. Besides tissue specific GRN, there are some common set of genes which are expressed in different magnitudes in different tissue with respect to seasonality. Expression of plectin gene was found to be most abundant in liver and gonads rather than brain and pituitary. Plectin proteins are abundant (75%) in cytoskeleton of ovary, which is a major phosphoacceptor, in phosphorylation, thus important for the protein’s association with the cytoskeleton [43]. Plectin is a hemidesmosomal protein that mediates hyperproliferation of gonads and liver. It also increases the expression of estrogen receptor in fish gonads which triggers seasonal response in ovary manifesting seasonality [44]. Claudin gene was found to be upregulated in all the 4 tissues but its highest expression was found in gonads. Claudin in teleost fish is reported to be tissue specific in expression. It functions as pore forming TJ protein increasing luminal fluid accumulation and volume expansion in ovary [45] as observed in histological picture (Fig. 1).

Above findings are reflective of neuro-endochrinological and physiological mechanism which operates in all oviparous vertebrates, including fish. For such events, egg-yolk and chorionic proteins like vitellogenin and choriogenin are synthesized heterologously in the liver [46]. Our findings are in consonance of the reports in other fishes. Fish liver is reported to be involved in synthesis, degradation, transportation, and storage of lipid. These major activities are well coordinated by BPGL axis along with oocyte growth and development in tune of seasonality. The axis synchronises the food availability and ingestion behaviour along with photoperiod, temperature and gonadal maturation leading to seasonal reproduction in temperate fishes [47]. Our observation in two phases of reproductive seasonality observed in rohu is also having contrasting GSI and HSI. Similar finding is also reported in other fish where vitellogenesis associated with GSI and HSI (hepato-somatic index) due to increase in total protein, glycogen and cholesterol content of the ovary and liver and stored fat reserve, which is used for gonadal maturation [48].

qPCR for validation of DEG

To validate the finding, qPCR analysis was performed. Relative gene expression of randomly selected 20 differentially expressed genes along with one housekeeping gene (beta actin) showed similar magnitude with the computed log fold change value except three genes (Supplementary file 11). This could be due to cross-reactivity in the designed primer [49].

Biological clock genes in carp fish for circannual breeding

Blast hit of four sets of DEGs revealed presence of at least 32 well known biological genes in carp fish. These sets of transcripts represent carp fish specific BCG sequences which is also a direct evidence for presence of these biological clock genes in this fish. These differentially expressed BCGs in BPGL tissues are involved in manifestation of circannual rhythm of breeding. In case of European fish reproductive axis gene expression in brain is known to be modulated by photoperiod which is regulated by biological clock genes and conserved transcription factors [3].

This is the first report having transcriptomic evidence of biological clock genes in rohu mediating seasonality over BPGL axis controlling reproductive behaviour. Interestingly, all these 32 biological clock genes are well conserved from lower vertebrate like cyclostome (Lamprey) to higher vertebrate like mammals. In pituitary, we found differential expression of Bmal1, Clock, Per1, Per2, Cry1 and Cry2 genes sensing duration of photoperiod for synchronization in BPGL axis. This culminates in manifestation of reproductive seasonality with circannual rhythm [50].

Few already reported genes of rohu fish involved in seasonality regulation were also found in our dataset. For example Kisspeptin and its receptor Kissr [51]. We found differential expression of known sex determining and seasonality genes for example aromatase cyp19a1a, estrogen receptor esr1a, and foxl2 in ovaries similar to Iberian cyprinid fish (Squalius pyrenaicus) [52] and GnRH3, similar to other teleost [53]. Similar gene family like Kiss and gnrh are reported to control reproductive axis under influence of contrasting day length as photoperiod in fish sea bass [54].

Though expression of biological clock genes are tissue specific as well as magnitude specific [55] but they are well coordinated and synchronized across various organs involved in sensing photoperiod, ambient temperature, food supplies, immunity, locomotory activity and body growth [56] enabling the fish for seasonal breeding. Isotocin-neurophysin gene was found to be upregulated in brain. Isotocin is teleost homologue of mammalian oxytocin. In case of zebra fish it has been reported that isotocin stimulates the proliferation of cells enhancing the functional activities of ionocytes. Isotocin stimulates the proliferation of cells in fish gonads [57]. Similarly, higher expression of synaptotagmin gene was also found in hypothalamus and pituitary which is involved in neurotransmitter secretion [58].

Reproductive seasonality in fish is reported to be controlled by hypothalamus/ SV organ which acts as a sensor of seasonal changes in day length [29]. Existence of this organ is already reported in Labeo rohita [59]. Present transcriptome based study in rohu also confirms the seasonality control by similar regulatory pathway of BPGL axis. Seasonality in tropical fishes is an adaptation to environment and its optimal growth and perpetuation of species in its specific ecological niches. Fish represents vertebrate circadian timing system operated by circadian biological clock and photoreceptor organ. This ability originates during embryogenesis and continues till reproductive age [60]. Seasonal reproductive behaviour is induced by environmental cues. Various sets of biological clock sense the environmental cues and acts as a transducer eliciting neuroendocrine mediated developmental and physiological events under coordination of BPGL axis [61].

Evolutionary conservation of BCGs from cyclostome to eutherian mammals

Extent of conservation in primitive jawless fish like lamprey, cartilaginous, teleost, lung fish; amphibian, bird, reptile and three mammalian classes (monotreme, marsupials and eutherians) is summarized in Table 11. The comparative analysis of carp fish BCGs sequence with different vertebrate phyla indicates that these biological clock genes are well conserved in all. Maximum percentage similarity of all the 32 biological clock genes of rohu was found with jawless fish, lamprey. However, query coverage was maximum up to 100% with cyprinid (common carp) and zebrafish up to 100% (Supplementary file 12). Circular plot shows further qualitative depiction of extent of conservation in selected six phyla namely, amphibian, fish, birds reptiles, eutherian mammals and ovoviviparous mammals (Fig. 5).

Table 11 Similarity of rohu biological clock gene across different classes of chordates showing extent of conservation
Fig. 5
figure5

Extent of rohu biological clock gene conservation across different classes of chordates

These findings infer that genes and neuroendocrine pathways sensing environmental cues are well conserved during evolutionary process of species diversification into different classes of vertebrates [62]. Such biological clock genes are reported to be conserved in other fishes also like Atlantic cod fish having 18 conserved biological clock genes [63]. Thus, investigation of present transcriptomic landscape revealed that at least 32 reported biological clock genes of other diverse species are also present in rohu fish. With seasonal changes these genes are differentially expressed synchronizing environment and reproduction.

In evolution, biological clock has been found universal, which originated 2.5 billion years ago [64]. This is because of universality in response of all evolving organisms against change in day length with rotation of earth having daily and seasonal variation in environment. Environmental signal entrains the BCGs to trigger seasonal behaviour for annual biological rhythm for reproduction/migration [65]. Magnitude of biological rhythm varies from 20 min (bacterial reproduction) to daily as circadian rhythm but such rhythm may have much longer duration like tidal rhythm (14 days), lunar rhythm (28 days), circannual rhythm one year, sometime it may be of 17 years (cicada species emergence from ground) to 60 years (bamboo flowering) [66].

Circadian rhythm is well connected to lunar, tidal, circannual rhythm in hibernating mammals. It is also present in invertebrates [67]. Circannual rhythm of breeding is well known in farm fishes, birds and mammals [68]. Artificial manipulation of day length affects rhythmicity of biological clock affecting reproductive activities like ovulation, spermatogenic activity, gamete quality and sexual behaviour in fish, bird and mammals [69]. Oscillation in biological rhythm has adaptive significance to reduce competition between species by occupying differential species specific spatial niches. This is done by differential phasing of biological rhythm which is entrained by both internal and external force [8]. Thus, we observe ubiquity of circadian rhythm phylogenetically having strong adaptive significance [70]. This evolutionary comparative genomic study endorses that fish biological clock can be used as a model for research and knowledge discovery in terms of depicting molecular and cellular mechanism in response to environment in applied research like drug efficacy for jetlag etc. [71].

Web genomic resources of evolutionary conserved BCGs and candidate genes of seasonality

Present investigation of environmentally triggered reproduction can be of much relevance in the endeavour of out of seasonal breeding leading to higher fish productivity [56]. The present genomic resource catalogues 75,554 differentially expressed transcripts, 115,534 SSRs, 65,584 SNPs from de novo assembly, 514 pathways, 5379 TFs, 187 mature miRNA regulating 1576 DEGs and 32 BCGs. It also catalogues ready to use primers. This is freely accessibility for academic purpose.

Fish transcriptome can be a powerful tool to establish relationship between genotype and phenotype using structural and functional annotation approach [72]. In rohu fish, SSRs and SNPs have been reported to be associated with important disease resistant trait against aeromoniasis bacterial disease [72]. Our transcriptome database can be used as genomic resources as it is having putative markers useful for similar association studies in future [73]. Such genomic resources having SNP markers has also been used in linkage mapping and QTL discovery in rohu fish [74].

Candidate genes involved in GRN of ovary can be used in SNP discovery and future association studies in fish selection. Liver GRN genes can be used for biomonitoring of fish breeding area with seasonality to predict magnitude of fish productivity [75]. GRN of liver having TOR signaling pathway which is associated with lipogenesis can be over-activated to utilize glucose more efficiently by genetic selection [76]. Similarly, genes expressed in liver associated with innate immune system can be used for biomonitoring of fish mortality [77].

Conclusion

This is the first report in carp fish rohu deciphering candidate genes of circannual biological rhythm of reproduction along with BCGs. Study reveals reproductive seasonality of rohu fish with elucidation of neuroendocrine physiological events at molecular level. We present the transcriptomic landscape of BPGL axis operating through various gene regulatory network associated with tissue involved in manifestation of circannual breeding having environmental synchronization by BCGs and BPGL axis. Differentially expressed genes in BPGL tissues along with putative molecular markers (SSRs and SNPs), pathways, transcription factors, mature miRNA and BCGs are available in the form of web genomic resources (LrSATDb: http://webtom.cabgrid.res.in/lrsatdb/). This study reveals atleast 32 BCGs in this carp fish for the first time which are also highly conserved across jawless, cartilaginous, teleost fish, amphibian, reptile, birds, prototheria, marsupials and placental mammals. This is because of universal molecular mechanism of rhythmicity in response to environment and earth rotation. Biological rhythm and its amplitude has different adaptive and reproductive significance during speciation and evolutionary divergence of different taxa.

Gene regulatory network operating in BPGL tissues mediating photo-periodism sensing, neuroendocrine secretion, metabolism and yolk synthesis in liver, gonadal maturation, muscular growth with sensory and auditory perception in this fish are elucidated. Study further reveals that fish can be a good model for research on biological clock besides its immense use in the endeavour of reproductive efficiency improvement.

Methods

Ethics statement

As per Chapter 4, Section 15(1) of the Prevention of Cruelty to Animals (PCA) Act 1960 of Govt of India, using guidelines recommended by CPCSEA (Committee for the Purpose of Controlling and Supervising Experiments in Animals), Institute Animal Ethics Committee of ICAR-CIFA approved the experiments on 1st March 2017.

Rearing and maintenance of the brooder in the pond

Adult male and female rohu of two-year-old (body weight 0.8 to 1.5 kg) was reared in earthen ponds of 0.05 ha in the carp breeding unit of ICAR-CIFA farm (Latitude 20°1′06″–20°11′45″N, Longitude 80°50′52″–85°51′35″E, 45 m MSL) as per standard procedure. For optimum condition, stocking density was maintained at 1500 kg/ha along with routine monitoring of water temperature, pH, dissolved oxygen, total alkalinity and total hardness at 24–30 °C, 7.5–8.5, 100–140 ppm, and 100–130 ppm, respectively. Fishes were fed with a commercial feed (CP-24%) @2% of their body weight once daily.

Reproductive phase identification by GSI estimation with ovarian morphology and histology

Length and body weight of each euthanized individual were recorded before opening of the abdominal cavity. Fishes were humanely sacrificed as per SOP for dissection permitted by the Institute. The gonad position in the posterior half of the abdominal cavity on both side below the swim bladder was observed. Testis in male or ovary in female was identified by using the features like position of the gonadal ducts, length, appearance and color etc. Only the ovaries in females were taken out intact from the peritoneal covering, photographed and weighed. GSI was calculated using the formula: GSI = weight of the ovary X 100/ weight of body [78]. For histology, middle portion of the right and left lobes of the gonads from the sampled fish were taken after dissection. The ovary samples (1–1.5 mm thickness) were immersed in Bouin’s fixative for 24 h and 5 μm thick sections were cut using a mechanical microtome (WESWOX Optik Rotary Microtome, Ambala Cantt, India). Sections were stained using Delafield’s haematoxylin, counterstained by eosin and were processed for examination with light microscopy using routine histological procedures [79]. The diameter of the oocytes was measured by the oculometer standardized against a stage micrometer on random sampling basis.

Tissue sample collection

The fishes were euthanized with MS-222 at 300 mg/L before dissection. Brain, pituitary, ovary and liver tissues of adult rohu female were collected during post spawning resting (PSR) phase in December and initiating gonad activity (IGA) phase in February, respectively. Tissues were collected from minimum ten fishes for each phase, quickly frozen in liquid nitrogen and stored at -80 °C, until used for RNA extraction. In order to avoid biological variability within each phase of fishes we used one ovary for RNA extraction and other one for histological and GSI confirmation.

RNA extraction and sequencing

Total RNA was extracted using Qiagen RNeasy mini kit using manufacturer’s protocol. RNA concentration and purity were estimated using Bioanalyzer 2100 (Agilent, Santa Clara, CA, USA) and samples having RIN (RNA Integrity Number) > 7.0 were used for sequence data generation. Library preparation was performed using Illumina TruSeq RNA library protocol and sequenced on Illumina HiSeq 2000 platform (Illumina, San Diego, CA) to generate 100 nucleotide pair-end reads.

Pre-processing, de novo assembly, abundance estimation and annotation

Data were pre-processed for quality check and low-quality reads, reads with ambiguous bases ‘N’ and adaptor sequences were trimmed using FASTQC [80] and trimmomatic [81] tool. The processed paired end reads of PSR and IGA phases were pooled together to make rohu specific de novo transcriptome assembly using Trinity [82] which was used to obtain differentially expressed genes. We preferred this transcriptome assembly rather than non-seasonal breeding zebrafish reference sequence to avoid exclusion of rohu specific novel transcripts/ genes. Further, abundance estimation was performed using RSEM (RNA-Seq by Expectation Maximization) [83] and DEGs were identified from the data set of four tissues. DEGs from comparative studies of IGA season vs. PSR season were identified in each of the BPGL tissues using edgeR package [84].

Homology search and functional characterization of candidate genes

Standalone NCBI blast was used for sequence similarity search of all novel transcripts and differential expressed genes against the NCBI non-redundant database (nr.32). Annotation was done for the DEGs along with functional categorization, interproscan, mapping and identification of KEGG pathways using Blast2GO PRO 3.1 software [11]. Another ciprinid fish, Danio rerio transcription factors available in AnimalTFDB 2.0 database [85] were retrieved to predict similar transcription factors but differential expression of these genes in BPGL tissues were done using de novo transcriptome assembly of rohu.

Discovery of putative markers of seasonality associated genes

Simple Sequence Repeats (SSRs) were detected from L. rohita by using MISA perl script [86]. Mononucleotides, dinucleotides, trinucleotides, tetranucleotides, pentanucleotides and hexanucleotides repeating units were identified with default parameters i.e. 10 repeating units for mono-, 6 repeating units for di-, 5 for tri-, tetra-, penta- and hexa-nucleotides.

SNPs and Indels were detected using two references i.e. Danio rerio reference genome data and generated de novo transcription reference assembly. The reference genome of Danio rerio GRCz10 was retrieved from ensemble database (http://asia.ensembl.org/Danio_rerio/Info/Index?redirect=no) as it is closely related species of L. rohita and both belongs to same family, Cyprinidae. We separately identified the SNPs and Indels from IGA and PSR tissue samples. Burrows-Wheeler Aligner (BWA) tool [87] was used to alignment and SAMtools package used for calling SNPs and Indels [88]. Annotation of SNPs and Indels were performed using SnpEff tool (http://snpeff.sourceforge.net/). Stringent filtering criteria was used for the detection of SNPs and Indels, i.e., minimum 15x [19] read depth and quality of 30 were used. Annotations of detected variants were performed against Danio rerio GRCz10 reference and identified the region where maximum number of SNPs were present.

miRNA target prediction of DEGs

A total of 1637 mature miRNA sequences of nine fish species such as Danio rerio, Fugu rubripes, Cyprinus carpio, Hippoglossus hippoglossus, Ictalurus punctatus, Oryzias latipes, Paralichthys olivaceus, Salmo salar and Tetraodon nigroviridis were retrieved from miRBase database release 21 [89]. MiRNA target predictions of L. rohita from DEGs of all stages were performed by using miRanda-3.3a tool [90].

Construction of gene regulatory network

Gene regulatory network (GRN) was constructed using up and down regulated DEGs having > 8 fold difference in expression. Visualization was done using open source tool, Cytoscape version 3.2.1 [91]. Network centrality was analysed using Network Analyzer plug-in. Analysis of degree, betweenness and stress were used to identify hub genes and its network. The “logical model” was constructed based on sample size [26].

qPCR for validation of DEG

For quantitative PCR, 20 transcripts were randomly selected. Primers were designed using Primer 3 software [92]. The primer were synthesized in 10 nm scale and purified by HPLC. The total RNA were converted into cDNA using affinity script qPCR cDNA synthesis kit (Agilent Technologies, USA) as per manufacturer’s protocol. Gene expressions were measured using SYBR Green chemistry (Brilliant II SYBR Green qPCR master mix, Agilent Technologies, USA) with standard 40 cycles in Stratagene mx3005P instrument (Agilent Technologies, USA). The dissociation curve analysis was performed to ensure specificity of amplification. PCR efficiency was calibrated by 10-fold serial dilution of PCR product of each gene. Specific gene amplification was cross-checked by melt curve analysis with each primer set. Each PCR was done with technical duplicate along with negative control (without template). The mean Ct value of technical replicates were calculated. Magnitude of differential gene expression were calculated in terms of ΔΔCT fold change value as described by Pfaffl [93]. Housekeeping beta actin gene was used as reference to normalize the qPCR data. This gene was selected based on literature reporting successful use in similar targeted tissues in rohu fish [73]. Fold change values were computed by log transformation (Log2).

Discovery of biological clock genes in carp fish

Four sets of DEG, each set obtained from each of BPGL tissue were subjected to homology search using BLAST. From literature, names of known BCGs [50, 63, 79] were used to identify carp fish specific BCG transcripts as an evidence of BCGs. These selected BCG transcripts were used for further analysis (supplementary file 12).

Evolutionary conservation of BCGs from cyclostome to eutherian mammals

To identify biological clock genes and their extent of conservation, comparative genomics approach was used. For this study, 32 genes were selected viz. Bmal2,clock, cry1,cry2, foxl2, GNRH, kisspeptin, LH (Luteinizing hormone), Per1, Per2, Per3, esr1a, PVR, Bty, wisp1, mdka, amh, ependymin, CaMK2, ehmt1, ehmt2, racgap1, fstl3, fstl4, fstl5, prdm, Otx, Myc, srebf1, pac, picalm and aanat and homology search was done using Blastn programme of NCBI of 10 different genomes representing various classes of chordates such as jawless fish Petromyzon marinus (GCA_000148955.1); cartilaginous fish whale shark Rhincodon typus (GCA_001642345.2); teleost Danio rerio (GCA_000002035.4); Common carp Cyprinus carpio (GCA_000951615.2); Amphibia tropical clawed frog Xenopus tropicalis (GCA_000004195.3); Bird Gallus gallus (GCA_000002315.3); reptile Green anole Anolis carolinensis (GCA_000090745.2); egg laying mammal Ornithorhynchus anatinus (GCF_000002275.2), ovoviviparous mammal Monodelphis domestica (GCF_000002295.2), eutherian mammals Ovis aries (sheep) (GCA_000298735.2).

Web genomic resources of evolutionary conserved BCGs and candidate genes of seasonality

An online relational database of rohu fish transcriptome was developed which catalogues tissue wise transcripts/contigs, putative SSRs, SNPs, Indels, transcription factors, miRNA targets representing two reproductive phases (IGA and PSR). The architecture is “three-tier architecture” viz., client-, middle- and database tier. This genomic resource is freely accessible for non-commercial use at http://webtom.cabgrid.res.in/lrsatdb/. In order to browse and query, user can go through the web pages in client tier. All the information are available in various tables corresponding to MySQL in the database tier. Server side scripting in PHP was done in the middle tier for database connectivity, query execution and fetching. In order to generate primers over selected markers, Primer3 executable was integrated at the backend.

Availability of data and materials

The transcriptome dataset of the study used in this article are available in the NCBI repository with following accessions and is kept at hold till the publication. These would be made public after publication. (Bioproject: PRJNA401304; BioSamples: SAMN07602341, SAMN07602342, SAMN07602343, SAMN07602344, SAMN07602345, SAMN07602346, SAMN07602347, SAMN07602348).

Abbreviations

BCGs:

Biological clock genes

BPGL:

Brain-Pituitary-Gonad-Liver

QPCR:

Quantitative polymerase chain reaction

GRN:

Gene regulatory networks

EST:

Expressed sequence tag

SSR:

Simple sequence repeats

SNP:

Single nucleotide polymorphism

DEG:

Differentially expressed genes

PSR:

Post spawning resting

IGA:

Initiation of gonadal activity

POF:

Primary ovarian follicles

SOF:

Secondary ovarian follicles

BR:

Brain

PIT:

Pituitary

LIV:

Liver

OVA:

Ovary

KEGG:

Kyoto Encyclopedia of Genes and Genomes

TF:

Transcription factors

Blast:

Basic local alignment search tool

References

  1. 1.

    Jhingran VG: Fish and fisheries of India. 1975, 1–954.

  2. 2.

    Mistlberger RE, Antle MC. Entrainment of circadian clocks in mammals by arousal and food. Essays Biochem. 2011;49(1):119–36. https://doi.org/10.1042/bse0490119.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Martins RS, Gomez A, Zanuy S, Carrillo M, Canário AV. Photoperiodic modulation of circadian clock and reproductive axis gene expression in the pre-pubertal European sea bass brain. PLoS One. 2015;10(12):e0144158. https://doi.org/10.1371/journal.pone.0144158.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Kullgren A, Jutfelt F, Fontanillas R, Sundell K, Samuelsson L, Wiklander K, et al. The impact of temperature on the metabolome and endocrine metabolic signals in Atlantic salmon (Salmo salar). Comp Biochem Physiol A Mol Integr Physiol. 2013;164(1):44–53. https://doi.org/10.1016/j.cbpa.2012.10.005.

  5. 5.

    Isorna E, De Pedro N, Valenciano AI, Alonso-Gómez ÁL, Delgado MJ. Interplay between the endocrine and circadian systems in fishes. J Endocrinol. 2017;232(3):R141–59. https://doi.org/10.1530/JOE-16-0330.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Menaker M. Biological clocks. Bioscience. 1969;19(8):681–92. https://doi.org/10.2307/1294894.

    Article  Google Scholar 

  7. 7.

    Sahu DK, Panda SP, Panda S, Das P, Meher PK, Hazra RK, et al. Identification of reproduction-related genes and SSR-markers through expressed sequence tags analysis of a monsoon breeding carp rohu, Labeo rohita (Hamilton). Gene. 2013;524(1):1–4. https://doi.org/10.1016/j.gene.2013.03.111.

  8. 8.

    Sahu DK, Panda SP, Meher PK, Das P, Routray P, Sundaray JK, et al. Construction, de-novo assembly and analysis of transcriptome for identification of reproduction-related genes and pathways from rohu, Labeo rohita (Hamilton). PLoS One. 2015;10(7):e0132450. https://doi.org/10.1371/journal.pone.0132450.

  9. 9.

    Koç ND, Aytekin Y, Yüce R. Ovary Maturatıon stages and histological investigation of ovary of the zebrafish (Danio rerio). Brazilian Arch Biol Technol-An Int J. 2008;51(3):513–22. https://doi.org/10.1590/S1516-89132008000300010.

    Article  Google Scholar 

  10. 10.

    Sundararaj BI, Vasal S. Photoperiod and temperature control in the regulation of reproduction in the female catfish Heteropneustes fossilis. J Fish Board Canada. 1976;33(4):959–73. https://doi.org/10.1139/f76-123.

    Article  Google Scholar 

  11. 11.

    Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6. https://doi.org/10.1093/bioinformatics/bti610.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Ishida N, Kaneko M, Allada R. Biological clocks. Proc Natl Acad Sci U S A. 2000;97(16):9347.

    CAS  Google Scholar 

  13. 13.

    Betancor MB, Mcstay E, Minghetti M, Migaud H, Tocher DR, Davie A. Daily Rhythms in Expression of Genes of Hepatic Lipid Metabolism in Atlantic Salmon (Salmo salar L.). PLoS ONE. 2014;9(9):e106739.

    Article  Google Scholar 

  14. 14.

    Ma H, Weber GM, Hostuttler MA, Wei H, Wang L, Yao J. MicroRNA expression profiles from eggs of different qualities associated with post-ovulatory ageing in rainbow trout (Oncorhynchus mykiss). BMC Genomics. 2015;16(1):201. https://doi.org/10.1186/s12864-015-1400-0.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Ramachandra RK, Salem M, Gahr S, Rexroad CE, Yao J. Cloning and characterization of microRNAs from rainbow trout (Oncorhynchus mykiss): their expression during early embryonic development. BMC Dev Biol. 2008;8(1):41. https://doi.org/10.1186/1471-213X-8-41.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Soni K, Choudhary A, Patowary A, Singh AR, Bhatia S, Sivasubbu S, et al. Pillai B: miR-34 is maternally inherited in Drosophila melanogaster and Danio rerio. Nucleic Acids Res. 2013;41(8):4470–80. https://doi.org/10.1093/nar/gkt139.

  17. 17.

    Juanchich A, Le Cam A, Montfort J, Guiguen Y, Bobe J. Identification of differentially expressed miRNAs and their potential targets during fish ovarian development. Biol Reprod. 2013;88(5):128–31. https://doi.org/10.1095/biolreprod.112.105361.

    Article  PubMed  Google Scholar 

  18. 18.

    Presslauer C, Bizuayehu TT, Kopp M, Fernandes JM, Babiak I. Dynamics of miRNA transcriptome during gonadal development of zebrafish. Sci Rep. 2017;7(1):43850. https://doi.org/10.1038/srep43850.

    Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Uitdewilligen JG, Wolters AM, Bjorn B, Borm TJ, Visser RG, van Eck HJ. A next-generation sequencing method for genotyping-by-sequencing of highly heterozygous autotetraploid potato. PLoS One. 2013;8(5):e62355. https://doi.org/10.1371/journal.pone.0062355.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Hur B, Chae H, Kim S. Combined analysis of gene regulatory network and SNV information enhances identification of potential gene markers in mouse knockout studies with small number of samples. BMC Med Genet. 2015;8(2):S10. https://doi.org/10.1186/1755-8794-8-S2-S10.

    CAS  Article  Google Scholar 

  21. 21.

    Kim DC, Wang J, Liu C, Gao J. Inference of SNP-gene regulatory networks by integrating gene expressions and genetic perturbations. Biomed Res Int. 2014;629697:1–9. https://doi.org/10.1155/2014/629697.

    Article  Google Scholar 

  22. 22.

    Karlebach G, Shamir R. Modelling and analysis of gene regulatory networks. Nat Rev Mol Cell Biol. 2008;9(10):770–80. https://doi.org/10.1038/nrm2503.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Zhu J, Alvarez S, Marsh EL, LeNoble ME, Cho IJ, Sivaguru M, et al. Cell wall proteome in the maize primary root elongation zone. II. Region-specific changes in water soluble and lightly ionically bound proteins under water deficit. Plant Physiol. 2007;145(4):1533–48. https://doi.org/10.1104/pp.107.107250.

  24. 24.

    Narang V, Ramli MA, Singhal A, Kumar P, de Libero G, Poidinger M, et al. Automated identification of core regulatory genes in human gene regulatory networks. PLoS Comput Biol. 2015;11(9):e1004504. https://doi.org/10.1371/journal.pcbi.1004504.

  25. 25.

    Cai X, Bazerque JA, Giannakis GB. Inference of gene regulatory networks with sparse structural equation models exploiting genetic perturbations. PLoS Comput Biol. 2013;9(5):e1003068. https://doi.org/10.1371/journal.pcbi.1003068.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Vijesh N, Chakrabarti SK, Sreekumar J. Modeling of gene regulatory networks: a review. J Biomed Sci Eng. 2013;6(02):223–31. https://doi.org/10.4236/jbise.2013.62A027.

    Article  Google Scholar 

  27. 27.

    Taranger GL, Carrillo M, Schulz RW, Fontaine P, Zanuy S, Felip A, et al. Control of puberty in farmed fish. Gen Comp Endocrinol. 2010;165(3):483–515. https://doi.org/10.1016/j.ygcen.2009.05.004.

  28. 28.

    Matsui H. Dopamine system, cerebellum, and nucleus ruber in fish and mammals. Develop Growth Differ. 2017;59(4):219–27. https://doi.org/10.1111/dgd.12357.

    CAS  Article  Google Scholar 

  29. 29.

    Nakane Y, Ikegami K, Iigo M, Ono H, Takeda K, Takahashi D, et al. The saccus vasculosus of fish is a sensor of seasonal changes in day length. Nat Commun. 2013;4(1):2108. https://doi.org/10.1038/ncomms3108.

  30. 30.

    Dardente H, Hazlerigg DG, Ebling FJ. Thyroid hormone and seasonal rhythmicity. Front Endocrinol. 2014;5:19. https://doi.org/10.3389/fendo.2014.00019.

    Article  Google Scholar 

  31. 31.

    Danzmann RG, Kocmarek AL, Norman JD, Rexroad CE, Palti Y. Transcriptome profiling in fast versus slow-growing rainbow trout across seasonal gradients. BMC Genomics. 2016;17(1):60. https://doi.org/10.1186/s12864-016-2363-5.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Basinou V, Park JS, Cederroth CR, Canlon B. Circadian regulation of auditory function. Hear Res. 2017;347:47–55. https://doi.org/10.1016/j.heares.2016.08.018.

    Article  PubMed  Google Scholar 

  33. 33.

    Trudeau VL. Neuroendocrine regulation of gonadotrophin II release and gonadal growth in the goldfish. Carassius auratus Reviews of reproduction. 1997;2(1):55–68. https://doi.org/10.1530/ror.0.0020055.

  34. 34.

    Habibi HR, De Leeuw R, Nahorniak CS, Goos HT, Peter RE. Pituitary gonadotropin-releasing hormone (GnRH) receptor activity in goldfish and catfish: seasonal and gonadal effects. Fish Physiol Biochem. 1989;7(1–6):109–18. https://doi.org/10.1007/BF00004696.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Melamed P, Rosenfeld H, Elizur A, Yaron Z. Endocrine regulation of gonadotropin and growth hormone gene transcription in fish. Comp Biochem Physiol C: Pharmacol Toxicol Endocrinol. 1998;119(3):325–38. https://doi.org/10.1016/S0742-8413(98)00021-8.

    CAS  Article  Google Scholar 

  36. 36.

    Sambroni E, Lareyre JJ, Le Gac F. Fsh controls gene expression in fish both independently of and through steroid mediation. PLoS One. 2013;8(10):e76684. https://doi.org/10.1371/journal.pone.0076684.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Guzmán JM, Luckenbach JA, Yamamoto Y, Swanson P. Expression profiles of Fsh-regulated ovarian genes during oogenesis in coho salmon. PLoS One. 2014;9(12):e114176. https://doi.org/10.1371/journal.pone.0114176.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  38. 38.

    Sen U, Mukherjee D, Bhattacharyya SP, Mukherjee D. Seasonal changes in plasma steroid levels in Indian major carp Labeo rohita: influence of homologous pituitary extract on steroid production and development of oocyte maturational competence. Gen Comp Endocrinol. 2002;128(2):123–34. https://doi.org/10.1016/S0016-6480(02)00060-6.

    CAS  Article  PubMed  Google Scholar 

  39. 39.

    Goetz FW, Garczynski M. The ovarian regulation of ovulation in teleost fish. Fish Physiol Biochem. 1997;17(1–6):33–8. https://doi.org/10.1023/A:1007765902327.

    CAS  Article  Google Scholar 

  40. 40.

    Prasad P, Ogawa S, Parhar IS. Role of serotonin in fish reproduction. Front Neurosci. 2015;9:195. https://doi.org/10.3389/fnins.2015.00195.

    Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Hemati N, Ross SE, Erickson RL, Groblewski GE, MacDougald OA. Signaling pathways through which insulin regulates CCAAT/enhancer binding protein α (C/EBPα) phosphorylation and gene expression in 3T3-L1 adipocytes correlation with GLUT4 gene expression. J Biol Chem. 1997;272(41):25913–9. https://doi.org/10.1074/jbc.272.41.25913.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Faber-Hammond J, Samanta MP, Whitchurch EA, Manning D, Sisneros JA, Coffin AB. Saccular transcriptome profiles of the seasonal breeding Plainfin midshipman fish (Porichthys notatus), a teleost with divergent sexual phenotypes. PLoS One. 2015;10(11):e0142814. https://doi.org/10.1371/journal.pone.0142814.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Herrmann H, Wiche G. Specific in situ phosphorylation of plectin in detergent-resistant cytoskeletons from cultured Chinese hamster ovary cells. J Biol Chem. 1983;258(23):14610–8. https://doi.org/10.1016/S0021-9258(17)43906-8.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Wada-Hiraike O, Imamov O, Hiraike H, Hultenby K, Schwend T, Omoto Y, et al. Role of estrogen receptor β in colonic epithelium. Proc Natl Acad Sci U S A. 2006;103(8):2959–64. https://doi.org/10.1073/pnas.0511271103.

  45. 45.

    Kolosov D, Bui P, Chasiotis H, Kelly SP. Claudins in teleost fishes. Tissue Barriers. 2013;1(3):e25391. https://doi.org/10.4161/tisb.25391.

    Article  PubMed  PubMed Central  Google Scholar 

  46. 46.

    Arukwe A, Goksøyr A. Eggshell and egg yolk proteins in fish: hepatic proteins for the next generation: oogenetic, population, and evolutionary implications of endocrine disruption. Comp Hepatol. 2003;2(1):4. https://doi.org/10.1186/1476-5926-2-4.

    Article  PubMed  PubMed Central  Google Scholar 

  47. 47.

    Ballesteros TM, Torres-Mejia M, Ramírez-Pinilla MP. How does diet influence the reproductive seasonality of tropical freshwater fish?: a case study of a characin in a tropical mountain river. Neotropical Ichthyology. 2009;7(4):693–700. https://doi.org/10.1590/S1679-62252009000400019.

    Article  Google Scholar 

  48. 48.

    Shaikh Abdullah S, Lohar PS. Biochemical composition and gonadosomatic index of three major carps in Hatnoor reservoir, Maharashtra. India J Ecobiotechnol. 2011;3(6):1–4.

    CAS  Google Scholar 

  49. 49.

    Spandidos A, Wang X, Wang H, Dragnev S, Thurber T, Seed B. A comprehensive collection of experimentally validated primers for polymerase chain reaction quantitation of murine transcript abundance. BMC Genomics. 2008;9(1):633. https://doi.org/10.1186/1471-2164-9-633.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  50. 50.

    Lincoln GA, Andersson H, Loudon A. Clock genes in calendar cells as the basis of annual timekeeping in mammals--a unifying hypothesis. J Endocrinol. 2003;179(1):1–13. https://doi.org/10.1677/joe.0.1790001.

    CAS  Article  PubMed  Google Scholar 

  51. 51.

    Chi L, Li X, Liu Q, Liu Y. Photoperiod regulate gonad development via kisspeptin/kissr in hypothalamus and saccus vasculosus of Atlantic salmon (Salmo salar). PLoS One. 2017;12(2):e0169569. https://doi.org/10.1371/journal.pone.0169569.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Machado MP, Matos I, Grosso AR, Schartl M, Coelho MM. Non-canonical expression patterns and evolutionary rates of sex-biased genes in a seasonal fish. Mol Reprod Dev. 2016;83(12):1102–15. https://doi.org/10.1002/mrd.22752.

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    Bar I, Cummins S, Elizur A. Transcriptome analysis reveals differentially expressed genes associated with germ cell and gonad development in the southern bluefin tuna (Thunnus maccoyii). BMC Genomics. 2016;17(1):217. https://doi.org/10.1186/s12864-016-2397-8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Espigares F, Rocha A, Gómez A, Carrillo M, Zanuy S. Photoperiod modulates the reproductive axis of European sea bass through regulation of kiss1 and gnrh2 neuronal expression. Gen Comp Endocrinol. 2017;240:35–45. https://doi.org/10.1016/j.ygcen.2016.09.007.

    CAS  Article  PubMed  Google Scholar 

  55. 55.

    Ko CH, Takahashi JS. Molecular components of the mammalian circadian clock. Hum Mol Genet. 2006;15(2):R271–7. https://doi.org/10.1093/hmg/ddl207.

    CAS  Article  PubMed  Google Scholar 

  56. 56.

    Bowden TJ, Thompson KD, Morgan AL, Gratacap RM, Nikoskelainen S. Seasonal variation and the immune response: a fish perspective. Fish Shellfish Immunol. 2007;22(6):695–706. https://doi.org/10.1016/j.fsi.2006.08.016.

    Article  PubMed  Google Scholar 

  57. 57.

    Chou MY, Hung JC, Wu LC, Hwang SP, Hwang PP. Isotocin controls ion regulation through regulating ionocyte progenitor differentiation and proliferation. Cell Mol Life Sci. 2011;68(16):2797–809. https://doi.org/10.1007/s00018-010-0593-2.

    CAS  Article  PubMed  Google Scholar 

  58. 58.

    Xi D, Chi H, Gainer H. Analysis of synaptotagmin I–IV messenger RNA expression and developmental regulation in the rat hypothalamus and pituitary. Neuroscience. 1999;88(2):425–35. https://doi.org/10.1016/S0306-4522(98)00234-6.

    CAS  Article  PubMed  Google Scholar 

  59. 59.

    Rao PP. Studies on the structural variations in the brain of teleosts and their significance. Cells Tissues Organs. 1967;68(3):379–99. https://doi.org/10.1159/000143042.

    Article  Google Scholar 

  60. 60.

    Vatine G, Vallone D, Gothilf Y, Foulkes NS. It's time to swim! Zebrafish and the circadian clock. FEBS Lett. 2011;585(10):1485–94. https://doi.org/10.1016/j.febslet.2011.04.007.

    CAS  Article  PubMed  Google Scholar 

  61. 61.

    Paibomesai MI, Moghadam HK, Ferguson MM, Danzmann RG. Clock genes and their genomic distributions in three species of salmonid fishes: associations with genes regulating sexual maturation and cell cycling. BMC Res Notes. 2010;3(1):215. https://doi.org/10.1186/1756-0500-3-215.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  62. 62.

    Reaume CJ, Sokolowski MB. Conservation of gene function in behaviour. Philosoph Tran Royal Soc B: Biol Sci. 2011;366(1574):2100–10. https://doi.org/10.1098/rstb.2011.0028.

    Article  Google Scholar 

  63. 63.

    Lazado CC, Kumaratunga HP, Nagasawa K, Babiak I, Giannetto A, Fernandes JM. Daily rhythmicity of clock gene transcripts in Atlantic cod fast skeletal muscle. PLoS One. 2014;9(6):e99172. https://doi.org/10.1371/journal.pone.0099172.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  64. 64.

    Callaway E. A biological clock to wind them all. Nat News (Nature Publishing Group). 2012. https://doi.org/10.1038/nature.2012.10654.

  65. 65.

    Kumar V, Wingfield JC, Dawson A, Ramenofsky M, Rani S, Bartell P. Biological clocks and regulation of seasonal reproduction and migration in birds. Physiol Biochem Zool. 2010;83(5):827–35. https://doi.org/10.1086/652243.

    Article  PubMed  Google Scholar 

  66. 66.

    Williams Sarah CP. Hacking the biological clock. Life time. 2015. Stanford Med.

  67. 67.

    Mary AB: Circannual Rhythmicity in Invertebrates In Circannual Clocks: Annual Biological Rhythms Ed Pengelley E T, 1974, Pages 11–53, Academic Press.

  68. 68.

    Chemineau P, Malpaux B, Brillard JP, Fostier A. Seasonality of reproduction and production in farm fishes, birds and mammals. Animal. 2007;1(3):419–32. https://doi.org/10.1017/S1751731107691873.

    CAS  Article  PubMed  Google Scholar 

  69. 69.

    Nishiwaki-Ohkawa T, Yoshimura T. Molecular basis for regulating seasonal reproduction in vertebrates. J Endocrinol. 2016;229(3):R117–27. https://doi.org/10.1530/JOE-16-0066.

    CAS  Article  PubMed  Google Scholar 

  70. 70.

    Pittendrigh CS: On temporal organization in living systems. In: Harvey lectures, 1961, p. 93-125. Academic press, Inc., N Y.

  71. 71.

    Foulkes NS, Whitmore D, Vallone D, Bertolucci C. Studying the evolution of the vertebrate circadian clock: the power of fish as comparative models. Adv Genet. 2016;95:1–30. https://doi.org/10.1016/bs.adgen.2016.05.002.

    CAS  Article  PubMed  Google Scholar 

  72. 72.

    Qian X, Ba Y, Zhuang Q, Zhong G. RNA-Seq technology and its application in fish transcriptomics. Omics: J Int Biol. 2014;18(2):98–110. https://doi.org/10.1089/omi.2013.0110.

    CAS  Article  Google Scholar 

  73. 73.

    Robinson N, Sahoo PK, Baranski M, Mahapatra KD, Saha JN, Das S, et al. Expressed sequences and polymorphisms in rohu carp (Labeo rohita, Hamilton) revealed by mRNA-seq. Mar Biotechnol. 2012;14(5):620–33. https://doi.org/10.1007/s10126-012-9433-8.

  74. 74.

    Robinson N, Baranski M, Mahapatra KD, Saha JN, Das S, Mishra J, et al. A linkage map of transcribed single nucleotide polymorphisms in rohu (Labeo rohita) and QTL associated with resistance to Aeromonas hydrophila. BMC Genomics. 2014;15(1):541. https://doi.org/10.1186/1471-2164-15-541.

  75. 75.

    Skiba-Cassy S, Lansard M, Panserat S, Médale F. Rainbow trout genetically selected for greater muscle fat content display increased activation of liver TOR signaling and lipogenic gene expression. Am J Phys Regul Integr Comp Phys. 2009;297(5):R1421–9. https://doi.org/10.1152/ajpregu.00312.2009.

    CAS  Article  Google Scholar 

  76. 76.

    Hu G, Gu W, Sun P, Bai Q. Wang: transcriptome analyses reveal lipid metabolic process in liver related to the difference of carcass fat content in rainbow trout (Oncorhynchus mykiss). Int J genomics. 2016;7281585.

  77. 77.

    Asker N, Kristiansson E, Albertsson E, Larsson DJ, Förlin L. Hepatic transcriptome profiling indicates differential mRNA expression of apoptosis and immune related genes in eelpout (Zoarces viviparus) caught at Göteborg harbor. Sweden Aquatic Toxicol. 2013;130:58–67. https://doi.org/10.1016/j.aquatox.2012.12.017.

    CAS  Article  Google Scholar 

  78. 78.

    Barber BJ, Blake NJ. Reproductive physiology. In: Shumway SE, Parsons GJ, editors. Scallops: biology, ecology, and aquaculture. Amsterdam: Elsevier; 2006. p. 357–416. https://doi.org/10.1016/S0167-9309(06)80033-5.

    Chapter  Google Scholar 

  79. 79.

    Luna LG: Manual of histologic staining methods of the Armed Forces Institute of Pathology. 1968, 1–258 (New York: McGraw-Hill).

  80. 80.

    Andrews S: FastQC: a quality control tool for high throughput sequence data. 2010.

  81. 81.

    Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. https://doi.org/10.1093/bioinformatics/btu170.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  82. 82.

    Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-Seq: reference generation and analysis with trinity. Nat Protoc. 2013;8(8):1494–512. https://doi.org/10.1038/nprot.2013.084.

  83. 83.

    Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC bioinformatics. 2011;12(1):323. https://doi.org/10.1186/1471-2105-12-323.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  84. 84.

    Robinson MD, McCarthy DJ. Smyth GK: edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40. https://doi.org/10.1093/bioinformatics/btp616.

    CAS  Article  Google Scholar 

  85. 85.

    Zhang HM, Liu T, Liu CJ, Song S, Zhang X, Liu W, et al. AnimalTFDB 2.0: a resource for expression, prediction and functional study of animal transcription factors. Nucleic Acids Res. 2014;43(D1):D76–81. https://doi.org/10.1093/nar/gku887.

  86. 86.

    Thiel T, Michalek W, Varshney R, Graner A: Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.). Theor Appl Gene. 2003, 106:411–422.

  87. 87.

    Li H. Durbin: fast and accurate short read alignment with burrows–wheeler transform. Bioinformatics. 2009;25(14):1754–60. https://doi.org/10.1093/bioinformatics/btp324.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  88. 88.

    Li H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics. 2011;27(21):2987–93. https://doi.org/10.1093/bioinformatics/btr509.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  89. 89.

    Griffiths-Jones S, Enright AJ, Grocock RJ, Van Dongen S, Bateman A. miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006;34(90001):D140–4. https://doi.org/10.1093/nar/gkj112.

    CAS  Article  PubMed  Google Scholar 

  90. 90.

    Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS. MicroRNA targets in Drosophila. Genome Biol. 2003;5(1):R1. https://doi.org/10.1186/gb-2003-5-1-r1.

    Article  PubMed  PubMed Central  Google Scholar 

  91. 91.

    Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504. https://doi.org/10.1101/gr.1239303.

  92. 92.

    Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3—new capabilities and interfaces. Nucleic Acids Res. 2012;40(15):e115–5. https://doi.org/10.1093/nar/gks596.

  93. 93.

    Pfaffl MW. A new mathematical model for relative quantification in real-time RT–PCR. Nucleic Acids Res. 2001;29(9):e45–445. https://doi.org/10.1093/nar/29.9.e45.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  94. 94.

    Bennett V, Lorenzo DN. Spectrin-and ankyrin-based membrane domains and the evolution of vertebrates. Curr Top Membr. 2013;72:1–37. https://doi.org/10.1016/B978-0-12-417027-8.00001-5.

    CAS  Article  PubMed  Google Scholar 

  95. 95.

    Ferrari R, Forabosco P, Vandrovcova J, Botía JA, Guelfi S, Warren JD, et al. Frontotemporal dementia: insights into the biological underpinnings of disease through gene co-expression network analysis. Mol Neurodegener. 2016;11(1):21. https://doi.org/10.1186/s13024-016-0085-4.

  96. 96.

    Herrero-Turrión MJ, Rodríguez-Martín I, López-Bellido R, Rodríguez RE. Whole-genome expression profile in zebrafish embryos after chronic exposure to morphine: identification of new genes associated with neuronal function and mu opioid receptor expression. BMC Gen. 2014;15(1):874. https://doi.org/10.1186/1471-2164-15-874.

    CAS  Article  Google Scholar 

  97. 97.

    Greenwood AK, Peichel CL. Social regulation of gene expression in threespine sticklebacks. PLoS One. 2015;10(9):e0137726. https://doi.org/10.1371/journal.pone.0137726.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  98. 98.

    Nonet ML. A window into domain amplification through piccolo in teleost fish. G3: genes, genomes. Genetics. 2012;2(11):1325–39. https://doi.org/10.1534/g3.112.003624.

    CAS  Article  Google Scholar 

  99. 99.

    Olave I, Wang W, Xue Y, Kuo A, Crabtree GR. Identification of a polymorphic, neuron-specific chromatin remodeling complex. Genes Dev. 2002;16(19):2509–17. https://doi.org/10.1101/gad.992102.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  100. 100.

    Heimel JA, Hermans JM, Sommeijer JP, Levelt CN. Genetic control of experience-dependent plasticity in the visual cortex. Genes Brain Behav. 2008;7(8):915–23. https://doi.org/10.1111/j.1601-183X.2008.00431.x.

    CAS  Article  PubMed  Google Scholar 

  101. 101.

    Teng H, Cai W, Zhou L, Zhang J, Liu Q, Wang Y, et al. Evolutionary mode and functional divergence of vertebrate NMDA receptor subunit 2 genes. PLoS One. 2010;5(10):e13342. https://doi.org/10.1371/journal.pone.0013342.

  102. 102.

    Drew RE, Rodnick KJ, Settles M, Wacyk J, Churchill E, Powell MS, et al. Effect of starvation on transcriptomes of brain and liver in adult female zebrafish (Danio rerio). Physiol Genomics. 2008;35(3):283–95. https://doi.org/10.1152/physiolgenomics.90213.2008.

  103. 103.

    Guillot R, Cortés R, Navarro S, Mischitelli M, García-Herranz V, Sánchez E, et al. Behind melanocortin antagonist overexpression in the zebrafish brain: a behavioral and transcriptomic approach. Hormones Behav. 2016;82:87–100. https://doi.org/10.1016/j.yhbeh.2016.04.011.

  104. 104.

    Boothroyd CE, Wijnen H, Naef F, Saez L, Young MW. Integration of light and temperature in the regulation of circadian gene expression in Drosophila. PLoS Genet. 2007;3(4):e54. https://doi.org/10.1371/journal.pgen.0030054.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  105. 105.

    Ciampi R, Nikiforov YE. Alterations of the BRAF gene in thyroid tumors. Endocr Pathol. 2005;16(3):163–71. https://doi.org/10.1385/EP:16:3:163.

    CAS  Article  PubMed  Google Scholar 

  106. 106.

    Li Z, Zheng W, Wang Z, Zeng Z, Zhan H, Li C, et al. A transgenic zebrafish liver tumor model with inducible Myc expression reveals conserved Myc signatures with mammalian liver tumors. Dis Model Mech. 2013;6(2):414–23. https://doi.org/10.1242/dmm.010462.

  107. 107.

    Ahlemann M, Zeidler R, Lang S, Mack B, Münz M, Gires O. Carcinoma-associated eIF3i overexpression facilitates mTOR-dependent growth transformation. Mol Carcinog. 2006;45(12):957–67. https://doi.org/10.1002/mc.20269.

    CAS  Article  PubMed  Google Scholar 

  108. 108.

    Choudhuri A, Evans T, Maitra U. Non-core subunit eIF3h of translation initiation factor eIF3 regulates zebrafish embryonic development. Dev Dyn. 2010;239(6):1632–44. https://doi.org/10.1002/dvdy.22289.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  109. 109.

    Wiseman SB, He Y, Gamal-El Din M, Martin JW, Jones PD, Hecker M, et al. Transcriptional responses of male fathead minnows exposed to oil sands process-affected water. Comp Biochem Physiol Part C: Toxicol Pharmacol. 2013;157(2):227–35. https://doi.org/10.1016/j.cbpc.2012.12.002.

  110. 110.

    Fujita M, Mitsuhashi H, Isogai S, Nakata T, Kawakami A, Nonaka I, et al. Filamin C plays an essential role in the maintenance of the structural integrity of cardiac and skeletal muscles, revealed by the medaka mutant zacro. Dev Biol. 2012;361(1):79–89. https://doi.org/10.1016/j.ydbio.2011.10.008.

  111. 111.

    Shiraishi E, Hosseini H, Kang DK, Kitano T, Akiyama H. Nanosecond pulsed electric field suppresses development of eyes and germ cells through blocking synthesis of retinoic acid in Medaka (Oryzias latipes). PLoS One. 2013;8(8):e70670. https://doi.org/10.1371/journal.pone.0070670.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  112. 112.

    Adolfi MC, Herpin A, Regensburger M, Sacquegno J, Waxman JS, Schartl M. Retinoic acid and meiosis induction in adult versus embryonic gonads of medaka. Sci Rep. 2016;6(1):34281. https://doi.org/10.1038/srep34281.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  113. 113.

    Cao LX, Le Bousse-Kerdilès MC, Clay D, Oshevski S, Jasmin C, Krief P. Implication of a new molecule IK in CD34+ hematopoietic progenitor cell proliferation and differentiation. Blood. 1997;89(10):3615–23. https://doi.org/10.1182/blood.V89.10.3615.

    CAS  Article  PubMed  Google Scholar 

  114. 114.

    Vacaru AM, Unlu G, Spitzner M, Mione M, Knapik EW, Sadler KC. In vivo cell biology in zebrafish–providing insights into vertebrate development and disease. J Cell Sci. 2014;127(3):485–95. https://doi.org/10.1242/jcs.140194.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  115. 115.

    Gu M, Lin G, Lai Q, Zhong B, Liu Y, Mi Y, et al. Ctenopharyngodon idella IRF2 plays an antagonistic role to IRF1 in transcriptional regulation of IFN and ISG genes. Dev Comp Immunol. 2015;49(1):103–12. https://doi.org/10.1016/j.dci.2014.11.014.

  116. 116.

    Klangnurak W, Tokumoto T. Fine selection of up-regulated genes during ovulation by in vivo induction of oocyte maturation and ovulation in zebrafish. Zool Lett. 2017;3(1):2. https://doi.org/10.1186/s40851-017-0065-8.

    Article  Google Scholar 

  117. 117.

    Ribas LP, Planas JV, Barton B, Monetti C, Bernadini G, Saroglia M, et al. A differentially expressed enolase gene isolated from the gilthead sea bream (Sparus aurata) under high-density conditions is up-regulated in brain after in vivo lipopolysaccharide challenge. Aquaculture. 2004;241(1–4):195–206. https://doi.org/10.1016/j.aquaculture.2004.07.015.

  118. 118.

    Fort P, Blangy A. The evolutionary landscape of Dbl-like RhoGEF families: adapting eukaryotic cells to environmental signals. Genome Biol Evol. 2017;9(6):1471–86. https://doi.org/10.1093/gbe/evx100.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  119. 119.

    Barat A, Kumar R, Goel C, Singh AK, Sahoo PK. De novo assembly and characterization of tissue-specific transcriptome in the endangered golden mahseer, Tor putitora. Meta Gene. 2016;7:28–33. https://doi.org/10.1016/j.mgene.2015.11.001.

    Article  PubMed  Google Scholar 

  120. 120.

    Birk OS, Casiano DE, Wassif CA, Cogliati T. The LIM homeobox gene Lhx9 is essential for mouse gonad formation. Nature. 2000;403(6772):909–13. https://doi.org/10.1038/35002622.

    CAS  Article  PubMed  Google Scholar 

  121. 121.

    Knoll-Gellida A, André M, Gattegno T, Forgue J, Admon A, Babin PJ. Molecular phenotype of zebrafish ovarian follicle by serial analysis of gene expression and proteomic profiling, and comparison with the transcriptomes of other animals. BMC Genomics. 2006;7(1):46. https://doi.org/10.1186/1471-2164-7-46.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  122. 122.

    Rhee JS, Kim BM, Kim RO, Choi BS, Choi IY, Lee YM, et al. Analysis of expressed sequence tags from the liver and ovary of the euryhaline hermaphroditic fish, Kryptolebias marmoratus. Comp Biochem Physiol Part D: Genom Proteomics. 2011;6(3):244–55. https://doi.org/10.1016/j.cbd.2011.05.001.

  123. 123.

    Gouttenoire J, Valcourt U, Bougault C, Aubert-Foucher E, Arnaud E, Giraud L, et al. Knockdown of the intraflagellar transport protein IFT46 stimulates selective gene expression in mouse chondrocytes and affects early development in zebrafish. J Biol Chem. 2007;282(42):30960–73. https://doi.org/10.1074/jbc.M705730200.

  124. 124.

    Baumgart M, Groth M, Priebe S, Savino A, Testa G, Dix A, et al. RNA-seq of the aging brain in the short-lived fish N. furzeri - conserved pathways and novel genes associated with neurogenesis. Aging Cell. 2014;13(6):965–74. https://doi.org/10.1111/acel.12257.

  125. 125.

    Stacey NE, Pandey S. Effects of indomethacin and prostaglandins on ovulation of goldfish. Prostaglandins. 1975;9(4):597–607. https://doi.org/10.1016/0090-6980(75)90065-9.

    CAS  Article  PubMed  Google Scholar 

  126. 126.

    Touhata K, Tanaka M, Toyohara H, Tanaka H, Sakaguchi M. Seasonal change in collagen content of red seabream muscle. Fish Sci. 2000;66(3):553–7. https://doi.org/10.1046/j.1444-2906.2000.00073.x.

    CAS  Article  Google Scholar 

  127. 127.

    Palstra AP, Planas JV. Fish under exercise. Fish Physiol Biochem. 2011;37(2):259–72. https://doi.org/10.1007/s10695-011-9505-0.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  128. 128.

    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. https://doi.org/10.1016/j.ygcen.2015.01.025.

  129. 129.

    Bayne S, Li H, Jones ME, Pinto AR, van Sinderen M, Drummond A, et al. Estrogen deficiency reversibly induces telomere shortening in mouse granulosa cells and ovarian aging in vivo. Protein Cell. 2011;2(4):333–46. https://doi.org/10.1007/s13238-011-1033-2.

  130. 130.

    Philip AM, Vijayan MM. Stress-immune-growth interactions: cortisol modulates suppressors of cytokine signaling and JAK/STAT pathway in rainbow trout liver. PLoS One. 2015;10(6):e0129299. https://doi.org/10.1371/journal.pone.0129299.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  131. 131.

    Fedorov VB, Goropashnaya AV, Tøien Ø, Stewart NC, Chang C, Wang H, et al. Modulation of gene expression in heart and liver of hibernating black bears (Ursus americanus). BMC Genomics. 2011;12(1):171. https://doi.org/10.1186/1471-2164-12-171.

  132. 132.

    Komoroske LM, Connon RE, Jeffries KM, Fangue NA. Linking transcriptional responses to organismal tolerance reveals mechanisms of thermal sensitivity in a mesothermal endangered fish. Mol Ecol. 2015;24(19):4960–81. https://doi.org/10.1111/mec.13373.

    CAS  Article  PubMed  Google Scholar 

  133. 133.

    Tingaud-Sequeira A, Forgue J, André M, Babin PJ. Epidermal transient down-regulation of retinol-binding protein 4 and mirror expression of apolipoprotein Eb and estrogen receptor 2a during zebrafish fin and scale development. Dev Dyn. 2006;235(11):3071–9. https://doi.org/10.1002/dvdy.20921.

    CAS  Article  PubMed  Google Scholar 

  134. 134.

    Yang D, Liu Q, Yang M, Wu H, Wang Q, Xiao J, et al. RNA-seq liver transcriptome analysis reveals an activated MHC-I pathway and an inhibited MHC-II pathway at the early stage of vaccine immunization in zebrafish. BMC Genomics. 2012;13(1):319. https://doi.org/10.1186/1471-2164-13-319.

  135. 135.

    Jackson B, Brocker C, Thompson DC, Black W, Vasiliou K, Nebert DW, et al. Update on the aldehyde dehydrogenase gene (ALDH) superfamily. Human Genomics. 2011;5(4):283–303. https://doi.org/10.1186/1479-7364-5-4-283.

  136. 136.

    Groh KJ, Nesatyy VJ, Segner H, Eggen RI, Suter MJF. Global proteomics analysis of testis and ovary in adult zebrafish (Danio rerio). Fish Physiol Biochem. 2011;37(3):619–47. https://doi.org/10.1007/s10695-010-9464-x.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  137. 137.

    Um HN, Han JM, Hwang JI, Hong SI, Vaudry H, Seong JY. Molecular coevolution of kisspeptins and their receptors from fish to mammals. Ann N Y Acad Sci. 2012;1200(1):67–74. https://doi.org/10.1111/j.1749-6632.2010.05508.x.

    CAS  Article  Google Scholar 

  138. 138.

    Chen YJ, Zhang TY, Chen HY, Lin SM, Luo L, Wang DS. Simultaneous stimulation of glycolysis and gluconeogenesis by feeding in the anterior intestine of the omnivorous GIFT tilapia, Oreochromis niloticus. Biol Open. 2017;6(6):818–24. https://doi.org/10.1242/bio.024836.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  139. 139.

    Liu C, Yao J, Wei Q, Xing G, Cao X. Spatial and temporal expression patterns of Osbpl2a and Osbpl2b during zebrafish embryonic development. Int J Pediatr Otorhinolaryngol. 2016;84:174–9. https://doi.org/10.1016/j.ijporl.2016.02.030.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by ICAR-CIFA, Bhubaneshwar (Institute Project Code: I-59(n), P1-2012/2-FA-L-10/8110) and CABin grant of ICAR-IASRI (F. no. Agril. Edn.4-1/2013-A&P). Authors are thankful to Director, ICAR-CIFA, Bhubaneshwar for sanctioning the project to SN as Principal Investigator for all financial support in genomic data generation. Authors are also thankful to Indian Council of Agricultural Research, Ministry of Agriculture and Farmers’ Welfare, Govt. of India for providing financial assistance in the form of CABin grant as well as Advanced Super Computing Hub for Omics Knowledge in Agriculture (ASHOKA) facility at ICAR-IASRI, New Delhi, India created under National Agricultural Innovation Project, funded by World Bank.

Funding

This work was supported by ICAR-CIFA, Bhubaneshwar (Institute Project Code: I-59(n), P1–2012/2-FA-L-10/8110) and CABin grant of ICAR-IASRI (F. no. Agril. Edn.4–1/2013-A&P). Authors are thankful to Director, ICAR-CIFA, Bhubaneshwar for sanctioning the project to SN as Principal Investigator for all financial support in genomic data generation.

Author information

Affiliations

Authors

Contributions

SN, DK and JKS conceived theme of the study. SN designed the experiment and generated the genomic data along with SP, GM, UKU, DKS, PKM, DKV, PR and PD. SJ, MAI, RSJ and UBA did computational analysis and developed genomic resource. SJ, MAI, RSJ, DK and SN drafted the manuscript. SN, JKS, PJ, DK and AR edited the manuscript. All co-authors read and approved the final manuscript.

Corresponding authors

Correspondence to Samiran Nandi or Dinesh Kumar.

Ethics declarations

Ethics statement

As per Chapter 4, Section 15(1) of the Prevention of Cruelty to Animals (PCA) Act 1960 of Govt of India, using guidelines recommended by CPCSEA (Committee for the Purpose of Controlling and Supervising Experiments in Animals), Institute Animal Ethics Committee of ICAR-CIFA approved the experiments on 1st March 2017.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Jaiswal, S., Nandi, S., Iquebal, M.A. et al. Revelation of candidate genes and molecular mechanism of reproductive seasonality in female rohu (Labeo rohita Ham.) by RNA sequencing. BMC Genomics 22, 685 (2021). https://doi.org/10.1186/s12864-021-08001-6

Download citation

Keywords

  • Labeo rohita
  • Circannual
  • Cyprinid
  • Biological clock genes
  • Marker
  • Reproductive seasonality