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.