Evaluation of potential reference genes for real time RT-PCR studies in Atlantic halibut (Hippoglossus Hippoglossus L.); during development, in tissues of healthy and NNV-injected fish, and in anterior kidney leucocytes

Cell death induced by Lepeophtheirus salmonis labial gland protein 3 in salmonid fish leukocytes: A mechanism for disabling host immune responses

The salmon louse (Lepeophtheirus salmonis) is an ectoparasite feeding on mucus, skin, and blood of salmonids. On parasitised fish erosions and, at later lice stages, ulcerations appear at the louse feeding site. In susceptible species like Atlantic salmon (Salmo salar) with a limited rejection of lice, only a mild inflammatory response with minor influx of immune cells is seen at these lesions, as the salmon louse secrete proteins that can dampen immune responses. In a previous study, Lepeophtheirus salmonis labial gland protein 3 (LsLGP3) was suggested to dampen cellular responses, and the present study aimed at increasing our understanding of its mode of action. LsLGP3 was found to be secreted on to the host skin, and both in vivo and in vitro experiments were performed to elucidate its function. Histological analysis of the louse attachment site revealed an epidermal and dermal influx of mainly macrophages and granulocytes after 5 days post infestation. The immune cell influx was deeper in the dermis throughout the louse infestation, and LsLGP3 may be involved in dampening this response. Enriched populations of Atlantic salmon B-cells, T-cells, granulocytes, and monocytes were exposed to recombinant LsLGP3 (recLGP3) in vitro, resulting in a significant decrease in cell viability compared to non-exposed controls. An apoptotic cell morphology with "beads-on-a-string" like protrusions was seen in all leukocyte cell fractions after recLGP3 exposure, but not in erythrocytes or keratocytes. A decreased viability was also detected in pink salmon leucocytes, which was not in leucocytes from non-salmonid species. These functional insights suggest that LsLGP3 specifically induces apoptosis of salmonid leukocytes and is likely a key protein secreted by the lice that disables the Atlantic salmon ability to mount an adequate immune response towards the salmon louse. In vivo LsLGP3 knock down studies indicated that the effect is localised primarily at the lice feeding site, without affecting immune cells that are not situated adjacent to the lice-inflicted lesion. The findings from this study could significantly aid in the development of new immune based anti-salmon louse prophylactic measures and treatments.

Salmon louse labial gland enzymes: implications for host settlement and immune modulation

Salmon louse (Lepeophtheirus salmonis) is a skin- and blood-feeding ectoparasite, infesting salmonids. While feeding, labial gland proteins from the salmon louse may be deposited on the Atlantic salmon (Salmo salar) skin. Previously characterized labial gland proteins are involved in anti-coagulation and may contribute to inhibiting Atlantic salmon from mounting a sufficient immune response against the ectoparasite. As labial gland proteins seem to be important in the host-parasite interaction, we have, therefore, identified and characterized ten enzymes localized to the labial gland. They are a large group of astacins named L. salmonis labial gland astacin 1-8 (LsLGA 1-8), one serine protease named L. salmonis labial gland serine protease 1 (LsLGSP1), and one apyrase named L. salmonis labial gland apyrase 1 (LsLGAp1). Protein domain predictions showed that LsLGA proteins all have N-terminal ShK domains, which may bind to potassium channels targeting the astacins to its substrate. LsLGA1 and -4 are, in addition, expressed in another gland type, whose secrete also meets the host-parasite interface. This suggests that LsLGA proteins may have an anti-microbial function and may prevent secondary infections in the wounds. LsLGAp1 is predicted to hydrolyze ATP or AMP and is, thereby, suggested to have an immune dampening function. In a knockdown study targeting LsLGSP1, a significant increase in IL-8 and MMP13 at the skin infestation site was seen under LsLGSP1 knockdown salmon louse compared to the control, suggesting that LsLGSP1 may have an anti-inflammatory effect. Moreover, most of the identified labial gland proteins are expressed in mature copepodids prior to host settlement, are not regulated by starvation, and are expressed at similar or higher levels in lice infesting the salmon louse-resistant pink salmon (Oncorhynchus gorbuscha). This study, thereby, emphasizes the importance of labial gland proteins for host settlement and their immune dampening function. This work can further contribute to anti-salmon louse treatment such as vaccine development, functional feed, or gene-edited salmon louse-resistant Atlantic salmon.

Identification and characterization of two salmon louse heme peroxidases and their potential as vaccine antigens

Salmon louse, Lepeophtheirus salmonis, represents major challenge for salmon farming. Current treatments impose welfare issues and are costly, whereas prophylactic measures are unavailable. Two salmon louse heme peroxidases (LsPxtl-1 and LsPxtl-2) were tested for their importance for parasite development and as potential vaccine candidates. LsPxtl-1 possesses two heme peroxidase domains and is expressed in ovaries and gut, whereas LsPxtl-2 encodes one domain and contains N-terminal signal peptide and an Eph receptor ligand-binding domain. LsPxtl-1, but not LsPxtl-2, knockdown in nauplius II stage caused poor swimming and death, indicating its importance for parasite development. Immunizations using single DNA plasmid injection encoding the peroxidases or heterologous prime (DNA) and boost (recombinant LsPxtl-2 protein) gave non-significant reduction in lice numbers. Single injection gave low specific antibody levels compared with the prime-boost. The findings suggest LsPxtl-1 is important for parasite development but formulations and vaccination modalities used did not significantly reduce lice infestation.

Transcriptomic and targeted immune transcript analyses confirm localized skin immune responses in Atlantic salmon towards the salmon louse

Atlantic salmon (Salmo salar) are highly susceptible to infestations with the ectoparasite Lepeophtheirus salmonis, the salmon louse. Infestations elicit an immune response in the fish, but the response does not lead to parasite clearance, nor does it protect against subsequent infestations. It is, however, not known why the immune response is not adequate, possibly because the local response directly underneath the louse has been poorly evaluated. The present study describes the transcriptomic response by RNA sequencing of skin at the site of copepodid attachment. Analysing differentially expressed genes, 2864 were higher and 1357 were lower expressed at the louse attachment site compared to uninfested sites in the louse infested fish, while gene expression at uninfested sites were similar to uninfested control fish. The transcriptional patterns of selected immune genes were further detailed in three skin compartments/types: Whole skin, scales only and fin tissue. The elevation of pro-inflammatory cytokines and immune cell marker transcripts observed in whole skin and scale samples were not induced in fin, and a higher cytokine transcript level in scale samples suggest it can be used as a nonlethal sampling method to enhance selective breeding trials. Furthermore, the immune response was followed in both skin and anterior kidney as the infestation developed. Here, newly moulted preadult 1 stage lice induced a higher immune response than chalimi and adult lice. Overall, infestation with salmon louse induce a modest but early immune response with an elevation of mainly innate immune transcripts, with the response primarily localized to the site of attachment.

Small, charged proteins in salmon louse (Lepeophtheirus salmonis) secretions modulate Atlantic salmon (Salmo salar) immune responses and coagulation

Little is known about glandular proteins secreted from the skin- and blood-feeding ectoparasite salmon louse (Lepeophtheirus salmonis). The labial gland has ducts extending into the oral cavity of the lice, and the present study aimed to identify novel genes expressed by this gland type and to investigate their role in modulation of host parameters at the lice feeding site. Five genes associated with labial gland function were identified and named Lepeophteirus salmonis labial gland protein (LsLGP) 1-4 and 1 like (LsLGP1L). All LsLGPs were predicted to be small charged secreted proteins not encoding any known protein domains. Functional studies revealed that LsLGP1 and/or LsLGP1L regulated the expression of other labial gland genes. Immune dampening functions were indicated for LsLGP2 and 3. Whereas LsLGP2 was expressed throughout the parasitic life cycle and found to dampen inflammatory cytokines, LsLGP3 displayed an increased expression in mobile stages and appeared to dampen adaptive immune responses. Expression of LsLGP4 coincided with moulting to the mobile pre-adult I stage where hematophagous feeding is initiated, and synthetic LsLGP4 decreased the clotting time of Atlantic salmon plasma. Results from the present study confirm that the salmon louse secretes immune modulating and anti-coagulative proteins with a potential application in new immune based anti-salmon louse treatments.

Re-evaluation of common carp (Cyprinus carpio L.) housekeeping genes for gene expression studies - considering duplicated genes

Quantitative real-time PCR is one of the most widely used techniques for measuring changes in the expression of target transcripts due to its sensitivity, specificity, and cost-effectiveness. However, the essential step that determines appropriate and correct data interpretation is the selection of proper endogenous control genes. Identifying useful reference genes with stable expression is critical for accurate normalization and precise results. Functional divergence of duplicated genes in tetraploid species, like common carp, can complicate the choice for a proper reference gene. In the present study, we determined the expression stability of duplicated genes of 40s, b2m, ef1α, gapdh, g6pd, and odc1 in different tissues of common carp (Cyprinus carpio L.). Gene expression analysis comprised healthy control fish, fish under bacterial and parasitic infections, and across the early stage of common carp development. Obtained data were compared with the actb gene, which is used widely as a reference in RT-qPCR analysis. The application of the three different algorithms - geNorm, NormFinder, BestKeeper, allowed comparative evaluation of the expression stability of the tested genes. Subsequently, the RefFinder, a web-based tool, was used to rank the examined housekeeping genes comprehensively. We demonstrate variable transcription stability levels in the examined mRNAs as well as differences in expression between paralog gene copies. The 40s, b2m, ef1α and actb genes showed the most stable expression across all physiological conditions and tissues. The gapdh, odc1, and g6pd gene variants demonstrated lower stability. Differences in expression patterns between duplicated genes underline the possibility of functional divergence between them. This aspect should be considered in polyploid species before selecting the reference gene(s). Our study also points on the importance of choice for a reference gene (paralog) when expressing newly identified genes and the spatial expression profile is performed. SUBJECTS: Aquaculture, Molecular Biology, Fish Science.

Evaluation of candidate reference genes for quantitative RTqPCR analysis in goldfish (Carassius auratus L.) in healthy and CyHV-2 infected fish

The accuracy of quantitative real time PCR (RTqPCR) can be attained only when a suitable reference gene is used. The gene expression for a particular gene may vary within different cells at different conditions. Hence, the suitability and stability of various potential reference genes have to be determined for expression studies. In this study, we have examined the potential of four different reference genes including β-Actin (ACTB), 18S ribosomal RNA (18S), glyceraldehyde-3P-dehydrogenase (GAPDH), and elongation factor 1 alpha (EF1A_A) in seven different tissues including gill, liver, kidney, spleen, heart, muscle and intestine of goldfish (Carassius auratus). The housekeeping genes were analyzed from healthy fish and in CyHV-2 challenged fish. Based upon the real time PCR results the gene expression varied among the genes and in tissues. The expression levels of the housekeeping genes were then compared and evaluated with the RefFinder web tool which analyses results using four different algorithms - BestKeeper, delta Ct, geNorm and NormFinder. EF1A_A was ranked to be the best gene in healthy fish by BestKeeper and geNorm analysis. The delta Ct and NormFinder algorithm have found 18S to be a stable gene in healthy fish but 18S was given to be least expressed in challenged fish. ACTB was also given as a stable gene by geNorm analysis in both healthy and challenged fish. Also, in CyHV-2 challenged fish, EF1A_A was identified as the best gene by all the three analysis except by BestKeeper analysis, where it has ranked GADPH as the best housekeeping gene. Expression of the four candidate reference genes differed across all tissue types tested, inferring that a thorough study of the reference genes is necessary for cross tissue comparison. These results can be further used in the immune gene response study of goldfish infected with any viral pathogen to develop better health strategies in the disease management of goldfish aquaculture.

Roles of three putative salmon louse (Lepeophtheirus salmonis) prostaglandin E2 synthases in physiology and host-parasite interactions

BACKGROUND

The salmon louse (Lepeophtheirus salmonis) is a parasite of salmonid fish. Atlantic salmon (Salmo salar) exhibit only a limited and ineffective immune response when infested with this parasite. Prostaglandins (PGs) have many biological functions in both invertebrates and vertebrates, one of which is the regulation of immune responses. This has led to the suggestion that prostaglandin E₂ (PGE₂) is important in the salmon louse host-parasite interaction, although studies of a salmon louse prostaglandin E₂ synthase (PGES) 2 gene have not enabled conformation of this hypothesis. The aim of the present study was, therefore, to characterize two additional PGES-like genes.

METHODS

Lepeophtheirus salmonis microsomal glutathione S-transferase 1 like (LsMGST1L) and LsPGES3L were investigated by sequencing, phylogenetics, transcript localization and expression studies. Moreover, the function of these putative PGES genes in addition to the previously identified LsPGES2 gene was analyzed in double stranded (ds) RNA-mediated knockdown (KD) salmon louse.

RESULTS

Analysis of the three putative LsPGES genes showed a rather constitutive transcript level throughout development from nauplius to the adult stages, and in a range of tissues, with the highest levels in the ovaries or gut. DsRNA-mediated KD of these transcripts did not produce any characteristic changes in phenotype, and KD animals displayed a normal reproductive output. The ability of the parasite to infect or modulate the immune response of the host fish was also not affected by KD.

CONCLUSIONS

Salmon louse prostaglandins may play endogenous roles in the management of reproduction and oxidative stress and may be a product of salmon louse blood digestions.

Distress Regulates Different Pathways in the Brain of Common Carp: A Preliminary Study

In this study, a stress trial was conducted with common carp, one of the most important species in aquaculture worldwide, to identify relevant gene regulation pathways in different areas of the brain. Acute distress due to exposure to air significantly activated the expression of the immediate early gene c-fos in the telencephalon. In addition, evidence for regulation of the two corticotropin-releasing factor (crf) genes in relation to their binding protein (corticotropin-releasing hormone-binding protein, crh-bp) is presented in this preliminary study. Inferences on the effects of due to exposure to air were obtained by using point estimation, which allows the prediction of a single value. This constitutes the best description to date of the previously generally unknown effects of stress in different brain regions in carp. Furthermore, principal component analyses were performed to reveal possible regulation patterns in the different regions of the fish brain. In conclusion, these preliminary studies on gene regulation in the carp brain that has been influenced by exposure to a stressor reveal that a number of genes may be successfully used as markers for exposure to unfavourable conditions.

Cold Adaptation in Antarctic Notothenioids: Comparative Transcriptomics Reveals Novel Insights in the Peculiar Role of Gills and Highlights Signatures of Cobalamin Deficiency

Far from being devoid of life, Antarctic waters are home to Cryonotothenioidea, which represent one of the fascinating cases of evolutionary adaptation to extreme environmental conditions in vertebrates. Thanks to a series of unique morphological and physiological peculiarities, which include the paradigmatic case of loss of hemoglobin in the family Channichthyidae, these fish survive and thrive at sub-zero temperatures. While some of the distinctive features of such adaptations have been known for decades, our knowledge of their genetic and molecular bases is still limited. We generated a reference de novo assembly of the icefish Chionodraco hamatus transcriptome and used this resource for a large-scale comparative analysis among five red-blooded Cryonotothenioidea, the sub-Antarctic notothenioid Eleginops maclovinus and seven temperate teleost species. Our investigation targeted the gills, a tissue of primary importance for gaseous exchange, osmoregulation, ammonia excretion, and its role in fish immunity. One hundred and twenty genes were identified as significantly up-regulated in Antarctic species and surprisingly shared by red- and white-blooded notothenioids, unveiling several previously unreported molecular players that might have contributed to the evolutionary success of Cryonotothenioidea in Antarctica. In particular, we detected cobalamin deficiency signatures and discussed the possible biological implications of this condition concerning hematological alterations and the heavy parasitic loads typically observed in all Cryonotothenioidea.

MicroRNA-29b modulates the innate immune response by suppressing IFNγs production in orange-spotted grouper (Epinephelus coioides)

Interferon-γ (IFNγ), a type II interferon, is essential to host resistance against various infections. Unlike other vertebrates, fish have two types of IFNγs, IFNγ1 (also named IFNγ-rel) and IFNγ2. MicroRNAs (miRNAs) regulate multiple biological processes by suppressing mRNA translation or inducing mRNA degradation. Among them, miR-29 can directly target IFNγ and affact innate and adaptive immune responses in mice. There are five members of the miR-29 family in orange-spotted grouper (Epinephelus coioides), which share the same miRNA seed region. However, whether miR-29 directly targets E. coioides IFNγs and regulate IFNγ production is still unknown. In the present study, the negative correlation between miR-29b and both IFNγs in immune tissues of healthy E. coioides and grouper spleen cells (GS cells) stimulated with LPS or poly I:C was demonstrated. Moreover, dual-luciferase reporter assays and western blotting were performed to demonstrate that miR-29b suppressed E. coioides IFNγ production. Studies of NO production in GS cells after miR-29b transfection revealed that miR-29b overexpression affected NO production through the downregulation of IFNγ expression. Taken together, our results suggest that miR-29b may directly target E. coioides IFNγs and modulate IFNγ-mediated innate immune responses by suppressing E. coioides IFNγs production.

The Role of APOSTART in Switching between Sexuality and Apomixis in Poa pratensis

The production of seeds without sex is considered the holy grail of plant biology. The transfer of apomixis to various crop species has the potential to transform plant breeding, since it will allow new varieties to retain valuable traits thorough asexual reproduction. Therefore, a greater molecular understanding of apomixis is fundamental. In a previous work we identified a gene, namely APOSTART, that seemed to be involved in this asexual mode of reproduction, which is very common in Poa pratensis L., and here we present a detailed work aimed at clarifying its role in apomixis. In situ hybridization showed that PpAPOSTART is expressed in reproductive tissues from pre-meiosis to embryo development. Interestingly, it is expressed early in few nucellar cells of apomictic individuals possibly switching from a somatic to a reproductive cell as in aposporic apomixis. Moreover, out of 13 APOSTART members, we identified one, APOSTART_6, as specifically expressed in flower tissue. APOSTART_6 also exhibited delayed expression in apomictic genotypes when compared with sexual types. Most importantly, the SCAR (Sequence Characterized Amplified Region) derived from the APOSTART_6 sequence completely co-segregated with apomixis.

Transcriptome-based evaluation and validation of suitable housekeeping gene for quantification real-time PCR under specific experiment condition in teleost fishes

Quantification real-time PCR (qRT-PCR) is a common method in analysis of gene expression, but the stable reference genes for the normalization analysis have not been appreciated before identifying expression pattern of genes in teleost fishes. In this study, we selected eight candidate reference genes (18S, Actin, EF-1α, 40S, B2M, TUBA, UBCE and GAPDH) basing on transcriptome analysis and the traditional housekeeping genes, and analyzed the stability of the reference genes in spleen, head kidney and head kidney leukocytes (HKL) after pathogen challenge in Schizothorax prenanti (S. prenanti). Three common programs (geNorm, NormFinder and Bestkeeper) were used to evaluate the stability of the candidate reference genes. Two reference genes, Actin and EF-1α presented higher stability, while 18S and GAPDH were the lower stable genes, both in in vitro and in vivo. An important immune gene, toll-like receptor 22a (TLR22a), was selected to validate the stability of the proposed reference genes (Actin and EF-1α) across different experiment treatments. The results reveal that Actin and EF-1α are quite suitable reference genes for the normalization analysis. Otherwise, using the most stable gene Actin to validate the reliable of transcriptome data showed the high correlation between the fold change of transcriptome data and qRT-PCR data. In conclusion, our study not only acquired the suitable reference gene for the qRT-PCR assay under specific experiment condition, but also provided a comprehensive method to evaluate and validate the reference gene based on transcriptome analysis in teleost fishes.

Reference gene validation for quantification of gene expression during ovarian development of turbot (Scophthalmus maximus)

Quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR) is a powerful and sensitive method used in gene expression analysis. Suitable reference genes, which are stable under all experimental circumstances and tissues significantly improve the accuracy of qRT-PCR data. In this study, the stability of six genes, namely, 18S ribosomal RNA (18s), beta-actin (actb), elongation factor 1-alpha (ef1α), glyceraldehyde-3-phosphate-dehydrogenase (gapdh), cathepsin D (ctsd), and beta-2-microglobulin (b2m) were evaluated as potential references for qRT-PCR analysis. The genes were examined in the hypothalamus-pituitary-ovary-liver (HPOL) axis throughout turbot ovarian development via using the geNorm, NormFinder and BestKeeper algorithms. Results showed that the most stable reference genes were ef1α, actb, and ctsd in the hypothalamus, pituitary, ovary and liver, respectively. The best-suited gene combinations for normalization were 18s, ef1α, and ctsd in the hypothalamus; actb, ctsd, and 18s in the pituitary; actb, and ctsd in the ovary; gapdh and ctsd in the liver. Moreover, the expression profile of estrogen receptor α (erα) manifested no significant difference normalization to the aforementioned best-suited gene during turbot ovarian development. However, no single gene or pair of genes is suitable as an internal control and account for the amplification differences among the four tissues during ovarian development. In summary, these results provide a basic data for the optimal reference gene selection and obtain highly accurate normalization of qRT-PCR data in HPOL axis-related gene expression analysis during turbot ovarian development.

Selection of reference genes for quantitative real-time PCR normalisation in largemouth bass Micropterus salmoides fed on alternative diets

The partial cDNA sequences of eight reference genes (actb, tuba1, gapdh58, gapdh59, eef1a1, RNA 18 s, pabpc1, ube2I) were cloned from largemouth bass Micropterus salmoides. The expression levels of these eight genes were compared in the various tissues (eye, spleen, kidney, gill, muscle, brain, liver, heart, gut and gonad) of M. salmoides fed on forage fish. The results showed that the candidate genes exhibited tissue-specific expression to various degrees and the stability ranking order was eef1a1 > tuba1 > RNA 18 s > pabpc1 > ube2I > actb > gapdh58 > gapdh59 among tissue types. Four candidate genes eef1a1, tuba1, RNA 18 s and actb were used to analyse the stability in liver tissues of largemouth bass between the forage-fish group and the formulated-feed group. The candidate genes also showed some changes in expression levels in the livers, while eef1a1 and tuba1 had the most stable expression in livers of fish fed on alternative diets within 10 candidates. So eef1a1 and tuba1 were recommended as optimal reference gene in quantitative real-time PCR analysis to normalise the expression levels of target genes in tissues and lives of the M. salmoides fed on alternative diets. In livers, the expression levels of gck normalised by eef1a1 and tuba1 showed the significant up-regulation in formulated feed group (P < 0.05) than those in forage-fish group. While sex difference has no significant effects on the expression levels of gck in both groups.

β-actin gene expression is variable among individuals and not suitable for normalizing mRNA levels in Portunus trituberculatus

The housekeeping gene encoding β-actin appears to be the most widely-used internal reference for gene expression studies in experimental animals or their cell lines. However, the effectiveness of β-actin to normalize mRNA levels expression in many crustacean species is still object of debate. To date, it is still unclear if β-actin is suitable to be utilized as the internal reference in qualitative real-time gene expression study in crab species. To address this concern, we evaluated 5 candidate reference genes encoding β-actin, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), cyclophilin A, elongation factor 1-α (EF1-α), and 18 S ribosomal RNA (18 S rRNA) in the swimming crabs (Portunus trituberculatus) models. Our data showed that the β-actin gene expression varied significantly across individual swimming crab individuals in gills or hemocytes and the expression of 18 S rRNA, EF1-α, cyclophilin or GAPDH gene were relatively stable compared to that of β-actin. Moreover, the expression stability of the reference genes among different tissues in normal crabs or after WSSV challenge was also tested by geNorm and NormFinder software. Among tissues, 18 S rRNA was most stably expressed in different tissues, followed by cyclophilin A and EF1-α, compared to β-actin and GAPDH. Upon to viral simulation, GAPDH was found to be the most stable internal control gene in gills and cyclophilin A was ranked as the most stable gene in hemocytes.

Salmon louse rhabdoviruses: Impact on louse development and transcription of selected Atlantic salmon immune genes

Recently, it has been shown that the salmon louse (Lepeophtheirus salmonis) is commonly infected by one or two vertically transmitted Lepeophtheirus salmonis rhabdoviruses (LsRVs). As shown in the present study, the viruses have limited effect on louse survival, developmental rate and fecundity. Since the LsRVs were confirmed to be present in the louse salivary glands, the salmon cutaneous immune response towards LsRV positive and negative lice was analyzed. In general, L. salmonis increased the expression of IL1β, IL8 and IL4/13A at the attachment site, in addition to the non-specific cytotoxic cell receptor protein 1 (NCCRP-1). Interestingly, LsRV free lice induced a higher skin expression of IL1β, IL8, and NCCRP-1 than the LsRV infected lice. The inflammatory response is important for louse clearance, and the present results suggest that the LsRVs can be beneficial for the louse by dampening inflammation. Further research is, however; needed to ascertain whether this is a direct modulatory effect of secreted virions, or if virus replication is altering the level of louse salivary gland proteins.

The brain-gut axis of fish: Rainbow trout with low and high cortisol response show innate differences in intestinal integrity and brain gene expression

In fish, the stress hormone cortisol is released through the action of the hypothalamic pituitary interrenal axis (HPI-axis). The reactivity of this axis differs between individuals and previous studies have linked this to different behavioural characteristics and stress coping styles. In the current study, low and high responding (LR and HR) rainbow trout in terms of cortisol release during stress were identified, using a repeated confinements stress test. The expression of stress related genes in the forebrain and the integrity of the stress sensitive primary barrier of the intestine was examined. The HR trout displayed higher expression levels of mineralocorticoid and serotonergic receptors and serotonergic re-uptake pumps in the telencephalon during both basal and stressed conditions. This confirms that HPI-axis reactivity is linked also to other neuronal behavioural modulators, as both the serotonergic and the corticoid system in the telencephalon are involved in behavioural reactivity and cognitive processes. Involvement of the HPI-axis in the brain-gut-axis was also found. LR trout displayed a lower integrity in the primary barrier of the intestine during basal conditions compared to the HR trout. However, following stress exposure, LR trout showed an unexpected increase in intestinal integrity whereas the HR trout instead suffered a reduction. This could make the LR individuals more susceptible to pathogens during basal conditions where instead HR individuals would be more vulnerable during stressed conditions. We hypothesize that these barrier differences are caused by regulation/effects on tight junction proteins possibly controlled by secondary effects of cortisol on the intestinal immune barrier or differences in parasympathetic reactivity.

Association of genetic variants and expression levels of porcine FABP4 and FABP5 genes

The FABP4 and FABP5 genes, coding for fatty acid transport proteins, have long been studied as positional candidate genes for SSC4 QTL affecting fat deposition and composition traits in pigs. Polymorphisms in these genes, FABP4:g.2634_2635insC and FABP5:g.3000T>G, have previously been associated with fatness traits in an Iberian by Landrace cross (IBMAP). The aim of the present work was to evaluate the functional implication of these genetic variants. For this purpose, FABP4 and FABP5 mRNA expression levels in 114 BC1_LD animals (25% Iberian × 75% Landrace) were analyzed using real-time quantitative PCR in backfat and muscle. FABP4 gene expression in backfat, but not in muscle, was associated with FABP4:g.2634_2635insC. In contrast, FABP5:g.3000T>G was not associated with gene expression levels. An expression-based genome-wide association study highlighted the FABP4:g.2634_2635insC polymorphism as the polymorphism most associated with FABP4 gene expression in backfat. Furthermore, other genomic regions associated in trans with the mRNA expression of FABP4 in backfat and FABP5 in muscle were also identified. Finally, two putative transcription binding sites for PPARG and NR4A2 may be affected by the FABP4:g.2634_2635insC polymorphism, modifying FABP4 gene expression. Our results reinforce FABP4 as a candidate gene for fatness traits on SSC4.

RNAi-mediated treatment of two vertically transmitted rhabdovirus infecting the salmon louse (Lepeophtheirus salmonis)

Rhabdoviruses are a family of enveloped negative-sense single-stranded RNA viruses infecting a variety of hosts. Recently, two vertically transmitted salmon louse (Lepeophtheirus salmonis) rhabdoviruses (LsRV) have been identified. The prevalence of these viruses was measured along the Norwegian coast and found to be close to 100%, and with the present lack of suitable cell lines to propagate these viruses, it is challenging to obtain material to study their host impact and infection routes. Thus, virus free lice strains were established from virus infected lice carrying one or both LsRVs by treating them with N protein dsRNA twice during development. The viral replication of the N protein was specifically down-regulated following introduction of virus-specific dsRNA, and virus-free lice strains were maintained for several generations. A preliminary study on infection routes suggested that the LsRV-No9 is maternally transmitted, and that the virus transmits from males to females horizontally. The ability to produce virus free strains allows for further studies on transmission modes and how these viruses influences on the L.salmonis interaction with its salmonid host. Moreover, this study provides a general fundament for future studies on how vertically transmitted rhabdoviruses influence the biology of their arthropod hosts.

Expression stability and selection of optimal reference genes for gene expression normalization in early life stage rainbow trout exposed to cadmium and copper

Gene expression analysis represents a powerful approach to characterize the specific mechanisms by which contaminants interact with organisms. One of the key considerations when conducting gene expression analyses using quantitative real-time reverse transcription-polymerase chain reaction (qPCR) is the selection of appropriate reference genes, which is often overlooked. Specifically, to reach meaningful conclusions when using relative quantification approaches, expression levels of reference genes must be highly stable and cannot vary as a function of experimental conditions. However, to date, information on the stability of commonly used reference genes across developmental stages, tissues and after exposure to contaminants such as metals is lacking for many vertebrate species including teleost fish. Therefore, in this study, we assessed the stability of expression of 8 reference gene candidates in the gills and skin of three different early life-stages of rainbow trout after acute exposure (24h) to two metals, cadmium (Cd) and copper (Cu) using qPCR. Candidate housekeeping genes were: beta actin (b-actin), DNA directed RNA polymerase II subunit I (DRP2), elongation factor-1 alpha (EF1a), glyceraldehyde 3-phosphate dehydrogenase (GAPDH), glucose-6-phosphate dehydrogenase (G6PD), hypoxanthine phosphoribosyltransferase (HPRT), ribosomal protein L8 (RPL8), and 18S ribosomal RNA (18S). Four algorithms, geNorm, NormFinder, BestKeeper, and the comparative ΔCt method were employed to systematically evaluate the expression stability of these candidate genes under control and exposed conditions as well as across three different life-stages. Finally, stability of genes was ranked by taking geometric means of the ranks established by the different methods. Stability of reference genes was ranked in the following order (from lower to higher stability): HPRT<GAPDH<EF1a<G6PD<RPL8<DRP2<b-actin in gills of fish exposed to Cd; b-actin<GAPDH<G6PD<DRP2<RPL8<HPRT<EF1a in gills of fish exposed to Cu; RPL8<HPRT<GAPDH<G6PD<EF1a<DRP2<b-actin in the skin of fish exposed to Cd; and EF1a<GAPDH<RPL8<HPRT<G6PD<b-actin<DRP2 in the skin of fish exposed to Cu. The results demonstrated that the stability of reference genes depended on the metal, life-stage and/or organ in question. Thus, attention should be paid to these factors before selection of reference gene for relative quantification of the gene expressions.

Appropriateness of reference genes for normalizing messenger RNA in mouse 2,4-dinitrobenzene sulfonic acid (DNBS)-induced colitis using quantitative real time PCR

2,4-Dinitrobenzene sulfonic acid (DNBS)-induced colitis is an experimental model that mimics Crohn's disease. Appropriateness of reference genes is crucial for RT-qPCR. This is the first study to determine the stability of reference gene expression (RGE) in mice treated with DNBS. DNBS experimental Colitis was induced in male C57BL/6 mice. RNA was extracted from colon tissue and comprehensive analysis of 13 RGE was performed according to predefined criteria. Relative colonic TNF-α and IL-1β mRNA levels were calculated. Colitis significantly altered the stability of mucosal RGE. Commonly used glyceraldehyde-3-phosphate dehydrogenase (Gapdh), β-actin (Actb), or β2-microglobulin (β2m) showed the highest fluctuation within the inflamed and control groups. Conversely, ribosomal protein large P0 (Rplp0), non-POU domain containing (Nono), TATA-box-binding protein (Tbp) and eukaryotic translation elongation factor 2 (Eef2) were not affected by inflammation and were the most stable genes. TNF-α and IL-1β mRNA levels was dependent on the reference gene used and varied from significant when the most stable genes were used to non-significant when the least stable genes were used. The appropriate choice of RGE is critical to guarantee satisfactory normalization of RT-qPCR data when using DNBS-Model. We recommend using Rplp0, Nono, Tbp, Hprt and Eef2 instead of common reference genes.

Selection and evaluation of new reference genes for RT-qPCR analysis in Epinephelus akaara based on transcriptome data

Groupers are an economically important fish species in world fishery markets. Because many studies using RT-qPCR have addressed gene expression in groupers, appropriate reference genes are required to obtain reliable and accurate results. In this study, the most suitable reference genes were identified from eleven candidate genes of one of the most valuable species, Epinephelus akaara, in a range of different experimental conditions. Using the software packages geNorm, NormFinder, BestKeeper and refFinder, three traditionally used reference genes, β-actin (β-ACT), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and beta-2-microglobulin (B2M), were identified as not suitable for E. akaara gene expression studies, whereas two newly identified reference genes, conserved oligomeric Golgi complex subunit 5 (Cog5) and brefeldin a-inhibited guanine nucleotide-exchange protein 1 (ARFGEF1), could be universally applied under all the tested conditions. These data provide the foundation for more precise results in RT-qPCR studies of gene expression in E. akaara and other Epinephelus species.

Pathogen-dependent role of turbot (Scophthalmus maximus) interferon-gamma

Interferon-gamma has been typically described as a pro-inflammatory cytokine playing an important role in the resolution of both viral and bacterial diseases. Nevertheless, some anti-inflammatory functions are also attributed to this molecule. In this work we have characterized for the first time the turbot (Scophthalmus maximus) interferon-gamma gene (ifng) and its expression pattern under basal conditions, after type I IFNs administration, and viral and bacterial infection. The intramuscular injection of an expression plasmid encoding turbot Ifng (pMCV1.4-ifng) was not able to affect the transcription of numerous immune genes directly related to the activity of IFN-gamma, with the exception of macrophage-colony stimulating factor (csf1). It was also unable to reduce the mortality caused by a Viral Hemorrhagic Septicemia Virus (VHSV) or Aeromonas salmonicida challenge. Interestingly, at 24 h post-infection, turbot previously inoculated with pMCV1.4-ifng and infected with VHSV showed an increase in the expression of pro-inflammatory cytokines and type I IFNs compared to those fish not receiving expression plasmid, indicating a synergic effect of Ifng and VHSV. On the other hand, some macrophage markers, such as the macrophage receptor with collagenous structure (marco), were down-regulated by Ifng during the viral infection. Ifng had the opposite effect in those turbot infected with the bacteria, showing a reduction in the transcription of pro-inflammatory and type I IFNs genes, and an increase in the expression of genes related to the activity of macrophages.

Sequence Evolution and Expression of the Androgen Receptor and Other Pathway-Related Genes in a Unisexual Fish, the Amazon Molly, Poecilia formosa, and Its Bisexual Ancestors

The all-female Amazon molly (Poecilia formosa) originated from a single hybridization of two bisexual ancestors, Atlantic molly (Poecilia mexicana) and sailfin molly (Poecilia latipinna). As a gynogenetic species, the Amazon molly needs to copulate with a heterospecific male, but the genetic information of the sperm-donor does not contribute to the next generation, as the sperm only acts as the trigger for the diploid eggs' embryogenesis. Here, we study the sequence evolution and gene expression of the duplicated genes coding for androgen receptors (ars) and other pathway-related genes, i.e., the estrogen receptors (ers) and cytochrome P450, family19, subfamily A, aromatase genes (cyp19as), in the Amazon molly, in comparison to its bisexual ancestors. Mollies possess-as most other teleost fish-two copies of the ar, er, and cyp19a genes, i.e., arα/arβ, erα/erβ1, and cyp19a1 (also referred as cyp19a1a)/cyp19a2 (also referred to as cyp19a1b), respectively. Non-synonymous single nucleotide polymorphisms (SNPs) among the ancestral bisexual species were generally predicted not to alter protein function. Some derived substitutions in the P. mexicana and one in P. formosa are predicted to impact protein function. We also describe the gene expression pattern of the ars and pathway-related genes in various tissues (i.e., brain, gill, and ovary) and provide SNP markers for allele specific expression research. As a general tendency, the levels of gene expression were lowest in gill and highest in ovarian tissues, while expression levels in the brain were intermediate in most cases. Expression levels in P. formosa were conserved where expression did not differ between the two bisexual ancestors. In those cases where gene expression levels significantly differed between the bisexual species, P. formosa expression was always comparable to the higher expression level among the two ancestors. Interestingly, erβ1 was expressed neither in brain nor in gill in the analyzed three molly species, which implies a more important role of erα in the estradiol synthesis pathway in these tissues. Furthermore, our data suggest that interactions of steroid-signaling pathway genes differ across tissues, in particular the interactions of ars and cyp19as.

Ribosomal protein L7 as a suitable reference gene for quantifying gene expression in gastropod Bellamya aeruginosa

Expression levels of eight candidate reference genes were quantified in tissues of gastropod Bellamya aeruginosa exposed for 10 d to various stressors, including fasting, 17β-estradiol, 17α-methyltestosterone, and Cd(2+). The results showed that 18s rRNA was the most highly expressed of the candidate reference genes, while H2A was the least expressed. There were no significant changes (p>0.05) in the expression of the eight genes in tissues among the different treatments. Using RefFinder to evaluate the expression stabilities of the eight candidate reference genes, ribosomal protein was shown to be the most stable reference gene, and no effects were observed among the different stressor treatments. These results indicate that RPL 7 is the most suitable reference gene for quantifying gene expression in B. aeruginosa under environmental stress, which was verified in B. aeruginosa exposed to high doses of E2 for 24 and 72h.

Developmental profiles and expression of the DNA methyltransferase genes in the fathead minnow (Pimephales promelas) following exposure to di-2-ethylhexyl phthalate

DNA methylation is an epigenetic regulator of gene expression, and this process has been shown to be disrupted by environmental contaminants. Di-2-(ethylhexyl) phthalate (DEHP) and related phthalate esters have been shown to affect development in early life stages of fish and can alter genomic methylation patterns in vertebrates. The objectives of this study were the following: (1) Describe the expression patterns of the DNA methyltransferase (dnmt) genes during early fathead minnow (FHM) development. These genes are critical for methylation and imprinting during development. (2) Determine the effects of DEHP on the development of FHM larvae [1 and 14 days post-hatch (dph)]. (3) Determine the effect of DEHP on dnmt expression and global methylation status in larval FHM. FHMs were first collected over a developmental time course [1, 3, 5, 6, and 14 days post-fertilization (dpf)] to investigate the expression patterns of five dnmt isoforms. The expression of dnmt1 and dnmt7 was relatively high in embryos at 1 dpf but was variable in expression, and these transcripts were later expressed at a lower level (>3 dpf); dnmt3 was significantly higher in embryos at 1 dpf compared to those at 3 dpf. Dnmt6 showed more of a constitutive pattern of expression during the first 2 weeks of development, and the mRNA levels of dnmt8 were higher in embryos at 5 and 6 dpf compared to those at 1 and 3 dpf, corresponding to the hatching period of the embryos. A waterborne exposure to three concentrations of DEHP (1, 10 and 100 µg/L) was conducted on 1-day FHM embryos for 24 h and on larval fish for 2 weeks, ending at 14 dpf. DEHP did not negatively affect survival, hatch rate, or the expression of dnmt isoforms in FHMs. There were no differences in global cytosine methylation following DEHP treatments in 14 dpf larvae, suggesting that environmentally relevant levels of DEHP may not affect global methylation at this stage of FHM development. However, additional targeted methylome studies are required to determine whether specific gene promoters are differently methylated following exposure to DEHP. This study offers new insight into the roles of the dnmt enzymes during FHM development.

Stability evaluation of reference genes for real-time PCR in zebrafish (Danio rerio) exposed to cadmium chloride and subsequently infected by bacteria Aeromonas hydrophila

Environmental and occupational cadmium (Cd) toxicity is a global concern, and the model organism zebrafish is an ideal species to investigate Cd toxicity. Among various detecting techniques, quantitative real-time PCR (qPCR) is a sensitive and efficient tool. Stable reference genes are critical for relative qPCR analysis. However, accumulated evidence shows that conventional reference genes can vary significantly under different experimental setups. Here we evaluated the stability of eight candidate reference genes of zebrafish with or without exposure to different concentrations of Cd. The results showed that the best four suitable reference genes in the five selected organs were: (1) spleen: β-actin>gapdh>ef1α>rpl13α; (2) kidney: rplp2>rpl7>β-actin>ef1α; (3) liver: rpl7>rpl13α>β-actin>ef1α; (4) gills: rplp2>gapdh>rnf7>ef1α; (5) intestine: ef1α>rnf7>rplp2>rpl13α. Moreover, we further assessed the expression stability of the four reference genes for Cd immunotoxicology studies in zebrafish. The expression profiles showed that ef1α in spleen and kidney, rpl13a in liver and rplp2 in intestine were the most suitable reference genes at 12h and 9 days after the injection with Aeromonas hydrophila following Cd exposure. In gills, the expression of gapdh was more stable than ef1α after 9 days of bacteria challenge while ef1α showed a higher stability than gapdh at 12h after bacteria injection. In conclusion, this study has demonstrated that different tissues of zebrafish have different suitable reference genes after Cd exposure and the subsequently pathogenic insults for qPCR. It emphasized the importance of reference gene evaluation for studies using qPCR, in particular when investigations involve factors not explored previously.

Evaluation of potential internal references for quantitative real-time RT-PCR normalization of gene expression in red drum (Sciaenops ocellatus)

Quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) has been used extensively for studying gene expression in diverse organisms including fish. In this study, with an aim to identify reliable reference genes for qRT-PCR in red drum (Sciaenops ocellatus), an economic fish species, we determined the expression stability of seven housekeeping genes in healthy and bacterium-infected red drum. Each of the selected candidate genes was amplified by qRT-PCR from the brain, gill, heart, intestine, kidney, liver, muscle, and spleen of red drum infected with or without a bacterial pathogen for 12 and 48 h. The mRNA levels of the genes were analyzed with the geNorm and NormFinder algorithms. The results showed that in the absence of bacterial infection, translation initiation factor 3, NADH dehydrogenase 1, and QM-like protein may be used together as internal references across the eight examined tissues. Bacterial infection caused variations in the rankings of the most stable genes in a tissue-dependent manner. For all tissues, two genes sufficed for reliable normalization at both 12 and 48 h post-infection. However, the optimal gene pairs differed among tissues and, for four of the examined eight tissues, between infection points. These results indicate that when studying gene expression in red drum under conditions of bacterial infection, the optimal reference genes should be selected on the basis of tissue type and, for accurate normalization, infection stage.

A new aquareovirus causing high mortality in farmed Atlantic halibut fry in Norway

A new aquareovirus was isolated from cultured Atlantic halibut (Hippoglossus hippoglossus) fry at a facility where massive mortalities had occurred during the start-feeding phase. The same virus was also detected in juveniles (about 10 grams) of the 2013 generation at two other production sites, but not in larger fish from generations 2007-2012. The virus replicated in BF-2 and CHSE-214 cell cultures and produced syncytia and plaque-like cytopathic effects. This Atlantic halibut reovirus (AHRV) was associated with necrosis of the liver and pancreas, syncytium formation in these tissues, and distinct viroplasm areas within the syncytium in halibut fry. Transmission electron microscopy revealed that the viroplasm contained virions, non-enveloped, icosahedral particles approximately 70 nm in diameter with a double capsid layer, amorphous material, and tubular structures. The RNA-dependent RNA polymerase (RdRp) gene from the AHRV isolates showed the highest amino acid sequence identity (80 %) to an isolate belonging to the species Aquareovirus A, Atlantic salmon reovirus TS (ASRV-TS). A partial sequence from the putative fusion-associated small transmembrane (FAST) protein of AHRV was obtained, and this sequence showed the highest amino acid sequence identity (46.8 %) to Green River Chinook virus which is an unassigned member of the genus Aquareovirus, while a comparison with isolates belonging to the species Aquareovirus A showed <33 % identity. A proper assessment of the relationship of AHRV to all members of the genus Aquareovirus, however, is hampered by the absence of genetic data from members of several Aquareovirus species. AHRV is the first aquareovirus isolated from a marine coldwater fish species and the second reovirus detected in farmed fish in Norway. A similar disease of halibut fry, as described in this paper, has also been described in halibut production facilities in Canada and Scotland.

Reference gene selection for real-time RT-PCR normalization in rice field eel (Monopterus albus) during gonad development

Real-time reverse transcriptase (RT) polymerase chain reaction (PCR) requires data normalization using an appropriate reference gene in order to obtain more reliable results with biological significance. We cloned a partial sequence of elongation factor-1-α (EF1α) and ribosomal protein L17 (RPL17) from Monopterus albus. We investigated the suitability of five commonly used reference genes [18S ribosomal RNA (18S), cytoskeletal protein (β-actin), glyceraldehyde phosphate dehydrogenase (GAPDH), EF1α and RPL17] as potential quantitative reference genes for normalizing real-time RT-PCR data generated in gonads of different developmental stages and in other tissues of M. albus. Analysis of the data indicated that 18S, β-actin and GAPDH are not suitable as reference genes because of their levels of variations of expression. EF1α and RPL17 might be suitable as reference genes in the gonads of different developmental stages as well as in other tissues of M. albus.

Cyclic LIPopeptides from Bacillus subtilis ABS-S14 elicit defense-related gene expression in citrus fruit

Effects of cyclic lipopeptides (CLPs) obtained from Bacillus subtilis ABS-S14 on eliciting defense-related gene transcription and activity of defense-related enzymes; glucanase (GLU), chitinase (CHI), peroxidase (POX) and lipoxygenase (LOX) in Citrus sinensis cv. Valencia fruit were determined. The maximum level of GLU transcripts induced in fruit treated with fengycin was significantly greatest among treatments at 48 h. Surfactin enhanced the LOX and POX transcripts. In parallel, corresponding enzyme activities were correlated with changes in gene expression observed in fruit inoculated with Penicillium digitatum following treatment with individual CLPs. Synergistic effects of fengycin and iturin A, fengycin and surfactin were shown in gene transcript of GLU and CHI, respectively, and surfactin induced POX and LOX gene expression of citrus flavedo without pathogen infection. These results suggest that fengycin and surfactin act as elicitors of defense-related gene expression in “Valencia” fruit following infection.

Identification of a novel RNA virus lethal to tilapia

Tilapines are important for the sustainability of ecological systems and serve as the second most important group of farmed fish worldwide. Significant mortality of wild and cultured tilapia has been observed recently in Israel. The etiological agent of this disease, a novel RNA virus, is described here, and procedures allowing its isolation and detection are revealed. The virus, denominated tilapia lake virus (TiLV), was propagated in primary tilapia brain cells or in an E-11 cell line, and it induced a cytopathic effect at 5 to 10 days postinfection. Electron microscopy revealed enveloped icosahedral particles of 55 to 75 nm. Low-passage TiLV, injected intraperitoneally in tilapia, induced a disease resembling the natural disease, which typically presents with lethargy, ocular alterations, and skin erosions, with >80% mortality. Histological changes included congestion of the internal organs (kidneys and brain) with foci of gliosis and perivascular cuffing of lymphocytes in the brain cortex; ocular inflammation included endophthalmitis and cataractous changes of the lens. The cohabitation of healthy and diseased fish demonstrated that the disease is contagious and that mortalities (80 to 100%) occur within a few days. Fish surviving the initial mortality were immune to further TiLV infections, suggesting the mounting of a protective immune response. Screening cDNA libraries identified a TiLV-specific sequence, allowing the design of a PCR-based diagnostic test. This test enables the specific identification of TiLV in tilapines and should help control the spread of this virus worldwide.

Analysis of qPCR reference gene stability determination methods and a practical approach for efficiency calculation on a turbot (Scophthalmus maximus) gonad dataset

Background

Gene expression analysis by reverse transcription quantitative PCR (qPCR) is the most widely used method for analyzing the expression of a moderate number of genes and also for the validation of microarray results. Several issues are crucial for a successful qPCR study, particularly the selection of internal reference genes for normalization and efficiency determination. There is no agreement on which method is the best to detect the most stable genes neither on how to perform efficiency determination. In this study we offer a comprehensive evaluation of the characteristics of reference gene selection methods and how to decide which one is more reliable when they show discordant outcomes. Also, we analyze the current efficiency calculation controversy. Our dataset is composed by gonad samples of turbot at different development times reared at different temperatures. Turbot (Scophthalmus maximus) is a relevant marine aquaculture European species with increasing production in the incoming years. Since females largely outgrow males, identification of genes related to sex determination, gonad development and reproductive behavior, and analysis of their expression profiles are of primary importance for turbot industry.

Results

We analyzed gene stability of six reference genes: RPS4, RPL17, GAPDH, ACTB, UBQ and B2M using the comparative delta-CT method, Bestkeeper, NormFinder and GeNorm approaches in gonad samples of turbot. Supported by descriptive statistics, we found NormFinder to be the best method, while on the other side, GeNorm results proved to be unreliable. According to our analysis, UBQ and RPS4 were the most stable genes, while B2M was the least stable gene. We also analyzed the efficiency calculation softwares LinRegPCR, LREanalyzer, DART and PCR-Miner and we recommend LinRegPCR for research purposes since it does not systematically overestimate efficiency.

Conclusion

Our results indicate that NormFinder and LinRegPCR are the best approaches for reference gene selection and efficiency determination, respectively. We also recommend the use of UBQ and RPS4 for normalization of gonad development samples in turbot.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-648) contains supplementary material, which is available to authorized users.

Selection of Reliable Reference Genes for Real-time qRT-PCR Analysis of Zi Geese (Anser anser domestica) Gene Expression

Zi geese (Anser anser domestica) belong to the white geese and are excellent layers with a superior feed-to-egg conversion ratio. Quantitative gene expression analysis, such as Real-time qRT-PCR, will provide a good understanding of ovarian function during egg-laying and consequently improve egg production. However, we still don’t know what reference genes in geese, which show stable expression, should be used for such quantitative analysis. In order to reveal such reference genes, the stability of seven genes were tested in five tissues of Zi geese. Methodology/Principal Findings: The relative transcription levels of genes encoding hypoxanthine guanine phosphoribosyl transferase 1 (HPRT1), β-actin (ACTB), β-tubulin (TUB), glyceraldehyde-3-phosphate-dehydrogenase (GADPH), succinate dehydrogenase flavoprotein (SDH), 28S rRNA (28S) and 18S rRNA (18S) have been quantified in heart, liver, kidney, muscle and ovary in Zi geese respectively at different developmental stages (1 d, 2, 4, 6 and 8 months). The expression stability of these genes was analyzed using geNorm, NormFinder and BestKeeper software. Conclusions: The expression of 28S in heart, GAPDH in liver and ovary, ACTB in kidney and HPRT1 in muscle are the most stable genes as identified by the three different analysis methods. Thus, these genes are recommended for use as candidate reference genes to compare mRNA transcription in various developmental stages of geese.

A multicopy Y-chromosomal SGNH hydrolase gene expressed in the testis of the platyfish has been captured and mobilized by a Helitron transposon

Background

Teleost fish present a high diversity of sex determination systems, with possible frequent evolutionary turnover of sex chromosomes and sex-determining genes. In order to identify genes involved in male sex determination and differentiation in the platyfish Xiphophorus maculatus, bacterial artificial chromosome contigs from the sex-determining region differentiating the Y from the X chromosome have been assembled and analyzed.

Results

A novel three-copy gene called teximY (for testis-expressed in Xiphophorus maculatus on the Y) was identified on the Y but not on the X chromosome. A highly related sequence called texim1, probably at the origin of the Y-linked genes, as well as three more divergent texim genes were detected in (pseudo)autosomal regions of the platyfish genome. Texim genes, for which no functional data are available so far in any organism, encode predicted esterases/lipases with a SGNH hydrolase domain. Texim proteins are related to proteins from very different origins, including proteins encoded by animal CR1 retrotransposons, animal platelet-activating factor acetylhydrolases (PAFah) and bacterial hydrolases. Texim gene distribution is patchy in animals. Texim sequences were detected in several fish species including killifish, medaka, pufferfish, sea bass, cod and gar, but not in zebrafish. Texim-like genes are also present in Oikopleura (urochordate), Amphioxus (cephalochordate) and sea urchin (echinoderm) but absent from mammals and other tetrapods. Interestingly, texim genes are associated with a Helitron transposon in different fish species but not in urochordates, cephalochordates and echinoderms, suggesting capture and mobilization of an ancestral texim gene in the bony fish lineage. RT-qPCR analyses showed that Y-linked teximY genes are preferentially expressed in testis, with expression at late stages of spermatogenesis (late spermatids and spermatozeugmata).

Conclusions

These observations suggest either that TeximY proteins play a role in Helitron transposition in the male germ line in fish, or that texim genes are spermatogenesis genes mobilized and spread by transposable elements in fish genomes.

Reference gene selection for quantitative real-time RT-PCR normalization in the half-smooth tongue sole (Cynoglossus semilaevis) at different developmental stages, in various tissue types and on exposure to chemicals

Quantitative real time RT-PCR has been described as the most sensitive method for the detection of low abundance mRNA. To date, no reference genes have been screened in the half-smooth tongue sole (Cynoglossus semilaevis). The aim of this study was to select the most stable genes for quantitative real-time RT-PCR. Eight housekeeping genes (18S, TUBA, B2M, ACTB, EF1A, GAPDH, RPL17 and UBCE) were tested at different developmental stages, in different tissues, and following exposure to the drug SB-431542. Using geNorm, BestKeeper and NormFinder software, GAPDH/B2M, GAPDH/18S and UBCE/GAPDH were identified as the most suitable genes from samples taken of different developmental stages while 18S/RPL17 were consistently ranked as the best reference genes for different tissue types. Furthermore, TUBA/B2M, TUBA/UBCE and B2M/TUBA were found to be the most suitable genes in samples treated with the drug, SB-431542 by geNorm, BestKeeper and NormFinder respectively. Across both different developmental stages and tissue types, the combination of 18S and GAPDH was the most stable reference gene analyzed by Ref-Finder. To test and verify the screened reference genes, the expression profiles of LEFTY-normalized to the combination of GAPDH/18S and ACTB were presented. These results will be useful for future gene-expression studies in the half-smooth tongue sole.

Quantitative real time RT-PCR study of pathogen-induced gene expression in rock bream (Oplegnathus fasciatus): internal controls for data normalization

Rock bream Oplegnathus fasciatus is an important economic fish species. In this study, we evaluated the appropriateness of six housekeeping genes as internal controls for quantitative real-time PCR (RT-qPCR) analysis of gene expression in rock bream before and after pathogen infection. The expression of the selected genes in eight tissues infected with Vibrio alginolyticus or megalocytivirus was determined by RT-qPCR, and the PCR data were analyzed with geNorm and NormFinder algorithms. The results showed that before pathogen infection, mediator of RNA polymerase II transcription subunit 8 and β-actin were ranked as the most stable genes across the examined tissues. After bacterial or viral infection, the stabilities of the housekeeping genes varied to significant extents in tissue-dependent manners, and no single pair of genes was identified as suitable references for all tissues for either of the pathogen stimuli. In addition, for the majority of tissues, the most stable genes during bacterial infection differed from those during viral infection. Nevertheless, optimum reference genes were identified for each tissue under different conditions. Taken together, these results indicate that tissue type and the nature of the infectious agent used in the study can all influence the choice of normalization factors, and that the optimum reference genes identified in this study will provide a useful guidance for the selection of internal controls in future RT-PCR study of gene expression in rock bream.

Evaluation of reference genes for quantitative real-time RT-PCR analysis of gene expression in Nile tilapia (Oreochromis niloticus)

Quantitative real-time reverse-transcriptase polymerase chain reaction (RT-qPCR) has been used frequently to study gene expression related to fish immunology. In such studies, a stable reference gene should be selected to correct the expression of the target gene. In this study, seven candidate reference genes (glyceraldehyde-3-phosphate dehydrogenase (GADPH), ubiquitin-conjugating enzyme (UBCE), 18S ribosomal RNA (18S rRNA), beta-2-microglobulin (B2M), elongation factor 1 alpha (EF1A), tubulin alpha chain-like (TUBA) and beta actin (ACTB)), were selected to analyze their stability and normalization in seven tissues (liver, spleen, kidney, brain, heart, muscle and intestine) of Nile tilapia (Oreochromis niloticus) challenged with Streptococcus agalactiae or Streptococcus iniae, respectively. The results showed that all the candidate reference genes exhibited tissue-dependent transcriptional variations. With PBS injection as a control, UBCE was the most stable and suitable single reference gene in the intestine, liver, brain, kidney, and spleen after S. iniae infection, and in the liver, kidney, and spleen after S. agalactiae infection. EF1A was the most suitable in heart and muscle after S. iniae or S. agalactiae infection. GADPH was the most suitable gene in intestine and brain after S. agalactiae infection. In normal conditions, UBCE and 18S rRNA were the most stably expressed genes across the various tissues. These results showed that for RT-qPCR analysis of tilapia, selecting two or more reference genes may be more suitable for cross-tissue analysis of gene expression.

Validation of reference genes for expression studies during craniofacial development in arctic charr

Arctic charr (Salvelinus alpinus) is a highly polymorphic species and in Lake Thingvallavatn, Iceland, four phenotypic morphs have evolved. These differences in morphology, especially in craniofacial structures are already apparent during embryonic development, indicating that genes important in the formation of the craniofacial features are expressed differentially between the morphs. In order to generate tools to examine these expression differences in Arctic charr, the aim of the present study was to identify reference genes for quantitative real-time PCR (qPCR). The specific aim was to select reference genes which are able to detect very small expression differences among different morphs. We selected twelve candidate reference genes from the literature, identified corresponding charr sequences using data derived from transcriptome sequencing (RNA-seq) and examined their expression using qPCR. Many of the candidate reference genes were found to be stably expressed, yet their quality-rank as reference genes varied considerably depending on the type of analysis used. In addition to commonly used software for reference gene validation, we used classical statistics to evaluate expression profiles avoiding a bias for reference genes with similar expression patterns (co-regulation). Based on these analyses we chose three reference genes, ACTB, UB2L3 and IF5A1 for further evaluation. Their consistency was assessed in an expression study of three known craniofacially expressed genes, sparc (or osteonectin), matrix metalloprotease 2 (mmp2) and sox9 (sex-determining region Y box 9 protein) using qPCR in embryo heads derived from four charr groups at three developmental time points. The three reference genes were found to be very suitable for studying expression differences between the morphotypes, enabling robust detection of small relative expression changes during charr development. Further, the results showed that sparc and mmp2 are differentially expressed in embryos of different Arctic charr morphotypes.

Atlantic halibut (Hippoglossus hippoglossus L.) T-cell and cytokine response after vaccination and challenge with nodavirus

Viral encephalopathy and retinopathy (VER), caused by nodavirus, is one of the major infectious diseases affecting the marine fish farming industry, yet no effective vaccine is available. In this study, we examined the halibut immune response following administration of an experimental vaccine comprising a recombinant nodavirus capsid protein in combination with an oil adjuvant (OA). Four groups of halibut were injected with either: PBS alone, PBS plus OA, 10μg recCP plus OA, or 50μg recCP plus OA. 15 weeks later, half the fish in each group were challenged with nodavirus and the immune response investigated by analysis of: serum levels of recCP-specific halibut immunoglobulins (Igs), and mRNA transcript levels of several T-cell markers (CD3ɛ, Lck, CD4, CD4-2, CD8α and CD8β) and cytokines (IL-1β, IL-6, IL-12βc and IFNγ). Additionally, the presence of nodaviral RNA2 transcripts in the brains of infected halibut was analysed. After vaccination, the level of IL-6 was consistently elevated in the spleens of fish given injections containing the OA. The combination of recCP and OA increased the expression of IL-1β and IFNγ, as well as the level of recCP-specific Igs in blood plasma. Following challenge with nodavirus, IL-1β and IFNγ were elevated in halibut spleens after 24h in all groups that had received OA with or without recCP antigen. In brain, a general increase in the expression levels of all T-cell markers and IFNγ was observed following challenge with nodavirus. The viral load at 8 weeks post-challenge was lower in the fish that received 50μg recCP, with 5 out of 8 individuals being negative for nodavirus. Additionally, a better correlation between these markers (apart from the CD8 markers), and the viral RNA2 was also observed in this group, suggesting that the activation of CD4+T-cells might be important in reducing the viral load. In conclusion, this study identifies recCP as a promising candidate antigen for the future development of a vaccine against nodavirus.

Polymorphism in the ELOVL6 gene is associated with a major QTL effect on fatty acid composition in pigs

Background

The ELOVL fatty acid elongase 6 (ELOVL6), the only elongase related to de novo lipogenesis, catalyzes the rate-limiting step in the elongation cycle by controlling the fatty acid balance in mammals. It is located on pig chromosome 8 (SSC8) in a region where a QTL affecting palmitic, and palmitoleic acid composition was previously detected, using an Iberian x Landrace intercross. The main goal of this work was to fine-map the QTL and to evaluate the ELOVL6 gene as a positional candidate gene affecting the percentages of palmitic and palmitoleic fatty acids in pigs.

Methodology and Principal Findings

The combination of a haplotype-based approach and single-marker analysis allowed us to identify the main, associated interval for the QTL, in which the ELOVL6 gene was identified and selected as a positional candidate gene. A polymorphism in the promoter region of ELOVL6, ELOVL6:c.-533C>T, was highly associated with the percentage of palmitic and palmitoleic acids in muscle and backfat. Significant differences in ELOVL6 gene expression were observed in backfat when animals were classified by the ELOVL6:c.-533C>T genotype. Accordingly, animals carrying the allele associated with a decrease in ELOVL6 gene expression presented an increase in C16:0 and C16:1(n-7) fatty acid content and a decrease of elongation activity ratios in muscle and backfat. Furthermore, a SNP genome-wide association study with ELOVL6 relative expression levels in backfat showed the strongest effect on the SSC8 region in which the ELOVL6 gene is located. Finally, different potential genomic regions associated with ELOVL6 gene expression were also identified by GWAS in liver and muscle, suggesting a differential tissue regulation of the ELOVL6 gene.

Conclusions and Significance

Our results suggest ELOVL6 as a potential causal gene for the QTL analyzed and, subsequently, for controlling the overall balance of fatty acid composition in pigs.

Selection of normalization factors for quantitative real time RT-PCR studies in Japanese flounder (Paralichthys olivaceus) and turbot (Scophthalmus maximus) under conditions of viral infection

Disease outbreaks caused by iridoviruses are known to affect farmed flounder (Paralichthys olivaceus) and turbot (Scophthalmus maximus). To facilitate quantitative real time RT-PCR (qRT-PCR) analysis of gene expression in flounder and turbot during viral infection, we in this study examined the potentials of 9 housekeeping genes of flounder and turbot as internal references for qRT-PCR under conditions of experimental infection with megalocytivirus, a member of the Iridoviridae family. The mRNA levels of the 9 housekeeping genes in the brain, gill, heart, intestine, kidney, liver, muscle, and spleen of flounder and turbot were determined by qRT-PCR at 24h and 72h post-viral infection, and the expression stabilities of the genes were determined with geNorm and NormFinder algorithms. The results showed that (i) viral infection induced significant changes in the mRNA levels of the all the examined genes in a manner that was dependent on both tissue type and infection stage; (ii) for a given time point of infection, stability predictions made by the two algorisms were highly consistent for most tissues; (iii) the optimum reference genes differed at different infection time points at least in some tissues; (iv) at both examined time points, no common reference genes were identified across all tissue types. These results indicate that when studying gene expression in flounder and turbot in relation to viral infection, different internal references may have to be used not only for different tissues but also for different infection stages.

Liver transcriptome profile in pigs with extreme phenotypes of intramuscular fatty acid composition

BACKGROUND

New advances in high-throughput technologies have allowed for the massive analysis of genomic data, providing new opportunities for the characterization of the transcriptome architectures. Recent studies in pigs have employed RNA-Seq to explore the transcriptome of different tissues in a reduced number of animals. The main goal of this study was the identification of differentially-expressed genes in the liver of Iberian x Landrace crossbred pigs showing extreme phenotypes for intramuscular fatty acid composition using RNA-Seq.

RESULTS

The liver transcriptomes of two female groups (H and L) with phenotypically extreme intramuscular fatty acid composition were sequenced using RNA-Seq. A total of 146 and 180 unannotated protein-coding genes were identified in intergenic regions for the L and H groups, respectively. In addition, a range of 5.8 to 7.3% of repetitive elements was found, with SINEs being the most abundant elements. The expression in liver of 186 (L) and 270 (H) lncRNAs was also detected. The higher reproducibility of the RNA-Seq data was validated by RT-qPCR and porcine expression microarrays, therefore showing a strong correlation between RT-qPCR and RNA-Seq data (ranking from 0.79 to 0.96), as well as between microarrays and RNA-Seq (r=0.72). A differential expression analysis between H and L animals identified 55 genes differentially-expressed between groups. Pathways analysis revealed that these genes belong to biological functions, canonical pathways and three gene networks related to lipid and fatty acid metabolism. In concordance with the phenotypic classification, the pathways analysis inferred that linolenic and arachidonic acids metabolism was altered between extreme individuals. In addition, a connection was observed among the top three networks, hence suggesting that these genes are interconnected and play an important role in lipid and fatty acid metabolism.

CONCLUSIONS

In the present study RNA-Seq was used as a tool to explore the liver transcriptome of pigs with extreme phenotypes for intramuscular fatty acid composition. The differential gene expression analysis showed potential gene networks which affect lipid and fatty acid metabolism. These results may help in the design of selection strategies to improve the sensorial and nutritional quality of pork meat.

Identification of housekeeping genes suitable for gene expression analysis in Jian carp (Cyprinus carpio var. jian)

Jian carp (Cyprinus carpio var. jian) is an important economic fish species cultured in China. In this report, we performed a systematic analysis to identify an appropriate housekeeping (HK) gene for the study of gene expression in Jian carp. For this purpose, partial DNA sequences of four potential candidate genes (elongation factor 1 alpha (EF-1α), glyceraldehyde-3-phosphate (GAPDH), beta-actin (ACTB), and 18S ribosomal RNA (18S rRNA) were isolated, and their expression levels were studied using RNA extracted from nine tissues (forebrain, hypothalamus, liver, fore-intestine, hind-intestine, ovary, muscle, heart, kidney) in juvenile and adult Jian carp. Gene expression levels were quantified by quantitative real time RT-PCR (qRT-PCR), and expression stability was evaluated by comparing the coefficients of variation (CV) of the Ct values. The results showed that EF-1α was the most suitable HK gene in all tissues of juvenile and adult Jian carp. However, at distinct juvenile and adult developmental stages, there was not a single optimal gene for normalization of expression levels in all tissues. EF-1α was the most stable gene only in forebrain, hypothalamus, liver, heart, and kidney. These results provide data that can be expected to aid gene expression analysis in Jian carp research, but underline the importance of identifying the optimal HK gene for each new experimental paradigm.

Reference genes to quantify gene expression during oogenesis in a teleost fish

Understanding the molecular events involved in the acquisition of competence during oogenesis is a key step to determine the secret of 'high quality' eggs for aquaculture. Quantitative real time polymerase chain reaction (qPCR) is the technique of election to determine changes in transcript abundance in such studies, but choosing reference genes for normalization, in particular during oogenesis, remains a challenge. In the present study, transcription of 6 functionally distinct genes, β actin (ACTB), cathepsin D (CTSD), cathepsin Z (CTSZ), elongation factor 1 α (EEF1A), TATA binding protein (TBP) and tubulin A (TUBA1A) was assessed as normalizers of bone morphogenetic protein (BMP) and activin membrane-bound inhibitor (BAMBI) gene expression in mRNA from Mozambique tilapia oocytes during oogenesis. Reverse transcription was equally efficient and varies little in all samples. Most of the genes considered for reference were stable during early stages of oogenesis but variations were observed during vitellogenesis. A single gene and up to 3 genes were shown to be insufficient for reliable normalization throughout the whole oogenesis. The combination of the genes ACTB, CTSD, EEF1A and CTSZ as reference was found to minimize variation and has the most stable expression pattern between maturation stages.

Atlantic halibut experimentally infected with nodavirus shows increased levels of T-cell marker and IFNγ transcripts

The transcript levels of viral RNAs, selected T-cell marker and cytokine genes, toll like receptor (TLR) 7, and two interferon stimulated genes (ISG) were analysed in sexually immature adult Atlantic halibut (Hippoglossus hippoglossus L.) experimentally infected with nodavirus. The expression of the T-cell markers, TLR7 and the cytokine genes was further explored in in vitro stimulated anterior kidney leucocytes (AK leucocytes) isolated from the experiment fish and from additional untreated non-injected fish. The levels of viral RNA1 and RNA2 were increasing in brain and eye at around 4 and 8weeks post injection (wpi), respectively, and still increasing at the end of the experiment, especially in eye. Immuno-positive cells and signs of vacuolisation in both brain and eye were seen at 14wpi. Increased transcript levels of TCRβ, CD4-2, CD4, CD8α, and Lck in brain and eye of the experimentally infected halibut suggested an involvement of halibut T-cells in the immune response against nodavirus. Interestingly, a similar expression pattern of TCRβ, CD4 and Lck was seen in both brain and eye. However, compared to brain that showed elevated transcript levels of TCRβ, CD4 and Lck mainly at 10 and 14wpi, the increase appeared earlier between 3 and 4wpi in the eye. Yet, an increase in the transcript level of IFNγ was seen at 10 and 14wpi in both organs. Moreover, elevated levels of TLR7, IL-1β, IL-6, ISG15 and Mx were detected in vivo. The in vitro experiments, stimulating AK leucocytes with ConA-PMA, imiquimod or nodavirus, further supported an involvement of IL-6 and IFNγ in the immune response against nodavirus and the involvement of CD8β(+) cells. Results from the present study thus indicate an importance of T-cells, IFNγ and the analysed ISGs in the immune response against nodavirus in Atlantic halibut, and would be of great help in future vaccination trials giving the possibility to monitor the immune response rather than mortality during post-vaccination challenge experiments.

Genetic resistance to rhabdovirus infection in teleost fish is paralleled to the derived cell resistance status

Genetic factors of resistance and predisposition to viral diseases explain a significant part of the clinical variability observed within host populations. Predisposition to viral diseases has been associated to MHC haplotypes and T cell immunity, but a growing repertoire of innate/intrinsic factors are implicated in the genetic determinism of the host susceptibility to viruses. In a long-term study of the genetics of host resistance to fish rhabdoviruses, we produced a collection of double-haploid rainbow trout clones showing a wide range of susceptibility to Viral Hemorrhagic Septicemia Virus (VHSV) waterborne infection. The susceptibility of fibroblastic cell lines derived from these clonal fish was fully consistent with the susceptibility of the parental fish clones. The mechanisms determining the host resistance therefore did not associate with specific host immunity, but rather with innate or intrinsic factors. One cell line was resistant to rhabdovirus infection due to the combination of an early interferon IFN induction - that was not observed in the susceptible cells - and of yet unknown factors that hamper the first steps of the viral cycle. The implication of IFN was well consistent with the wide range of resistance of this genetic background to VSHV and IHNV, to the birnavirus IPNV and the orthomyxovirus ISAV. Another cell line was even more refractory to the VHSV infection through different antiviral mechanisms. This collection of clonal fish and isogenic cell lines provides an interesting model to analyze the relative contribution of antiviral pathways to the resistance to different viruses.