Identification and Validation of the Reference Genes in the Echiuran Worm Urechis Unicinctus Based on Transcriptome Data

Overview

Journal BMC Genomics

Publisher Biomed Central

Specialty Genetics

Date 2023 May 10

PMID 37165306

Authors

Jiao Chen

Yunjian Wang

Zhi Yang

Danwen Liu

Yao Jin

Xixi Li

Yuhang Deng

Boya Wang

Zhifeng Zhang

Yubin Ma

Affiliations

Soon will be listed here.

Abstract

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.

Citing Articles

Identification and validation of stable reference genes for RT-qPCR analyses of Kobresia littledalei seedlings.

Sun H, Li C, Li S, Ma J, Li S, Li X BMC Plant Biol. 2024; 24(1):389.

PMID: 38730341 PMC: 11088182. DOI: 10.1186/s12870-024-04924-w.

References

Li Z, Li X, Zhang Q, Yuan L, Zhou X . Reference gene selection for transcriptional profiling in Cryptocercus punctulatus, an evolutionary link between Isoptera and Blattodea. Sci Rep. 2020; 10(1):22169. PMC: 7746730. DOI: 10.1038/s41598-020-79030-6. View

Kerchner K, Mou T, Sun Y, Rusnac D, Sprang S, Briknarova K . The structure of the cysteine-rich region from human histone-lysine N-methyltransferase EHMT2 (G9a). J Struct Biol X. 2021; 5:100050. PMC: 8261083. DOI: 10.1016/j.yjsbx.2021.100050. View

Tang B, Tan A, Lim L, Lee S, Low D, Hong W . Syntaxin 12, a member of the syntaxin family localized to the endosome. J Biol Chem. 1998; 273(12):6944-50. DOI: 10.1074/jbc.273.12.6944. View

Zhang Y, Zhang Z, Ren M, Liu X, Zhou X, Yang J . Selection of Reference Genes for RT-qPCR Analysis in the Hawthorn Spider Mite, Amphitetranychus viennensis (Acarina: Tetranychidae), Under Acaricide Treatments. J Econ Entomol. 2022; 115(2):662-670. DOI: 10.1093/jee/toac019. View

Liu X, Zhang L, Zhang Z, Ma X, Liu J . Transcriptional response to sulfide in the Echiuran Worm Urechis unicinctus by digital gene expression analysis. BMC Genomics. 2015; 16:829. PMC: 4618349. DOI: 10.1186/s12864-015-2094-z. View

Wei M, Qin Z, Kong D, Liu D, Zheng Q, Bai S . Echiuran Hox genes provide new insights into the correspondence between Hox subcluster organization and collinearity pattern. Proc Biol Sci. 2022; 289(1982):20220705. PMC: 9449475. DOI: 10.1098/rspb.2022.0705. View

Song J, Cho J, Park J, Hwang J . Identification and validation of stable reference genes for quantitative real time PCR in different minipig tissues at developmental stages. BMC Genomics. 2022; 23(1):585. PMC: 9374586. DOI: 10.1186/s12864-022-08830-z. View

Cagin U, Puzzo F, Gomez M, Moya-Nilges M, Sellier P, Abad C . Rescue of Advanced Pompe Disease in Mice with Hepatic Expression of Secretable Acid α-Glucosidase. Mol Ther. 2020; 28(9):2056-2072. PMC: 7474269. DOI: 10.1016/j.ymthe.2020.05.025. View

Mughal B, Leemans M, Spirhanzlova P, Demeneix B, Fini J . Reference gene identification and validation for quantitative real-time PCR studies in developing Xenopus laevis. Sci Rep. 2018; 8(1):496. PMC: 5764961. DOI: 10.1038/s41598-017-18684-1. View

10.

Zhang L, Liu X, Liu J, Zhang Z . Characteristics and function of sulfur dioxygenase in Echiuran worm Urechis unicinctus. PLoS One. 2013; 8(12):e81885. PMC: 3846777. DOI: 10.1371/journal.pone.0081885. View

11.

Andersen C, Jensen J, Orntoft T . Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004; 64(15):5245-50. DOI: 10.1158/0008-5472.CAN-04-0496. View

12.

Bustin S, Benes V, Garson J, Hellemans J, Huggett J, Kubista M . The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009; 55(4):611-22. DOI: 10.1373/clinchem.2008.112797. View

13.

Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A . Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002; 3(7):RESEARCH0034. PMC: 126239. DOI: 10.1186/gb-2002-3-7-research0034. View

14.

Liu X, Zhang Z, Ma X, Li X, Zhou D, Gao B . Sulfide exposure results in enhanced sqr transcription through upregulating the expression and activation of HSF1 in echiuran worm Urechis unicinctus. Aquat Toxicol. 2015; 170:229-239. DOI: 10.1016/j.aquatox.2015.11.021. View

15.

Guo H, Jiang L, Xia Q . Selection of reference genes for analysis of stress-responsive genes after challenge with viruses and temperature changes in the silkworm Bombyx mori. Mol Genet Genomics. 2015; 291(2):999-1004. DOI: 10.1007/s00438-015-1125-4. View

16.

Fujiwara A, Tazawa E, Hino A, Asami K, Yasumasu I . Respiration in Eggs of the Echiuroid, Urechis unicinctus, Before and After Fertilization: echiuroid eggs/fertilization/respiration/redox dyes/uncoupler of oxidative phosphorylation. Dev Growth Differ. 2023; 28(5):431-442. DOI: 10.1111/j.1440-169X.1986.00431.x. View

17.

Wang J, Zhang L, Lian S, Qin Z, Zhu X, Dai X . Evolutionary transcriptomics of metazoan biphasic life cycle supports a single intercalation origin of metazoan larvae. Nat Ecol Evol. 2020; 4(5):725-736. DOI: 10.1038/s41559-020-1138-1. View

18.

Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M . KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2022; 51(D1):D587-D592. PMC: 9825424. DOI: 10.1093/nar/gkac963. View

19.

Liu D, Qin Z, Wei M, Kong D, Zheng Q, Bai S . Genome-Wide Analyses of Heat Shock Protein Superfamily Provide New Insights on Adaptation to Sulfide-Rich Environments in (Annelida, Echiura). Int J Mol Sci. 2022; 23(5). PMC: 8910992. DOI: 10.3390/ijms23052715. View

20.

Ma X, Liu X, Zhou D, Bai Y, Gao B, Zhang Z . The NF-κB pathway participates in the response to sulfide stress in Urechis unicinctus. Fish Shellfish Immunol. 2016; 58:229-238. DOI: 10.1016/j.fsi.2016.09.020. View