Differential gene expression and protein-protein interaction network profiling of sulfur mustard-exposed rabbit corneas employing RNA-seq data and bioinformatics tools

Abstract

Sulfur mustard (SM) ocular exposure severely damages the cornea and causes vision impairment. At present, no specific therapy exists to mitigate SM-induced corneal injury and vision loss. This study performed transcriptome profiling of naïve, SM-damaged, and SM-undamaged rabbit corneas using RNA-seq analysis and bioinformatic tools to gain a better mechanistic understanding and develop SM-specific medical countermeasures. The mRNA profiles of rabbit corneas 4 weeks post SM vapor exposure were generated using Illumina-NextSeq deep sequencing (Gene Expression Omnibus accession # GSE127708). The RNA sequences of naïve (n = 4), SM-damaged (n = 5), and SM-undamaged (n = 5) corneas were subjected to differential expression (DE) analysis after quality control profiling with FastQC. DE analysis was performed using HISAT2, StringTie, and DESeq2. The log2(FC)±2 and adjusted p˂0.05 were chosen to identify the most relevant genes. A total of 5930 differentially expressed genes (DEGs) (upregulated: 3196, downregulated: 2734) were found in SM-damaged corneas compared to naïve corneas, whereas SM-undamaged corneas showed 1884 DEGs (upregulated: 1029, downregulated: 855) compared to naïve corneas. DE profiling of SM-damaged corneas to SM-undamaged corneas revealed 985 genes (upregulated: 308, downregulated: 677). The DE profiles were subsequently subjected to signaling pathway enrichment, and protein‒protein interactions (PPIs) were analyzed. Pathway enrichment was performed for the genes associated with cellular apoptosis, death, adhesion, migration, differentiation, proliferation, extracellular matrix, and tumor necrosis factor production. To identify novel targets, we narrowed the pathway analysis to upregulated and downregulated genes associated with cell proliferation and differentiation, and PPI networks were developed. Furthermore, protein targets associated with cell differentiation and proliferation that may play vital roles in corneal fibrosis and wound healing post SM injury were identified.