Arachidonic Acid Metabolism Zebrafish Mouse

Arachidonic Acid Metabolism Zebrafish Eicosanoids A. Nodes Edges 5,6-EET ephx2 5,6-DHET node 8 hpl 16hpl compound

cyp2p10 8,9-EET ephx2 8,9-DHET node complex

14,15-EET ephx2 C14775 enrichment score -1 0 1 Linoleic acid metaolism 20-OH-Leukotriene Arachidonic 12(S)-HPETE 5(S)-HETE Leukotriene B4 cyp4f3 B4 acid alox12 gpx1a lta4h pla2g6 Leukotriene alox5a 5-HPETE alox5a A4 ltc4s Leukotriene C4 ggt5a Leukotriene D4

ptgs2a Prostaglandin Phosphatidylcholine G2 Prostaglandin I2 ptgis ptgs2a K17687 Prostaglandin E2 Prostaglandin ptges cyp4f3 20-HETE H2 Thromboxane cbr1 A2 cyp2u1 tbxas1 Prostaglandin F2 cbr1 alpha ptgdsb Prostaglandin D2 K07415 fam213b K17728

K17721 19(S)-HETE K17687 cyp2p10 cyp2u1

K07425 11,12-EET ephx2 11,12-DHET cyp2p10

Mouse B. Eicosanoids 5,6-EET Ephx2 5,6-DHET Nodes Edges rd1/rd1 -/- compound 8,9-EET Ephx2 8,9-DHET node Pde6b Rho Cyp2j6 node 14,15-EET Ephx2 C14775 complex

Cyp2c29 16(R)-HETE enrichment score -1 0 1 Alox8 8(S)-HPETE

Alox12b 12(R)-HPETE C14781 Arachidonic Cyp2j6 Cyp2j6 C14782 acid Linoleic acid Pla2g5 metaolism 12(S)-HPETE Cyp2j6 C14813 Cyp2j6 C14814 Alox12 15(S)-HPETE Gpx1 15(S)-HETE Phosphatidylcholine Alox15 20-OH-Leukotriene Alox5 5-HPETE Gpx1 5(S)-HETE Leukotriene B4 Cyp4f14 B4 Lta4h Ptgs1 Prostaglandin Leukotriene G2 Alox5 A4 Ltc4s Leukotriene C4 Ggt1 Leukotriene D4

Cyp4a10 Prostaglandin I2 Ptgis Cyp4f14 20-HETE Ptgs1 Prostaglandin E2 Cyp2u1 Prostaglandin Ptges H2 Thromboxane Cbr1 Tbxas1 A2 Cyp2e1 Cbr1 Prostaglandin F2 alpha K17728 Ptgdsb Prostaglandin D2 K17721 Fam213b 19(S)-HETE Cyp4a10 Cyp2j6 Cyp2u1

Cyp4a10 11,12-EET Ephx2 11,12-DHET Cyp2j6

Downloaded from iovs.arvojournals.org on 10/03/2021 Figure S1. Gene expression changes in the Arachidonic Acid Metabolism pathway. A) A schematic of the KEGG Arachidonic Acid Metabolism pathway with differential gene expression plotted

for genes with FDR ≤ 0.05 in zebrafish Müller glia. |Log2FC ≥ 1| is represented by the maximum and minimum color value according to the color bar. Genes without color (white) either do not have any data, or did not meet the FDR criterion. For each gene, the left half corresponds to expression changes at 8 hpl and the right half corresponds to changes at 16 hpl. B) A schematic of the KEGG Arachi- donic Acid Metabolism pathway with differential gene expression plotted for all genes with expression

data (no FDR threshold applied) in mouse Müller glia from retinal degeneration models. |Log2FC ≥ 1| is represented by the maximum and minimum color value according to the color bar. Genes without color (white) do not have any data. For each gene, the left half corresponds to expression changes in Pde6brd1/rd1 and the right half corresponds to changes in Rho-/-.

Downloaded from iovs.arvojournals.org on 10/03/2021 Circadian Rhythm A.

Creb1 Zebrafish 3’,5’-Cyclic AMP Nodes

Csnk1e node 8 hpl 16 hpl +p Per node Fbxw11 +u complex

Cul1 enrichment score Prkaa1 -1 0 1

+p Edges compound

Fbxl3 Cry1 expression Nr1d1 +u repression

activation Per Clock Bhlhe40 inhibition Cry1 +p Arntl phophorylation Rora +u Csnk1e ubiquitination

B. Creb1 3’,5’-Cyclic AMP Mouse Nodes Csnk1d rd1/rd1 -/- +p node Pde6b Rho Per Btrc +u node complex

Skp1a enrichment score Prkab1 -1 0 1 +p Edges compound Fbxl3 Cry1 Nr1d1 +u expression

repression

Per Clock activation Bhlhe40 Cry1 inhibition Arntl +p phophorylation Csnk1d Rora +u ubiquitination

Downloaded from iovs.arvojournals.org on 10/03/2021 Figure S2. Expression changes in the Circadian Rhythm pathway. A) A schematic of the KEGG Circadian Rhythm pathway with differential gene expression plotted for genes with a FDR ≤ 0.05 in

zebrafish Müller glia. |Log2FC ≥ 1| is represented by the maximum and minimum color value according to the color bar. Genes without color (white) either do not have any data, or did not meet the FDR criterion. For each gene, the left half corresponds to gene expression changes at 8 hpl and the right half corresponds to changes at 16 hpl. B) A schematic of the KEGG Circadian Rhythm pathway with differential gene expression plotted for all genes with expression data (no FDR threshold applied)

in mouse Müller glia retinal degeneration models. |Log2FC ≥ 1| is represented by the maximum and minimum color value according to the color bar. Genes without color (white) do not have any data. For each gene, the left half corresponds to expression changes in Pde6brd1/rd1 and the right half corresponds to changes in Rho-/-. Note that the KEGG pathway does not distinguish between Per1 and Per2. In zebrafish, per2 is positively enriched at both 8 and 16 hpl (5.36, 2.29) (Table S2; Suppl. File 1); in the two mouse degeneration models, Per1 (-2,54, -2.26) and Per2 (-1.14, -0.32) are negatively enriched (Suppl. File 1).

Downloaded from iovs.arvojournals.org on 10/03/2021 A. D.

B. E.

C. F.

Downloaded from iovs.arvojournals.org on 10/03/2021 Figure S3. Comparative analysis of gene expression in FACS-isolated Müller glia. Histogram plots of the frequency distribution of gene expression levels in FACS-isolated, GFP+ Müller glia as determined by microarray profiling or RNA-seq analysis. Colored vertical lines represent the expression levels of selected, cell-type specific genes in the associated key. A) Data from the present study (RNA-seq). (B,C) Published datasets from zebrafish Müller glia (microarray). Note: the microarray fromAgilent that was used in Ram- achandran et al., (B) did not contain probes for rho, gnat1, or rlbp1a. D,E, F) Published datasets from mouse Müller glia. Information about the published datasets and cell-type specificity of the selected genes can be found in SupplementaryTable S4 and Supplemen- tary Table S5, respectively.

Downloaded from iovs.arvojournals.org on 10/03/2021 A. D. G.

B. E. H.

C. F. I.

Downloaded from iovs.arvojournals.org on 10/03/2021 Figure S4. Comparative analysis of gene expression in single mouse Müller glia. A-H) Histogram plots of the frequency distribution of gene expression levels determined by microarray profiling from a published dataset of hand-picked, individual Müller glia. Colored vertical lines represent the expression levels of selected, cell-type specific genes in the associated key. A-E) 8 wk wild-type

(WT) mouse Müller glia. F-H) P13 WT mouse Müller glia. I) Log2 Rlbp1 expression values plotted (y-axis) vs. log2 Rho expression values (x-axis) for each of the individually, hand-picked WT mouse Müller glia from P13 (blue) or 8wk (red) mice. Information about the published datasets and cell-type specificity of the selected genes can be found in SupplementaryTable S4 and Supplementary Table S5, respectively.

Downloaded from iovs.arvojournals.org on 10/03/2021 A. D.

B. E.

C. F.

Downloaded from iovs.arvojournals.org on 10/03/2021 Figure S5. Comparative analysis of gene expression in FACS-isolated mouse rod and cone photoreceptors. A, B, D, E) Histo- gram plots of the frequency distribution of gene expression levels determined by RNA-seq from published datasets of FACS-isolated, GFP+ mouse photoreceptors. Colored vertical lines represent the expression levels of selected, cell-type specific genes in the asso-

ciated key. A) Adult, wild-type (WT) rods. B) Adult, WT cones. C) Log2 gene expression levels in rods vs. cones of selected cell-type specific genes, as indicated in the color code in the associated key. D) P28, WT Nrl:GFP mouse rod photoreceptors. E) P28, Nrl-/-

;Nrl:GFP mouse S-cone-like photoreceptors. Note: this mutant retina does not contain rods. F) Log2 gene expression levels in rods vs. S-cone-like photoreceptors of selected cell-type specific genes, as indicated in the color code in the associated key. Information about the published datasets and cell-type specificity of the selected genes can be found in SupplementaryTable S4 and Supplementary Table S5, respectively.

Downloaded from iovs.arvojournals.org on 10/03/2021