Determination of reference microRNAs for relative quantification in grass carp (Ctenopharyngodon idella)

Selection of references for quantitative real-time PCR analysis of microRNAs in Nile tilapia (Oreochromis niloticus) under osmotic stress

MicroRNAs play crucial regulatory roles in various aspects of development and physiology, including environmental adaptation and stress responses in teleosts. RT-qPCR is the most commonly used method for studying microRNA expression, with the accuracy and reliability of results depending on the use of an appropriate reference gene for normalization. This study aimed to evaluate seven miRNAs (U6, Let-7a, miR-23a, miR-25-3, miR-103, miR-99-5, and miR-455) expression stability in different tissues of Nile tilapia subjected to osmotic stress. Fish were divided into two groups: a control and an experimental group, raised in 0 and 12 ppt salinity water respectively. After 21 days, brain, gills, liver, and posterior intestine were collected for analysis. Different mathematical algorithms (geNorm, NormFinder, BestKeeper, and the comparative ΔCt method) were employed to identify the most suitable reference miRNAs. The results indicate that the miR-455/miR-23a combination is a robust reference for normalizing miRNA expression levels in studies of osmotic stress responses in Nile tilapia. The stability of miRNA expression can vary depending on specific stress conditions and biological processes, underscoring the necessity of selecting appropriate normalizing miRNAs for each experimental context. This study identifies reliable reference genes for future RT-qPCR analyses of miRNA expression, thereby enhancing our understanding of molecular responses in fish to environmental challenges. These insights are fundamental to the development of new technologies for the improved management and sustainability of aquaculture practices.

LncRNA-adm2 targets adm2 via cid-miR-n3 and negatively regulates the inflammatory response in grass carp (Ctenopharyngodon idella)

Long non-coding RNAs (lncRNAs), which impact gene expression following pathogen infections, have garnered significant attention in recent years. Recent discoveries have revealed that lncRNAs play a crucial role in fish immune responses to pathogen infections. We investigated the influence of lncRNA-adm2 on the antibacterial immune response generated by Aeromonas hydrophila in grass carp (Ctenopharyngodon idella) through the adsorption of cid-miR-n3. Furthermore, we found that cid-miR-n3 interacts with lncRNA-adm2 and targets the 3' UTR of adm2. The upregulation of lncRNA-adm2 expression led to the suppression of pro-inflammatory cytokines (il-1β and il-6) in CIK cells, while anti-inflammatory cytokines (il-10) increased. Our research provides evidence that lncRNAs are involved in the antibacterial immune response of fish, expanding our understanding of the function of lncRNAs in teleosts.

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.

miR-130a targets CiGadd45bb to modulate the inflammatory response to bacterial infection in Ctenopharyngodon idella kidney (CIK) cells

Septicemia is a systemic inflammatory response to bacterial infection that results in a hyper-inflammatory state, which could lead to septic shock and death in grass carp (Ctenopharyngodon idella). The aim of this study was to determine the underlying mechanism of microRNA (miR-130a) in bacteria-infected grass carp. Expression levels of miR-130a against Aeromonas hydrophila (A. hydrophila) infection in Ctenopharyngodon idella kidney cells (CIK) were analyzed. Luciferase reporter assay, quantitative reverse transcription-polymerase chain reaction were performed to explore whether Ctenopharyngodon idella growth arrest and DNA damage-inducible 45 (CiGadd45bb) was a target of miR-130a. MiR-130a mimic, inhibitor and miR-control were transfected to CIK respectively. After transfection, the expression levels of proinflammatory genes were determined. Here we show that CiGadd45bb as a target of miR-130a. We also confirmed that miR-130a levels were significantly higher after being stimulated for 4 h and lower after 12 h (P < 0.01). Overexpressing miR-130a strikingly inhibited p38, JNK, ERK and TNF-a genes (P < 0.01) and silencing miR-130a activated p38, JNK, ERK, TNF-a, IFN and IL-8 (P < 0.01). Our results provide a theoretical basis for studying the molecular mechanism underlying the regulation of inflammation by miR-130a in grass carp.

Multiomics analysis revealed miRNAs as potential regulators of the immune response in Carassius auratus gills to Aeromonas hydrophila infection

The gill of fish is an important immune organ for pathogen defense, but its microRNA (miRNA) expression and regulatory mechanism remain unclear. In this study, we report on the histopathological and immunohistochemical features of the gills of the crucian carp Carassius auratus challenged with Aeromonas hydrophila. Small RNA libraries of the gills were constructed and sequenced on the Illumina HiSeq 2000 platform. A total of 1,165 differentially expressed miRNAs (DEMs) were identified in gills, of which 539 known and 7 unknown DEMs were significantly screened (p < 0.05). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that the potential target genes/proteins were primarily involved in 33 immune-related pathways, in which the inflammatory responses were focused on the Toll-like receptor (TLR), mitogen-activated protein kinase (MAPK), and nuclear factor kappa B (NF-κB) signaling pathways. Moreover, the expression levels of 14 key miRNAs (e.g., miR-10, miR-17, miR-26a, miR-144, miR-145, and miR-146a) and their target genes (e.g., TNFα, TLR4, NF-κB, TAB1, PI3K, and IRAK1) were verified. In addition, the protein levels based on isobaric tags for relative and absolute quantification (iTRAQ) were significantly associated with the results of the quantitative real-time PCR (qRT-PCR) analysis (p < 0.01). miR-17/pre-miR-17 were identified in the regulation expression of the NF-κB target gene, and the phylogenetic tree analysis showed that the pre-miR-17 of C. auratus with the closest similarity to the zebrafish Danio rerio is highly conserved in teleosts. This is the first report of the multi-omics analysis of the miRNAs and proteins in the gills of C. auratus infected with A. hydrophila, thus enriching knowledge on the regulation mechanism of the local immune response in Cyprinidae fish.

elk1/miR-462-731 Feedback Loop Regulates Macrophages Polarization and Phagocytosis in Grass Carp (Ctenopharyngodon idella)

MicroRNA clusters are microRNAs (miRNAs) that are distributed in close proximity on chromosomes. In this study, we report a miRNA cluster identified from grass carp (Ctenopharyngodon idella), miR-462-731, which plays a positive role in host antibacterial immunity. The expression of miR-462-731 was disrupted after infection by Aeromonas hydrophila. Transcription factor ETS transcription factor ELK1 was identified to bind to the promoter of the miR-462-731 cluster and suppress its expression. In addition, miR-731 negatively regulates the expression of elk1, forms an elk1/miR-462-731 double negative feedback loop. In addition, we found that miR-731 directly targets ezrin a (ezra), participates in inducing PI3K/AKT signaling in macrophage, to induce macrophage polarization to the M1 phenotype with stronger phagocytosis. Our results demonstrate a novel elk1/miR-462-731 feedback loop. The data deepen our understanding of the relationship between macrophage polarization and phagocytosis in teleost fish.

Selection of References for microRNA Quantification in Japanese Flounder (Paralichthys olivaceus) Normal Tissues and Edwardsiella tarda-Infected Livers

MicroRNA (miRNA) plays essential roles in post-transcriptional regulation of protein coding genes, and the quantitative real-time polymerase chain reaction (qRT-PCR) is the powerful and broadly employed tool to conduct studies of miRNA expression. Identifying appropriate references to normalize quantitative data is a prerequisite to ensure the qRT-PCR accuracy. Until now, there has been no report about miRNA reference for qRT-PCR in Japanese flounder (Paralichthys olivaceus), one important marine cultured fish along the coast of Northern Asia. In this study, combined with miRNA-Seq analysis and literature search, 10 candidates (miR-34a-5p, miR-205-5p, miR-101a-3p, miR-22-3p, miR-23a-3p, miR-210-5p, miR-30c-5p, U6, 5S rRNA, and 18S rRNA) were chosen as potential references to test their expression stability among P. olivaceus tissues, and in livers of P. olivaceus infected with Edwardsiella tarda at different time points. The expression stability of these candidates was analyzed by qRT-PCR and evaluated with Delta CT, BestKeeper, geNorm, as well as NormFinder methods, and RefFinder was employed to estimate the comprehensive ranking according to the four methods. As the result, miR-22-3p and miR-23a-3p were proved to be the suitable combination as reference miRNAs for both P. olivaceus normal tissues and livers infected with E. tarda, and they were successfully applied to normalize miR-7a and miR-221-5p expression in P. olivaceus livers in response to E. tarda infection. All these results provide valuable information for P. olivaceus miRNA quantitative expression analysis in the future.

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.

Functional miR-142a-3p Induces Apoptosis and Macrophage Polarization by Targeting tnfaip2 and glut3 in Grass Carp (Ctenopharyngodon idella)

In the process of microbial invasion, the inflammation reaction is induced to eliminate the pathogen. However, un-controlled or un-resolved inflammation can lead to tissue damage and death of the host. MicroRNAs (miRNAs) are the signaling regulators that prevent the uncontrolled progress of an inflammatory response. Our previous work strongly indicated that miR-142a-3p is related to the immune regulation in grass carp. In the present study, we found that the expression of miR-142a-3p was down-regulated after infection by Aeromonas hydrophila. tnfaip2 and glut3 were confirmed as be the target genes of miR-142a-3p, which were confirmed by expression correlation analysis, gene overexpression, and dual luciferase reporter assay. The miR-142a-3p can reduce cell viability and stimulate cell apoptosis by targeting tnfaip2 and glut3. In addition, miR-142a-3p also regulates macrophage polarization induced by A. hydrophila. Our results suggest that miR-142a-3p has multiple functions in host antibacterial immune response. Our research provides further understanding of the molecular mechanisms between miRNAs and their target genes, and provides a new insights for the development of pro-resolution strategies for the treatment of complex inflammatory diseases in fish.

LncRNA-WAS and lncRNA-C8807 interact with miR-142a-3p to regulate the inflammatory response in grass carp

Septicemia of grass carp is a systemic inflammatory reaction caused by bacterial infection. More and more evidences show that long non-coding RNAs (lncRNAs) can participate in the regulation of inflammatory response. In the present study, lncRNA-WAS and lncRNA-C8807 were confirmed to be involved in the inflammatory response following infection with Aeromonas hydrophila. LncRNA-WAS and lncRNA-C8807 could interact with miR-142a-3p. LncRNA-WAS and lncRNA-C8807 interact with miR-142a-3p to effect pro-inflammatory genes and NF-κB pathway. Our results provide a theoretical basis for studying the molecular mechanism underlying the regulation of inflammation by lncRNA in grass carp.

miR-148 targets CiGadd45ba and CiGadd45bb to modulate the inflammatory response to bacterial infection in grass carp

In grass carp (Ctenopharyngodon idella), septicemia is a systemic inflammatory response to bacterial infection. Once infected bacteria, a hyper-inflammatory state that could lead to septic shock and death. There is increasing evidence that microRNAs are involved in the regulation of the inflammatory response. Ctenopharyngodon idella growth arrest and DNA damage-inducible 45 ba and bb (CiGadd45ba and CiGadd45bb) are two subtypes of Gadd45b. In the present study, miR-148 was confirmed to be involved in the inflammatory response after infection with Aeromonas hydrophila. Dual-luciferase reporter assays and miRNA expression profiling confirmed that miR-148 targeted both CiGadd45ba and CiGadd45bb. Transfection with miR148 mimics and inhibitors altered the expression levels of proinflammatory genes, suggesting that miRNAs regulate the immune response in grass carp. Our results provide a theoretical basis for studying the molecular mechanism underlying the regulation of inflammation by miR-148 in grass carp.

miR-23a-3p and miR-23a-5p target CiGadd45ab to modulate inflammatory response and apoptosis in grass carp

Ctenopharyngodon idella growth arrest and DNA damage-inducible 45 ab (CiGadd45ab) is a subtype of the Gadd45a gene of the Gadd45 family in grass carp. There is increasing evidence that microRNAs (miRNAs) are involved in the regulation of inflammatory and apoptotic responses. However, little is known about the regulatory effects of miRNAs on CiGadd45ab expression. In the present study, CiGadd45ab was identified as a target gene of miR-23a-3p and miR-23a-5p, based on miRNA expression profiling and a dual-luciferase reporter assay. In addition, miR-23a-3p and miR-23a-5p were both confirmed to be involved in the inflammatory response following infection with Aeromonas hydrophila by targeting CiGadd45ab. Transfection with miR-23a-3p and miR-23a-5p mimics and inhibitor altered proinflammatory gene expression and apoptosis rate, thereby suggesting that miRNAs regulate immune response and anti-apoptosis by targeting CiGadd45ab in grass carp. Our results provide a theoretical basis for exploring the molecular mechanisms by which miR-23a-3p and miR-23a-5p target CiGadd45ab to regulate inflammation and apoptosis against bacterial infection in grass carp.

miR-21 targets jnk and ccr7 to modulate the inflammatory response of grass carp following bacterial infection

Grass carp septicemia is a systemic inflammatory response that develops following a bacterial infection. The hyperinflammatory state develops could lead to septic shock and lethality. There is increasing evidence that microRNAs are involved in the regulation of the inflammatory response. In the present study, miR-21 was confirmed to be involved in the inflammatory response following infection with Aeromonas hydrophila and LPS stimulation. Both jnk and ccr7 were identified as target gene of miR-21 by overexpression, inhibition, and dual luciferase reporter assays experiments. Meanwhile, miR-21 targets the jnk and ccr7 to modulate downstream pro-inflammatory factors tnf-α, il-1β, il-6, and il-12. Our results provide a theoretical basis for exploring the molecular mechanism of grass carp miR-21 regulating inflammation.

Cloning, functional analysis, and microRNA-induced negative regulation of growth arrest and DNA damage-inducible 45 γ (Gadd45g) in grass carp (Ctenopharyngodon idella)

Growth arrest and DNA damage-inducible 45 gamma (Gadd45g) is a member of Gadd45 gene family of immunological proteins in mammals. Herein, we identified and characterised Gadd45g from grass carp. The cDNA spans over 1189 bp, with an open reading frame of 480 bp encoding a 159 amino acid protein. CiGadd45g mRNAs were expressed in all tissues investigated, with abundant expression in liver, kidney, heart, brain, blood and skin. Following infection with Aeromonas hydrophila, CiGadd45g expression was upregulated in these immune-related tissues (gill, liver, spleen, intestine, kidney and head kidney). Immune-related cytokines (p38 and JNK) and proinflammatory cytokines (IL-8, IFN-1 and TNF-α) were activated by CiGadd45g. CiGadd45g and downstream genes were regulated by microRNA miR-429b. These results indicate that CiGadd45g plays an important immune role in the response to A. hydrophila infection in grass carp.

Selection of reference genes for miRNA quantitative PCR and its application in miR-34a/Sirtuin-1 mediated energy metabolism in Megalobrama amblycephala

MiRNAs are small, non-coding RNAs that downregulate gene expression at post-transcriptional levels. They have emerged as important regulators involved in metabolism, immunity, and cancer. Real-time quantitative PCR (RT-qPCR) is an effective and main method for quantifying target miRNA. For robust RT-qPCR method, suitable reference genes play crucial roles in data normalization. Blunt snout bream (Megalobrama amblycephala) is an economically important aquaculture species; however, no reference genes dedicated for qPCR method has been identified for this species so far. The objective of this study was to screen stable reference genes for miRNA RT-qPCR and demonstrated its application in energy metabolism in blunt snout bream. The stabilities of ten potential reference genes (miR-21-1-5p, miR-107a-3p, miR-222a-3p, miR-146a-5p, miR-101a-3p, miR-22a-3p, miR-103-3p, miR-456-3p, miR-221-3p, and U6 (RNU6A)) were evaluated in nine tissues (brain, muscle, liver, skin, spleen, heart, gill, intestine, and eye) under normal condition and in three tissues (liver, intestine, and spleen) under four stresses (heat stress, ammonia stress, bacterial challenge, and glycolipid stress). Using GeNorm, NormFinder, and RefFinder softwares, we discovered that different tissues and stresses are both important variability factors for the expression stability of miRNAs. After verifying miR-34a/Sirtuin-1 expressions in high-carbohydrate diet-induced blunt snout bream, we eventually identified that the most stable reference gene in this species was miR-221-3p, and the best combination of reference genes were miR-221-3p and miR-103-3p.

Selection of suitable candidate genes for miRNA expression normalization in Yellow River Carp (Cyprinus carpio. var)

Yellow River carp is widely cultivated in the world due to its economic value in aquaculture, and the faster growth of females compared to males. It is believed that microRNAs (miRNA) are involved in gonadal differentiation and development. qPCR is the most preferred method for miRNA functional analysis. Reliable reference genes for normalization in qRT-PCR are the key to ensuring the accuracy of this method. The aim of present research was to evaluate as well as identify the efficacy of reference genes for miRNA expression using qRT-PCR in Yellow River carp. Nine ncRNAs (miR-101, miR-23a, let7a, miR-26a, miR-146a, miR-451, U6, 5S, and 18S) were chosen and tested in four sample sets: (1) different tissues in adult carp, (2) different tissues in juvenile carp, (3) different early developmental stages of carp, and (4) different developmental stages of carp gonads. The stability and suitability values were calculated using NormFinder, geNorm, and BestKeeper software. The results showed that 5S was a suitable reference gene in different tissues of adult and juvenile carp. The genes 5S, 18S, and U6 were the most stable reference genes in the early developmental stages of carp. Let-7a and miR-23a were considered as the suitable reference genes in the development of gonads. All these reference genes were subsequently validated using miR-430. The results showed that genes 5S and 18S were the most suitable reference genes to normalize miRNA expression under normal growth conditions in early different developmental stages. The genes Let-7a, and miR-23a were the most suitable in different developmental stages. The present study is the first comprehensive study of the stability of miRNA reference genes in Yellow River carp, providing valuable as well as basic data for investigating more accurate miRNA expression during gonadal differentiation and development of carp.

Identification and characterization of skin color microRNAs in Koi carp (Cyprinus carpio L.) by Illumina sequencing

Background

MicroRNAs (miRNAs) are endogenous, small (21–25 nucleotide), non-coding RNAs that play important roles in numerous biological processes. Koi carp exhibit diverse color patterns, making it an ideal subject for studying the genetics of pigmentation. However, the influence of miRNAs on skin color regulation and variation in Koi carp is poorly understood.

Results

Herein, we performed small RNA (sRNA) analysis of the three main skin colors in Koi carp by Illumina sequencing. The results revealed 330, 397, and 335 conserved miRNAs (belonging to 81 families) and 340, 353, and 351 candidate miRNAs in black, red, and white libraries, respectively. A total of 164 differentially expressed miRNAs (DEMs) and 14 overlapping DEMs were identified, including miR-196a, miR-125b, miR-202, miR-205-5p, miR-200b, and etc. Target prediction and functional analysis of color-related miRNAs such as miR-200b, miR-206, and miR-196a highlighted putative target genes, including Mitf, Mc1r, Foxd3, and Sox10 that are potentially related to pigmentation. Determination of reference miRNAs for relative quantification showed that let-7a was the most abundant single reference gene, and let-7a and miR-26b was the most abundant combination.

Conclusions

The findings provide novel insight into the molecular mechanisms determining skin color differentiation in Koi carp, and serve as a valuable reference for future studies on tissue-specific miRNA abundance in Koi carp.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5189-5) contains supplementary material, which is available to authorized users.

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

Some mammalian reference genes, which are widely used to normalize the qRT-PCR, could not be used for this purpose due to its high expression variation. The normalization with false reference genes leads to misinterpretation of results. The silversides (Odontesthes spp.) has been used as models for evolutionary, osmoregulatory and environmental pollution studies but, up to now, there are no studies about reference genes in any Odontesthes species. Furthermore, many studies on silversides have used reference genes without previous validations. Thus, present study aimed to was to clone and sequence potential reference genes, thereby identifying the best ones in Odontesthes humensis considering different tissues, ages and conditions. For this purpose, animals belonging to three ages (adults, juveniles, and immature) were exposed to control, Roundup®, and seawater treatments for 24 h. Blood samples were subjected to flow-cytometry and other collected tissues to RNA extraction; cDNA synthesis; molecular cloning; DNA sequencing; and qRT-PCR. The candidate genes tested included 18s, actb, ef1a, eif3g, gapdh, h3a, atp1a, and tuba. Gene expression results were analyzed using five algorithms that ranked the candidate genes. The flow-cytometry data showed that the environmental challenges could trigger a systemic response in the treated fish. Even during this systemic physiological disorder, the consensus analysis of gene expression revealed h3a to be the most stable gene expression when only the treatments were considered. On the other hand, tuba was the least stable gene in the control and gapdh was the least stable in both Roundup® and seawater groups. In conclusion, the consensus analyses of different tissues, ages, and treatments groups revealed that h3a is the most stable gene whereas gapdh and tuba are the least stable genes, even being considered two constitutive genes.

Selection of reference genes for expression analysis in mouse models of acute alcoholic liver injury

Investigations of hepatic gene expression are crucial for determining the molecular factors involved in acute alcoholic liver injury. The results of liver molecular investigations may reveal etiologically important genomic alterations. Therefore, it is necessary to normalize gene expression data to identify stable genes, which may be used as a reference under different experimental conditions. The aim of the present study was to apply reverse transcription‑quantitative polymerase chain reaction analysis and use analysis software to investigate the expression stability of candidate reference genes in hepatic tissues from mice with acute alcoholic liver injury. The acute alcoholic liver injury models were established by the intragastric administration of alcohol (5 mg/kg) in Imprinting Control Region mice. Total RNA was isolated from the mouse livers, following which the expression levels of seven reference genes, β-actin, glyceraldehyde 3-phosphate dehydrogenase (Gadph), glucuronidase β, hypoxanthine phosphoribosyltransferase 1 (Hprt1), 18S ribosomal RNA, TATA binding protein and β‑2 microglobulin, were examined, and gene expression stability was assessed using the geNorm, NormFinder and BestKeeper tools. The geNorm analysis revealed that the gene with the lowest variability was Hprt1. Hprt1 and Gapdh were validated as the optimal reference gene pair in all samples from all groups. The NormFinder and BestKeeper results showed that Hprt1 was the most stable gene in all samples. Alcohol induces endoplasmic reticulum (ER) stress, causing changes in the expression levels of ER stress‑associated genes. The stability of Hprt1 was verified by the expression analysis of ER stress‑associated genes, and gene expression levels in the ethanol groups were upregulated, with a significant difference in expression, compared with those in the control group. Therefore, Hprt1 was selected as the most stable gene, and Hprt1 and Gapdh were determined to be the optimum gene pair in mouse models of acute alcoholic liver injury. The reliability of the Hprt1 gene was confirmed by expression analysis of ER stress‑associated genes.

Discovery and characterization of conserved and novel microRNAs from blunt snout bream (Megalobrama amblycephala) by deep sequencing

MicroRNAs (miRNAs) are short, single stranded RNA molecules with approximately 22 nts in length, which regulate the stability and translation of messenger RNAs in several organisms. To increase the repertoire of miRNAs characterized in M. amblycephala, we used the deep sequencing technology to sequence a small RNA library using pooled RNA sample isolated from the 4 different tissues of M. amblycephala. A total of 309 conserved miRNAs that originated from 131 miRNA families were detected. 15 novel candidates miRNA were identified. Randomly selected 6 miRNAs were analyzed by stem-loop qRT-PCR and differential expression patterns were observed in 6 different tissues of M. amblycephala. Furthermore, the potential targets were predicted. GO analysis showed that most of the targets were involved in a broad range of physiological functions including fish growth, development, metabolism, stress responses and so on. Overall, our results significantly increased the number of novel miRNAs in M. amblycephala, which should be useful for further investigation into the role of miRNAs in regulating diverse biological processes.

Integrated mRNA and microRNA transcriptome analyses reveal regulation of thermal acclimation in Gymnocypris przewalskii: A case study in Tibetan Schizothoracine fish

Environmental acclimation is important episode in wildlife occupation of the high-altitude Tibetan Plateau (TP). Transcriptome-wide studies on thermal acclimation mechanism in fish species are rarely revealed in Tibetan Plateau fish at high altitude. Thus, we used mRNA and miRNA transcriptome sequencing to investigate regulation of thermal acclimation in larval Tibetan naked carp, Gymnocypris przewalskii. We first remodeled the regulation network of mRNA and miRNA in thermal acclimation, and then identified differential expression of miRNAs and target mRNAs enriched in metabolic and digestive pathways. Interestingly, we identified two candidate genes contributed to normal skeletal development. The altered expression of these gene groups could potentially be associated with the developmental issues of deformity and induced larval death. Our results have three important implications: first, these findings provide strong evidences to support our hypothesis that G. przewalskii possess ability to build heat-tolerance against the controversial issue. Second, this study shows that transcriptional and post-transcriptional regulations are extensively involved in thermal acclimation. Third, the integrated mRNA and microRNA transcriptome analyses provide a large number of valuable genetic resources for future studies on environmental stress response in G. przewalskii and as a case study in Tibetan Schizothoracine fish.

MicroRNA-induced negative regulation of TLR-5 in grass carp, Ctenopharyngodon idella

MicroRNAs (miRNAs) are endogenous small non-coding RNAs that play crucial roles in numerous biological processes. However, the role of miRNAs in antibacterial defence in fish has not been fully determined. Here, we identified that nine miRNAs are differentially expressed in kidney between susceptible and resistant grass carp strains. Analysis of spatial and temporal miRNA expression patterns suggests that cid-miRn-115 and miR-142a-3p are potential regulators of anti-bacterial activity. Overexpressing of cid-miRn-115 and miR-142a-3p results in a visible change in Ctenopharyngodon idella kidney (CIK) cells immune effector activity. Bioinformatics analysis and overexpressing assay shows that cid-miRn-115 and miR-142a-3p directly regulate tlr5 expression. cid-miRn-115 and miR-142a-3p overexpressing leads to a significant decrease in tlr5 expression in CIK, thereby repressing its downstream genes, such as il-1β, il-8 and tnf-α. These findings provide a novel insight into the determination of anti-bacterial compounds in grass carp.

Next-generation sequencing identified microRNAs that associate with motile aeromonad septicemia in grass carp

MicroRNAs (miRNAs) modulate many biological processes and can be up-or down-regulated in a variety of diseases states, including in bacterial infection. In this study, we characterized miRNAs associated with susceptibility or resistance to motile Aeromonas hydrophila in grass carp (Ctenopharyngodon idella), a commercially important farmed fish species in China. Using Illumina next-generation sequencing, we detected 185 miRNAs in A. hydrophila-susceptible (SGC) and -resistant grass carp (RGC) during immune activation. Twenty-one miRNAs were found to be differentially expressed between SGC and RGC, with 23.8% (5 out 21) exhibiting elevated expression in SGC. Further investigation of a selected five miRNAs indicated differences in the timing and tissue sources of expression. Target genes likely regulated by these miRNAs were identified using computational prediction combined with transcriptome expression data. MiRNA targets of differentially expressed mRNAs included those with functions important for bacterial resistance, such as tlr4 and nfil3-6. In addition, we demonstrated that nfil3-6 and tlr4 genes are direct targets of let-7i and cid-miRn-118. The present study suggested that the multiple miRNAs, displaying diverse kinetics and tissue sources on antibacterial immune processes.

Selection of reference genes for microRNA quantitative expression analysis in Chinese perch, Siniperca chuatsi

Real-time quantitative reverse transcription PCR (RT-qPCR) is one of the most effective and sensitive techniques in gene expression assay, for which selection of reference genes is a prerequisite. In teleost species, such as Chinese perch, the expression profiling of miRNAs as reference genes for RT-qPCR has not been intensively studied. In the present study, the expression profiles of six miRNAs (miR-101a, miR-146a, miR-22a, miR-23a, miR-26a and let-7a) and one small nuclear RNA (U6) were assayed with RT-qPCR in different adult tissues, developmental stages and growth conditions of Chinese perch, Siniperca chuatsi. The analyses revealed that embryonic developmental stage is an important variability factor in the expression stability of miRNAs. All six miRNAs exhibited better expression consistency than U6 in most of the conditions examined, and therefore, they may be more suitable as a reference gene for miRNA quantification. When different tissues and developmental stages were considered, miR-22a demonstrated the most consistent expression pattern, and the best combination of reference genes was miR-22a and miR-23a. Our study offers useful data for selecting miRNAs as reference genes for RT-qPCR analysis of miRNAs in teleost fishes under different conditions.

De novo assembly of the grass carp Ctenopharyngodon idella transcriptome to identify miRNA targets associated with motile aeromonad septicemia

Background

De novo transcriptome sequencing is a robust method of predicting miRNA target genes, especially for organisms without reference genomes. Differentially expressed miRNAs had been identified previously in kidney samples collected from susceptible and resistant grass carp (Ctenopharyngodon idella) affected by Aeromonas hydrophila. Target identification for these differentially expressed miRNAs poses a major challenge in this non-model organism.

Results

Two cDNA libraries constructed from mRNAs of susceptible and resistant C. idella were sequenced by Illumina Hiseq 2000 technology. A total of more than 100 million reads were generated and de novo assembled into 199,593 transcripts which were further extensively annotated by comparing their sequences to different protein databases. Biochemical pathways were predicted from these transcript sequences. A BLASTx analysis against a non-redundant protein database revealed that 61,373 unigenes coded for 28,311 annotated proteins. Two cDNA libraries from susceptible and resistant samples showed that 721 unigenes were expressed at significantly different levels; 475 were significantly up-regulated and 246 were significantly down-regulated in the SG samples compared to the RG samples. The computational prediction of miRNA targets from these differentially expressed genes identified 188 unigenes as the targets of 5 conserved and 4 putative novel miRNA families.

Conclusion

This study demonstrates the feasibility of identifying miRNA targets by transcriptome analysis. The transcriptome assembly data represent a substantial increase in the genomic resources available for C. idella and will provide insights into the gene expression profile analysis and the miRNA function annotations in further studies.