Transcriptomic Study of Late Defocus Recognition in Chicks Based on RNA-seq Technology
Liyuan Sun, Li Zhu, Sitong Chen, Kai Wang, Mingwei Zhao
Institute of Medical Technology, Peking University Health Science Center, Department of Ophthalmology, Peking University People's Hospital, Eye Diseases and Optometry Institute, Beijing Key Laboratory of Diagnosis and Therapy of Retinal and Choroid Diseases, College of Optometry, Peking University Health Science Center, Beijing 100044, China
Abstract:Objective: To investigate the molecular mechanism of late defocus recognition. Methods: In this experimental study, 30 7-day-old White Leghorn chicks were randomly divided into a minus-lens group, a plus lens group and a control group, and the chicks were induced with -10 D/+10 D/0 D lens in the right eye of the chicks respectively, with 6 days' intervention time. Biological parameters such as diopter, axial length, and thickness of each tissue in the posterior pole were measured in all animals before and after the optical intervention, and the data were analyzed by one-way ANOVA. After the measurement of the biological parameters of the eyeball, the posterior pole tissues of the eyes were isolated, RNA was extracted respectively, and differentially expressed genes (DEGs) were screened by RNA-seq technology. GO analysis, KEGG pathway analysis, and PPI network analysis was performed to observe the effect of defocus intervention on the transcriptome of each tissue in the posterior pole of chicks. Results: The diopter, axial length, and choroidal thickness of the chicks in the minus-lens group and the plus lens group were significantly different from those in the control group (P<0.05); with P<0.05 and |log2FC|>1 as the screening criteria, compared with the control group and the minus-lens group, there were 203 DEGs in the retina, 757 DEGs in the choroid, and 1 509 DEGs in the sclera; compared with the control group and plus-lens group, the retina had 191 DEGs, the choroid had 378 DEGs, and the sclera had 1 918 DEGs. Compared with the control group, the overlapping genes of minus-lens group and plus-lens group all showed the same differential expression trend except for LOC121112941. There was more than 50% overlap in the GO analysis pathways. The common pathways in the KEGG analysis at the retinal level were arachidonic acid metabolism, tyrosine metabolism, and retinol metabolism; overlapping pathways at the choroidal level were protein processing in the endoplasmic reticulum, arginine and proline metabolism, and retinol metabolism; at the scleral level were cell adhesion molecules, ECM-receptor interactions, focal adhesions, cytokines - cytokine receptor interaction, sulfur metabolism, glycosaminoglycan biosynthesis - chondroitin sulfate/dermatan sulfate, Toll-like receptor signaling, TGF-β signaling, and MAPK signaling. Key genes in each tissue during out-of-focus recognition were identified by protein-protein interaction network analysis. Conclusion: This study screened out most of the overlapping genes in the late defocus recognition process at the transcriptome level, maintaining the same up- or down-regulation trend, and defocus signals in different directions may induce similar changes in some pathways.
孙丽媛, 朱莉, 陈思童, 等.. 基于RNA-seq技术的雏鸡晚期离焦识别转录组学研究[J]. 中华眼视光学与视觉科学杂志, 2022, 24(7): 494-505.
Liyuan Sun, Li Zhu, Sitong Chen,et al. Transcriptomic Study of Late Defocus Recognition in Chicks Based on RNA-seq Technology. Chinese Journal of Optometry Ophthalmology and Visual science, 2022, 24(7): 494-505. DOI: 1.3760/cma.j.cn115909-20220303-00080
Troilo D, Smith EL 3rd, Nickla DL, et al. IMI - report on experimental models of emmetropization and myopia. Invest Ophthalmol Vis Sci, 2019, 60(3): M31-M88. DOI: 10.1167/ iovs.18-25967.
[2]
Rucker F. Monochromatic and white light and the regulation of eye growth. Exp Eye Res, 2019, 184: 172-182. DOI: 10.1016/j.exer.2019.04.020.
[3]
Williams RA, Boothe RG. Effects of defocus on monkey (Macaca nemestrina) contrast sensitivity: behavioral measurements and predictions. Am J Optom Physiol Opt, 1983, 60(2): 106-111. DOI: 10.1097/00006324-198302000-00003.
[4]
Tkatchenko TV, Tkatchenko AV. Pharmacogenomic approach to antimyopia drug development: pathways lead the way. Trends Pharmacol Sci, 2019, 40(11): 833-852. DOI: 10.1016/ j.tips.2019.09.009.
[5]
Richardson R, Tracey-White D, Webster A, et al. The zebrafish eye-a paradigm for investigating human ocular genetics. Eye (Lond), 2017, 31(1):68-86. DOI: 10.1038/eye.2016.198.
[6]
Irving EL, Yakobchuk-Stanger C. Myopia progression control lens reverses induced myopia in chicks. Ophthalmic Physiol Opt, 2017, 37(5): 576-584. DOI: 10.1111/opo.12400.
[7]
Jiang X, Kurihara T, Kunimi H, et al. A highly efficient murine model of experimental myopia. Sci Rep, 2018, 8(1): 2026. DOI: 10.1038/s41598-018-20272-w.
[8]
Grytz R, Siegwart JT Jr. Changing material properties of the tree shrew sclera during minus lens compensation and recovery. Invest Ophthalmol Vis Sci, 2015, 56(3): 2065-2078. DOI: 10.1167/iovs.14-15352.
[9]
Graham B, Judge SJ. The effects of spectacle wear in infancy on eye growth and refractive error in the marmoset (Callithrix jacchus). Vision Res, 1999, 39(2): 189-206. DOI: 10.1016/ s0042-6989(98)00189-8.
[10]
Smith EL 3rd, Hung LF, Harwerth RS. Effects of optically induced blur on the refractive status of young monkeys. Vision Res, 1994, 34(3): 293-301. DOI: 10.1016/0042-6989(94)90088- 4.
[11]
Howlett MH, McFadden SA. Spectacle lens compensation in the pigmented guinea pig. Vision Res, 2009, 49(2): 219-227. DOI: 10.1016/j.visres.2008.10.008.
[12]
Dong Z, Chen Y. Transcriptomics: advances and approaches. Sci China Life Sci, 2013, 56(10): 960-967. DOI: 10.1007/s11427- 013-4557-2.
[13]
Li J, Du W, Xu N, et al. RNA-Seq analysis for exploring the pathogenesis of retinitis pigmentosa in P23H knockin mice. Ophthalmic Res, 2021, 64(5): 798-810. DOI: 10.1159/000515727.
[14]
Weiss S, Schaeffel F. Diurnal growth rhythms in the chicken eye: relation to myopia development and retinal dopamine levels. J Comp Physiol A, 1993, 172(3): 263-270. DOI: 10.1007/ BF00216608.
[15]
Guggenheim JA, Chen YP, Yip E, et al. Pre-treatment choroidal thickness is not predictive of susceptibility to form-deprivation myopia in chickens. Ophthalmic Physiol Opt, 2011, 31(5): 516- 528. DOI: 10.1111/j.1475-1313.2011.00827.x.
[16]
Wisely CE, Sayed JA, Tamez H, et al. The chick eye in vision research: an excellent model for the study of ocular disease. Prog Retin Eye Res, 2017, 61: 72-97. DOI: 10.1016/ j.preteyeres.2017.06.004.
[17]
Li R, Li Y, Kristiansen K, et al. SOAP: short oligonucleotide alignment program. Bioinformatics, 2008, 24(5): 713-714. DOI: 10.1093/bioinformatics/btn025.
[18]
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods, 2015,12(4):357- 360. DOI: 10.1038/nmeth.3317.
[19]
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods, 2012, 9(4): 357-359. DOI: 10.1038/ nmeth.1923.
[20]
Pertea M, Pertea GM, Antonescu CM, et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol, 2015, 33(3): 290-295. DOI: 10.1038/ nbt.3122.
[21]
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol, 2014, 15(12): 550. DOI: 10.1186/s13059-014-0550-8.
[22]
Li C, McFadden SA, Morgan I, et al. All-trans retinoic acid regulates the expression of the extracellular matrix protein fibulin-1 in the guinea pig sclera and human scleral fibroblasts. Mol Vis, 2010, 16:689-697.
[23]
Zuehlke AD, Beebe K, Neckers L, et al. Regulation and function of the human HSP90AA1 gene. Gene, 2015, 570(1): 8-16. DOI: 10.1016/j.gene.2015.06.018.
[24]
Wang J, Lee J, Liem D, et al. HSPA5 Gene encoding Hsp70 chaperone BiP in the endoplasmic reticulum. Gene, 2017, 618: 14-23. DOI: 10.1016/j.gene.2017.03.005.
[25]
Huang J, Li Y, Liu Z, et al. Transcriptomic responses to heat stress in rainbow trout Oncorhynchus mykiss head kidney. Fish Shellfish Immunol, 2018, 82: 32-40. DOI: 10.1016/ j.fsi.2018.08.002.
[26]
Frost MR, Norton TT. Differential protein expression in tree shrew sclera during development of lens-induced myopia and recovery. Mol Vis, 2007, 13:1580-1588.
[27]
Ma WQ, Sun XJ, Zhu Y, et al. PDK4 promotes vascular calcification by interfering with autophagic activity and metabolic reprogramming. Cell Death Dis, 2020, 11(11): 991. DOI: 10.1038/s41419-020-03162-w.
[28]
Zhang S, Hulver MW, McMillan RP, et al. The pivotal role of pyruvate dehydrogenase kinases in metabolic flexibility. Nutr Metab (Lond), 2014, 11(1): 10. DOI: 10.1186/1743-7075-11-10.
[29]
Tkatchenko TV, Troilo D, Benavente-Perez A, et al. Gene expression in response to optical defocus of opposite signs reveals bidirectional mechanism of visually guided eye growth. PLoS Biol, 2018, 16(10): e2006021. DOI: 10.1371/journal. pbio.2006021.
[30]
Simon P, Feldkaemper M, Bitzer M, et al. Early transcriptional changes of retinal and choroidal TGFbeta-2, RALDH-2, and ZENK following imposed positive and negative defocus in chickens. Mol Vis, 2004, 10: 588-597.
[31]
Schippert R, Schaeffel F, Feldkaemper MP. Microarray analysis of retinal gene expression in chicks during imposed myopic defocus. Mol Vis, 2008, 14:1589-1599.
[32]
Ashby RS, Feldkaemper MP. Gene expression within the amacrine cell layer of chicks after myopic and hyperopic defocus. Invest Ophthalmol Vis Sci, 2010, 51(7): 3726-3735. DOI: 10.1167/iovs.09-4615.
[33]
Yang J, Reinach PS, Zhang S, et al. Changes in retinal metabolic profiles associated with form deprivation myopia development in guinea pigs. Sci Rep, 2017, 7(1): 2777. DOI: 10.1038/s41598- 017-03075-3.
[34]
Zeng L, Li X, Liu J, et al. RNA-Seq analysis reveals an essential role of the tyrosine metabolic pathway and inflammation in myopia-induced retinal degeneration in Guinea pigs. Int J Mol Sci, 2021, 22(22): 12598. DOI: 10.3390/ijms222212598.
[35]
Mertz JR, Wallman J. Choroidal retinoic acid synthesis: a possible mediator between refractive error and compensatory eye growth. Exp Eye Res, 2000, 70(4): 519-527. DOI: 10.1006/ exer.1999.0813.
[36]
Bertrand E, Fritsch C, Diether S, et al. Identification of apolipoprotein A-I as a "STOP" signal for myopia. Mol Cell Proteomics, 2006, 5(11): 2158-2166. DOI: 10.1074/mcp. M600073-MCP200.
[37]
Yu FJ, Lam TC, Sze AY, et al. Alteration of retinal metabolism and oxidative stress may implicate myopic eye growth: evidence from discovery and targeted proteomics in an animal model. J Proteomics, 2020, 221: 103684. DOI: 10.1016/ j.jprot.2020.103684.
[38]
Liu K, Fang J, Jin J, et al. Serum metabolomics reveals personalized metabolic patterns for macular neovascular disease patient stratification. J Proteome Res, 2020, 19(2): 699-707. DOI: 10.1021/acs.jproteome.9b00574.
[39]
Flitcroft DI, Loughman J, Wildsoet CF, et al. Novel myopia genes and pathways identified from syndromic forms of myopia. Invest Ophthalmol Vis Sci, 2018, 59(1): 338-348. DOI: 10.1167/ iovs.17-22173.
[40]
Guo D, Ding M, Song X, et al. Regulatory roles of differentially expressed microRNAs in metabolic processes in negative Lensinduced myopia Guinea pigs. BMC Genomics, 2020, 21(1): 13. DOI: 10.1186/s12864-020-6447-x.
[41]
Zhu Y, Li W, Zhu D, et al. MicroRNA profiling in the aqueous humor of highly myopic eyes using next generation sequencing. Exp Eye Res, 2020, 195: 108034. DOI: 10.1016/ j.exer.2020.108034.
[42]
Moring AG, Baker JR, Norton T T. Modulation o f glycosaminoglycan levels in tree shrew sclera during lensinduced myopia development and recovery. Invest Ophthalmol Vis Sci, 2007, 48(7): 2947-2956. DOI: 10.1167/iovs.06-0906.
[43]
Kumar S, Kim J. PLK-1 Targeted inhibitors and their potential against tumorigenesis. Biomed Res Int, 2015, 2015: 705745. DOI: 10.1155/2015/705745. [44] Zhang X, Zhang X, Li X, et al. Connection between CDC20 expression and hepatocellular carcinoma prognosis. Med Sci Monit, 2021, 27:e926760. DOI: 10.12659/MSM.926760.
[45]
Ustinov NB, Korshunova AV, Gudimchuk NB. Protein complex NDC80: properties, functions, and possible role in pathophysiology of cell division. Biochemistry (Mosc), 2020, 85(4): 448-462. DOI: 10.1134/S0006297920040057.
[46]
Zhang Y, Zhou L, Wang S, et al. Exploration of retinoblastoma pathogenesis with bioinformatics. Transl Cancer Res, 2021, 10(7): 3527-3537. DOI: 10.21037/tcr-21-1034.
[47]
Kim T, Gartner A. Bub1 kinase in the regulation of mitosis. Anim Cells Syst (Seoul), 2021, 25(1): 1-10. DOI: 10.1080/1976 8354.2021.1884599.