Normalizing RT-qPCR data: are we getting the right answers? An appraisal of normalization approaches and internal reference genes from a case study in the finfish Lates calcarifer

Deciphering the impact of stickwater hydrolysate on growth performance, immune response, and IGF-1/PI3K/AKT/mTOR signaling pathway in Siberian sturgeon (Acipenser baerii) fingerlings

A feeding trial lasting 56 days was carried out to assess how the inclusion of stickwater hydrolysate (SWH) in the diet of Siberian sturgeon (Acipenser baerii) fingerlings affected their growth performance, immunity, digestive enzyme activity, and gene expression linked to the IGF-1/PI3K/AKT/mTOR signaling pathway. Siberian sturgeon fingerlings were acclimatized and fed isonitrogenous, isoenergetic diets with varying SWH concentrations (0%, 0.5%, 1.5%, and 2.5%). Growth parameters, serum proteins, immunological and digestive enzyme activities, and gene expression levels were assessed post-trial. Results demonstrated that 0.5%, and 1.5% SWH treatments significantly improved weight gain, specific growth rate, feed conversion ratio, and protein efficiency ratio. Notably, these diets also elevated serum protein and plasma globulin levels, reduced albumin-to-globulin ratios, and enhanced lysozyme, myeloperoxidase (MPO) activities, and immunoglobulin (Ig) M levels, indicating an immunostimulatory effect. Digestive enzyme activities were markedly increased in the SWH groups, particularly at 1.5%. Gene expression analyses revealed upregulation of mtorc1, s6K, akt, pi3k, and igf1, with concurrent downregulation of 4e-bp1 in the muscle of fish, signifying activation of the IGF-1/PI3K/AKT/mTOR pathway, which is central to protein synthesis and muscle growth. In conclusion, SWH at appropriate levels significantly enhances growth, digestive efficiency, and immune function in Siberian sturgeon fingerlings, while also activating key metabolic pathways.

Harnessing paraprobiotics and postbiotics for enhanced immune function in Asian seabass (Lates calcarifer): Insights into pattern recognition receptor signaling

The Asian seabass, Lates calcarifer, is a key species in Asian aquaculture due to its nutritional value and adaptability. However, disease outbreaks, particularly viral and bacterial infections, pose significant challenges to its production. Immunostimulants offer promising solutions but raise safety concerns. Paraprobiotics and postbiotics (CPP) emerge as safer alternatives, exerting health benefits without live microorganisms. This study investigated the potential of probiotic paraprobiotic and postbiotic supplements derived from Bacillus subtilis to enhance the immune response and antioxidant capacity of Asian seabass and improve their resistance to Streptococcus iniae infection. Analysis of antioxidant activity and lipid peroxidation revealed significant improvements in fish supplemented with CPP, indicating their effectiveness in mitigating oxidative stress. Immunological assays demonstrated enhanced growth performance and serum immunity, including increased alternative complement activity, immunoglobulin levels, and phagocytic activity, in supplemented fish. Furthermore, upregulated expression of proinflammatory cytokines (TNF-α, IL-6, IL-1β) and pattern recognition receptors (NLRC3, TLR22, MDA5) in immune tissues. Fish supplemented with CPP exhibited higher resistance and survival rates against S. iniae infection challenge compared to control groups. The study elucidates the mechanisms underlying the immunomodulatory effects of CPP, shedding light on their potential applications in aquaculture.

Paralog switching facilitates diadromy: ontogenetic, microevolutionary and macroevolutionary evidence

Identifying how the demands of migration are met at the level of gene expression is critical for understanding migratory physiology and can potentially reveal how migratory forms evolve from nonmigratory forms and vice versa. Among fishes, migration between freshwater and seawater (diadromy) requires considerable osmoregulatory adjustments, powered by the ion pump Na+, K+-ATPase (NKA) in the gills. Paralogs of the catalytic α-subunit of the pump (NKA α1a and α1b) are reciprocally upregulated in fresh- and seawater, a response known as paralog-switching, in gills of some diadromous species. We tested ontogenetic changes in NKA α-subunit paralog expression patterns, comparing pre-migrant and migrant alewife (Alosa pseudoharengus) sampled in their natal freshwater environment and after 24 h in seawater. In comparison to pre-migrants, juvenile out-migrants exhibited stronger paralog switching via greater downregulation of NKA α1a in seawater. We also tested microevolutionary changes in the response, exposing juvenile diadromous and landlocked alewife to freshwater (0 ppt) and seawater (30 ppt) for 2, 5, and 15 days. Diadromous and landlocked alewife exhibited salinity-dependent paralog switching, but levels of NKA α1b transcription were higher and the decrease in NKA α1a was greater after seawater exposure in diadromous alewife. Finally, we placed alewife α-subunit NKA paralogs in a macroevolutionary context. Molecular phylogenies show alewife paralogs originated independently of paralogs in salmonids and other teleosts. This study demonstrated that NKA paralog switching is tied to halohabitat profile and that duplications of the NKA gene provided the substrate for multiple, independent molecular solutions that support a diadromous life history.

Effects of atrazine on movement, metabolism and gene expression in Pelophylax nigromaculatus larvae under global warming

Global warming and environmental pollutants both pose a threat to the behavior and physiology of animals, but research on the combined effects of the two is limited. Atrazine, a widely used herbicide, has toxic effects on organisms. In this study, the effects of environmental concentrations of atrazine exposure (100 μg/L) for seven days on the movement, metabolism and gene expression related to motility of Pelophylax nigromaculatus larvae (GS8) were investigated under global warming. The results showed that compared to the optimal growth temperature (18 °C), atrazine treatment under global warming (21 °C) significantly increased the average speed (about 11.2 times) and maximum acceleration (about 1.98 times) of P. nigromaculatus larvae, altered the relative abundance of 539 metabolites, including Formyl-5-hydroxykynurenamine, 2,4-Dihydroxybenzophenone, and FAPy-adenine, and changed the nucleotide metabolism, pyrimidine metabolism, glycerophospholipid metabolism, and purine metabolism, as well as increased the gene expression of SPLA2 (about 6.46 times) and CHK (about 3.25 times). In summary, atrazine treatment under global warming caused metabolic disorders in amphibian larvae and increased the expression of some movement-related genes in the brain, resulting in abnormally active.

Characterization of repetitive DNA on the genome of the marsh rat Holochilus nanus (Cricetidae: Sigmodontinae)

Repetitive DNA are sequences repeated hundreds or thousands of times and an abundant part of eukaryotic genomes. SatDNA represents the majority of the repetitive sequences, followed by transposable elements. The species Holochilus nanus (HNA) belongs to the rodent tribe Oryzomyini, the most taxonomically diverse of Sigmodontinae subfamily. Cytogenetic studies on Oryzomyini reflect such diversity by revealing an exceptional range of karyotype variability. However, little is known about the repetitive DNA content and its involvement in chromosomal diversification of these species. In the search for a more detailed understanding about the composition of repetitive DNA on the genome of HNA and other species of Oryzomyini, we employed a combination of bioinformatic, cytogenetic and molecular techniques to characterize the repetitive DNA content of these species. RepeatExplorer analysis showed that almost half of repetitive content of HNA genome are composed by Long Terminal Repeats and a less significant portion are composed by Short Interspersed Nuclear Elements and Long Interspersed Nuclear Elements. RepeatMasker showed that more than 30% of HNA genome are composed by repetitive sequences, with two main waves of repetitive element insertion. It was also possible to identify a satellite DNA sequence present in the centromeric region of Oryzomyini species, and a repetitive sequence enriched on the long arm of HNA X chromosome. Also, comparative analysis between HNA genome with and without B chromosome did not evidence any repeat element enriched on the supernumerary, suggesting that B chromosome of HNA is composed by a fraction of repeats from all the genome.

Freshwater transfer affected intestinal microbiota with correlation to cytokine gene expression in Asian sea bass

As a catadromous fish, Asian sea bass (Lates calcarifer) juveniles migrate from seawater (SW) to freshwater (FW) for growth and development. During migration, they undergo physiological changes to acclimate to environmental salinity. Thus, it is crucial to understand how SW-to-FW migration affects the gut microbiota of catadromous fish. To the best of our knowledge, no study has revealed the effects of transfer to hypotonic environments on a catadromous fish microbiota. In this study, we aimed to determine the effects of FW transfer on the microbiota and cytokine gene expression in the intestines of juvenile catadromous Asian sea bass. The relationship between the water and the gut microbiota of this euryhaline species was also examined. We found that FW transfer affected both mucosa- and digesta-associated microbiota of Asian sea bass. Plesiomonas and Cetobacterium were dominant in both the mucosa- and digesta-associated microbiota of FW-acclimated sea bass. The pathogenic genera Vibrio, Staphylococcus, and Acinetobacter were dominant in the SW group. Although dominant fish microbes were present in the water, fish had their own unique microbes. Vitamin B6 metabolism was highly expressed in the FW fish microbiota, whereas arginine, proline, and lipid metabolism were highly expressed in the SW fish microbiota. Additionally, the correlation between cytokine gene expression and microbiota was found to be affected by FW transfer. Taken together, our results demonstrated that FW transfer altered the composition and functions of mucosa- and digesta-associated microbiota of catadromous Asian sea bass intestines, which correlated with cytokine gene expression.

First characterization of PIWI-interacting RNA clusters in a cichlid fish with a B chromosome

BACKGROUND

B chromosomes are extra elements found in several eukaryote species. Usually, they do not express a phenotype in the host. However, advances in bioinformatics over the last decades have allowed us to describe several genes and molecular functions related to B chromosomes. These advances enable investigations of the relationship between the B chromosome and the host to understand how this element has been preserved in genomes. However, considering that transposable elements (TEs) are highly abundant in this supernumerary chromosome, there is a lack of knowledge concerning the dynamics of TE control in B-carrying cells. Thus, the present study characterized PIWI-interacting RNA (piRNA) clusters and pathways responsible for silencing the mobilization of TEs in gonads of the cichlid fish Astatotilapia latifasciata carrying the B chromosome.

RESULTS

Through small RNA-seq and genome assembly, we predicted and annotated piRNA clusters in the A. latifasciata genome for the first time. We observed that these clusters had biased expression related to sex and the presence of the B chromosome. Furthermore, three piRNA clusters, named curupira, were identified in the B chromosome. Two of them were expressed exclusively in gonads of samples with the B chromosome. The composition of these curupira sequences was derived from LTR, LINE, and DNA elements, representing old and recent transposition events in the A. latifasciata genome and the B chromosome. The presence of the B chromosome also affected the expression of piRNA pathway genes. The mitochondrial cardiolipin hydrolase-like (pld6) gene is present in the B chromosome, as previously reported, and an increase in its expression was detected in gonads with the B chromosome.

CONCLUSIONS

Due to the high abundance of TEs in the B chromosome, it was possible to investigate the origin of piRNA from these jumping genes. We hypothesize that the B chromosome has evolved its own genomic guardians to prevent uncontrolled TE mobilization. Furthermore, we also detected an expression bias in the presence of the B chromosome over A. latifasciata piRNA clusters and pathway genes.

The mechanism of Megalobrama amblycephala muscle injury repair based on RNA-seq

Skeletal muscle myogenesis and injury-induced muscle regeneration contribute to muscle formation. Skeletal muscle stem cells, termed satellite cells (SCs), proliferate to repair injured muscle. To identify the molecular mechanism of regeneration after muscle injury as well as the genes related to muscle development in fish, in this study, the immunohistochemistry and the high-throughput RNA sequencing (RNA-seq) analysis were performed after Megalobrama amblycephala muscle was injured by needle stab. The results showed that paired box7-positive (Pax7⁺) SCs increased, and peaked at 96 to 144 h-post injury (hpi). The 6729 differentially expressed genes (DEGs), including 2125 up-regulated and 4604 down-regulated genes were found. GO terms significantly enriched by DEGs contained intercellular connections, signaling transduction and enzyme activity. KEGG enrichment analysis showed that most of the pathways were related to immunity, metabolism and cells related molecules, including actin skeleton regulation, Epstein Barr virus infection and plaque adhesion. The WGCNA results revealed that actin cytoskeleton and lipid metabolism related genes probably played crucial roles during repair after muscle injury. Collectively, all these results will provide new insights into the molecular mechanisms underlying muscle injury repair of fish.

Screening of reference genes in tiger puffer (Takifugu rubripes) across tissues and under different nutritional conditions

The present study was aimed at screening suitable reference genes for quantitative real-time polymerase chain reaction (qRT-PCR) in tiger puffer (Takifugu rubripes), an important aquaculture species in Asia and also a good model species for lipid research. Specifically, this reference gene screening was targeted at standardization of gene expression in different tissues (liver, muscle, brain, intestine, heart, eye, skin, and spleen) or under different nutritional conditions (starvation and different dietary lipid levels). Eight candidate reference genes (ribosomal protein L19 and L13 (RPL19 and RPL13), elongation factor-1 alpha (EF1α), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), hypoxanthine guanine phosphoribosyl transferase1 (HPRT1), beta-2-Microglobulin (B2M), 18S ribosomal RNA (18SrRNA), and beta actin (ACTB)) were evaluated with four algorithms (geNorm, NormFinder, BestKeeper, and comparative ΔCt method). The results showed that different algorithms generated inconsistent results. Based on these findings, RPL19, EF1α, 18SrRNA, and RPL13 were relatively stable in different tissues of tiger puffer. During starvation conditions, ACTB/RPL19 was the best reference gene combination. Under different dietary lipid levels, ACTB/RPL13 was the most suitable reference gene combination. The present results will help researchers to obtain more accurate results in future qRT-PCR analysis in tiger puffer.

Influence of bamboo vinegar powder (BVP) enriched diet on antioxidant status, immunity level, and pro-anti-inflammatory cytokines modulation in Asian sea bass, Lates calcarifer (Bloch 1790) against Vibrio anguillarum

Effect of bamboo vinegar powder (BVP) on growth, immunity, disease resistance, and immune-related gene expressions in juvenile Asian sea bass (barramundi), Lates calcarifer against Vibrio anguillarum was investigated. V. anguillarum infected fish fed by 2g BVP kg-1 diet exhibited significant growth, albumin (AB), serum total protein (TP), and globulin (GB) levels. Similarly, enhanced lysozyme (LZM), phagocytic (PC), anti-protease (AP), respiratory burst (RB), bactericidal (BC) activities and total immunoglobulin (Ig) levels were observed in fish fed 2g BVP kg-1. In addition, use of 2g BVP kg-1 in diet modulated several immune related gene expressions against Vibrio disease. This study exhibited a positive effect of dietary 2g BVP kg-1 diet on hemato-biochemical and immunological responses, which enhance the immunocompetence and health status of L. calcarifer against V. anguillarum infection.

Identification of stable pollen development related reference genes for accurate qRT-PCR analysis and morphological variations in autotetraploid and diploid rice

Autotetraploid rice exhibited hybrid vigor and greater genetic variation compared to diploid rice, but low pollen fertility is a major hindrance for its utilization. Our previous analysis revealed that large number of pollen fertility genes were exhibited down-regulation in autotetraploid rice. Hence, it is of utmost importance to reveal the expression patterns of pollen fertility genes with high accuracy. To find stable reference genes for autotetraploid rice, we compared the pollen development stages between diploid and autotetraploid rice, and 14 candidate genes were selected based on transcriptome analysis to evaluate their expression levels. Autotetraploid rice (i.e. Taichung65-4x) displayed lower seed set (40.40%) and higher percentage of abnormalities during the pollen development process than its diploid counterpart. To detect the candidate reference genes for pollen development of autotetraploid and diploid rice, we used five different algorithms, including NormFinder, BestKeeper, ΔCt method, geNorm and Re-Finder to evaluate their expression patterns stability. Consequently, we identified two genes, Cytochrome b5 and CPI, as the best candidate reference genes for qRT-PCR normalization in autotetraploid and diploid rice during pre-meiosis, meiosis, single microspore and bicellular pollen development stages. However, Cytochrome b5 was found to be the most stably expressed gene during different pollen development stages in autotetraploid rice. The results of our study provide a platform for subsequent gene expression analyses in autotetraploid rice, which could also be used in other polyploid plants.

Differential expression of miRNAs in the presence of B chromosome in the cichlid fish Astatotilapia latifasciata

BACKGROUND

B chromosomes (Bs) are extra elements observed in diverse eukaryotes, including animals, plants and fungi. Although Bs were first identified a century ago and have been studied in hundreds of species, their biology is still enigmatic. Recent advances in omics and big data technologies are revolutionizing the B biology field. These advances allow analyses of DNA, RNA, proteins and the construction of interactive networks for understanding the B composition and behavior in the cell. Several genes have been detected on the B chromosomes, although the interaction of B sequences and the normal genome remains poorly understood.

RESULTS

We identified 727 miRNA precursors in the A. latifasciata genome, 66% which were novel predicted sequences that had not been identified before. We were able to report the A. latifasciata-specific miRNAs and common miRNAs identified in other fish species. For the samples carrying the B chromosome (B⁺), we identified 104 differentially expressed (DE) miRNAs that are down or upregulated compared to samples without B chromosome (B^-) (p < 0.05). These miRNAs share common targets in the brain, muscle and gonads. These targets were used to construct a protein-protein-miRNA network showing the high interaction between the targets of differentially expressed miRNAs in the B⁺ chromosome samples. Among the DE-miRNA targets there are protein-coding genes reported for the B chromosome that are present in the protein-protein-miRNA network. Additionally, Gene Ontology (GO) terms related to nuclear matrix organization and response to stimulus are exclusive to DE miRNA targets of B⁺ samples.

CONCLUSIONS

This study is the first to report the connection of B chromosomes and miRNAs in a vertebrate species. We observed that the B chromosome impacts the miRNAs expression in several tissues and these miRNAs target several mRNAs involved with important biological processes.

Adequate supply of dietary taurine stimulates expression of molecular markers of growth and protein turnover in juvenile barramundi (Lates calcarifer)

A trial was conducted to investigate the effect of dietary taurine (Tau) supply on the plasma amino acid composition and hepatic expression of several genes in juvenile barramundi (Lates calcarifer) after feeding. Triplicate tanks of fish (average weight, 89.3 g) were fed diets containing either a deficient (1 g kg^-1), adequate (8 g kg^-1) or excessive (19 g kg^-1) level of dietary Tau. Liver tissues collected before feeding, and at 2- and 4-h post-feeding, were analysed for expression of genes involved in pathways of sulphur amino acid turnover, Tau biosynthesis and transport, target of rapamycin (TOR) signalling, the somatotropic axis and protein turnover. The treatment had no significant effect on the profiles of any amino acid in plasma collected over time after feeding, other than Tau and glycine. The expression profile of cystine and Tau synthetic genes suggested an effect of Tau excess on the metabolism of cystine. Markers of two pathways of Tau biosynthesis appear to be active in this species, providing proof that this species possesses the ability to synthesise Tau from SAA precursors. A marker for the regulation of Tau transport and homeostasis was shown to be directly regulated by Tau availability, whilst a link between adequate supply of Tau and TOR pathway-mediated growth stimulation was also apparent. An observed depression in expression of genes of the somatotropic axis, coupled with upregulation of the proteolytic and TOR-suppressing genes, in response to excessive Tau supply in the diet, signalled that excessive Tau may not be conducive to optimal growth in this species.

Sex-specific dmrt1 and cyp19a1 methylation and alternative splicing in gonads of the protandrous hermaphrodite barramundi

Epigenetics is involved in sex differentiation of gonochoristic and hermaphroditic fish species, whereby two genes dmrt1 (pro-male) and cyp19a1 (pro-female) are known to play major roles. Barramundi, Lates calcarifer, is an important tropical aquaculture species that undergo natural and permanent male to female sex change, a process for which the exact underlying molecular mechanisms are still unknown. To elucidate whether DNA methylation is involved in sex control of barramundi, a next-generation bisulfite amplicon sequencing approach was used to target 146 CpG sites within proximal promoters and first exons of seven sex-related genes (dmrt1, cyp19a1, amh, foxl2, nr5a2, sox8 and sox9) of 24 testis and 18 ovaries of captive and wild adult barramundi. Moreover, comparative expression profiles of the key dmrt1 and cyp19a1 genes were further investigated using RT-qPCR and Sanger sequencing approaches, whereas expression levels of remaining targeted genes were based on available literature for the species. Results showed that cyp19a1 and amh were more methylated in males, whereas dmrt1 and nr5a2 were more methylated in females (P < 0.001), with no gender differences found for foxl2, sox8 or sox9 genes (P > 0.05). Sex-biased promoter DNA methylation was inversely related to gene expression only for dmrt1 and nr5a2, and directly related to amh expression, whereas no differences in cyp19a1 expression were found between testes and ovaries. Notably, unique sex-specific alternative splicing of dmrt1 and cyp19a1 were discovered, whereby males lacked the full-length aromatase coding cyp19a1 mRNA due to partial or total exon splicing, and females lacked the dmrt1 exon containing the DM-domain sequence. This study advances the current knowledge aiming to elucidate the genetic mechanisms within male and female gonads of this large protandrous hermaphrodite by providing the first evidence of epigenetics and alternative splicing simultaneously affecting key genes (cyp19a1 and dmrt1) central to sex differentiation pathways.

Characterization of lipid metabolism genes and the influence of fatty acid supplementation in the hepatic lipid metabolism of dusky grouper (Epinephelus marginatus)

Dusky grouper is an important commercial fish species in many countries, but some factors such as overfishing has significantly reduced their natural stocks. Aquaculture emerges as a unique way to conserve this species, but very little biological information is available, limiting the production of this endangered species. To understand and generate more knowledge about this species, liver transcriptome sequencing and de novo assembly was performed for E. marginatus by Next Generation Sequencing (NGS). Sequences obtained were used as a tool to validate the presence of key genes relevant to lipid metabolism, and their expression was quantified by qPCR. Moreover, we investigated the influence of supplementing different dietary fatty acids on hepatic lipid metabolism. The results showed that the different fatty acids added to the diet dramatically changed the gene expression of some key enzymes associated with lipid metabolism as well as hepatic fatty acid profiles. Elongase 5 gene expression was shown to influence intermediate hepatic fatty acid elongation in all experimental groups. Hepatic triglycerides reflected the diet composition more than hepatic phospholipids, and were characterized mainly by the high percentage of 18:3n3 in animals fed with a linseed oil rich diet. Results for the saturated and monounsaturated fatty acids suggest a self-regulatory potential for retention and oxidation processes in liver, since in general the tissues did not directly reflect these fatty acid diet compositions. These results indicated that genes involved in lipid metabolism pathways might be potential biomarkers to assess lipid requirements in the formulated diet for this species.

Usual normalization strategies for gene expression studies impair the detection and analysis of circadian patterns

Recent studies have shown that transcriptomes from different tissues present circadian oscillations. Therefore, the endogenous variation of total RNA should be considered as a potential bias in circadian studies of gene expression. However, normalization strategies generally include the equalization of total RNA concentration between samples prior to cDNA synthesis. Moreover, endogenous housekeeping genes (HKGs) frequently used for data normalization may exhibit circadian variation and distort experimental results if not detected or considered. In this study, we controlled experimental conditions from the amount of initial brain tissue samples through extraction steps, cDNA synthesis, and quantitative real time PCR (qPCR) to demonstrate a circadian oscillation of total RNA concentration. We also identified that the normalization of the RNA's yield affected the rhythmic profiles of different genes, including Per1-2 and Bmal1. Five widely used HKGs (Actb, Eif2a, Gapdh, Hprt1, and B2m) also presented rhythmic variations not detected by geNorm algorithm. In addition, the analysis of exogenous microRNAs (Cel-miR-54 and Cel-miR-39) spiked during RNA extraction suggests that the yield was affected by total RNA concentration, which may impact circadian studies of small RNAs. The results indicate that the approach of tissue normalization without total RNA equalization prior to cDNA synthesis can avoid bias from endogenous broad variations in transcript levels. Also, the circadian analysis of 2^{-Cycle threshold (Ct)} data, without HKGs, may be an alternative for chronobiological studies under controlled experimental conditions.

Analysis of Different Approaches for the Selection of Reference Genes in RT-qPCR Experiments: A Case Study in Skeletal Muscle of Growing Mice

The reliability of reverse transcription-quantitative PCR (RT-qPCR) results in gene expression studies depends on the approaches used to account for non-biological variations. In order to find a proper normalization strategy for the study of genes related to growth hormone signaling in skeletal muscle of growing mice, nine unrelated genes were evaluated as internal controls. According to the most used algorithms–geNorm, the Comparative ΔCq method, NormFinder and BestKeeper–GSK3B, YWHAZ, RPL13A and RN18S were found as the most stable. However, the relative expression levels of eight of the potential reference genes assessed decreased with age in cDNA samples obtained from the same amount of total RNA. In a different approach to analyze this apparent discrepancy, experiments were performed with cDNA obtained from equal amounts of purified mRNA. Since the decline was still observed, the hypothesis of an age-related change in mRNA to total RNA ratio that could account for the systematic decrease was rejected. Differences among experimental groups could be due to a substantial increase with age in highly expressed mRNAs, which would bias the quantitation of the remaining genes. Consequently, those reference genes reflecting this dilution effect, which would have been discarded considering their variable relative expression levels, arose as suitable internal controls.

Identification of Sex and Female's Reproductive Stage in Commercial Fish Species through the Quantification of Ribosomal Transcripts in Gonads

The estimation of maturity and sex of fish stocks in European waters is a requirement of the EU Data Collection Framework as part of the policy to improve fisheries management. On the other hand, research on fish biology is increasingly focused in molecular approaches, researchers needing correct identification of fish sex and reproductive stage without necessarily having in house the histological know-how necessary for the task. Taking advantage of the differential gene transcription occurring during fish sex differentiation and gametogenesis, the utility of 5S ribosomal RNA (5S rRNA) and General transcription factor IIIA (gtf3a) in the molecular identification of sex and gametogenic stage was tested in different economically-relevant fish species from the Bay of Biscay. Gonads of 9 fish species (, Atlantic, Atlantic-chub and horse mackerel, blue whiting, bogue, European anchovy, hake and pilchard and megrim), collected from local commercial fishing vessels were histologically sexed and 5S and 18S rRNA concentrations were quantified by capillary electrophoresis to calculate a 5S/18S rRNA index. Degenerate primers permitted cloning and sequencing of gtf3a fragments in 7 of the studied species. 5S rRNA and gtf3a transcript levels, together with 5S/18S rRNA index, distinguished clearly ovaries from testis in all of the studied species. The values were always higher in females than in males. 5S/18S rRNA index values in females were always highest when fish were captured in early phases of ovary development whilst, in later vitellogenic stages, the values decreased significantly. In megrim and European anchovy, where gonads in different oogenesis stages were obtained, the 5S/18S rRNA index identified clearly gametogenic stage. This approach, to the sexing and the quantitative non-subjective identification of the maturity stage of female fish, could have multiple applications in the study of fish stock dynamics, fish reproduction and fecundity and fish biology in general.

Identification of Novel Reference Genes Suitable for qRT-PCR Normalization with Respect to the Zebrafish Developmental Stage

Reference genes used in normalizing qRT-PCR data are critical for the accuracy of gene expression analysis. However, many traditional reference genes used in zebrafish early development are not appropriate because of their variable expression levels during embryogenesis. In the present study, we used our previous RNA-Seq dataset to identify novel reference genes suitable for gene expression analysis during zebrafish early developmental stages. We first selected 197 most stably expressed genes from an RNA-Seq dataset (29,291 genes in total), according to the ratio of their maximum to minimum RPKM values. Among the 197 genes, 4 genes with moderate expression levels and the least variation throughout 9 developmental stages were identified as candidate reference genes. Using four independent statistical algorithms (delta-CT, geNorm, BestKeeper and NormFinder), the stability of qRT-PCR expression of these candidates was then evaluated and compared to that of actb1 and actb2, two commonly used zebrafish reference genes. Stability rankings showed that two genes, namely mobk13 (mob4) and lsm12b, were more stable than actb1 and actb2 in most cases. To further test the suitability of mobk13 and lsm12b as novel reference genes, they were used to normalize three well-studied target genes. The results showed that mobk13 and lsm12b were more suitable than actb1 and actb2 with respect to zebrafish early development. We recommend mobk13 and lsm12b as new optimal reference genes for zebrafish qRT-PCR analysis during embryogenesis and early larval stages.

Rapid effects of essential fatty acid deficiency on growth and development parameters and transcription of key fatty acid metabolism genes in juvenile barramundi (Lates calcarifer)

Barramundi (Lates calcarifer), a catadromous teleost of significant and growing commercial importance, are reported to have limited fatty acid bioconversion capability and therefore require preformed long-chain PUFA (LC-PUFA) as dietary essential fatty acid (EFA). In this study, the response of juvenile barramundi (47·0 g/fish initial weight) fed isolipidic and isoenergetic diets with 8·2% added oil was tested. The experimental test diets were either devoid of fish oil (FO), and thus with no n-3 LC-PUFA (FO FREE diet), or with a low inclusion of FO (FO LOW diet). These were compared against a control diet containing only FO (FO CTRL diet) as the added lipid source, over an 8-week period. Interim samples and measurements were taken fortnightly during the trial in order to define the aetiology of the onset and progression of EFA deficiency. After 2 weeks, the fish fed the FO FREE and FO LOW diets had significantly lower live-weights, and after 8 weeks significant differences were detected for all performance parameters. The fish fed the FO FREE diet also had a significantly higher incidence of external abnormalities. The transcription of several genes involved in fatty acid metabolism was affected after 2 weeks of feeding, showing a rapid nutritional regulation. This experiment documents the aetiology of the onset and the progression of EFA deficiency in juvenile barramundi and demonstrates that such deficiencies can be detected within 2 weeks in juvenile fish.

Identification of Barramundi (Lates calcarifer) DC-SCRIPT, a Specific Molecular Marker for Dendritic Cells in Fish

Antigen presentation is a critical step bridging innate immune recognition and specific immune memory. In mammals, the process is orchestrated by dendritic cells (DCs) in the lymphatic system, which initiate clonal proliferation of antigen-specific lymphocytes. However, fish lack a classical lymphatic system and there are currently no cellular markers for DCs in fish, thus antigen-presentation in fish is poorly understood. Recently, antigen-presenting cells similar in structure and function to mammalian DCs were identified in various fish, including rainbow trout (Oncorhynchus mykiss) and zebrafish (Danio rerio). The present study aimed to identify a potential molecular marker for DCs in fish and therefore targeted DC-SCRIPT, a well-conserved zinc finger protein that is preferentially expressed in all sub-types of human DCs. Putative dendritic cells were obtained in culture by maturation of spleen and pronephros-derived monocytes. DC-SCRIPT was identified in barramundi by homology using RACE PCR and genome walking. Specific expression of DC-SCRIPT was detected in barramundi cells by Stellaris mRNA FISH, in combination with MHCII expression when exposed to bacterial derived peptidoglycan, suggesting the presence of DCs in L. calcarifer. Moreover, morphological identification was achieved by light microscopy of cytospins prepared from these cultures. The cultured cells were morphologically similar to mammalian and trout DCs. Migration assays determined that these cells have the ability to move towards pathogens and pathogen associated molecular patterns, with a preference for peptidoglycans over lipopolysaccharides. The cells were also strongly phagocytic, engulfing bacteria and rapidly breaking them down. Barramundi DCs induced significant proliferation of responder populations of T-lymphocytes, supporting their role as antigen presenting cells. DC-SCRIPT expression in head kidney was higher 6 and 24 h following intraperitoneal challenge with peptidoglycan and lipopolysaccharide and declined after 3 days relative to PBS-injected controls. Relative expression was also lower in the spleen at 3 days post challenge but increased again at 7 days. As DC-SCRIPT is a constitutively expressed nuclear receptor, independent of immune activation, this may indicate initial migration of immature DCs from head kidney and spleen to the injection site, followed by return to the spleen for maturation and antigen presentation. DC-SCRIPT may be a valuable tool in the investigation of antigen presentation in fish and facilitate optimisation of vaccines and adjuvants for aquaculture.

Selection of reference genes as internal controls for gene expression in tissues of red abalone Haliotis rufescens (Mollusca, Vetigastropoda; Swainson, 1822)

The red abalone Haliotis rufescens is one of the most important species for aquaculture in Baja California, México, and despite this, few gene expression studies have been done in tissues such as gill, head and gonad. For this purpose, reverse transcription and quantitative real time PCR (RT-qPCR) is a powerful tool for gene expression evaluation. For a reliable analysis, however, it is necessary to select and validate housekeeping genes that allow proper transcription quantification. Stability of nine housekeeping genes (ACTB, BGLU, TUBB, CY, GAPDH, HPRTI, RPL5, SDHA and UBC) was evaluated in different tissues of red abalone (gill, head and gonad/digestive gland). Four-fold serial dilutions of cDNA (from 25 ngμL(-1) to 0.39 ngμL(-1)) were used to prepare the standard curve, and it showed gene efficiencies between 0.95 and 0.99, with R(2)=0.99. geNorm and NormFinder analysis showed that RPL5 and CY were the most stable genes considering all tissues, whereas in gill HPRTI and BGLU were most stable. In gonad/digestive gland, RPL5 and TUBB were the most stable genes with geNorm, while SDHA and HPRTI were the best using NormFinder. Similarly, in head the best genes were RPL5 and UBC with geNorm, and GAPDH and CY with NormFinder. The technical variability analysis with RPL5 and abalone gonad/digestive gland tissue indicated a high repeatability with a variation coefficient within groups ≤ 0.56% and between groups ≤ 1.89%. These results will help us for further research in reproduction, thermoregulation and endocrinology in red abalone.

Comprehensive selection of reference genes for gene expression normalization in sugarcane by real time quantitative rt-PCR

The increasingly used real time quantitative reverse transcription-PCR (qRT-PCR) method for gene expression analysis requires one or several reference gene(s) acting as normalization factor(s). In order to facilitate gene expression studies in sugarcane (Saccharum officinarum), a non-model plant with limited genome information, the stability of 13 candidate reference genes was evaluated. The geNorm, NormFinder and deltaCt methods were used for selecting stably expressed internal controls across different tissues and under various experimental treatments. These results revealed that, among these 13 candidate reference genes, GAPDH, eEF-1a and eIF-4α were the most stable and suitable for use as normalization factors across all various experimental samples. In addition, APRT could be a candidate for examining the relationship between gene copy number and transcript levels in sugarcane tissue samples. According to the results evaluated by geNorm, combining CUL and eEF-1α in hormone treatment experiments; CAC and CUL in abiotic stress tests; GAPDH, eEF-1a and CUL in all treatment samples plus CAC, CUL, APRT and TIPS-41 in cultivar tissues as groups for normalization would lead to more accurate and reliable expression quantification in sugarcane. This is the first systematic validation of reference genes for quantification of transcript expression profiles in sugarcane. This study should provide useful information for selecting reference genes for more accurate quantification of gene expression in sugarcane and other plant species.

Postprandial molecular responses in the liver of the barramundi, Lates calcarifer

The regulation of gene expression by nutrients is an important mechanism governing energy storage and growth in most animals, including fish. At present, very few genes that regulate intermediary metabolism have been identified in barramundi, nor is there any understanding of their nutritional regulation. In this study, a partial barramundi liver transcriptome was assembled from next-generation sequencing data and published barramundi EST sequences. A large number of putative metabolism genes were identified in barramundi, and the changes in the expression of 24 key metabolic regulators of nutritional pathways were investigated in barramundi liver over a time series immediately after a meal of a nutritionally optimised diet for this species. Plasma glucose and free amino acid levels showed a mild postprandial elevation which peaked 2 h after feeding, and had returned to basal levels within 4 or 8 h, respectively. Significant activation or repression of metabolic nuclear receptor regulator genes were observed, in combination with activation of glycolytic and lipogenic pathways, repression of the final step of gluconeogenesis and activation of the Akt-mTOR pathway. Strong correlations were identified between a number of different metabolic genes, and the coordinated co-regulation of these genes may underlie the ability of this fish to utilise dietary nutrients. Overall, these data clearly demonstrate a number of unique postprandial responses in barramundi compared with other fish species and provide a critical step in defining the response to different dietary nutrient sources.

Persian sturgeon insulin-like growth factor I: molecular cloning and expression during various nutritional conditions

The effects of different periods of starvation (1, 2, 3, and 4 weeks) and subsequent re-feeding (over a 4 week) on the compensatory growth performance and insulin-like growth factor I (IGF-I) mRNA expression in liver and white muscle were investigated in juvenile Persian sturgeon (Acipenser persicus). First, a fragment of 617 nucleotides coding for IGF-I was cloned from liver, which included an open reading frame of 486 nucleotides, encoding a 162 amino acid preproIGF-I. This is composed of a 45 aa for signal peptide, a 117 aa for the mature peptide comprising the B, C, A, and D domains, and a 47 aa for E domain. The mature Persian sturgeon IGF-I exhibits high sequence identities with other sturgeon species and teleost, ranging between 68 and 95 %. The pattern of IGF-I mRNA expression in the liver and white muscle was measured in response to different periods of starvation and subsequent re-feeding. Nutritional status influenced IGF-I mRNA expression pattern in both liver and muscle. IGF-I mRNA expression in the liver increased during starvation, before decreasing after re-feeding. Furthermore, white muscle IGF-I mRNA expression showed better responses to nutritional status and decreased following starvation and increased by re-feeding. However, changes in the expression of IGF-I mRNA were not significantly different between any of the treatments in both tissues. These data suggest that muscle and liver IGF-I mRNA expression do not have a regulatory role for somatic growth induced by compensatory growth in Persain sturgeon.

Otx2 expression and implications for olfactory imprinting in the anemonefish, Amphiprion percula

The otx2 gene encodes a transcription factor (OTX2) essential in the formation of the brain and sensory systems. Specifically, OTX2-positive cells are associated with axons in the olfactory system of mice and otx2 is upregulated in odour-exposed zebrafish, indicating a possible role in olfactory imprinting. In this study, otx2 was used as a candidate gene to investigate the molecular mechanisms of olfactory imprinting to settlement cues in the coral reef anemonefish, Amphiprion percula. The A. percula otx2 (Ap-otx2) gene was elucidated, validated, and its expression tested in settlement-stage A. percula by exposing them to behaviourally relevant olfactory settlement cues in the first 24 hours post-hatching, or daily throughout the larval phase. In-situ hybridisation revealed expression of Ap-otx2 throughout the olfactory epithelium with increased transcript staining in odour-exposed settlement-stage larval fish compared to no-odour controls, in all scenarios. This suggests that Ap-otx2 may be involved in olfactory imprinting to behaviourally relevant settlement odours in A. percula.

Isolation and selection of suitable reference genes for real-time PCR analyses in the skeletal muscle of the fine flounder in response to nutritional status: assessment and normalization of gene expression of growth-related genes

In the present study, different reference genes were isolated, and their stability in the skeletal muscle of fine flounder subjected to different nutritional states was assessed using geNorm and NormFinder. The combinations between 18S and ActB; Fau and 18S; and Fau and Tubb were chosen as the most stable gene combinations in feeding, long-term fasting and refeeding, and short-term refeeding conditions, respectively. In all periods, ActB was identified as the single least stable gene. Subsequently, the expression of the myosin heavy chain (MYH) and the insulin-like growth factor-I receptor (IGF-IR) was assessed. A large variation in MYH and IGF-IR expression was found depending on the reference gene that was chosen for normalizing the expression of both genes. Using the most stable reference genes, mRNA levels of MYH decreased and IGF-IR increased during fasting, with both returning to basal levels during refeeding. However, the drop in mRNA levels for IGF-IR occurred during short-term refeeding, in contrast with the observed events in the expression of MYH, which occurred during long-term refeeding. The present study highlights the vast differences incurred when using unsuitable versus suitable reference genes for normalizing gene expression, pointing out that normalization without proper validation could result in a bias of gene expression.

How to perform RT-qPCR accurately in plant species? A case study on flower colour gene expression in an azalea (Rhododendron simsii hybrids) mapping population

Background

Flower colour variation is one of the most crucial selection criteria in the breeding of a flowering pot plant, as is also the case for azalea (Rhododendron simsii hybrids). Flavonoid biosynthesis was studied intensively in several species. In azalea, flower colour can be described by means of a 3-gene model. However, this model does not clarify pink-coloration. The last decade gene expression studies have been implemented widely for studying flower colour. However, the methods used were often only semi-quantitative or quantification was not done according to the MIQE-guidelines. We aimed to develop an accurate protocol for RT-qPCR and to validate the protocol to study flower colour in an azalea mapping population.

Results

An accurate RT-qPCR protocol had to be established. RNA quality was evaluated in a combined approach by means of different techniques e.g. SPUD-assay and Experion-analysis. We demonstrated the importance of testing noRT-samples for all genes under study to detect contaminating DNA. In spite of the limited sequence information available, we prepared a set of 11 reference genes which was validated in flower petals; a combination of three reference genes was most optimal. Finally we also used plasmids for the construction of standard curves. This allowed us to calculate gene-specific PCR efficiencies for every gene to assure an accurate quantification. The validity of the protocol was demonstrated by means of the study of six genes of the flavonoid biosynthesis pathway. No correlations were found between flower colour and the individual expression profiles. However, the combination of early pathway genes (CHS, F3H, F3'H and FLS) is clearly related to co-pigmentation with flavonols. The late pathway genes DFR and ANS are to a minor extent involved in differentiating between coloured and white flowers. Concerning pink coloration, we could demonstrate that the lower intensity in this type of flowers is correlated to the expression of F3'H.

Conclusions

Currently in plant research, validated and qualitative RT-qPCR protocols are still rare. The protocol in this study can be implemented on all plant species to assure accurate quantification of gene expression. We have been able to correlate flower colour to the combined regulation of structural genes, both in the early and late branch of the pathway. This allowed us to differentiate between flower colours in a broader genetic background as was done so far in flower colour studies. These data will now be used for eQTL mapping to comprehend even more the regulation of this pathway.

Cross-platform microarray meta-analysis for the mouse jejunum selects novel reference genes with highly uniform levels of expression

Reference genes (RGs) with uniform expression are used for normalization of reverse transcription quantitative PCR (RT-qPCR) data. Their optimization for a specific biological context, e.g. a specific tissue, has been increasingly considered. In this article, we compare RGs identified by expression data meta-analysis restricted to the context tissue, the jejunum of Mus musculus domesticus, i) to traditional RGs, ii) to expressed interspersed repeated DNA elements, and iii) to RGs identified by meta-analysis of expression data from diverse tissues and conditions. To select the set of candidate RGs, we developed a novel protocol for the cross-platform meta-analysis of microarray data. The expression stability of twenty-four putative RGs was analysed by RT-qPCR in at least 14 jejunum samples of the mouse strains C57Bl/6N, CD1, and OF1. Across strains, the levels of expression of the novel RGs Plekha7, Zfx, and Ube2v1 as well as of Oaz1 varied less than two-fold irrespective of genotype, sex or their combination. The gene set consisting of Plekha7 and Oaz1 showed superior expression stability analysed with the tool RefFinder. The novel RGs are functionally diverse. This facilitates expression studies over a wide range of conditions. The highly uniform expression of the optimized RGs in the jejunum points towards their involvement in tightly regulated pathways in this tissue. We also applied our novel protocol of cross-microarray platform meta-analysis to the identification of RGs in the duodenum, the ileum and the entire small intestine. The selection of RGs with improved expression stability in a specific biological context can reduce the number of RGs for the normalization step of RT-qPCR expression analysis, thus reducing the number of samples and experimental costs.

Gelatinases impart susceptibility to high-fat diet-induced obesity in mice

Gelatinases play a role in adipose and muscle hypertrophy and could be involved in tissue remodeling in response to high-fat diet (HFD) intake. This study tested potential roles of gelatinases (matrix metalloproteinses-2 and -9 [MMP-2 and -9]) in relationship to an antigrowth factor [myostatin (MSTN)] known to be dysregulated in relation to HFD-induced obesity (HFDIO) propensity. In vitro and ex vivo analyses demonstrated that MMP-9 increased mature MSTN levels, indicating a potential role of gelatinases in MSTN activation in vivo. HFD intake resulted in increased body weight and circulating blood glucose values in C57BL/6J and MMP-9 null mice, with no changes observed in SWR/J mice. HFD intake attenuated MMP-9 and MMP-2 mRNA levels in SWR/J mice while elevating MMP-2 levels in skeletal muscle in C57BL/6J mice. In MMP-9 null mice, the effects of HFD intake were muted. Consistent with changes in mRNA levels, HFD intake increased MMP-9 activity in muscle tissue of C57BL/6J mice, demonstrating a strong relationship between HFDIO susceptibility and local MMP regulation. Overall, resistance to HFDIO appears to correspond to low MMP-9 and MSTN levels, suggesting a role of MMP-9 in MSTN activation in local tissue responses to HFD intake.

Expression profile of immune-related genes in Lates calcarifer infected by Cryptocaryon irritans

Cryptocaryon irritans causes Cyptocaryonosis or white spot disease in a wide range of marine fish including Lates calcarifer (Asian seabass). However, the immune response of this fish to the parasite is still poorly understood. In this study, quantitative polymerase chain reaction (qPCR) was performed to assess the expression profile of immune-related genes in L. calcarifer infected by C. irritans. A total of 21 immune-related genes encoding various functions in the fish immune system were utilized for the qPCR analysis. The experiment was initiated with the infection of juvenile fish by exposure to theronts from 200 C. irritans cysts, and non-infected juvenile fish were used as controls. Spleen, liver, gills and kidney tissues were harvested at three days post-infection from control and infected fish. In addition, organs were also harvested on day-10 post-infection from fish that had been allowed to recover from day-4 up to day-10 post-infection. L. calcarifer exhibited pathological changes on day-3 post-infection with the characteristic presence of white spots on the entire fish body, excessive mucus production and formation of a flap over the fish eye. High quality total RNA was extracted from all tissues and qPCR was performed. The qPCR analysis on the cohort of 21 immune-related genes of the various organs harvested on day-3 post-infection demonstrated that most genes were induced significantly (p < 0.05) in all tissues, particularly liver (11/21 genes) and kidney (11/21). The expression profile demonstrated that induction of the MHC Class IIα gene was the highest compared to the other genes followed by serum amyloid A, CC chemokine and hepcidin-2 precursor genes. In fish that were allowed to recover from the C. irritans infection (10 days post-infection), expression of the immune-related genes was down-regulated to levels similar to the control fish. These results provide insights into the interaction between C. irritans and L. calcarifer and suggest that the innate immune system plays an important role in early defence against parasite infection allowing the fish to eventually recover from the infection.

Gene expression and activity of digestive enzymes during the larval development of Asian seabass (Lates calcarifer)

Knowledge of digestive enzyme development during larval stages provides a better understanding of the digestive and nutritional physiology of fish larvae. This study characterized the ontogeny of key digestive enzymes in Asian seabass larvae from hatching to juvenile stage (30 days post hatch, dph) using molecular and biochemical approaches. Gene expression and activity of pepsinogen (pg), trypsinogen (try), chymotrypsinogen (ctr), bile salt-activated lipase (bal), α-amylase (amy), leucine aminopeptidase (lap) and alkaline phosphatase (alp) were determined. Gene expression and enzyme activity of all digestive enzymes were detectable from hatching. Pepsinogen mRNA levels were close to detection limit during 0-15 dph, but were highly expressed from 18 dph and onwards. This coincided with a marked increase in specific and individual pepsin activity, indicating complete development of digestive function. Expression levels of try, ctr, amy and bal were high between 3 and 15 dph and thereafter a decreasing trend was observed. Intestinal enzymes, lap and alp, showed highest expression levels during the yolk sac stage, and thereafter levels decreased. Activity of all digestive enzymes increased from around 18 dph and onwards. In conclusion, the development of main digestive enzymes in Asian seabass larvae shows a similar pattern to that of other marine fish species.

Inbred strains of zebrafish exhibit variation in growth performance and myostatin expression following fasting

Although the zebrafish (Danio rerio) has been widely utilized as a model organism for several decades, there is little information available on physiological variation underlying genetic variation among the most commonly used inbred strains. This study evaluated growth performance using physiological and molecular markers of growth in response to fasting in six commonly used zebrafish strains [AB, TU, TL, SJA, WIK, and petstore (PET) zebrafish]. Fasting resulted in a standard decrease in whole blood glucose levels, a typical vertebrate glucose metabolism pattern, in AB, PET, TL, and TU zebrafish strains. Alternatively, fasting did not affect glucose levels in SJA and WIK zebrafish strains. Similarly, fasting had no effect on myostatin mRNA levels in AB, PET, TU, and WIK zebrafish strains, but decreased myostatin-1 and -2 mRNA levels in SJA zebrafish. Consistent with previous work, fasting increased myostatin-2 mRNA levels in TL zebrafish. These data demonstrate that variation is present in growth performance between commonly used inbred strains of zebrafish. These data can help future research endeavors by highlighting the attributes of each strain with regard to growth performance so that the most fitting strain may be utilized.

The gene expression response of the catadromous perciform barramundi Lates calcarifer to an acute heat stress

The acute heat-shock response of the tropical estuarine fish species barramundi Lates calcarifer as indicated by the expression of genes within stress (hsp 90AA, hsp 90AB, hsp 70 and hsc 70), metabolic (cisy, cco II and ldh) and growth (igf1 and mstn 1) related pathways was examined following an increase in water temperature from 28 to 36° C over 30 min. Lates calcarifer were maintained at the acute stress temperature of 36° C for 1 h before being returned to 28° C and allowed to recover at this temperature for a further 2 weeks. Muscle tissue sampling over the experimental period allowed for the expression quantification of stress, metabolic and growth-related genes via quantitative real-time polymerase chain reaction (qrt-PCR) where a robust and reliable normalization approach identified both α-tub and Rpl8 as appropriate genes for the analysis of gene expression in response to an acute heat stress. hsp90AA and hsp70 of the inducible heat-shock response pathway showed a massive up-regulation of gene expression in response to heat stress, whilst the constitutive heat-shock genes hsp90AB and hsp70 showed no change over the course of the experiment and a small increase after 2 weeks of recovery, respectively. Of the three genes representing the metabolic pathway (cisy, cco II and ldh) only cco II changed significantly showing a decrease in gene expression, which may suggest a small suppression of aerobic metabolism. igf1 of the growth pathway showed no significant differences in response to an acute heat stress, whilst mstn1 increased at the beginning of the heat stress but returned to basal levels soon after. Overall, the results demonstrate that an acute heat stress in L. calcarifer caused a significant increase in the expression of genes from the stress response pathway and a possible decrease in aerobic metabolism with only relatively minor changes to the growth pathway highlighting the hardy nature of L. calcarifer and its resilience in coping with sudden temperature changes routinely encountered within its natural environment.

Identification of a tubulin-α gene specifically expressed in testis and adductor muscle during stable reference gene selection in the hermaphrodite gonad of the lion's paw scallop Nodipecten subnodosus

For non-model species, as many used for aquaculture, with minimal or no genomic information, relative quantification of gene expression studies requires preliminary research including the isolation of potential reference genes and the identification of those stably expressed under the biological conditions of interest. Here we report on the isolation of five partial gene sequences from gonad tissue cDNA in the functional hermaphrodite scallop Nodipecten subnodosus to be evaluated as reference genes: 18S-rRNA, riboprotein l8 (rp-l8), actin-β (act-β), elongation factor 1α (ef-1α) and alpha-tubulin-α (tub-α). We found that 18S-rRNA was stably expressed independently of the priming method used to reverse transcribe RNA to cDNA, oligo-dT or random hexamer. Stability analysis for the five putative reference genes with geNorm and NormFinder indicated that 18S together with rp-l8 were the most stable genes for normalization of gene expression during gonad development in both, male and female sexual regions of the hermaphrodite N. subnodosus. The least stable gene was tub-α, showing a biased expression profile between sexual regions of the gonad, therefore this gene was analyzed thereafter as a target gene together with vitellogenin (vit) and a DEAD-box RNA helicase (dbx) gene. Relative expression, estimated by normalization with the combination of 18S and rp-l8 as reference genes, indicated that as gonad development advanced two of the target genes were up-regulated, tub-α in the male region and vit in the female region. Whereas an increased expression was expected during development for vit for its known role in vitellogenesis, the increased expression of tub-α in the male sexual region was unexpected, and pointed toward this gene being a testis-specific α-tubulin isotype. Further analyses of gene expression among tissues indicated that tub-α is specifically and highly expressed in the male gonad, although expression in adductor muscle was also observed at significantly lower levels. The existence of testis specific α- and β-tubulins has been previously reported in other taxa, relating their function to sperm axoneme formation. Tissue-specific tubulin genes, particularly their promoters, have recently found an application as native promoters for transgene tissue-specific expression in research and reproductive control of insect plagues. The third target gene, a putative member of the DEAD-box RNA helicase family (dbx), showed no changes in expression during gonad development or between sexual regions, therefore it was chosen to discuss the different statistical inferences resulting from the arbitrary use of 'randomly chosen' reference genes when normalizing gene expression.

Normalizing for biology: accounting for technical and biological variation in levels of reference gene and insulin-like growth factor 1 (igf1) transcripts in fish livers

Feeding, fasting and re-feeding is a common experimental paradigm for studying growth endocrinology. Herein we demonstrate dynamic changes in the livers of coho salmon under these conditions and how changes in liver composition can influence quantification and interpretation of liver gene expression data. A three-week fast resulted in decreases in hepatosomatic index (liver size), liver glycogen content, and liver DNA concentration. In addition, significant differences were found in liver transcript levels from fed and fasted fish for the reference genes, arp and ef1a, when these were normalized to total RNA. We took the additional step of normalizing reference gene transcript levels to the liver homogenate RNA/DNA ratio to account for differences in RNA yield/cell and the number of cells sampled, normalizing to transcript number per cell rather than transcript number per unit RNA. After this additional step no significant differences in liver transcript levels of reference genes were found. The significance of these results was illustrated by normalizing liver transcript levels of insulin-like growth factor 1 (igf1) to ef1a transcript levels or ef1a transcript levels by RNA/DNA. The different normalization strategies resulted in differing patterns of change in igf1 transcript levels between fed and fasted fish. The novelty of this work rests in a two-step normalization process, attempting to account for both 1) technical errors in reverse transcription and qPCR reactions, and 2) biological variance in liver samples.

Mechanisms of colour adaptation in the prawn Penaeus monodon

Exposure of prawns to dark- or light-coloured substrates is known to trigger a strong colour adaptation response through expansion or contraction of the colouration structures in the prawn hypodermis. Despite the difference in colour triggered by this adaptive response, total levels of the predominant carotenoid pigment, astaxanthin, are not modified, suggesting that another mechanism is regulating this phenomenon. Astaxanthin binds to a specific protein called crustacyanin (CRCN), and it is the interaction between the quantities of each of these compounds that produces the diverse range of colours seen in crustacean shells. In this study, we investigated the protein changes and genetic regulatory processes that occur in prawn hypodermal tissues during adaptation to black or white substrates. The amount of free astaxanthin was higher in animals adapted to dark substrate compared with those adapted to light substrate, and this difference was matched by a strong elevation of CRCN protein. However, there was no difference in the expression of CRCN genes either across the moult cycle or in response to background substrate colour. These results indicate that exposure to a dark-coloured substrate causes an accumulation of CRCN protein, bound with free astaxanthin, in the prawn hypodermis without modification of CRCN gene expression. On light-coloured substrates, levels of CRCN protein in the hypodermis are reduced, but the carotenoid is retained, undispersed in the hypodermal tissue, in an esterified form. Therefore, the abundance of CRCN protein affects the distribution of pigment in prawn hypodermal tissues, and is a crucial regulator of the colour adaptation response in prawns.

Expression of zona pellucida B proteins in juvenile rare minnow (Gobiocypris rarus) exposed to 17α-ethinylestradiol, 4-nonylphenol and bisphenol A

Zona pellucida (ZP) containing proteins are glycoproteins in teleost chorion and are encoded by several gene subfamilies, mainly including ZPA, ZPB, ZPC and ZPX genes. In teleost species, ZP genes are expressed either in liver under regulation of estrogen or in ovary. In the present study, five ZP gene isoforms were isolated and characterized in Gobiocypris rarus. The putative amino acid sequences of these ZP gene isoforms contain the typical trefoil motif and a ZP domain. These five G. rarus ZP gene isoforms were named as grZPB.1, grZPB.2, grZPB.3, grZPB.4 and grZPB.5. Real-time quantitative reverse transcription PCR (RT-qPCR) analysis indicated that all these ZP mRNA isoforms were exclusively expressed in ovary. G. rarus juveniles at the age of 21 days postfertilization were exposed to 17α-ethinylestradiol (EE2; 0.01, 0.1 and 1 nM), 4-nonylphenol (4-NP; 10, 100 and 1000 nM) or bisphenol A (BPA; 0.1, 1 and 10nM) for 3 days. mRNA expressions of ZPB isoforms following the exposure to xenoestrogen were detected by RT-qPCR. Data were analyzed by the 2(-△△Cq) method. The results indicate that induction by 0.1-1nM EE2 on mRNA expression of the grZPB isoforms is weaker than for vitellogenin. 4-NP exposures at three concentrations had differential effects on the grZPBs. BPA at three concentrations weakly induced mRNA expression of the grZPB isoforms.

Gene expression specificity of the mussel antifungal mytimycin (MytM)

We previously reported the nucleotide sequences and diversity of mytimycin (MytM) from the Mediterranean mussel, Mytilus galloprovincialis. Using real-time PCR (q-PCR), we observed that the MytM gene was mainly expressed in circulating hemocytes and to a less extent in the mantle. In vivo challenge with bacteria or with the yeast, Candida albicans, did not increase the expression as measured by q-PCR in hemocytes. By contrast, injection of the filamentous fungus, Fusarium oxysporum, induced a sudden and strong increase of expression at 9h p.i. (stimulation index of 25.7 ± 2.1). Optimum stimulating dose was 10(4) spores of F. oxysporum per mussel. In the same samples, AMP mytilin and myticin showed no stimulation. Consequently, we hypothesized the existence of 2 different signal transduction pathways, one activated by bacteria and yeast, the other triggered by filamentous fungi. A second challenge performed with F. oxysporum 24 h after the first challenge induced an increase of MytM gene expression (stimulation index of 3.5 ± 1.7). However, this second increase was significantly lower than the first, suggesting less efficient response rather than significant protection.

Growing backwards: an inverted role for the shrimp ortholog of vertebrate myostatin and GDF11

Myostatin (MSTN) and growth differentiation factor-11 (GDF11) are closely related proteins involved in muscle cell growth and differentiation as well as neurogenesis of vertebrates. Both MSTN and GDF11 negatively regulate their functions. Invertebrates possess a single ortholog of the MSTN/GDF11 family. In order to understand the role of MSTN/GDF11 in crustaceans, the gene ortholog was identified and characterized in the penaeid shrimp Penaeus monodon. The overall protein sequence and specific functional sites were highly conserved with other members of the MSTN/GDF11 family. Gene transcripts of pmMstn/Gdf11, assessed by real-time PCR, were detected in a variety of tissue types and were actively regulated in muscle across the moult cycle. To assess phenotypic function in shrimp, pmMstn/Gdf11 gene expression was downregulated by tail-muscle injection of sequence-specific double-stranded RNA. Shrimp with reduced levels of pmMstn/Gdf11 transcripts displayed a dramatic slowing in growth rate compared with control groups. Findings from this study place the MSTN/GDF11 gene at the centre of growth regulation in shrimp, but suggest that, compared with higher vertebrates, this gene has an opposite role in invertebrates such as shrimp, where levels of gene expression may positively regulate growth.