Integrative RNA-seq and ATAC-seq Analysis of Retinitis Pigmentosa Caused by PDE6 Mutation

Integrative RNA-seq and ATAC-seq Analysis of Retinitis Pigmentosa Caused by PDE6 Mutation | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Integrative RNA-seq and ATAC-seq Analysis of Retinitis Pigmentosa Caused by PDE6 Mutation Wenrong Xu, Yan Li, Xiaoxiao Feng, Libo Xiao, Lan Li, Kangwei Jiao This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-735055/v1 This work is licensed under a CC BY 4.0 License Status: Under Review Version 1 posted 4 You are reading this latest preprint version Abstract Purpose : Inhibition of PARP1 could relieve PDE6 mutation-induced Retinitis pigmentosa (RP). However, the mechanism related with PARP1 overexpression in the RP has not been clarified. We attempted to explore the potential regulatory mechanism related with PARP1 underlying RP. Methods : ATAC-seq and RNA-seq were preformed for retina tissues of C3H and rd1 mice. The differential expressed genes (DEGs) were identified, followed by PARP1-DEG coexpression network, and PPI network construction. GO-BP and pathway enrichment of DEGs of interest were performed by clusterprofiler software. The overlapped genes that might play regulatory role in PARP1 expression were mined by integrated analysis of RNA-seq and ATAC-seq data. Results : Total 1061 DEGs were identified between C3H and rd1 group. Co-expression network was constructed with 313 PARP1-gene coexpression pairs. The down-regulated DEGs were closely related with visual perception, light stimulus related biological process, while the up-regulated DEGs were significantly enriched in phototransduction and PPAR signaling pathway. PPI network was constructed with 202 nodes and 375 edges, which was clustered into 3 modules. Module 1 genes were closely related with detection of light stimulus, visual perception related biological process and phototransduction pathway (involved with Gnat1/Guca1b/Gnat2/Sag/Pde6g). By integrated analysis of the RNA-seq and ATAC-seq, the overlapped up-regulated genes were Asxl3 and Nyap2, while the down-regulated genes were Tmem136 and Susd3. Conclusion : Gnat1 may play a key role in RP development by interacting with PARP1. Susd3 may play a regulatory role in PARP1 expression and affect RP formation. Ophthalmology ATAC-seq RNA-seq Retinitis pigmentosa PARP1 retina tissues Figures Figure 1 Figure 2 Figure 3 Figure 4 Highlights ATAC-seq and RNA-seq were preformed for retina tissues of C3H and rd1 mice. PARP1-gene coexpression network was constructed. Susd3 was a overlapped down-regulated gene by ATAC-seq and RNA-seq. Gnat1 was a module 1 gene with highest nod degree. Introduction Retinitis pigmentosa (RP) is a genetic retinal neurodegenerative disease, characterized by the pathogenic gene mutation and progressive apoptosis of retina photoreceptor cells. RP is one of the most common causes of blindness and severe visual decline [ 1 ]. It is estimated that the incidence rate of RP is about 1/4000 worldwide [ 2 ] and it mostly common affects people aged 20 to 60 years [ 3 ]. Despite the advances in the treatment measures for RP, such as gene therapy [ 4 ], stem cell therapy [ 5 ], there is no cure for RP currently. RP gene mutation is responsible for majority of RP cases, which is associated with retinal molecular pathway defects, leading to the final photoreceptor cell apoptosis. Thus, targeting the prevention of photoreceptor cell apoptosis may shed light on the effective therapy for RP. Many efforts have been made to the deep understanding of the photoreceptor cell apoptosis underlying RP. It is reported that PDE (Phosphodiesterases) and PARP1 (poly-ADP-ribose polymerase) play a key role in the apoptosis process of photoreceptor cell [ 6 ]. PDE plays a regulatory role in signal transduction. The PDE defects-induced by RP gene mutation could cause the accumulation of cGMP, which posed toxicity to photoreceptor cell [ 7 ]. PDE6 mutation is responsible for 4–8% RP in human homologous animal models [ 7 , 8 ]. In addition, PARP1 exerts various biological functions in cell apoptosis and gene expression regulation [ 9 ]. PARP-1 is involved in the chromatin function regulation by interacting with transcriptional factor (TF), DNA methyltransferas 1, and histone modifying enzymes [ 10 , 11 ]. Previous evidence has revealed that PARP1 is proposed as the therapeutic target for intervening retinal photoreceptor function and survival in RP induced by PDE6A mutation [ 12 ], indicating the causal relationship between PARP1 dysregulation and PDE6 gene mutations. The previous study for the progression of retinal degeneration in RP mice model with PDE6A mutation suggested that the apoptotic photoreceptors cells caused by PDE6 mutation peaked at postnatal 15 day (P15) [ 13 ]. RP caused by PDE6 dysfunction is primarily studied in rd1 mice with Pde6b gene mutations [ 14 ]. Therefore, in this study, we performed RNA-seq and ATAC-seq for retina tissues from PDE6B mutant mice (Rd1) and C3H wild type (wt) mice. The differential accessible peaks and differential expressed genes were analyzed. We aimed to mine the regulatory mechanism underlying PARP1 overexpression and photoreceptors cells apoptosis. Methods Animals and sampling The C3H wild type (wt) mice during the first 31 postnatal days and age-paired Rd1 (Pde6b rd1) mice with Pde6b mutation were presented by Tubingen university. All the specific pathogen free (SPF) animals were maintained in the IVC cages with food and water freely. The mice were acclimatized to experimental environment under a 12:12 light/dark cycle at 22°C. The C3H mice and Rd1 mice at P15 were sacrificed, respectively. Then, the mice retinal tissues were extracted and the single cells were isolated according the previous description [12]. Our study was approved by the Ethics Committee of Yunnan University (No. YNUCARE20210024) and Yunnan University affiliated Hospital (No. 2021045). All the animal experimental procedures conformed to Animal Care guidelines. ATAC sequencing and processing For ATAT-seq, 50,000 single cells were prepared and collected by centrifugation at 500 g for 5 min at 4°C. Cell nuclei extraction was achieved with cold lysis buffer followed by transposase reaction with Nextera kit. After purified, DNA was amplified by PCR analysis in a 50 µl reaction mixture with the reaction condition of 1 cycle of 72°C for 5 min, 98°C for 30 s, followed by 5 cycles of 98°C for 10s, 63°C fro 30s and 2°C for 2 min. After purified, the generated DNA libraries was pooled and subjected to Illumina HiSeq. The quality of raw data was controlled by FastQC. The pair-end reads were mapped to reference genome sequences by the bwa program. Peak calling was performed by MACS2 software with q value 1 and P value < 0.05. The differential Peak adjacent genes were annotated by annotatePeaks.pl tool in HOMER. RNA sequencing and processing Total RNA was extracted from retinal tissues and the quality control was achieved by Nanodrop and Agilent 2100. After rRNA was removed, RNA samples were fragmented by fragmentation buffer. The first strand of cDNA was synthesized with random hexamers, followed by the second strand synthesis with dNTPs, DNA polymerase and RNase H. The cDNAs were subjected to end repair, A tail addition, and the second strand (containing U) degradation. The cDNA library was constructed by PCR amplification and sequenced based on HiSeq platform. The raw data were filtered by removing the reads with adapter, >10% N bases and low quality. The clean reads were mapped to reference genome based on modified BWT algorithm by Hisat2 software. The gene expression levels were analyzed by reads count located at genomic regions or exon regions of genes according to FPKM method with the application of HTSeq software. The differentially expressed genes (DEGs) between C3H and Rd1 samples were screened by DESeq2 software with cutoff value of logFC > 0.585 and adjust P value < 0.05. Coexpression of PARP1 with DEGs The expression correlation between PARP1 and DEGs was estimated by Pearson correlation coefficient (r). The absolute value of r ranged from 0 to 1. When the absolute of r was closed to 1, the most significant correlation between the PARP1 and DEG pair was considered. The coexpression gene pairs were collected with |r| >0.9 and p value <0.05, followed by coexpression network construction. Protein-protein interaction (PPI) network and module analysis The protein interaction pairs were retrieved from STRING[ 15 ] (Version:10.0, http://www.string-db.org/). The DEGs were mapped to protein pairs and those with PPI score ≥0.4 were collected for PPI network construction by Cytoscape (version:3.2.0, http://www.cytoscape.org/). The topological properties of protein interaction network were analyzed by CytoNCA plugin [ 16 ] (Version 2.1.6, http://apps.cytoscape.org/apps/cytonca). The inteconnections of one gene with others were scored and the hub proteins in network was screened with highest node degrees. The functional modules in PPI network were classified by MCODE Cytoscape plugin [ 17 ] (Version1.4.2, http://apps.cytoscape.org/apps/MCODE). The threshold value was set as score＞5. GO and pathway enrichment analysis The DEGs in coexpression network and PPI network were subjected to GO-BP[ 18 ] and KEGG pathway [ 19 ] enrichment analysis with the extensively used tool of clusterprofiler [ 20 ] (Version3.18.1,http://www.bioconductor.org/packages/release/bioc/html/clusterProfiler.html). P value <0.05 was set as the cutoff value. Integrated analysis of RNA-seq and ATAC-seq data The analysis of ATAC-seq data unravels transcription factor binding sites. To explore the regulatory role of TFs in gene expression, A TAC-seq and RNA-seq data were integrated. The down-regulated genes obtained by RNA-seq data were overlapped with those around Peaks with depressed reads signals analyzed by A TAC-seq data, otherwise the up-regulated genes were overlapped. In the present study, the differential Peak adjacent genes analyzed by ATAC-seq data were overlapped with nodes in co-expression network. The overlaps were subjected to functional enrichment analysis. Results Identification of DEGs The genes that may be perturbed by Pde6b mutation and play a key role in the development of photoreceptor cells in the retina were identified by RNA-seqbn n . Total 1061 genes were found to be differentially expressed in Rd1 samples compared with C3H ones at P15, including 340 up-regulated genes and 721 down-regulated ones. The gene differential expression was illustrated by heatmap, which showed differential expression profile of DEGs between Rd1 and C3H samples (Figure 1A). The volcano plot represented differential expressed genes based on the defined p value and fold change cutoff values (Figure 1B). Coexpression network of PARP1 with DEGs Based on Pearson correlation analysis, 313 PARP1-DEG coexpressed pairs were obtained, including 102 pairs with positive correlation and 211 ones with negative expression. The coexpression network was illustrated in Figure 2, which was comprised with 102 up-regulated genes and 211 down-regulated genes. Function enrichment analysis was performed for up- and down-regulated genes, respectively. Results showed that the up-regulated genes were closely related with tube formation related biological process, PPAR signaling pathway and B cell receptor signaling pathway (Table 1). The down-regulated genes were significantly enriched in visual perception and sensory perception of light stimulus related biological process. The most significant pathway for down-regulated genes was phototransduction (mmu04744) (Table 1). PPI network and the functional modules The PPI network was constructed with 202 nodes connected with 375 edges (Figure 3A). Three significant network modules were divided by MCODE with PPI score ≥ 5, including Module1(score=8.222, 10 nodes and 37 edges), Module 2(score=6.333, 7 nodes and 19 edges) and Module 3 (7 nodes and 18 edges) (Figure 3B). After topological property analysis, the node degrees of modular genes were obtained (Table 2). Functional analysis showed that the genes in Module 1 were closely related with detection of light stimulus (GO:0009583) and visual perception (GO:0007601) related biological process and involved with Phototransduction pathway (mmu04744). The Module 2 were correlated with regulation of mitotic cell cycle (GO:0007346) and chromosome segregation (GO:0007059) related biological process and DNA replication (mmu03030) and cell cycle (mmu04110) related pathways. Additionally, the genes in module 3 were significantly enriched in behavioral fear response (GO:0001662), regulation of amino acid transport (GO:0051955), Neuroactive ligand-receptor interaction (mmu04080), PI3K-Akt signaling pathway (mmu04151) (Table 3). ATAC-seq and integrated analysis Based on ATAC-seq, we obtained 1394 peaks up-regulated signals and 891 differential peaks with down-regulated signal. The differential accessible peaks covered genes were overlapped with DEGs in co-expression network. The Venn diagram depicted that the overlapped up-regulated genes were Asxl3 and Nyap2 (Figure 4A, while the down-regulated overlaps were Tmem136 and Susd3 (Figure 4B). Discussion RP is one of the leading cause of inherited blindness, which has affected 1/3000 to 1/5000 individuals globally. PDE6B mutation is responsible for 5–8% RP patients, which affects rod-specific cyclic guanosine monophosphate (cGMP) level and results in rod degeneration [ 21 , 22 ]. Recent evidence showed that PARP1 was a promising therapeutic target for PDE6 mutation-induced RP. However, the regulatory mechanism of PARP1 in photoreceptor cell apoptosis underlying RP has not been clarified. In this study, we attempted to explore the potential mechanism related with photoreceptor cell apoptosis in PDE6 mutation-induced RP with ATAC-seq and RNA-seq. Our results of RNA-seq showed that there were 1061 DEGs between Rd1 and C3H mice at P15. The expression profile of DEGs could distinguish the samples from Rd1 and C3H clearly, which revealed that DEGs identified were significant and could play key roles in the pathogenesis of RP. In order to explore the mechanism related with PARP1 in RP, PARP1-DEG coexpression network was constructed. Functional analysis showed that the down-regulated genes with coexpression correlation with PARP1 were closely related with visual perception, sensory perception of light stimulus related biological process and phototransduction pathway. The manifestations of RP were characterized by loss of depth perception, night blindness, and loss of central vision [ 23 ]. The variety of eye disorders in RP patients was related with rod-opsin protein disruption, in parallel with the defects of phototransduction cascade of central nervous [ 24 ]. Our findings were consistent with the reports mentioned above and revealed that the visual perception and phototransduction pathway were dysregulated in RP. Additionally, our data showed that PPAR signaling pathway was a significant pathway enriched by up-regulated genes. It is reported that PPAR plays a regulatory role in the transcription of genes involved in intracellular signaling pathways [ 25 ]. PPAR is involved in lipid metabolism and homeostasis [ 26 ]. Accumulating evidence indicates that abnormal lipid metabolism is associated with RP development [ 26 , 27 ]. In this study, the PPAR pathway was found to be enhanced by up-regulated genes. Thus, we speculated that PARP1 could be associated with abnormal lipid metabolism induced photoreceptor cell apoptosis in RP. Besides, the PPI network indicated that the Module 1 genes were also enriched in phototransduction pathway, which was in agree with GO-BP analysis of DEGs in coexpression network. Gnat 1 was a down-regulated gene in retinal tissues of Rd1 mice with highest degree of 19 in Module 1 PPI network. As described in previous study, Gnat1 encodes the α-subunit of transducin in rod cells. Gnat1 variants were previously reported to be associated with congenital stationary night-blindness [ 28 ]. New evidence showed that Gnat1 mutation was responsible for the slow and late onset of retinal degeneration [ 29 ]. Similarly, our GO-BP enrichment analysis showed that Gnat1 was closely related with visual perception and light stimulus related biological process. Gnat1 showed negative expression correlation with PARP1 (r=-0.90 and p = 0.013). All these above suggested that the down-regulation of Gnat1 was correlated with PARP1 overexpression, which caused the dysfunction of visual perception, light stimulus, and phototransduction, leading to the onset of RP. Furthermore, the RNA-seq and ATAC-seq integrated analysis in our study revealed that Susd3 (sushi domain-containing protein 3) was a down-regulated overlap. Susd3 is cell surface protein and has been suggested as the potential biomarker for breast cancer [ 30 ]. Besides, Susd3 was found to be one of the genes encoding protein at the surface of young rod photoreceptors [ 31 ]. Susd2 was identified to be specially expressed for photoreceptors in the retina. However, the significant role of Susd3 in RP onset has not been clarified. ATAC-seq is a simple protocol for open chromatin analysis, which provides possibility of discovering TF binding sites. Susd3 was the gene covered by differential peaks by ATAC-seq, which indicated that Susd3 could exert TF function. Susd3 showed coexpression correlation with PARP1. Currently, little is known about whether Susd3 elicited transcriptional regulation on PARP1 and affected PDE6-mutation induced RP, which should be verified in the following studies. In summary, the down-regulated genes with coexpression correlation of PARP1 were significantly enriched in visual perception, sensory perception of light stimulus related biological process and phototransduction pathway. The PPAR signaling pathway was enhanced by up-regulated genes. Gnat1 could play a key role in visual perception, light stimulus, and phototransduction of RP by interacting with PARP1. Susd3 might play a regulatory role in PARP1 expression and affect RP formation. Our study may provide new insights in understanding the mechanism of PDE6-mutation induced RP. References Campochiaro, P.A. and T.A. Mir, The mechanism of cone cell death in Retinitis Pigmentosa. Prog Retin Eye Res, 2018. 62 : p. 24-37. Zhang, Q., Retinitis Pigmentosa: Progress and Perspective. Asia Pac J Ophthalmol (Phila), 2016. 5 (4): p. 265-71. Parmeggiani, F., Clinics, epidemiology and genetics of retinitis pigmentosa. Curr Genomics, 2011. 12 (4): p. 236-7. Schon, C., et al., Gene therapy successfully delays degeneration in a mouse model of PDE6A-linked retinitis pigmentosa (RP 43). Hum Gene Ther, 2017. Kim, J.Y., et al., Epiretinal Membrane Formation after Intravitreal Autologous Stem Cell Implantation in a Retinitis Pigmentosa Patient. Retin Cases Brief Rep, 2017. 11 (3): p. 227-231. Kaur, J., et al., Calpain and PARP activation during photoreceptor cell death in P23H and S334ter rhodopsin mutant rats. PLoS One, 2011. 6 (7): p. e22181. Gopalakrishna, K.N., K. Boyd, and N.O. Artemyev, Mechanisms of mutant PDE6 proteins underlying retinal diseases. Cell Signal, 2017. 37 : p. 74-80. Dvir, L., et al., Autosomal-recessive early-onset retinitis pigmentosa caused by a mutation in PDE6G, the gene encoding the gamma subunit of rod cGMP phosphodiesterase. Am J Hum Genet, 2010. 87 (2): p. 258-64. Andrabi, S.A., et al., Poly(ADP-ribose) (PAR) polymer is a death signal. Proc Natl Acad Sci U S A, 2006. 103 (48): p. 18308-13. Hottiger, M.O., ADP-ribosylation of histones by ARTD1: an additional module of the histone code? FEBS Lett, 2011. 585 (11): p. 1595-9. Thomas, C. and A.V. Tulin, Poly-ADP-ribose polymerase: machinery for nuclear processes. Mol Aspects Med, 2013. 34 (6): p. 1124-37. Jiao, K., et al., Efficacy of PARP inhibition in Pde6a mutant mouse models for retinitis pigmentosa depends on the quality and composition of individual human mutations. Cell Death Discov, 2016. 2 : p. 16040. Sothilingam, V., et al., Retinitis pigmentosa: impact of different Pde6a point mutations on the disease phenotype. Hum Mol Genet, 2015. 24 (19): p. 5486-99. Keeler, C.E., The Inheritance of a Retinal Abnormality in White Mice. Proc Natl Acad Sci U S A, 1924. 10 (7): p. 329-33. Szklarczyk, D., et al., STRING v10: protein–protein interaction networks, integrated over the tree of life. Nucleic acids research, 2014: p. gku1003. Tang, Y., et al., CytoNCA: A cytoscape plugin for centrality analysis and evaluation of protein interaction networks. BioSystems, 2015. 127 (Complete): p. 67-72. Bandettini, W.P., et al., MultiContrast Delayed Enhancement (MCODE) improves detection of subendocardial myocardial infarction by late gadolinium enhancement cardiovascular magnetic resonance: a clinical validation study. J Cardiovasc Magn Reson, 2012. 14 : p. 83. Ashburner, M., et al., Gene Ontology: tool for the unification of biology. Nat Genet, 2000. 25 (1): p. 25-29. Kanehisa, M. and S. Goto, KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Research, 2000. 28 (1): p. 27-30. Yu, G., et al., clusterProfiler: an R package for comparing biological themes among gene clusters. Omics, 2012. 16 (5): p. 284-7. Kim, M.S. and K. Joo, Genetic Mutation Profiles in Korean Patients with Inherited Retinal Diseases. 2019. 34 (21): p. e161. Yeo, J.H., et al., Development of a Pde6b Gene Knockout Rat Model for Studies of Degenerative Retinal Diseases. Invest Ophthalmol Vis Sci, 2019. 60 (5): p. 1519-1526. Prem Senthil, M., J. Khadka, and K. Pesudovs, Seeing through their eyes: lived experiences of people with retinitis pigmentosa. Eye (Lond), 2017. 31 (5): p. 741-748. Mendes, H.F., et al., Mechanisms of cell death in rhodopsin retinitis pigmentosa: implications for therapy. Trends Mol Med, 2005. 11 (4): p. 177-85. Yanagi, Y., Role of Peoxisome Proliferator Activator Receptor gamma on Blood Retinal Barrier Breakdown. PPAR Res, 2008. 2008 : p. 679237. Donato, L., et al., Transcriptome Analyses of lncRNAs in A2E-Stressed Retinal Epithelial Cells Unveil Advanced Links between Metabolic Impairments Related to Oxidative Stress and Retinitis Pigmentosa. Antioxidants (Basel), 2020. 9 (4). Holman, R.T., et al., Abnormal plasma lipids of patients with Retinitis pigmentosa. Lipids, 1994. 29 (1): p. 61-5. Szabo, V., et al., p.Gln200Glu, a putative constitutively active mutant of rod alpha-transducin (GNAT1) in autosomal dominant congenital stationary night blindness. Hum Mutat, 2007. 28 (7): p. 741-2. Carrigan, M. and E. Duignan, A novel homozygous truncating GNAT1 mutation implicated in retinal degeneration. 2016. 100 (4): p. 495-500. Yu, Z., et al., Sushi Domain-Containing Protein 3: A Potential Target for Breast Cancer. Cell Biochemistry and Biophysics, 2015. 72 (2): p. 321-324. Postel, K., et al., Analysis of cell surface markers specific for transplantable rod photoreceptors. Mol Vis, 2013. 19 : p. 2058-67. Cite Share Download PDF Status: Under Review Version 1 posted Reviews received at journal 27 Aug, 2021 Reviewers invited by journal 27 Aug, 2021 Editor assigned by journal 20 Jul, 2021 First submitted to journal 19 Jul, 2021 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-735055","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":48739085,"identity":"5513c45f-54e2-44ec-a5b4-c699fdc9f279","order_by":0,"name":"Wenrong Xu","email":"","orcid":"","institution":"Kunming Medical University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Wenrong","middleName":"","lastName":"Xu","suffix":""},{"id":48739086,"identity":"34e5afeb-683d-436a-8738-cfbab91f2301","order_by":1,"name":"Yan Li","email":"","orcid":"","institution":"Affiliated Hospital of Yunnan University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Yan","middleName":"","lastName":"Li","suffix":""},{"id":48739087,"identity":"d348e843-0255-4e94-b645-948814e4bd64","order_by":2,"name":"Xiaoxiao Feng","email":"","orcid":"","institution":"Affiliated Hospital of Yunnan University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Xiaoxiao","middleName":"","lastName":"Feng","suffix":""},{"id":48739088,"identity":"30dc69c1-88a7-4011-b504-145e88e0ef24","order_by":3,"name":"Libo Xiao","email":"","orcid":"","institution":"Affiliated Hospital of Yunnan University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Libo","middleName":"","lastName":"Xiao","suffix":""},{"id":48739089,"identity":"148a7190-6854-434d-9646-2b20ab6a3f4b","order_by":4,"name":"Lan Li","email":"","orcid":"","institution":"The Affiliated Calmette Hospital of Kunming Medical University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Lan","middleName":"","lastName":"Li","suffix":""},{"id":48739090,"identity":"9619dc65-99e5-480a-a4d8-e55e019eb972","order_by":5,"name":"Kangwei Jiao","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA/0lEQVRIiWNgGAWjYBACgwNgSoKBsRlIfYAK4tVi2ICkhXEGMVqMkTnMPMRoMWM/e/g1T4UFA3M788PHtm13EhvYm7dJMNTcwanFhicvzZrnDMhhbMbGuW3PEht4jpVJMBx7hlsLQ44ZUCXYL2bSuW2HExskcswkGBsO43YY/xugln8gLezfpC1BWuTf4NdiLJFj/Di3AaSFx0yaEWwLD34thjPemDH/OSbBA9RSbNhz7rBxG09asUXCMdxaDM7nGH+cUVMnZ9h/fOODH2WHZfvZD2+88aEGtxYgYJMAEjzgOGVkA3JBYgn4NAAjEJxM5MHsP/iVjoJRMApGwcgEAACPTwK2/QBeAAAAAElFTkSuQmCC","orcid":"","institution":"Affiliated Hospital of Yunnan University","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"Kangwei","middleName":"","lastName":"Jiao","suffix":""}],"badges":[],"createdAt":"2021-07-20 05:48:47","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-735055/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-735055/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":12908748,"identity":"13f15944-4da1-4c05-9ea5-6fb8ca4234dd","added_by":"auto","created_at":"2021-08-30 20:36:22","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":343397,"visible":true,"origin":"","legend":"Identification of differentially expressed genes by RNA-seq\nA, Heatmap for differentially expressed genes; B, volcano-plot for differentially expressed genes. Top 5 genes with highest fold changes were listed.\n","description":"","filename":"Figure1.png","url":"https://assets-eu.researchsquare.com/files/rs-735055/v1/a5f8d8e6376d03f1350b79ef.png"},{"id":12908749,"identity":"737b3f31-b2f1-4274-882c-c3220aeb14d5","added_by":"auto","created_at":"2021-08-30 20:36:22","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":1678080,"visible":true,"origin":"","legend":"PARP1-DEGs coexpression network \nRed triangle represents as PARP1, yellow circle indicates up-regulated genes, green prismatic indicates down-regulated genes, blue lines mean negative correlation and red lines mean positive correlation. \n","description":"","filename":"Figure2.png","url":"https://assets-eu.researchsquare.com/files/rs-735055/v1/1c627927e393db3dfa0b7f50.png"},{"id":12908751,"identity":"1effe054-9606-4629-adf3-a5bafeaac7e3","added_by":"auto","created_at":"2021-08-30 20:36:22","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":1160777,"visible":true,"origin":"","legend":"Protein-protein interaction network and its modules \nA, PPI network of DEGs were constructed with 202 nodes and 375 edges. B, PPI network was clustered into 3 modules, such as module-1(score=8.222), module-2(score=6.333) and module-3(score=6). \nRed triangle represents as PARP1, yellow circle indicates up-regulated genes, green prismatic indicates down-regulated genes. The node size is based on degree value.\n","description":"","filename":"Figure3.png","url":"https://assets-eu.researchsquare.com/files/rs-735055/v1/833fca61857ac4486404ac73.png"},{"id":12908750,"identity":"e8873da7-13a6-4a94-b475-1a201504a64e","added_by":"auto","created_at":"2021-08-30 20:36:22","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":222811,"visible":true,"origin":"","legend":"Integrated analysis of ATAC-seq and RNA-seq data \nThe overlapped genes by ATAC-seq and RNA-seq data were analyzed by Venn diagram. A, the overlapped genes for up-regulated genes; B, the overlaps for down-regulated genes. \n","description":"","filename":"Figure4.png","url":"https://assets-eu.researchsquare.com/files/rs-735055/v1/609616a91f18f6e309e2bc24.png"},{"id":13712596,"identity":"2da1f9e7-31a7-486b-accd-6bfe1a8e8da1","added_by":"auto","created_at":"2021-09-17 14:28:58","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":2190746,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-735055/v1/4652bce3-8502-46db-9eea-63f43ff53f5a.pdf"}],"financialInterests":"","formattedTitle":"\u003cp\u003eIntegrative RNA-seq and ATAC-seq Analysis of Retinitis Pigmentosa Caused by PDE6 Mutation\u003c/p\u003e","fulltext":[{"header":"Highlights","content":"\u003cul\u003e\n \u003cli\u003eATAC-seq and RNA-seq were preformed for retina tissues of C3H and rd1 mice.\u003c/li\u003e\n \u003cli\u003ePARP1-gene coexpression network was constructed.\u003c/li\u003e\n \u003cli\u003eSusd3 was a overlapped down-regulated gene by\u0026nbsp;ATAC-seq and RNA-seq.\u003c/li\u003e\n \u003cli\u003eGnat1 was a module 1 gene with highest nod degree. \u003c/li\u003e\n\u003c/ul\u003e"},{"header":"Introduction","content":"\u003cp\u003eRetinitis pigmentosa (RP) is a genetic retinal neurodegenerative disease, characterized by the pathogenic gene mutation and progressive apoptosis of retina photoreceptor cells. RP is one of the most common causes of blindness and severe visual decline [\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e]. It is estimated that the incidence rate of RP is about 1/4000 worldwide [\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e] and it mostly common affects people aged 20 to 60 years [\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e]. Despite the advances in the treatment measures for RP, such as gene therapy [\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e], stem cell therapy [\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e], there is no cure for RP currently. RP gene mutation is responsible for majority of RP cases, which is associated with retinal molecular pathway defects, leading to the final photoreceptor cell apoptosis. Thus, targeting the prevention of photoreceptor cell apoptosis may shed light on the effective therapy for RP.\u003c/p\u003e \u003cp\u003eMany efforts have been made to the deep understanding of the photoreceptor cell apoptosis underlying RP. It is reported that PDE (Phosphodiesterases) and PARP1 (poly-ADP-ribose polymerase) play a key role in the apoptosis process of photoreceptor cell [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. PDE plays a regulatory role in signal transduction. The PDE defects-induced by RP gene mutation could cause the accumulation of cGMP, which posed toxicity to photoreceptor cell [\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e]. PDE6 mutation is responsible for 4\u0026ndash;8% RP in human homologous animal models [\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e, \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e]. In addition, PARP1 exerts various biological functions in cell apoptosis and gene expression regulation [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e]. PARP-1 is involved in the chromatin function regulation by interacting with transcriptional factor (TF), DNA methyltransferas 1, and histone modifying enzymes [\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e, \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e]. Previous evidence has revealed that PARP1 is proposed as the therapeutic target for intervening retinal photoreceptor function and survival in RP induced by PDE6A mutation [\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e], indicating the causal relationship between PARP1 dysregulation and PDE6 gene mutations.\u003c/p\u003e \u003cp\u003eThe previous study for the progression of retinal degeneration in RP mice model with PDE6A mutation suggested that the apoptotic photoreceptors cells caused by PDE6 mutation peaked at postnatal 15 day (P15) [\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e]. RP caused by PDE6 dysfunction is primarily studied in rd1 mice with Pde6b gene mutations [\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e]. Therefore, in this study, we performed RNA-seq and ATAC-seq for retina tissues from PDE6B mutant mice (Rd1) and C3H wild type (wt) mice. The differential accessible peaks and differential expressed genes were analyzed. We aimed to mine the regulatory mechanism underlying PARP1 overexpression and photoreceptors cells apoptosis.\u003c/p\u003e"},{"header":"Methods","content":"\u003cp\u003e\u003cstrong\u003eAnimals and sampling\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe C3H wild type (wt) mice during the first 31 postnatal days and age-paired Rd1 (Pde6b rd1) mice with Pde6b mutation were presented by Tubingen university. All the specific pathogen free (SPF) animals were maintained in the IVC cages with food and water freely. The mice were acclimatized to experimental environment under a 12:12 light/dark cycle at 22\u0026deg;C. The C3H mice and Rd1 mice at P15 were sacrificed, respectively. Then, the mice retinal tissues were extracted and the single cells were isolated according the previous description [12]. Our study was approved by the Ethics Committee of Yunnan University (No. YNUCARE20210024) and Yunnan University affiliated Hospital (No. 2021045). All the animal experimental procedures conformed to Animal Care guidelines.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eATAC sequencing and processing\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eFor ATAT-seq, 50,000 single cells were prepared and collected by centrifugation at 500 g for 5 min at 4\u0026deg;C. Cell nuclei extraction was achieved with cold lysis buffer followed by transposase reaction with Nextera kit. After purified, DNA was amplified by PCR analysis in a 50 \u0026micro;l reaction mixture with the reaction condition of 1 cycle of 72\u0026deg;C for 5 min, 98\u0026deg;C for 30 s, followed by 5 cycles of 98\u0026deg;C for 10s, 63\u0026deg;C fro 30s and 2\u0026deg;C for 2 min. After purified, the generated DNA libraries was pooled and subjected to Illumina HiSeq. The quality of raw data was controlled by FastQC. The pair-end reads were mapped to reference genome sequences by the bwa program. Peak calling was performed by MACS2 software with q value \u0026lt;0.05. The differential peaks (open chromatin regions) between C3H and Rd1 samples were analyzed by DEGseq with the cutoff value of |log2FC| \u0026gt; 1 and P value \u0026lt; 0.05. The differential Peak adjacent genes were annotated by annotatePeaks.pl tool in HOMER.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eRNA sequencing and processing\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTotal RNA was extracted from retinal tissues and the quality control was achieved by Nanodrop and Agilent 2100. After rRNA was removed, RNA samples were fragmented by fragmentation buffer. The first strand of cDNA was synthesized \u0026nbsp;with random hexamers, followed by the second strand synthesis with dNTPs, DNA polymerase and RNase H. The cDNAs were subjected to end repair, A tail addition, and the second strand (containing U) degradation. The cDNA library was constructed by PCR amplification and sequenced based on HiSeq platform. The raw data were filtered by removing the reads with adapter, \u0026gt;10% N bases and low quality. The clean reads were mapped to reference genome based on modified BWT algorithm by Hisat2 software. The gene expression levels were analyzed by reads count located at genomic regions or exon regions of genes according to FPKM method with the application of HTSeq software. The differentially expressed genes (DEGs) between C3H and Rd1 samples were screened by DESeq2 software with cutoff value of logFC \u0026gt; 0.585 and adjust P value \u0026lt; 0.05.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCoexpression of PARP1 with DEGs\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe expression correlation between PARP1 and DEGs was estimated by Pearson correlation coefficient (r). The absolute value of r ranged from 0 to 1. When the absolute of r was closed to 1, the most significant correlation between the PARP1 and DEG pair was considered. The coexpression gene pairs were collected with |r| \u0026gt;0.9 and p value \u0026lt;0.05, followed by coexpression network construction.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eProtein-protein interaction (PPI) network and module analysis \u0026nbsp;\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe protein interaction pairs were retrieved from\u0026nbsp;STRING[\u003ca href=\"#_ENREF_15\" title=\"Szklarczyk, 2014 #34\"\u003e15\u003c/a\u003e]\u0026nbsp;(Version:10.0,\u0026nbsp;http://www.string-db.org/). The DEGs were mapped to protein pairs and those with PPI score \u0026ge;0.4 were collected for PPI network construction by\u0026nbsp;Cytoscape (version:3.2.0, http://www.cytoscape.org/). The topological properties of protein interaction network were analyzed by CytoNCA plugin\u0026nbsp;[\u003ca href=\"#_ENREF_16\" title=\"Tang, 2015 #46\"\u003e16\u003c/a\u003e]\u0026nbsp;(Version 2.1.6, http://apps.cytoscape.org/apps/cytonca). The inteconnections of one gene with others were scored and the hub proteins in network was screened with highest node degrees. The functional modules in PPI network were classified by MCODE Cytoscape plugin\u0026nbsp;[\u003ca href=\"#_ENREF_17\" title=\"Bandettini, 2012 #36\"\u003e17\u003c/a\u003e]\u0026nbsp;(Version1.4.2, http://apps.cytoscape.org/apps/MCODE). The threshold value was set as score＞5.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eGO and pathway enrichment analysis\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe DEGs in coexpression network and PPI network were subjected to GO-BP[\u003ca href=\"#_ENREF_18\" title=\"Ashburner, 2000 #4\"\u003e18\u003c/a\u003e]\u0026nbsp;and KEGG pathway\u0026nbsp;[\u003ca href=\"#_ENREF_19\" title=\"Kanehisa, 2000 #5\"\u003e19\u003c/a\u003e]\u0026nbsp;enrichment analysis with the extensively used tool of\u0026nbsp;clusterprofiler\u0026nbsp;[\u003ca href=\"#_ENREF_20\" title=\"Yu, 2012 #56\"\u003e20\u003c/a\u003e]\u0026nbsp;(Version3.18.1,http://www.bioconductor.org/packages/release/bioc/html/clusterProfiler.html). P value \u0026lt;0.05 was set as the cutoff value.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eIntegrated analysis of RNA-seq and ATAC-seq data\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe analysis of ATAC-seq data unravels transcription factor binding sites. To explore the regulatory role of TFs in gene expression,\u003cstrong\u003e\u0026nbsp;A\u003c/strong\u003eTAC-seq and RNA-seq data were integrated. The down-regulated genes obtained by RNA-seq data were overlapped with those around Peaks with depressed reads signals analyzed by \u003cstrong\u003eA\u003c/strong\u003eTAC-seq data, otherwise the up-regulated genes were overlapped. In the present study, the \u0026nbsp;differential Peak adjacent genes analyzed by ATAC-seq data were overlapped with nodes in co-expression network. The overlaps were subjected to functional enrichment analysis.\u0026nbsp;\u003c/p\u003e"},{"header":"Results","content":"\u003cp\u003e\u003cstrong\u003eIdentification of DEGs\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe genes that may be perturbed by Pde6b mutation and play a key role in the development of photoreceptor cells in the retina were identified by RNA-seqbn \u0026nbsp; n . \u0026nbsp;Total 1061 genes were found to be differentially expressed in Rd1 samples compared with C3H ones at P15, including 340 up-regulated genes and 721 down-regulated ones. The gene differential expression was illustrated by heatmap, which showed differential expression profile of DEGs between Rd1 and C3H samples (Figure 1A). The volcano plot represented differential expressed genes based on the defined p value and fold change cutoff values (Figure 1B).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCoexpression network of PARP1 with DEGs\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eBased on Pearson correlation analysis, 313 PARP1-DEG coexpressed pairs were obtained, including 102 pairs with positive correlation and 211 ones with negative expression. The coexpression network was illustrated in Figure 2, which was comprised with 102 up-regulated genes and 211 down-regulated genes. Function enrichment analysis was performed for up- and down-regulated genes, respectively. Results showed that the up-regulated genes were closely related with tube formation related biological process, PPAR signaling pathway and B cell receptor signaling pathway (Table 1). The down-regulated genes were significantly enriched in visual perception and sensory perception of light stimulus related biological process. The most significant pathway for down-regulated genes was phototransduction (mmu04744) (Table 1).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003ePPI network and the functional modules\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe PPI network was constructed with 202 nodes connected with 375 edges (Figure 3A). Three significant network modules were divided by MCODE with PPI score \u0026ge; 5, including Module1(score=8.222, 10 nodes and 37 edges), Module 2(score=6.333, 7 nodes and 19 edges) and Module 3 (7 nodes and 18 edges) (Figure 3B). After topological property analysis, the node degrees of modular genes were obtained (Table 2). Functional analysis showed that the genes in Module 1 were closely related with detection of light stimulus (GO:0009583) and visual perception (GO:0007601) related biological process and involved with Phototransduction pathway (mmu04744). The Module 2 were correlated with regulation of mitotic cell cycle (GO:0007346) and chromosome segregation (GO:0007059) related biological process and DNA replication (mmu03030) and cell cycle (mmu04110) related pathways. Additionally, the genes in module 3 were significantly enriched in behavioral fear response (GO:0001662), regulation of amino acid transport (GO:0051955), Neuroactive ligand-receptor interaction (mmu04080), PI3K-Akt signaling pathway (mmu04151) (Table 3).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eATAC-seq and integrated analysis\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eBased on ATAC-seq, we obtained 1394 peaks up-regulated signals and 891 differential peaks with down-regulated signal. The differential accessible peaks covered genes were overlapped with DEGs in co-expression network. The Venn diagram depicted that the overlapped up-regulated genes were Asxl3 and Nyap2 (Figure 4A, while the down-regulated overlaps were Tmem136 and Susd3 (Figure 4B).\u0026nbsp;\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003eRP is one of the leading cause of inherited blindness, which has affected 1/3000 to 1/5000 individuals globally. PDE6B mutation is responsible for 5\u0026ndash;8% RP patients, which affects rod-specific cyclic guanosine monophosphate (cGMP) level and results in rod degeneration [\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e, \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e]. Recent evidence showed that PARP1 was a promising therapeutic target for PDE6 mutation-induced RP. However, the regulatory mechanism of PARP1 in photoreceptor cell apoptosis underlying RP has not been clarified. In this study, we attempted to explore the potential mechanism related with photoreceptor cell apoptosis in PDE6 mutation-induced RP with ATAC-seq and RNA-seq.\u003c/p\u003e \u003cp\u003eOur results of RNA-seq showed that there were 1061 DEGs between Rd1 and C3H mice at P15. The expression profile of DEGs could distinguish the samples from Rd1 and C3H clearly, which revealed that DEGs identified were significant and could play key roles in the pathogenesis of RP. In order to explore the mechanism related with PARP1 in RP, PARP1-DEG coexpression network was constructed. Functional analysis showed that the down-regulated genes with coexpression correlation with PARP1 were closely related with visual perception, sensory perception of light stimulus related biological process and phototransduction pathway. The manifestations of RP were characterized by loss of depth perception, night blindness, and loss of central vision [\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e]. The variety of eye disorders in RP patients was related with rod-opsin protein disruption, in parallel with the defects of phototransduction cascade of central nervous [\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e]. Our findings were consistent with the reports mentioned above and revealed that the visual perception and phototransduction pathway were dysregulated in RP. Additionally, our data showed that PPAR signaling pathway was a significant pathway enriched by up-regulated genes. It is reported that PPAR plays a regulatory role in the transcription of genes involved in intracellular signaling pathways [\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e]. PPAR is involved in lipid metabolism and homeostasis [\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e]. Accumulating evidence indicates that abnormal lipid metabolism is associated with RP development [\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e, \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e]. In this study, the PPAR pathway was found to be enhanced by up-regulated genes. Thus, we speculated that PARP1 could be associated with abnormal lipid metabolism induced photoreceptor cell apoptosis in RP.\u003c/p\u003e \u003cp\u003eBesides, the PPI network indicated that the Module 1 genes were also enriched in phototransduction pathway, which was in agree with GO-BP analysis of DEGs in coexpression network. Gnat 1 was a down-regulated gene in retinal tissues of Rd1 mice with highest degree of 19 in Module 1 PPI network. As described in previous study, \u003cem\u003eGnat1\u003c/em\u003e encodes the α-subunit of transducin in rod cells. \u003cem\u003eGnat1\u003c/em\u003e variants were previously reported to be associated with congenital stationary night-blindness [\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e]. New evidence showed that \u003cem\u003eGnat1\u003c/em\u003e mutation was responsible for the slow and late onset of retinal degeneration [\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e]. Similarly, our GO-BP enrichment analysis showed that \u003cem\u003eGnat1\u003c/em\u003e was closely related with visual perception and light stimulus related biological process. Gnat1 showed negative expression correlation with PARP1 (r=-0.90 and p\u0026thinsp;=\u0026thinsp;0.013). All these above suggested that the down-regulation of Gnat1 was correlated with PARP1 overexpression, which caused the dysfunction of visual perception, light stimulus, and phototransduction, leading to the onset of RP.\u003c/p\u003e \u003cp\u003eFurthermore, the RNA-seq and ATAC-seq integrated analysis in our study revealed that Susd3 (sushi domain-containing protein 3) was a down-regulated overlap. Susd3 is cell surface protein and has been suggested as the potential biomarker for breast cancer [\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e]. Besides, Susd3 was found to be one of the genes encoding protein at the surface of young rod photoreceptors [\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e]. Susd2 was identified to be specially expressed for photoreceptors in the retina. However, the significant role of Susd3 in RP onset has not been clarified. ATAC-seq is a simple protocol for open chromatin analysis, which provides possibility of discovering TF binding sites. Susd3 was the gene covered by differential peaks by ATAC-seq, which indicated that Susd3 could exert TF function. Susd3 showed coexpression correlation with PARP1. Currently, little is known about whether Susd3 elicited transcriptional regulation on PARP1 and affected PDE6-mutation induced RP, which should be verified in the following studies.\u003c/p\u003e \u003cp\u003eIn summary, the down-regulated genes with coexpression correlation of PARP1 were significantly enriched in visual perception, sensory perception of light stimulus related biological process and phototransduction pathway. The PPAR signaling pathway was enhanced by up-regulated genes. Gnat1 could play a key role in visual perception, light stimulus, and phototransduction of RP by interacting with PARP1. Susd3 might play a regulatory role in PARP1 expression and affect RP formation. Our study may provide new insights in understanding the mechanism of PDE6-mutation induced RP.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eCampochiaro, P.A. and T.A. Mir, \u003cem\u003eThe mechanism of cone cell death in Retinitis Pigmentosa.\u003c/em\u003e Prog Retin Eye Res, 2018. \u003cstrong\u003e62\u003c/strong\u003e: p. 24-37.\u003c/li\u003e\n\u003cli\u003eZhang, Q., \u003cem\u003eRetinitis Pigmentosa: Progress and Perspective.\u003c/em\u003e Asia Pac J Ophthalmol (Phila), 2016. \u003cstrong\u003e5\u003c/strong\u003e(4): p. 265-71.\u003c/li\u003e\n\u003cli\u003eParmeggiani, F., \u003cem\u003eClinics, epidemiology and genetics of retinitis pigmentosa.\u003c/em\u003e Curr Genomics, 2011. \u003cstrong\u003e12\u003c/strong\u003e(4): p. 236-7.\u003c/li\u003e\n\u003cli\u003eSchon, C., et al., \u003cem\u003eGene therapy successfully delays degeneration in a mouse model of PDE6A-linked retinitis pigmentosa (RP 43).\u003c/em\u003e Hum Gene Ther, 2017.\u003c/li\u003e\n\u003cli\u003eKim, J.Y., et al., \u003cem\u003eEpiretinal Membrane Formation after Intravitreal Autologous Stem Cell Implantation in a Retinitis Pigmentosa Patient.\u003c/em\u003e Retin Cases Brief Rep, 2017. \u003cstrong\u003e11\u003c/strong\u003e(3): p. 227-231.\u003c/li\u003e\n\u003cli\u003eKaur, J., et al., \u003cem\u003eCalpain and PARP activation during photoreceptor cell death in P23H and S334ter rhodopsin mutant rats.\u003c/em\u003e PLoS One, 2011. \u003cstrong\u003e6\u003c/strong\u003e(7): p. e22181.\u003c/li\u003e\n\u003cli\u003eGopalakrishna, K.N., K. Boyd, and N.O. Artemyev, \u003cem\u003eMechanisms of mutant PDE6 proteins underlying retinal diseases.\u003c/em\u003e Cell Signal, 2017. \u003cstrong\u003e37\u003c/strong\u003e: p. 74-80.\u003c/li\u003e\n\u003cli\u003eDvir, L., et al., \u003cem\u003eAutosomal-recessive early-onset retinitis pigmentosa caused by a mutation in PDE6G, the gene encoding the gamma subunit of rod cGMP phosphodiesterase.\u003c/em\u003e Am J Hum Genet, 2010. \u003cstrong\u003e87\u003c/strong\u003e(2): p. 258-64.\u003c/li\u003e\n\u003cli\u003eAndrabi, S.A., et al., \u003cem\u003ePoly(ADP-ribose) (PAR) polymer is a death signal.\u003c/em\u003e Proc Natl Acad Sci U S A, 2006. \u003cstrong\u003e103\u003c/strong\u003e(48): p. 18308-13.\u003c/li\u003e\n\u003cli\u003eHottiger, M.O., \u003cem\u003eADP-ribosylation of histones by ARTD1: an additional module of the histone code?\u003c/em\u003e FEBS Lett, 2011. \u003cstrong\u003e585\u003c/strong\u003e(11): p. 1595-9.\u003c/li\u003e\n\u003cli\u003eThomas, C. and A.V. Tulin, \u003cem\u003ePoly-ADP-ribose polymerase: machinery for nuclear processes.\u003c/em\u003e Mol Aspects Med, 2013. \u003cstrong\u003e34\u003c/strong\u003e(6): p. 1124-37.\u003c/li\u003e\n\u003cli\u003eJiao, K., et al., \u003cem\u003eEfficacy of PARP inhibition in Pde6a mutant mouse models for retinitis pigmentosa depends on the quality and composition of individual human mutations.\u003c/em\u003e Cell Death Discov, 2016. \u003cstrong\u003e2\u003c/strong\u003e: p. 16040.\u003c/li\u003e\n\u003cli\u003eSothilingam, V., et al., \u003cem\u003eRetinitis pigmentosa: impact of different Pde6a point mutations on the disease phenotype.\u003c/em\u003e Hum Mol Genet, 2015. \u003cstrong\u003e24\u003c/strong\u003e(19): p. 5486-99.\u003c/li\u003e\n\u003cli\u003eKeeler, C.E., \u003cem\u003eThe Inheritance of a Retinal Abnormality in White Mice.\u003c/em\u003e Proc Natl Acad Sci U S A, 1924. \u003cstrong\u003e10\u003c/strong\u003e(7): p. 329-33.\u003c/li\u003e\n\u003cli\u003eSzklarczyk, D., et al., \u003cem\u003eSTRING v10: protein\u0026ndash;protein interaction networks, integrated over the tree of life.\u003c/em\u003e Nucleic acids research, 2014: p. gku1003.\u003c/li\u003e\n\u003cli\u003eTang, Y., et al., \u003cem\u003eCytoNCA: A cytoscape plugin for centrality analysis and evaluation of protein interaction networks.\u003c/em\u003e BioSystems, 2015. \u003cstrong\u003e127\u003c/strong\u003e(Complete): p. 67-72.\u003c/li\u003e\n\u003cli\u003eBandettini, W.P., et al., \u003cem\u003eMultiContrast Delayed Enhancement (MCODE) improves detection of subendocardial myocardial infarction by late gadolinium enhancement cardiovascular magnetic resonance: a clinical validation study.\u003c/em\u003e J Cardiovasc Magn Reson, 2012. \u003cstrong\u003e14\u003c/strong\u003e: p. 83.\u003c/li\u003e\n\u003cli\u003eAshburner, M., et al., \u003cem\u003eGene Ontology: tool for the unification of biology.\u003c/em\u003e Nat Genet, 2000. \u003cstrong\u003e25\u003c/strong\u003e(1): p. 25-29.\u003c/li\u003e\n\u003cli\u003eKanehisa, M. and S. Goto, \u003cem\u003eKEGG: Kyoto Encyclopedia of Genes and Genomes.\u003c/em\u003e Nucleic Acids Research, 2000. \u003cstrong\u003e28\u003c/strong\u003e(1): p. 27-30.\u003c/li\u003e\n\u003cli\u003eYu, G., et al., \u003cem\u003eclusterProfiler: an R package for comparing biological themes among gene clusters.\u003c/em\u003e Omics, 2012. \u003cstrong\u003e16\u003c/strong\u003e(5): p. 284-7.\u003c/li\u003e\n\u003cli\u003eKim, M.S. and K. Joo, \u003cem\u003eGenetic Mutation Profiles in Korean Patients with Inherited Retinal Diseases.\u003c/em\u003e 2019. \u003cstrong\u003e34\u003c/strong\u003e(21): p. e161.\u003c/li\u003e\n\u003cli\u003eYeo, J.H., et al., \u003cem\u003eDevelopment of a Pde6b Gene Knockout Rat Model for Studies of Degenerative Retinal Diseases.\u003c/em\u003e Invest Ophthalmol Vis Sci, 2019. \u003cstrong\u003e60\u003c/strong\u003e(5): p. 1519-1526.\u003c/li\u003e\n\u003cli\u003ePrem Senthil, M., J. Khadka, and K. Pesudovs, \u003cem\u003eSeeing through their eyes: lived experiences of people with retinitis pigmentosa.\u003c/em\u003e Eye (Lond), 2017. \u003cstrong\u003e31\u003c/strong\u003e(5): p. 741-748.\u003c/li\u003e\n\u003cli\u003eMendes, H.F., et al., \u003cem\u003eMechanisms of cell death in rhodopsin retinitis pigmentosa: implications for therapy.\u003c/em\u003e Trends Mol Med, 2005. \u003cstrong\u003e11\u003c/strong\u003e(4): p. 177-85.\u003c/li\u003e\n\u003cli\u003eYanagi, Y., \u003cem\u003eRole of Peoxisome Proliferator Activator Receptor gamma on Blood Retinal Barrier Breakdown.\u003c/em\u003e PPAR Res, 2008. \u003cstrong\u003e2008\u003c/strong\u003e: p. 679237.\u003c/li\u003e\n\u003cli\u003eDonato, L., et al., \u003cem\u003eTranscriptome Analyses of lncRNAs in A2E-Stressed Retinal Epithelial Cells Unveil Advanced Links between Metabolic Impairments Related to Oxidative Stress and Retinitis Pigmentosa.\u003c/em\u003e Antioxidants (Basel), 2020. \u003cstrong\u003e9\u003c/strong\u003e(4).\u003c/li\u003e\n\u003cli\u003eHolman, R.T., et al., \u003cem\u003eAbnormal plasma lipids of patients with Retinitis pigmentosa.\u003c/em\u003e Lipids, 1994. \u003cstrong\u003e29\u003c/strong\u003e(1): p. 61-5.\u003c/li\u003e\n\u003cli\u003eSzabo, V., et al., \u003cem\u003ep.Gln200Glu, a putative constitutively active mutant of rod alpha-transducin (GNAT1) in autosomal dominant congenital stationary night blindness.\u003c/em\u003e Hum Mutat, 2007. \u003cstrong\u003e28\u003c/strong\u003e(7): p. 741-2.\u003c/li\u003e\n\u003cli\u003eCarrigan, M. and E. Duignan, \u003cem\u003eA novel homozygous truncating GNAT1 mutation implicated in retinal degeneration.\u003c/em\u003e 2016. \u003cstrong\u003e100\u003c/strong\u003e(4): p. 495-500.\u003c/li\u003e\n\u003cli\u003eYu, Z., et al., \u003cem\u003eSushi Domain-Containing Protein 3: A Potential Target for Breast Cancer.\u003c/em\u003e Cell Biochemistry and Biophysics, 2015. \u003cstrong\u003e72\u003c/strong\u003e(2): p. 321-324.\u003c/li\u003e\n\u003cli\u003ePostel, K., et al., \u003cem\u003eAnalysis of cell surface markers specific for transplantable rod photoreceptors.\u003c/em\u003e Mol Vis, 2013. \u003cstrong\u003e19\u003c/strong\u003e: p. 2058-67.\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":true,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"international-ophthalmology","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"inte","sideBox":"Learn more about [International Ophthalmology](https://www.springer.com/journal/10792)","snPcode":"10792","submissionUrl":"https://submission.nature.com/new-submission/10792/3","title":"International Ophthalmology","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"ATAC-seq, RNA-seq, Retinitis pigmentosa, PARP1, retina tissues ","lastPublishedDoi":"10.21203/rs.3.rs-735055/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-735055/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e\u003cstrong\u003ePurpose\u003c/strong\u003e: Inhibition of PARP1 could relieve PDE6 mutation-induced Retinitis pigmentosa (RP). However, the mechanism related with PARP1 overexpression in the RP has not been clarified. We attempted to explore the potential regulatory mechanism related with PARP1 underlying RP. \u003c/p\u003e\u003cp\u003e\u003cstrong\u003eMethods\u003c/strong\u003e: ATAC-seq and RNA-seq were preformed for retina tissues of C3H and rd1 mice. The differential expressed genes (DEGs) were identified, followed by PARP1-DEG coexpression network, and PPI network construction. GO-BP and pathway enrichment of DEGs of interest were performed by clusterprofiler software. The overlapped genes that might play regulatory role in PARP1 expression were mined by integrated analysis of RNA-seq and ATAC-seq data. \u003c/p\u003e\u003cp\u003e\u003cstrong\u003eResults\u003c/strong\u003e: Total 1061 DEGs were identified between C3H and rd1 group. Co-expression network was constructed with 313 PARP1-gene coexpression pairs. The down-regulated DEGs were closely related with visual perception, light stimulus related biological process, while the up-regulated DEGs were significantly enriched in phototransduction and PPAR signaling pathway. PPI network was constructed with 202 nodes and 375 edges, which was clustered into 3 modules. Module 1 genes were closely related with detection of light stimulus, visual perception related biological process and phototransduction pathway (involved with Gnat1/Guca1b/Gnat2/Sag/Pde6g). By integrated analysis of the RNA-seq and ATAC-seq, the overlapped up-regulated genes were Asxl3 and Nyap2, while the down-regulated genes were Tmem136 and Susd3. \u003c/p\u003e\u003cp\u003e\u003cstrong\u003eConclusion\u003c/strong\u003e: Gnat1 may play a key role in RP development by interacting with PARP1. Susd3 may play a regulatory role in PARP1 expression and affect RP formation.\u003c/p\u003e","manuscriptTitle":"Integrative RNA-seq and ATAC-seq Analysis of Retinitis Pigmentosa Caused by PDE6 Mutation","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2021-08-30 20:36:20","doi":"10.21203/rs.3.rs-735055/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"editorInvitedReview","content":"","date":"2021-08-27T18:31:54+00:00","index":0,"fulltext":""},{"type":"reviewersInvited","content":"","date":"2021-08-27T16:49:24+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2021-07-20T09:41:40+00:00","index":"","fulltext":""},{"type":"submitted","content":"International Ophthalmology","date":"2021-07-20T01:48:36+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"international-ophthalmology","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"inte","sideBox":"Learn more about [International Ophthalmology](https://www.springer.com/journal/10792)","snPcode":"10792","submissionUrl":"https://submission.nature.com/new-submission/10792/3","title":"International Ophthalmology","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"c95a1293-14bd-4a30-a4eb-46036eaeaaeb","owner":[],"postedDate":"August 30th, 2021","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"under-review","subjectAreas":[{"id":6817827,"name":"Ophthalmology"}],"tags":[],"updatedAt":"2022-01-15T17:27:34+00:00","versionOfRecord":[],"versionCreatedAt":"2021-08-30 20:36:20","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-735055","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-735055","identity":"rs-735055","version":["v1"]},"buildId":"7rjqhiLT3MXkJMwkYKINL","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}