Evaluation of Reference Genes to Analyze Gene Expression in Silverside Odontesthes humensis Under Different Environmental Conditions

Validation of qPCR reference genes in the endangered annual killifish Austrolebias charrua considering different tissues, gender and environmental conditions

The annual killifish Austrolebias charrua is an endangered species, endemic to the southern region of South America, which inhabits temporary ponds that emerges in the rainy season. The main anthropogenic threat driving the extinction of A. charrua stems from extensive agriculture, primarily due to the widrespread use of glyphosate-based herbicides near their habitats. Annual killifishes have been used as models for ecotoxicological studies but, up to now, there are no studies about reference genes in any Austrolebias species. This represents an obstacle to the use of qPCR-based technologies, the standard method for gene expression quantification. The present study aimed to select and validate potential reference genes for qPCR normalization in the annual killifish Austrolebias charrua considering different tissues, gender and environmental conditions. The candidate reference genes 18 s, actb, gapdh, ef1a, shox, eif3g, and the control gene atp1a1 were evaluated in male and female individuals in three different tissues (brain, liver, and gills) under two experimental conditions (control and acute exposition to Roundup Transorb®). The collected tissues were submitted to RNA extraction, followed by cDNA synthesis, cloning, sequencing, and qPCR. Overall, 18 s was the most stable reference gene, and 18 s and ef1a were the most stable combination. Otherwise, considering all variables, gapdh and shox were the least stable candidate genes. Foremost, suitable reference genes were validated in A. charrua, facilitating accurate mRNA quantification in this species, which might be useful for developing molecular tools of ecotoxicological assessment based on gene expression analysis for environmental monitoring of annual killifish.

Assessment of oxidative stress biomarkers in the threatened annual killifish Austrolebias charrua exposed to Roundup

This study aimed to analyze the toxic effects of Roundup Transorb® on the endangered Neotropical annual killifish Austrolebias charrua through the assessment of molecular and biochemical biomarkers. The fish were collected in temporary ponds and exposed to environmentally realistic concentrations of the herbicide (5 mg.L-1 for 96 h). The production of ROS, lipid peroxidation, DNA damage, and membrane fluidity were evaluated in the blood cells by flow cytometry. The mRNA expression of the antioxidant-related genes sod2, cat, gstα, atp1a1, gclc, and ucp1 across the brain, liver, and gills was quantified. The acute exposure of annual killifish to Roundup significantly increased ROS production, lipid peroxidation, and DNA damage in their erythrocytes. Likewise, Roundup Transorb® decreased membrane fluidity in the blood cells of the exposed fish. Gene expression analysis revealed that Roundup exposure alters the relative expression of genes associated with oxidative stress and antioxidant defense. Our results give rise to new insights into adaptive mechanisms of A. charrua in response to Roundup. Since Brazilian annual killifishes strongly risk extinction, this study paves the way for developing novel biotechnologies applied to environmental monitoring and aquatic toxicology assessment.

MicroRNA qPCR normalization in Nile tilapia (Oreochromis niloticus): Effects of acute cold stress on potential reference targets

The Nile tilapia (Oreochromis niloticus) is one of the most important cultured fish worldwide, but tilapia culture is largely affected by low temperatures. Recent studies suggest that microRNAs (miRNAs) regulate cold tolerance traits in fish. In general, qPCR-based methods are the simplest and most accurate forms of miRNA quantification. However, qPCR data heavily depends on appropriate normalization. Therefore, the aim of the present study is to determine whether the expression of previously tested, stably expressed miRNAs are affected by acute cold stress in Nile tilapia. For this purpose, one small nuclear RNA (U6) and six candidate reference miRNAs (miR-23a, miR-25-3, Let-7a, miR-103, miR-99-5, and miR-455) were evaluated in four tissues (blood, brain, liver, and gills) under two experimental conditions (acute cold stress and control) in O. niloticus. The stability of the expression of each candidate reference miRNA was analyzed by four independent methods (the delta Ct method, geNorm, NormFinder, and BestKeeper). Further, consensual comprehensive ranking of stability was built with RefFinder. Overall, miR-103 was the most stable reference miRNA in this study, and miR-103 and Let-7a were the best combination of reference targets. Equally important, Let-7a, miR-23a, and miR-25-3 remained consistently stable across different tissues and experimental groups. Considering all variables, U6, miR-99-5, and miR-455 were the least stable candidates under acute cold stress. Most important, suitable reference miRNAs were validated in O. niloticus, facilitating further accurate miRNA quantification in this species.

Stability evaluation and validation of appropriate reference genes for real-time PCR expression analysis of immune genes in the rohu (Labeo rohita) skin following argulosis

Argulosis is one of the most unrestrained economically significant freshwater fish ectoparasitic diseases. Proper selection or normalization of the best reference gene governs the accuracy of results of gene expression studies using real-time PCR. Earlier studies in rohu carp (Labeo rohita) have used reference genes without proper validation. Here, seven candidate reference genes viz., acidic ribosomal protein (ARP0), glyceraldehyde 3-phosphate dehydrogenase, RNA polymerase II (RPo), elongation factor1α (EF1α), α- tubulin (AT), ribosomal protein L 10, and β-actin were evaluated using four algorithms (geNorm, BestKeeper, NormFinder and ∆Ct) followed by a comprehensive gene expression analysis using skin tissue of rohu at varied time points of experimental Argulus siamensis infection. ARP0 and EF1α were found to be the most stable whereas RPo and AT were considered as least stable genes based on basal expression level and variation in expression levels. Validation of candidate reference genes was undertaken by looking into the expression of six immune-related genes using the two most stable and two least stable genes as housekeeping genes in Argulus-infected rohu skin at different time points of infection. An increased expression of immune genes indicated the role of inflammation and the immune modulation process at the site of attachment of parasites in governing infection.

Sublethal ammonia induces alterations of emotions, cognition, and social behaviors in zebrafish (Danio rerio)

Ammonia pollutants were usually found in aquatic environments is due to urban sewage, industrial wastewater discharge, and agricultural runoff and concentrations as high as 180 mg/L (NH₄⁺) have been reported in rivers. High ammonia levels are known to impair multiple tissue and cell functions and cause fish death. Although ammonia is a potent neurotoxin, how sublethal concentrations of ammonia influence the central nervous system (CNS) and the complex behaviors of fish is still unclear. In the present study, we demonstrated that acute sublethal ammonia exposure can change social behavior of adult zebrafish. The exposure to 90 mg /L of (NH₄⁺) for 4 h induced a strong fear response and lower shoaling cohesion; exposure to 180 mg /L of (NH₄⁺) for 4 h reduced the aggressiveness, and social recognition, while the anxiety, social preference, learning, and short-term memory were not affected. Messenger RNA expressions of glutaminase and glutamate dehydrogenase in the brain were induced, suggesting that ammonia exposure altered glutamate neurotransmitters in the CNS. Our findings in zebrafish provided delicate information of ammonia neurotoxicity in complex higher-order social behaviors, which has not been revealed previously. In conclusion, sublethal and acute ammonia exposure can change specific behaviors of fish, which might lead to reductions in individual and population fitness levels.

Modulation of miR-429 during osmotic stress in the silverside Odontesthes humensis

Silverside fish inhabit marine coastal waters, coastal lagoons, and estuarine regions in southern South America. Although silversides are not fully adapted to freshwater, they can tolerate a wide range of salinity variations. MicroRNAs (miRNAs) are a class of ∼22 nucleotide noncoding RNAs, which are crucial regulators of gene expression at post-transcriptional level. Current data indicate that miRNAs biogenesis is altered by situations of environmental stress, thereby altering the expression of target mRNAs. Foremost, the silversides were acutely exposed to 30 g.L−1 of salt to reveal in which tissue miR-429 could be differentially expressed. Thus, fish were acclimated to freshwater (0 g.L−1) and to brackish water (10 g.L−1), and then exposed to opposite salinity treatment. Here, we reveal that miR-429, a gill-enriched miRNA, emerges as a prime osmoregulator in silversides. Taken together, our findings suggest that miR-429 is an endogenous regulator of osmotic stress, which may be developed as a biomarker to assist silverside aquaculture.

Acute exposition to Roundup Transorb® induces systemic oxidative stress and alterations in the expression of newly sequenced genes in silverside fish (Odontesthes humensis)

Roundup Transorb® (RDT) is a glyphosate-based herbicide commonly used in agricultural practices worldwide. This herbicide exerts negative effects on the aquatic ecosystem and affects bioenergetic and detoxification pathways, oxidative stress, and cell damage in marine organisms. These effects might also occur at the transcriptional level; however, the expression of genes associated with oxidative stress has not been studied well. Odontesthes humensis is a native Brazilian aquatic species naturally distributed in the habitats affected by pesticides, including Roundup Transorb® (RDT). This study evaluated the toxic effects of short-term exposure to RDT on O. humensis. Moreover, the genes related to oxidative stress were sequenced and characterized, and their expressions in the gills, hepatopancreas, kidneys, and brain of the fish were quantified by quantitative reverse transcription-polymerase chain reaction. The animals were exposed to two environmentally relevant concentrations of RDT (2.07 and 3.68 mg L^-1) for 24 h. Lipid peroxidation, reactive oxygen species (ROS), DNA damage, and apoptosis in erythrocytes were quantified by flow cytometry. The expression of the target genes was modulated in most tissues in the presence of the highest tested concentration of RDT. In erythrocytes, the levels of lipid peroxidation, ROS, and DNA damage were increased in the presence of both the concentrations of RDT, whereas cell apoptosis was increased in the group exposed to 3.68 mg L^-1 RDT. In conclusion, acute exposure to RDT caused oxidative stress in the fish, induced negative effects on cells, and modulated the expression of genes related to the enzymatic antioxidant system in O. humensis.

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.

Evaluation of qPCR reference genes in GH-overexpressing transgenic zebrafish (Danio rerio)

Reference genes (RGs) must have a stable expression in tissues in all experimental conditions to normalize real-time quantitative reverse transcription PCR (qRT-PCR) data. F0104 is a highly studied lineage of zebrafish developed to overexpress the growth hormone (GH). It is assumed that the transgenic process may influence the expression levels of commonly used RGs. The objective of the present study was to make a comprehensive analysis of stability of canditade RGs actb1, actb2, b2m, eif2s2, eef1a1, gapdh, rplp2, rpl7, rpl13α, tuba1, and rps18, in gh-transgenic and non-transgenic zebrafish. Liver, brain, intestine and muscle samples from both groups had qRT-PCR results analyzed by dCt, geNorm, NormFinder, BestKeeper, and RefFinder softwares. Consensus analyses among software concluded that rpl13α, rpl7, and eef1a1 are the most stable genes for zebrafish, considering the studied groups and tissues. Gapdh, rps18, and tuba1 suffered variations in stability among different tissues of both groups, and so, they were listed as the genes with lowest stability. Results from an average pairwise variations test indicated that the use of two RGs would generate reliable results for gene expression analysis in the studied tissues. We conclude that genes that are commonly used in mammals for qRT-PCR assays have low stability in both non-transgenic and gh-transgenic zebrafish reinforcing the importance of using species-specific RGs.

Genome-Wide Identification and Evaluation of New Reference Genes for Gene Expression Analysis Under Temperature and Salinity Stresses in Ciona savignyi

Rapid adaptation/accommodation to changing environments largely contributes to maximal survival of invaders during biological invasions, usually leading to success in crossing multiple barriers and finally in varied environments in recipient habitats. Gene expression is one of the most important and rapid ways during responses to environmental stresses. Selection of proper reference genes is the crucial prerequisite for gene expression analysis using the common approach, real-time quantitative PCR (RT-qPCR). Here we identified eight candidate novel reference genes from the RNA-Seq data in an invasive model ascidian Ciona savignyi under temperature and salinity stresses. Subsequently, the expression stability of these eight novel reference genes, as well as other six traditionally used reference genes, was evaluated using RT-qPCR and comprehensive tool RefFinder. Under the temperature stress, two traditional reference genes, ribosomal proteins S15 and L17 (RPS15, RPL17), and one novel gene Ras homolog A (RhoA), were recommended as the top three stable genes, which can be used to normalize target genes with a high and moderate expression level, respectively. Under the salinity stress, transmembrane 9 superfamily member (TMN), MOB kinase activator 1A-like gene (MOB) and ubiquitin-conjugating enzyme (UBQ2) were suggested as the top three stable genes. On the other hand, several commonly used reference genes such as α-tubulin (TubA), β-tubulin (TubB) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) showed unstable expressions, thus these genes should not be used as internal controls for gene expression analysis. We also tested the expression level of an important stress response gene, large proline-rich protein bag6-like gene (BAG) using different reference genes. As expected, we observed different results and conclusions when using different normalization methods, thus suggesting the importance of selection of proper reference genes and associated normalization methods. Our results provide a valuable reference gene resource for the normalization of gene expression in the study of environmental adaptation/accommodation during biological invasions using C. savignyi as a model.

Identification of suitable reference genes for real-time qPCR in homocysteine-treated human umbilical vein endothelial cells

The imbalance in homocysteine (Hcy) metabolism has been implicated in the pathogenesis of human diseases, including cardiovascular and neurodegenerative disorders. When attempting to identify gene expression profiles using quantitative real-time reverse transcription polymerase chain reaction (RT-qPCR), the selection of suitable reference genes is important. Here, the expression levels of 10 commonly used reference genes were assessed for normalization of RT-qPCR in Hcy-treated human umbilical vein endothelial cells (HUVECs) and control cells. The suitability of eight selected candidate genes was comparatively analyzed across the tested samples and separately ranked by four programs, geNorm, NormFinder, BestKeeper, and the ΔCt method. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was the most stable gene in the final ranking using the RankAggreg package. Surprisingly, the β-actin (ACTB) levels decreased significantly in Hcy-treated HUVECs compared with control HUVECs (P<0.05), and further study indicated that Hcy suppressed the expression of ACTB by upregulating the miR-145-5p level in Hcy-treated HUVECs. Our data suggest that GAPDH can be used as a reliable reference gene, while ACTB cannot; normalization of gene expression in RT-qPCR experiments in Hcy-treated HUVECs. The data, which identifies a suitable reference gene in Hcy-treated HUVECs, will contribute to the design of an effective and accurate method for quantitation of gene expression.