{'fraction_non_alphanumeric': 0.04257179836988841, 'fraction_numerical': 0.030013642564802184, 'mean_word_length': 4.754428341384863, 'pattern_counts': {'":': 0, '<': 4, '<?xml version=': 0, '>': 0, 'https://': 2, 'lorem ipsum': 0, 'www.': 13, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Novel RvD6 stereoisomer induces corneal nerve regeneration and wound healing post-injury by modulating trigeminal transcriptomic signature

The high-density corneal innervation plays a pivotal role in sustaining the integrity of the ocular surface. We have previously demonstrated that pigment epithelium-derived factor (PEDF) plus docosahexaenoic acid (DHA) promotes corneal nerve regeneration; here, we report the mechanism involved and the discovery of a stereospecific Resolvin D6-isomer (RvD6si) that drives the process. RvD6si promotes corneal wound healing and functional recovery by restoring corneal innervation after injury. RvD6si applied to the eye surface elicits a specific transcriptome signature in the trigeminal ganglion (TG) that includes Rictor, the rapamycin-insensitive complex-2 of mTOR (mTORC2), and genes involved in axon growth, whereas genes related to neuropathic pain are decreased. As a result, attenuation of ocular neuropathic pain and dry eye will take place. Thus, RvD6si opens up new therapeutic avenues for pathologies that affect corneal innervation.

Dry eye perturbs vision mainly during aging. It also occurs in rheumatoid arthritis, diabetes, thyroid gland pathologies, environmental conditions (e.g., exposure to smoke or pollutants), long-term use of contact lenses, and after refractive surgery. This pathology is triggered by a shortage in tears that lubricate, arrest infections, and nourish and sustain a clear eye surface. Corneal innervation is required to maintain the integrity of the ocular surface 1 , and nerve damage decreases tear production, blinking reflex, and perturbs epithelial wound healing, resulting in loss of transparency and vision [2][3][4][5] . For this reason, there is a strong relationship between dry eye and corneal nerve damage.
Axons from sensory nerves from the ophthalmic branch of the trigeminal ganglion (TG) neurons penetrate the corneal stroma surrounding the limbal area and branch out as the subepithelial plexus before reaching the corneal epithelium, finalizing as free nerve endings [6][7][8] .
After nerve damage occurs from refractive surgeries (e.g., situ keratomileusis/LASIK or photorefractive keratectomy/PRK), it takes between 3-15 years to recover corneal nerve integrity [9][10][11] . As a consequence, corneal sensitivity decreases, and dry-eye disease can develop, causing neuropathic pain, corneal ulcers, and in severe cases, the necessity for corneal transplants [12][13][14] . In addition, dry eye is linked to cold receptor function, mainly the transient receptor potential melastatin 8 (TRPM8) channels 15 that control the corneal surface rate of cooling and maintain normal tear secretion [16][17][18] . In fact, a decrease in TRPM8 terminals takes place, even long after experimental corneal surgery, suggesting that these changes contribute to post-surgery neuropathic pain 19 .
Topical treatment of the neurotrophin pigment epithelium-derived factor (PEDF) plus the ω−3 fatty acid family member docosahexaenoic acid (DHA) enhances nerve regeneration and stimulates nerve regrowth in rabbit and mouse corneas after experimental surgery, as well as in pathologies like diabetes and herpes virus simplex (HSV1) infection [20][21][22][23][24] . Moreover, PEDF activates the Ca 2+ -independent phospholipase A2 (iPLA2ζ) activity of the PEDF receptor (PEDF-R) and releases DHA from membrane phospholipids that can be converted into bioactive docosanoids 25 , including neuroprotectin D1 (NPD1) that induces corneal nerve regeneration in a rabbit RvD6si is derived from DHA. To investigate whether the new RvD6si originated from the added DHA, an ex vivo corneal organ culture model was employed (16 corneas/sample). The injured corneas were cultured for 4 h in the presence of DHA or deuterium-labeled DHA (DHA-d5) plus PEDF, and the lipids from the media were extracted and analyzed. Since 5 atoms of deuterium (D) are attached to the end of the DHA backbone (at the 21 st and 22 nd C), the total mass of RvD6si-d5 was shifted to 365 Da (the [M-H] m/z is 364 in MS results) while some of its product ions were not changed after fragmentation ( Fig. 2A). The MRM detection method was designed to capture the DHA-d5 total structure. The RvD6si-d5 was detected in the media with a similar RT to the RvD6si produced by PEDF + DHA (Fig. 2B). The full fragmentation of RvD6si-d5 confirmed the structure as well (Fig. 2C). In addition, the origin of the RvD6si was validated at three different concentrations of added DHA (Fig. 2D) with an enhanced synthesis as a function of increased DHA concentration. When analyzing possible arachidonic acid (AA)-and DHA-hydroxy derivatives (HDHA), the results showed a proportional increase of DHA products such as 14-and 17-HDHAs while the amount of AA and its hydroxyl derivatives 12-and 15-HETEs were not changed. These data strongly support that the new RvD6si originates from exogenous DHA.

RvD6si enhances corneal wound healing and recovery of corneal sensitivity after injury.
Although the 2D structure of the new RvD6si matched RvD6, the different RT could make them distinct in their biological activities. To obtain enough RvD6si for testing, 60 mice were injured and treated with PEDF + DHA every 30 min for 4 h, and the tears collected. The next day, the mice were euthanized, and the corneas isolated and incubated in media with PEDF + DHA for 4 h. The lipids extracted from tears and corneal media were combined and run in UPLC employing a C18 column, and fractions were collected every 30 sec from 6 to 12 min. All fractions were subject to lipidomic analysis to detect the availability of the new RvD6si. Fractions 6 to 8 with clear detectable amounts of RvD6si were pooled ( Supplementary Fig. S1). The purity of our targeted lipid mediator was determined by lipidomic analysis before being tested in vivo. The isolated RvD6si showed very low contamination of other DHA derivatives as well as AA, eicosapentaenoic acid (EPA), and their derivatives ( Supplementary  Fig. S1).
Previous studies have shown that PEDF + DHA promotes corneal wound healing in rabbit 20,21 , and in normal and diabetic mice 24,25 after experimental surgery. We tested the ability of RvD6s (either RvD6 or RvD6si) in stimulating corneal wound healing. The right mouse eyes were injured, and the animals were divided into four groups: vehicle, PEDF + DHA, RvD6, and RvD6si (Fig. 3A). Twenty h after injury, all drug-treated mice had faster corneal wound healing than vehicle; however, the greatest increase was found in the animals treated with the RvD6si (Fig. 3B,C).
Corneal sensitivity was evaluated at days 3, 6, 9, and 12 after corneal injury and treatment (Fig. 3D). A new methodology to measure the sensation in mouse cornea was introduced using the Belmonte non-contact aesthesiometer. Figure 3E shows a Gaussian-curve of distribution from basal corneal sensation recorded-values (n = 40 corneas) at a flow rate of 100.45 to 110.05 ml of air/minute (α = 0.05). It is important to note that this range of normal corneal sensitivity is critical to evaluate corneal sensation since the Belmonte non-contact aesthesiometer working flow rate is from 20 to 200 ml/min. The range of normal corneal sensitivity from 100.45 to 110.05 ml in the mouse was regarded as successful recovery after injury and treatment, and it was used to normalize the measurements. There was a faster recovery of corneal sensation in the animals treated with the RvD6si at 3 and 6 days after injury (Fig. 3F). By 9 days, the three treatments increased the sensitivity compared to vehicle, and at 12 days, there was no significant difference in any of the studied groups.
RvD6si activates corneal nerve regeneration. PEDF + DHA stimulates corneal nerve regeneration in injury animal models [20][21][22][23][24][25] . It was important to confirm the biological activity of RvD6si as a lipid mediator  www.nature.com/scientificreports www.nature.com/scientificreports/ underlying the action of PEDF + DHA. To test this hypothesis, mice were injured and treated (as described in Fig. 3D). Isolated corneas were stained with PGP 9.5, a pan-neuronal marker, and with SP neuropeptide antibodies. The density of non-injured corneal nerves positive to PGP 9.5 and SP, respectively, was used to normalize the values (Fig. 4A). Substance P is a main neuropeptide in mammalian corneas [26][27][28] . Moreover, a previous study from our group has demonstrated that there is a correlation between corneal sensitivity and SP-positive nerves 29 . RvD6si-treated mice have significantly higher sensitivity at day 3, 6, and 9, while PEDF + DHA-and RvD6treated groups showed higher cornea sensation only at day 9. At day 12, there was no difference between the tested compounds and vehicle. The statistical p-value is derived from one-way ANOVA, followed by Tukey's honest significant difference (HSD) multiple pairwise comparisons. (2020) 10:4582 | https://doi.org/10.1038/s41598-020-61390-8 www.nature.com/scientificreports www.nature.com/scientificreports/ At 12 days after injury and treatment, total corneal nerve density was 45.9 ± 6.8% of the normal cornea in the vehicle-treated group and significantly higher in the RvD6si treated corneas 62.6 ± 4.2% (p < 0.05) (Fig. 4B). PEDF + DHA and RvD6 treatment also increased nerve density to 59.9 ± 63% and 59.7 ± 11.2%, respectively. There were no significant differences between RvD6si, RvD6, and PEDF + DHA. Similarly, the density of SP-positive nerves at 12 days after injury was higher with RvD6si, RvD6 and PEDF + DHA treatments compared to the vehicle-treated group (Fig. 4C). This result confirmed a faster recovery of corneal sensitivity in treated corneas (Fig. 3F) and strengthened the biological function of the RvD6si as the main mediator in the mechanism of PEDF + DHA to enhance corneal nerve regeneration.
RvD6si modulates a transcriptome signature in the trigeminal ganglion. Because corneal sensory nerves originate in TG neurons, we explored if corneal injury could be sensed in the TG and t, in turn, would elicit a gene expression response. Thus, TG were harvested 12 days after injury and treatment with RvD6si or RvD6 or vehicle treatment used as control (Fig. 3D), and then RNA-seq analysis was performed. Quality controls showed that mapped reads range from 84.63 to 93.00%, with about 20,000 expressed genes / sample. Principal component analysis (PCA) showed good separation of vehicle-treated from the RvD6si-or RvD6-treated groups (Fig. 5A). The two RvD6s shared 58 upregulated genes and 36 downregulated genes compared to controls (Fig. 5B). To classify upregulated genes of RvD6si_vs_vehicle and RvD6_vs_vehicle, gene enrichment analysis was used to demonstrate that the RvD6si showed a difference in cellular comparted locations ( Supplementary Fig. S2) and that activate axonal growth cone genes (gene ontology number 0044295). The box plots depict two activated genes by RvD6si in this class: C9orf72 and Gpm6a (Fig. 5C). We also detected specific genes related to neuropeptides and ion channel receptors in the cornea that are stimulated by the addition of PEDF + DHA 19,21,24 (Fig. 5D). The RNA-seq established that RvD6 or RvD6si reduced gene expression of two major neuropeptides, tachykinin precursor 1 (Tac1) that encodes Substance P (SP) and calcitonin-related polypeptide beta (Calcb). It is important to note that these neuropeptides, especially Calcb, are major pain-induced mediators in migraine and other primary headaches 30 . In contrast, the RvD6si selectively enhanced the expression of transient receptor potential melastatin www.nature.com/scientificreports www.nature.com/scientificreports/ 8 (Trpm8) channel, and neuropilin 1 (Nrp1), the co-receptor for several factors including class III/IV semaphorins, certain isoforms of vascular endothelial growth factor, and transforming growth factor beta 31 .
Further analysis revealed a strong induction by RvD6si of the transcriptional factor Rictor (Fig. 5E) that is a part of the rapamycin-insensitive mammalian target complex-2 (mTORC2) (Fig. 5E,G). There were 39 genes www.nature.com/scientificreports www.nature.com/scientificreports/ modulated by RICTOR modified by RvD6si. Among those, 37 (95%) genes matched the IPA knowledge collected from published data, while only two genes, Egr1 and Psme3, did not fit with the prediction (yellow arrows) (Fig. 5F). It is important to note that all genes subjected to IPA analysis are significantly different (FDR < 0.05 in the DESeq2 analysis) in comparison to the vehicle-treated group. For this reason, 95% of downstream genes matched to IPA knowledge; the Rictor signaling is clearly stimulated in TG by RvD6si (Fig. 5G).

Discussion
Previous studies from our laboratory have demonstrated the potential use of PEDF + DHA for corneal wound healing and nerve regeneration in post-surgical models of rabbits and mice [20][21][22][23][24][25] . This included the observation that activation of the iPLA2ζ activity of the PEDF-R releases DHA from phospholipids, suggesting that docosanoids could be synthesized in the cornea 25 . Here, we report the finding identification and characterization of its bioactivity of a novel Resolvin D6si in tears that is derived from DHA upon activation of PEDF on its receptor. The full MS/MS fragmentation of the RvD6si matches six characteristic ions with the RvD6 as well as the UV diode array profile (Fig. 1). The biological activity revealed that it enhances corneal wound healing and sensitivity recovery, more potently than PEDF + DHA after corneal PRK-mimic surgery (Fig. 3). These results indicate that RvD6si is the main lipid mediator that contributes to the signaling mechanism of the action of PEDF + DHA. Moreover, the RvD6s and PEDF + DHA treatments show similar enhancement in corneal innervation at 12 days after injury and treatment (Fig. 5).
Resolvin D6 was described in human polymorphonuclear neutrophils 32 , skin 33 , brain 34 , cerebrospinal fluid 35 , and plasma 36 . However, this is the first report demonstrating a biological function of RvD6 and of a novel stereoisomer. The formation of potent bioactive mediators from DHA was proposed when mono-, di-, and tri-hydroxy DHA-derivatives were detected as enzyme-mediated products of oxygenated metabolites of DHA in the retina 37 . Unlike the retina, where photoreceptor membranes have high DHA content esterified at the sn-2 position of phospholipids 38 , the cornea contains more AA at that position 25,39 . For this reason, the addition of exogenous DHA is required to synthesize docosanoids rather than eicosanoids. Further, the RvD6si was not detected when corneas were treated with DHA or PEDF alone (data not shown), supporting the premise that novel RvD6si is only detected when corneas are treated with PEDF + DHA. This observation is in agreement with a previous study showing that neither RvD6 nor its stereoisomers were detected in human tear samples 40 . Since the RvD6si was found primarily in the tears or media of corneas in organ culture, this suggests that the RvD6si needs to be secreted into the extracellular compartment to become functional. The biological activity may be elicited through a receptor and, in turn, modulates cell signaling and transcription factors, upregulating, as a consequence, neurotrophic genes in the cornea 25 . RvD6si may act in autocrine fashion and/or may diffuse through tears and act as a paracrine signal on other ocular surface cells.
Most of the corneal nerves originate from neurons localized in the TG 6 . Therefore, using unbiased RNA sequencing, we have deciphered here that RVD6 and RvD6si shared a small number of upregulated genes in the TG, implicating that the signaling mechanism of their biological activities have differences. The RNA-seq data reveal a strong activation by the RvD6si of two genes, C9orf72 and Gpm6a, that stimulate neurogenesis and growth cone formation 41,42 . We also found genes related to pain since corneal neuropathic pain can occur after nerve damage 43 . The expression of two genes involved in pain was decreased in corneas treated with the RvD6si: Tac1 that encodes SP, which is one of the most abundant neuropeptides expressed in corneal nerves [26][27][28] . SP exerts proinflammatory effects, and preclinical studies linked their action to chronic pain 44 . Besides, Tac1 is upregulated in the TG after corneal injury supporting the finding that corneal nerves "sense" and transmit the inflammatory signal to the TG 45 . The other is Calcb, which encodes Calcitonin gene-related peptide (CGRP). CGRP is also abundant in corneal nerves 21 and plays an essential role in neurogenic inflammation and pain 30 . Another important gene in this category is Trmp8. TRPM8 channels regulate the wetting of the ocular surface and have an analgesic effect on chronic pain 17,[46][47][48][49] . Previous studies in mice where the nerves had been damaged showed that TRPM8-positive nerve fibers only reach 50% of their normal density by 3 months after the injury, suggesting that the decrease in TRPM8 nerve terminals can contribute to dry eye-like pain 19 . The increased expression of Trpm8 after injury and treatment with RvD6si suggests that the new lipid could protect corneas from pain. In addition, the selective increase of Nrp1 is also interesting, since it is the co-receptor of SEMA3A that has been shown to attenuate mechanical allodynia in a rat model of sciatic nerve injury 50 .
Our results disclose that the RvD6si potently and selectively induces Rictor gene expression in the TG. As a regulator of PI3K/Akt pathways, RICTOR is a key component of mTORC2 and is clearly involved in cell proliferation and repair. In agreement with this concept, the deletion of Rictor or mTORC2 inhibited the sensory-axonal regeneration in mice after dorsal root ganglion injury 51 .
In conclusion, our data demonstrate that a new RvD6si produced by the injured cornea after PEDF + DHA treatment is necessary for corneal wound healing and nerve regeneration. This novel lipid mediator activates signaling that communicates from the cornea to TG neurons, and as a response, modulates specific gene signatures that enhance axon growth, decrease neuropathic pain, and foster containment of dry eye. Our findings open potential therapeutic avenues using RvD6sifor impaired-corneal nerve diseases, including dry eye, corneal neurotrophic ulcers, neurotrophic keratitis, and neuropathic pain.

Methods
Animals. Ten-week-old male CD1 mice were purchased from Charles River (Wilmington, MA, USA) and www.nature.com/scientificreports www.nature.com/scientificreports/ Corneal injury and treatment. Mice were anesthetized with a mix of ketamine (200 mg/kg) and xylazine (10 mg/kg) that was injected intraperitoneally, and one drop of proparacaine hydrochloride solution (0.5%) was applied to the right eye subjected to injury. As previously described 19,29 , the center of the cornea was demarcated with a 2 mm trephine, and the epithelium and the anterior stroma were gently removed under a surgical microscope using a corneal rust ring remover (Algerbrush II; Alger Equipment Co., Lago Vista, TX, USA). One drop of 0.3% of tobramycin ophthalmic solution (Henry Schein, Melville, NY, USA) was applied to the eye to prevent postoperative infection. The same investigator (J.H.) performed all surgeries. Afterward, 10 µl of PEDF (50 ng/ml) plus DHA (50 nM) or DHA-derived lipid mediators were applied topically, as explained in each experimental design.
Lipidomic analysis. Five microliters of sterile PBS was instilled in the inferior cul-de-sac of the mouse eye, and 30 s later, tears were collected in 1 mL of ice-cold MeOH containing 1 g/L Butylated hydroxytoluene followed by the addition of 2 ml of CHCl 3 and 5 μl of an internal standard mixture of deuterium-labeled lipids AA-d8 (5 ng/μl), PGD 2 -d4 (1 ng/μl), EPA-d5 (1 ng/μl), 15-HETE-d8 (1 ng/μl), and LTB 4 -d4 (1 ng/μl). The samples were sonicated in a water bath for 30 min and stored at −80 °C overnight. The next day, the samples were centrifuged, the supernatant was collected, and the pellet was washed with 1 ml of CHCl 3 /MeOH (2:1) and centrifuged, and then the supernatants were combined. Water at pH 3.5 was added to the supernatant at the ratio 1:5, vortexed, and centrifuged, and the pH of the upper phase was adjusted to 3.5-4.0 with 1 N HCl. The lower phase was collected, dried under N 2 , and then resuspended in 1 ml of MeOH and stored at −80 °C.
For corneal organ culture experiments, 2 mL of media was collected and centrifugated at 14,000 rpm for 15 min at 4 °C to remove cellular debris. Lipids were extracted by the Blight and Dyer method 52 . Briefly, 3.75 ml of a mixture of CHCl 3 : MeOH (1:2) was added to 1 ml of sample and 5 μl of the deuterium-labeled internal standard mixture of lipids. The samples were vortexed and stored at −80 °C overnight. Next, to make two phases, 2.5 ml of CHCl 3 was added and vortexed, and then 2.5 mL of water (pH 3.5) was added, vortexed, and the pH of the upper phase adjusted to 3.5-4.0 with 1 N HCl. The lower phase was dried down under N 2 , resuspended in 1 ml of MeOH, and stored at −80 °C.
LC-MS/MS analysis was performed in a Xevo TQ equipped with Acquity I class ultra-performance liquid chromatography (UPLC) with a flow-through needle (Waters Corporation, Milford, MA). As previously described 25,53 , samples were dried under N 2 , resuspended in 20 μl of MeOH/H 2 O (2:1), and injected into a CORTECS C18 2.7 μm 4.6 × 100 mm column (Water, MA). The column temperature was set at 45 °C with a flow of 0.6 ml/min. The initial mobile phase consisted of 45% solvent A (H 2 O + 0.01% acetic acid) and 55% solvent B (MeOH + 0.01% acetic acid) and then a gradient to 15% solvent A for the first 10 min followed by a gradient to 2% solvent A for 18 min, 2% solvent A run isocratically until 25 min, and then a gradient back to 45% solvent A for re-equilibration until 30 min. Lipid standards (Cayman, Ann Arbor, MI) were used for tuning and optimization, as well as to create calibration curves for each compound. RvD6 [4S,17S-dihydroxy-5E,7Z,10Z,13Z,15E,19Z-doc osahexaenoic acid] standard was a generous gift of Dr. Charles Serhan (Harvard University, Boston, MA, USA).

Production of Resolvin D6si from mouse tears and cornea.
Mouse corneas (n = 60) were injured and treated topically with PEDF + DHA for 4 h. Tears were collected in MeOH and stored at −80 °C. After 24 h, mice were euthanized, and injured corneas were excised and cultured with PEDF + DHA in DMEM/F12 media for 4 h. The medium was collected, and lipids were extracted as described above. Lipids from pooled tears and cornea-cultured media were subjected to UPLC separation using a C18 column (Water, MA). Twelve fractions (30 sec/fraction) between 6-12 min after injection were collected. The procedure was repeated at least 8 times with 25 µl of sample/run until all the sample was fractionated. Each fraction was dried under N 2 and resuspended in 1 mL of MeOH. The presence of RvD6si in 10 µl of each fraction was confirm using the described LC-MS/ MS system. The fractions with high purity and concentration of RvD6si were pooled and stored at −80 °C until needed for the in vivo experiments.
Corneal wound healing. As previously described 29 , mice were euthanized 20 h after injury and treatment, and corneas were stained with 0.5% methylene blue for 20sec and then washed with PBS for 2 min. Photographs were taken with a dissecting microscope (SMZ 1500; Nikon, Tokyo, Japan) through an attached digital camera (DXM 1200; Nikon). The images corresponding to the wounded area were quantified using Photoshop CC 2014 software (Adobe, San Jose, CA, USA).
Corneal sensitivity measurement. The non-contact corneal aesthesiometer has been described as a more reliable method than the standard Cochet-Bonnet aesthesiometer to determine the corneal sensation threshold 54 . Therefore, for corneal sensation measurement, the Belmonte non-contact corneal aesthesiometer 55 was used with some modification. Briefly, one researcher held the mouse and kept the air output needle at a distance of 3 mm from the cornea. Another researcher controlled the air flow rate. The measurements started at an air flow rate of 80 ml per minute and then increased gradually by ten units until the mouse started blinking. When the mouse blinked, the air flow rate was recorded as the final corneal sensitivity index.
Corneal nerve analysis. Twelve days after injury and treatment, mice were euthanized, and the eyes enucleated and fixed with Zamboni's fixative (American Master Tech Scientific, Lodi, CA, USA) for 45 min at room temperature. As previously described 19,24,29 , the corneas were then excised and fixed for an additional 15 min, followed by 3 washes with PBS. To block nonspecific binding, corneas were incubated with 10% normal goat serum plus 0.5% Triton X-100 in PBS for 1 h at room temperature. Afterward, corneas were incubated with the primary antibodies, rabbit monoclonal anti-PGP9.5 (1:500), (ab108986; Abcam, Cambridge, MA, USA), and rat monoclonal anti-substance-P (SP; 1:100) (sc-21715; Santa Cruz Biotechnology, Dallas, TX, USA) for 24 h at room temperature with constant shaking. After being washed with PBS, the corneas were incubated with the corresponding secondary antibodies goat anti-rabbit Alexa-Fluor 488 (1:1000) and goat anti-rat Alexa-Fluor 488 www.nature.com/scientificreports www.nature.com/scientificreports/ (1:1000) (Thermo Fisher Scientific, Waltham, MA, USA) for 24 h at 4 °C. Four radial cuts were performed on each cornea that was flatly mounted on a slide with the endothelium side up and examined with a fluorescent microscope (Deconvolution microscope DP80; Olympus, Tokyo, Japan). The images were merged together to build the entire view of the corneal nerve network. The corneal nerve density was measured using Photoshop CC 2014 (Adobe) as previously described 26,29 . Trigeminal ganglion RNA sequencing. TG corresponding to the injury eye side (n = 5) were harvested and kept in RNAlater solution (Thermo Fisher Scientific) until homogenized on ice using a Dounce homogenizer. Total mRNA was extracted using an RNeasy mini kit (Qiagen, Germantown, MD, USA) as described by the manufacturer. The purity and concentration of RNA were determined with a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific), and the samples were stored at −80 °C until used. RNA sequencing was performed using the adapted Smart-seq2 protocol 56 . Briefly, one ng of total RNA was reverse transcribed with the Oligo-dT 30 VN and template-switching oligo (TSO) primers. The total cDNAs were amplified using ISPCR primer, and the library was made using the Nextera XT DNA library preparation kit (Illumina, San Diego, CA, USA). The libraries were pooled using the same molarity and sequenced using the NextSeq 500/550 High Output Kit v2 (75 cycles, Illumina). After demultiplexing, RNA-seq data were aligned to the GENCODE GRCm38 mouse primary genome assembly (Release M22, gencodegenes.org/mouse/) using the RSubread package v1.34.6 for R v3.6.1 57 . The outputted BAM files for sequencing data alignment were counted using featureCounts function (Subread v1.6.5 in Ubuntu LTS 16.4 operating system) 58 . Next, the raw count data were subjected to differential gene expression analysis using DESeq2 package for R 59 . The adjusted p-values were regarded as the false discover rate (FDR). Significantly changed genes (FDR < 0.05) between RvD6si_vs_vehicle and RvD6_vs_vehicle were subjected to the enrichment analysis using Enrichr 60 and pathway analysis using the IPA (QIAGEN Inc., https:// www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis).
Statistical analysis. Data are expressed as mean ± SD of ≥3 independent experiments. The data were analyzed by 1-way ANOVA followed by Tukey honest significant difference post hoc test at 95% confidence level to compare the different groups and considered significant when p < 0.05. All statistical analyses were performed using the Stata 14 (StataCorp, College Station, TX, USA). Graphs were made using Prism 7 software (GraphPad Software, La Jolla, CA, USA) and Bio Vinci (BioTuring, La Jolla, CA, USA). For the sequencing data, since the DE-Seq2 analysis does not provide the multi-samples comparison, the normalized counts from DE-Seq2 were used as the input of ANOVA test.

Data availability
Completed RNA-Seq data that support the findings of this study have been deposited in Gene Expression Omnibus with the accession code GSE138685.