Reference gene selection for transcriptional profiling in Cryptocercus punctulatus, an evolutionary link between Isoptera and Blattodea

Identification and validation of the reference genes in the echiuran worm Urechis unicinctus based on transcriptome data

BACKGROUND

Real-time quantitative PCR (RT-qPCR) is a crucial and widely used method for gene expression analysis. Selecting suitable reference genes is extremely important for the accuracy of RT-qPCR results. Commonly used reference genes are not always stable in various organisms or under different environmental conditions. With the increasing application of high-throughput sequencing, transcriptome analysis has become an effective method for identifying novel stable reference genes.

RESULTS

In this study, we identified candidate reference genes based on transcriptome data covering embryos and larvae of early development, normal adult tissues, and the hindgut under sulfide stress using the coefficient of variation (CV) method in the echiuran Urechis unicinctus, resulting in 6834 (15.82%), 7110 (16.85%) and 13880 (35.87%) candidate reference genes, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that the candidate reference genes were significantly enriched in cellular metabolic process, protein metabolic process and ribosome in early development and normal adult tissues as well as in cellular localization and endocytosis in the hindgut under sulfide stress. Subsequently, ten genes including five new candidate reference genes and five commonly used reference genes, were validated by RT-qPCR. The expression stability of the ten genes was analyzed using four methods (geNorm, NormFinder, BestKeeper, and ∆Ct). The comprehensive results indicated that the new candidate reference genes were more stable than most commonly used reference genes. The commonly used ACTB was the most unstable gene. The candidate reference genes STX12, EHMT1, and LYAG were the most stable genes in early development, normal adult tissues, and hindgut under sulfide stress, respectively. The log₂(TPM) of the transcriptome data was significantly negatively correlated with the Ct values of RT-qPCR (Ct =  - 0.5405 log₂(TPM) + 34.51), which made it possible to estimate the Ct value before RT-qPCR using transcriptome data.

CONCLUSION

Our study is the first to select reference genes for RT-qPCR from transcriptome data in Echiura and provides important information for future gene expression studies in U. unicinctus.

NR4A1 inhibits the epithelial–mesenchymal transition of hepatic stellate cells: Involvement of TGF-β–Smad2/3/4–ZEB signaling

This study aimed to examine whether nuclear receptor 4a1 (NR4A1) is involved in inhibiting hepatic stellate cell (HSC) activation and liver fibrosis through the epithelial–mesenchymal transition (EMT). HSC-T6 cells were divided into the control group, the acetaldehyde (200 μM, an EMT activator) group, and the NR4A1 activation group (Cytosporone B; 1 μM). The expression levels of the epithelial marker E-cadherin, the mesenchymal markers fibronectin (FN), vimentin, smooth muscle alpha-actin (α-SMA), and fibroblast-specific protein 1 (FSP-1), and the components of the transforming growth factor (TGF)-β pathway were detected by real-time polymerase chain reaction and western blotting. Compared with the control group, E-cadherin in the acetaldehyde group was downregulated, whereas FN, FSP-1, vimentin, α-SMA, and COL1A1/COL1A2 were upregulated (P < 0.05). Compared with the acetaldehyde group, NR4A1 agonist upregulated E-cadherin and downregulated FN, FSP-1, vimentin, α-SMA, and COL1A1/COL1A2 (P < 0.05). After acetaldehyde stimulation, TGF-β, Smad2/3/4, and zinc finger E-box-binding homeobox (ZEB) were upregulated, while Smad7 mRNA levels were downregulated (all P < 0.05). Compared with acetaldehyde alone, NR4A1 agonist increased Smad7 mRNA levels and reduced TGF-β, Smad2/3/4, and ZEB mRNA levels (all P < 0.05). NR4A1 activation suppresses acetaldehyde-induced EMT, as shown by epithelial and mesenchymal marker expression. The inhibition of the TGF-β–Smad2/3/4–ZEB signaling during HSC activation might be involved.