Abstract
Retinal photoreceptor cells contain the highest concentration of docosahexaenoic acid (DHA) in our bodies, and it has been long assumed that this is critical for supporting normal vision. Indeed, early studies using DHA dietary restriction documented reduced light sensitivity by DHA-deprived retinas. Recently, it has been demonstrated that a major route of DHA entry in the retina is the delivery across the blood–retina barrier by the sodium-dependent lipid transporter, Mfsd2a. This discovery opened a unique opportunity to analyze photoreceptor health and function in DHA-deprived retinas using the Mfsd2a knock-out mouse as animal model. Our lipidome analyses of Mfsd2a−/− retinas and outer segment membranes corroborated the previously reported decrease in the fraction of DHA-containing phospholipids and a compensatory increase in phospholipids containing arachidonic acid. We also revealed an increase in the retinal content of monounsaturated fatty acids and a reduction in very long chain fatty acids. These changes could be explained by a combination of reduced DHA supply to the retina and a concomitant upregulation of several fatty acid desaturases controlled by sterol regulatory element-binding transcription factors, which are upregulated in Mfsd2a−/− retinas. Mfsd2a−/− retinas undergo slow progressive degeneration, with ∼30% of photoreceptor cells lost by the age of 6 months. Despite this pathology, the ultrastructure Mfsd2a−/− photoreceptors and their ability to produce light responses were essentially normal. These data demonstrate that, whereas maintaining the lysophosphatidylcholine route of DHA supply to the retina is essential for long-term photoreceptor survival, it is not important for supporting normal phototransduction.
SIGNIFICANCE STATEMENT Phospholipids containing docosahexaenoic acid (DHA) are greatly enriched in the nervous system, with the highest concentration found in the light-sensitive membranes of photoreceptor cells. In this study, we analyzed the consequences of impaired DHA transport across the blood–retina barrier. We have found that, in addition to a predictable reduction in the DHA level, the affected retinas undergo a complex, transcriptionally-driven rebuilding of their membrane lipidome in a pattern preserving the overall saturation/desaturation balance of retinal phospholipids. Remarkably, these changes do not affect the ability of photoreceptors to produce responses to light but are detrimental for the long-term survival of these cells.
Introduction
There is ample evidence supporting the health benefits of consuming ω-3 fats, i.e., phospholipids containing docosahexaenoic acid (DHA) moieties. DHA is highly concentrated in the nervous system despite the fact that DHA cannot be synthetized by neurons and has to be transported to the brain across the blood–brain barrier (BBB). A major breakthrough in understanding this transport was reported by Nguyen et al. (2014) who discovered that facilitative transporter Mfsd2a, expressed in the endothelium of the BBB microvessels, is responsible for the uptake of plasma lysophosphatidylcholines (LPCs) across the BBB, thereby serving as a primary mechanism of DHA delivery to the brain. Mfsd2a knock-out, as well as its mutations, results in microcephaly and hypomyelination in animal models and human patients (Nguyen et al., 2014; Alakbarzade et al., 2015; Guemez-Gamboa et al., 2015; Chan et al., 2018; Harel et al., 2018).
Remarkably, the retina contains even more DHA than the brain, with the highest concentration found in the light-sensitive outer segment membranes of rod and cone photoreceptor cells (Fliesler and Anderson, 1983). It is generally agreed that, like the brain, the retina cannot synthesize its own DHA de novo and relies on uptake from extraretinal sources, such as blood-borne lipids, for maintaining its DHA content (Scott and Bazan, 1989; Bazan et al., 2011). Multiple studies analyzing the consequences of DHA dietary restriction documented that DHA-deprived retinas produced light responses of reduced sensitivity (Benolken et al., 1973; Wheeler et al., 1975; Senapati et al., 2018), suggesting that DHA may be critical for supporting visual signaling. Recently, Silver and colleagues showed that Mfsd2a is expressed in retinal pigment epithelium (RPE) and retinal vasculature and that the Mfsd2a-mediated lipid transport in RPE is particularly important for maintaining the high DHA concentration in the retina (Wong et al., 2016). They reported that whole eyes of Mfsd2a−/− mice contain ∼40% less DHA than WT controls and that this decrease is compensated by an increase in the content of arachidonic acid (AA). They also reported that Mfsd2a−/− photoreceptors appear to have a morphological defect in their outer segment disc structure. Yet, their analysis of visual function by electroretinography yielded nearly normal light responses.
Here, we analyzed the retinal phenotype of an alternative Mfsd2a−/− mouse strain, which corroborated and expanded several key observations from the previous report (Wong et al., 2016), disagreed with two and revealed a number of novel findings (for a detailed comparison between the two studies, see Discussion). We have found that Mfsd2a−/− retinas undergo a global remodeling of their fatty acid composition, including a 45% reduction in the fraction of DHA lipids and a compensatory 66% increase in the fraction of AA. We also observed a 56% increase in the fraction of phospholipids containing monounsaturated fatty acid moieties and a 57% loss in phospholipids containing very long chain fatty acid moieties. Similar changes were observed in an independently analyzed lipidome of isolated photoreceptor outer segment membranes. Mfsd2a−/− photoreceptors underwent relatively slow progressive degeneration, although surviving photoreceptors retained an entirely normal ultrastructure. Remarkably, single-cell recording from Mfsd2a−/− rods did not reveal any major phototransduction phenotype. This indicates that impaired LPC supply to the retina and the corresponding ∼twofold reduction in the content of retinal DHA are detrimental for long-term viability of photoreceptor cells, but do not impact the ability of surviving cells to produce normal responses to light.
Materials and Methods
Animals.
Frozen embryos of Mfsd2a heterozygote mice were purchased from the KOMP repository (032467-UCD) and recovered by the Duke Transgenics Core. One of the recovered heterozygous female mice was used as a founder. This is the same line as used by Ben-Zvi et al. (2014) and Chow and Gu (2017), except that we backcrossed the line used in this study to the C57BL/6J background (Jackson Laboratories, stock #000664) for eight generations. The line was negative for rd1 and rd8 mutations. Mice were maintained on the Teklad global 18% protein rodent diet 2918 (Envigo). At 1 month of age Mfsd2a−/− mice were on average ∼25% smaller than their WT littermates [14.8 ± 1.1 g vs 19.7 ± 0.9 g for Mfsd2a−/− vs Mfsd2a+/+ males (mean ± SD; n = 4); 12.7 ± 1.3 g vs 17 ± 1.2 g for Mfsd2a−/− vs Mfsd2a+/+ females (n = 5)]. Knock-out and WT littermate mice used for experiments were obtained by in-crossing the Mfsd2a+/− mice. Note, that Mfsd2a−/− mice in our hands could produce progeny as well. Based on monitoring three Mfsd2a−/− breeding pairs, we noted that one breeding pair produced no litters, whereas others produced one to two litters with one or two surviving mice. RhoP23H/P23H mice were purchased from Jackson Laboratories (stock #017628). Animals for experiments were obtained by mating RhoP23H/WT mice with C57BL/6J WT mice (Jackson Laboratories, stock #000664). Mouse genotypes were determined using real time PCR with specific probes designed for each gene (Transnetyx).
Animals were reared under a normal day/night cycle and handled according to the protocols approved by the Institutional Animal Care and Use Committees of Duke University and the University of Florida.
Antibodies.
For Western blotting, mouse monoclonal antibody 4D2 against rhodopsin (1:5000) was a gift from R.S. Molday (University of British Columbia). Rabbit anti-Gαt (sc-389; 1:10,000) and anti-phosducin (sc-23774; 1:5000) antibodies were from Santa Cruz Biotechnology. Rabbit anti-arrestin antibody (PA1–7; 1:5000) was from Affinity Bioreagents. Rabbit antibodies against retinal guanylate cyclase isoforms 1 (1:2000) and 2 (1:2000) were a gift from A. M. Dizhoor (Pennsylvania College of Optometry). Rabbit anti-heat shock 70 protein antibodies (ADI-SPA-816; 1:5000) were from Enzo Life Sciences. Secondary goat or donkey antibodies for Western blotting conjugated with AlexaFluor 680 (A-21057) and 800 (SA5–10044; both 1:10,000) were from Invitrogen.
For immunohistochemistry, we used mouse monoclonal antibody 1D4 against rhodopsin (ab5417; 1:1000) from Abcam. Secondary donkey antibodies conjugated with AlexaFluor 488 (A-21202; 1:1000) were from Life Technologies.
Histology and microscopy.
Mice were killed by CO2 and transcardially perfused with 3 ml of 2% paraformaldehyde and 2.5% glutaraldehyde solution in PBS at a flow rate of 1–2 ml/min. The superior–inferior line was marked with a burnmark using low temperature cautery (Bovie, AA90). Next, eyes were enucleated, exterior tissue was removed, and eyes were fixed overnight at 4°C in 2% paraformaldehyde and 2.5% glutaraldehyde in PBS. The next day, eyes were rinsed with PBS, cornea and lens were removed, treated with 2% osmium tetroxide (Petters et al., 1997) and embedded in EPON 812 (Electron Microscopy Sciences). Plastic-embedded retinal cross-sections (0.5 μm) were stained with toluidine blue for light microscopy (Lobanova et al., 2013). Nuclear count was performed in sections cut through the optic nerve in 100 μm segments of the outer nuclear layer at 500 μm steps from the optic nerve head as described by Lobanova et al. (2018). For electron microcopy, 60–80 nm ultrathin sections were cut from the same blocks with an ultra-microtome (Leica), collected on copper grids (Electron Microscopy Sciences) and post-stained with 1% uranyl acetate and Sato's lead (Sato, 1968).
For immunofluorescent localization of rhodopsin, enucleated eyes were cleaned from extraocular tissue, drop fixed in 4% paraformaldehyde in PBS overnight at 4°C, washed with PBS and embedded into paraffin at Histology Tech Services. Four-micrometer-thick paraffin sections were cut through the superior–inferior line through the optic nerve, collected on slides, and stored at room temperature until used. Immediately before immunolocalization analysis, sections were deparaffinized using Histoclear II Reagent (HS-202, National Diagnostics), rehydrated, permeabilized with 0.3% Triton-100 in PBS, blocked with 3% donkey serum in PBS, probed with primary anti-rhodopsin antibody, rinsed with PBS, stained with secondary antibody and Hoechst (1:1000), rinsed with PBS, mounted with mounting media (Immu-Mount, ThermoFisher Scientific) under coverslips (Fisherbrand, 12541B) and visualized with an A1 confocal microscope (Nikon).
The TUNEL (terminal deoxynucleotidyl transferase dUTP nick end labeling) assay was performed on deparaffinized retinal sections using the Fluorescein In Situ Cell Death Detection Kit (11684795910, Roche) following the manufacturer's instructions. Results were expressed as an average number of TUNEL-positive nuclei per retinal section. Sections prepared from both retinas of two animals (1 male and female) of each genotype were analyzed.
Western blotting.
Retinas from 5-week-old mice were collected between 12:00 and 2:00 P.M., carefully dissected, snap-frozen in liquid nitrogen and stored at −80°C until used. Samples for Western blotting were obtained by solubilizing one mouse retina by sonication in 120 μl of 1% Triton X-100 in PBS supplemented with protease inhibitor cocktail (Roche). Total protein concentration was measured using the DC Protein Assay kit (Bio-Rad) and samples were diluted with SDS-PAGE sample buffer to achieve a protein concentration of 1 mg/ml. A small portion of each sample was used to analyze rhodopsin, whereas the rest of the sample was boiled for 5 min before analysis. 0.5 μg of total protein was loaded by SDS-PAGE gel for rhodopsin and 20 μg for other proteins. Protein bands were visualized and quantified using the Odyssey Infrared Imaging System (LI-COR Biosciences).
Rod outer segment isolation.
Retinas from six 5-week-old mice of each genotype were extracted and flash-frozen. Outer segments were isolated from pooled retinas using a procedure described by Tsang et al. (1998) with slight modifications. Combined retinas were thawed 200 μl of buffer containing 20 mm HEPES, pH 7.4, 100 mm KCl, 2 mm MgCl2, and 1 mm EDTA supplemented with 27% sucrose. Samples were vortexed for 30 s at maximum speed and spun down at for 30 s at 200 rcf on a tabletop Eppendorf centrifuge. The supernatant was removed and placed on ice in a 11 × 34 mm centrifuge tube (Beckman Coulter). Another 200 μl of the same buffer with 27% sucrose was added to the remaining pellet and the vortexing/spinning sequence was repeated. This second supernatant was combined with the first and the total volume was then brought to 800 μl. This solution containing 27% sucrose was then underlaid with 1 ml of the same buffer containing 32% sucrose. Tubes were transferred to a TLS55 swinging-bucket rotor and spun at 25,000 rpm for 20 min at 4°C in an Optima Max-e ultracentrifuge (Beckman Coulter). Approximately 200 μl material representing the outer segment band formed at the 27–32% sucrose interface was removed and transferred to a 1.5 ml specialized low-retention Eppendorf tube (catalog #022363204) on ice. The fraction was diluted by 1 ml of the same buffer without sucrose and spun down at 19,000 rpm in a Beckman-Coulter Avanti J25 ultracentrifuge for 30 min at 4°C. The supernatant was then discarded and the pellet resuspended solubilized in 30 μl buffer without sucrose. Protein concentration was measured using the Bradford assay and the rest used for lipid analysis.
RNA-seq.
For whole-retina sequencing, total RNA was prepared from the retinas of three Mfsd2a−/− (all females) and Mfsd2a+/+ (2 males and 1 female) mice of 5–7 weeks of age using the RNeasy Mini Kit (Qiagen) with DNase digestion (RNase-Free DNase Set, Qiagen). Retinas were collected between 1:00 and 2:00 P.M. from mice killed with isoflurane, carefully dissected under a microscope and snap-frozen in liquid nitrogen. The same procedure was used to prepare retinas from two RhoP23H/WT mice (1 male and 1 female) and their WT littermates (1 male and 1 female) at P33. Libraries were prepared and sequenced at GENEWIZ. FASTQ files were aligned to mouse genome GRCm38.p3 using DNAStar software v15.3.0 and normalized by either reads per kilobase of transcript (RPKM) or using the DeSeq2 method. Genes expressed at the level exceeding 1 RPKM or exceeding the linear threshold of 1 in DeSeq2 normalization in all samples were filtered out in ArrayStar. Genes that changed for >15% with p <0.05 (as determined using either two-tailed homoscedastic t test for RPKM normalization or p value for DeSeq2 normalization) were considered differentially expressed. Fold-changes in the expression of genes were calculated based on the mean values for Mfsd2a−/− and WT littermates. Changes in expression levels of representative genes were plotted in GraphPad as percentage of average values for WT littermate mice calculated from datasets normalized by RPKM. A list of differentially expressed genes in Mfsd2a−/− mice and selected genes from RhoP23H/WT mice is presented in Fig. 10-1.
To establish cellular pathways affected by Mfsd2a loss, differentially expressed genes listed in Figure 10-1 were analyzed with the Ingenuity Pathway Analysis Software (v47547484, Qiagen) using the default parameters. Enriched pathways with p < 0.05 are listed in Fig. 10-2.
Single-cell sequencing analysis was performed using 8-week-old albino mice on the BALB/cJ, background (Jackson Laboratories, stock #000651) exposed to 2000 lux white light for 4 h. The mice were killed by CO2 followed by cervical dislocation. The retinas were quickly collected for immediate dissociation using reagents and protocols of the Papain Dissociation System (Worthington Biochemical). Retinas were placed in the pre-equilibrated papain solution [20 U/ml papain in 1.0 mm[scap] l-cysteine with 0.5 mm EDTA and 0.005% DNase I in Earle's balanced salt solution (EBSS)] and incubated in a 37°C water bath for 30 min, gently mixing every 5–10 min, to break down the tissue into a single-cell suspension. The suspension was triturated by gentle pipetting and any remaining tissue clumps were allowed to settle at the bottom of the tube. The suspended cells were then transferred to a new tube and pelleted using a swinging bucket centrifuge at 300 rcf for 5 min at room temperature. The papain solution was removed, and the cells were gently resuspended in pre-equilibrated ovomucoid-albumin protease inhibitor solution and DNase (2000 U/ml in EBSS). To remove cell fragments, a discontinuous density gradient was prepared by carefully layering the cell suspension on top of pre-equilibrated EBSS containing papain inhibitors (10 mg/ml ovomucoid and 10 mg/ml albumin) and centrifuging at 70 rcf for 6 min at room temperature. The supernatants were removed and the cells were washed by gently resuspending in cold Dulbecco's PBS (DPBS) and centrifuging at 300 rcf for 5 min. For the second wash, the cells were resuspended in DPBS and gently passed through a 40 μm filter before centrifuging at 300 rcf for 5 min. The pelleted cells were fixed by gently resuspending in DPBS, followed by the addition of cold methanol (90% final concentration) dropwise to the cells, vortexing at the lowest speed setting to avoid cell clumping. The cells were incubated on ice until fixation of all cells was complete (∼15 min). Cell viability and counts were assessed before and after methanol fixation using Trypan blue on both the Countess II FL Automated Cell Counter (Life Technologies) and manually using a hemocytometer. Samples were stored at −80°C overnight before proceeding with the rehydration process. Libraries were prepared using the 10X Genomics Single Cell 3′/5′ RNAseq Library Prep kit (v2 Chemistry) and sequenced on the Illumina NextSeq 500 at UF Interdisciplinary Center for Biotechnology Research core facility. Datasets were processed with the 10X Chromium Cell Ranger pipeline using default parameters, and further analyzed using the 10X Chromium Loupe Browser.
Annotation of lipid classes and fatty acids.
PAs, Phosphatidic acids; PEs, phosphatidylethanolamines; PE O-s, 1-O-alkyl-2-acylglycerophosphoethanolamines; LPEs, lysophosphatidylethanolamines; PSs, phosphatidylserines; LPSs, lysophosphatidylserines; PCs, phosphatidylcholines; PGs, phosphatidylglycerols; PC O-s, 1-O-alkyl-2-acylglycerophosphocholines; SMs, sphingomyelins; PIs, phosphatidylinositols; LPIs, lysophosphatidylinositols; TGs, triacylglycerols; DGs, diacylglycerols; Cers, ceramides; CLs, cardiolipins; GalCers, galactosylceramides; CFs or Chols, cholesterols; CEs, cholesterol esters. Lyso- and phospholipid species were annotated using the total number of carbon atoms and total number of double bonds in fatty acid moieties. DHA (22:6n-3); DPAn-6, docosapentaenoic acid (22:5n-6); AA (20:4n-6); VLCFA, very long fatty acid, defined in this study as one containing 24 and more carbons; MUFA, monounsaturated fatty acid (e.g., 16:1, 18:1, 20:1); SFA, saturated fatty acid (e.g., 16:0, 18:0, 20:0).
Quantitative shotgun lipidomics.
Eyes were collected from mice transcardially perfused with 30 ml of PBS. Retinas were extracted in Hank's solution, collected in 2.0 ml Eppendorf tubes (#0030 120.094), topped off with nitrogen, flash frozen in liquid nitrogen, stored at −80°C, and transported on dry ice. Retinas were collected from three Mfsd2a−/− (2 females and 1 male) and three Mfsd2+/+ (2 males and 1 female) mice at P38–P40.
Before lipid extraction, two retinas were homogenized at 4°C in 250 μl of isopropanol using 1 mm zirconium beads (BioSpec, 11079110z) on a TissueLyser (Qiagen, 85300). Retina and outer segment membrane preparations were dried under vacuum and lipids were extracted with 1.4 ml or 0.7 ml, respectively, of methyl tert-butyl ether (MTBE)/MeOH (10:3; v:v) spiked with internal lipid standards. The internal standard mixture for retinal samples contained the following: 1.59 nmol CF D7, 1.81 nmol CE D7 16:0, 1.49 nmol TG D5 50:0, 1.04 nmol DG D5 34:0, 1.47 nmol PC 25:0, 1.94 nmol LPC 13:0, 0.48 nmol PS 25:0, 1.05 nmol PE 25:0, 1.52 nmol LPE 13:0, 0.86 nmol PI 25:0, 0.11 nmol LPI 13:0, 0.49 nmol PG 25:0, 1.3 nmol Cer 30:1, 0.48 nmol SM 30:1, 0.49 nmol GalCer 30:1, 0.55 nmol PA 25:0. The internal standard mixture for outer segment membranes contained 1.78 nmol CF D7, 2.21 nmol CE D7 16:0, 1.04 nmol TG D5 50:0, 0.6 nmol DG D5 34:0, 1.38 nmol PC 25:0, 0.39 nmol LPC 13:0, 0.14 nmol PS 25:0, 0.59 nmol PE 25:0, 0.09 nmol LPE 13:0, 0.48 nmol PI 25:0, 0.06 nmol LPI 13:0, 0.14 nmol PG 25:0, 0.07 nmol Cer 30:1, 0.48 nmol SM 30:1, 0.27 nmol GalCer 30:1, 0.13 nmol PA D7 33:0. Next, 280 μl (retina) or 140 μl (outer segments) of water was added and the mixture incubated in an Eppendorf mixer for 1.5 h at 4°C.
Phase separation was induced by centrifugation at 13,400 rpm for 30 min at 4°C on a benchtop centrifuge. The organic phase was transferred into glass vial and stored at −20°C until analyzed. Three hundred microliter aliquots of the extract were dried under vacuum, resuspended in 50 μl of MTBE/ methanol/ water (10:3:1; v:v:v) and subjected to thin-layer chromatography (Macherey-Nagel, Nano-SILGUR-20; 10 × 10 cm) using a mobile phase chloroform/ethanol/water/trimethylamine (TEA; 30:35:7:35; v:v:v:v). Three fractions were collected by scraping the silica layer and lipids re-extracted with MTBE/methanol as described above, but without adding internal standards. Finally, 1.2 ml of organic extracts from each fraction were collected and stored in glass vials at −20°C. For mass spectrometric analysis, aliquots of collected fractions were dried and re-dissolved in isopropanol/methanol/chloroform mixture (4:2:1; v:v:v) with 7.5 mm ammonium formate or in ethanol/chloroform 5:1 (v:v) mixture with 0.1% TEA. 170 μl of the fraction 1 (Rf = 0.5–1.0) were dissolved in 180 μl ammonium formate or the TEA mix, 80 μl fraction 2 (Rf = 0.24–0.5), in 200 μl of the TEA mix, and 80 μl of fraction 3 (Rf = 0–0.24) in 200 μl NH4HCO2 or the TEA mix.
The fractions were analyzed on a hybrid quadrupole, Orbitrap tandem mass spectrometer Q Exactive (ThermoFisher Scientific), equipped with a robotic nanoflow ion source TriVersa NanoMate (Advion BioSciences) using 4.1 μm nanoelectrospray chips. The ion source was controlled by the Chipsoft 8.3.1 software (Advion BioSciences). Ionization voltage and back pressure were set to ±0.96 kV and 1.25 psi for mixtures containing ammonium formate and −2.0 kV and 1.0 psi for mixtures containing TEA. The temperature of the ion transfer capillary was 200°C; S-lens RF level was set at 50%. For lipidomic measurements FT MS spectra in positive or negative mode were acquired within the m/z range of 400–1200 at the target mass resolution of Rm/z 200 = 140,000; automated gain control (AGC) was set at 3 × 106 and the maximum injection time at 3000 ms. Subsequently, HCD FT MS/MS spectra were acquired at the same mass resolution using an inclusion list of m/z of lipids recognized in FT MS spectra. Precursor isolation window was 1.0 Da; AGC was set at 2 × 104 and maximum injection time at 650 ms. Normalized collision energies (15% in positive mode and 25% in negative mode) were used to confirm identity of lipids. Cholesterol was quantified as described by Sales et al. (2016). Acquired spectra were filtered by PeakStrainer (https://git.mpi-cbg.de/labShevchenko/PeakStrainer/wikis/home; Schuhmann et al., 2017) and lipids were identified and quantified by LipidXplorer software (Herzog et al., 2011) considering better than 5 ppm mass accuracy and signal-to-noise ratio above the value of 3.0. Cardiolipins were quantified using the PG standard. CE, TG, GalCer, and PA were detected at trace amounts and were not quantified.
To estimate changes in the amounts of SFA, MUFA, DHA, AA, and VLCFA (defined as having 24 or more carbons) in whole retinas shown in Figure 7D, we selected lysophospholipids and phospholipids containing these fatty acids based on previous fragmentation studies performed with rat retinas (Schuhmann et al., 2011). The relative amounts of fatty acids of each type were determined by taking into consideration the presence of two fatty acids in phospholipids and one fatty acid in lysophospholipids. Results were expressed as percentage of total fatty acids incorporated into all measured lysophospholipids and phospholipids.
Changes in the fraction of outer segment phospholipids containing MUFA, DHA, and AA (see Fig. 9B) were calculated from molar fractions of major PC, PE, PS, and PI species containing fatty acid of each type and shown as percentage of total phospholipids. Phospholipids used in these calculations are highlighted on Fig. 7-1 (tabs “Major OS PL %” and “Major retinal PL %”).
Fundus fluorescein angiography.
Mice were anesthetized with a mixture of ketamine-xylazine (95–9 mg/kg, respectively) and their eyes were dilated with ophthalmic solutions of atropine (1%) and phenylephrine (2.5%). Mice were administered an intraperitoneal injection of 100 mg/kg fluorescein (∼0.02 ml of 100 mg/ml fluorescein). Images were acquired starting 1 min post-injection using a Micron IV fundus camera. Four Mfsd2a−/− and their WT littermates were analyzed at 6 weeks of age. Experiments were performed between 10:00 A.M. and 1:00 P.M.
Suction electrode recordings.
Suction electrode recordings from the outer segments of intact mouse rods were performed as previously described (Gross and Burns, 2010). Dark-adapted retinas were stored on ice in L-15 supplemented with 10 mm glucose. Recordings were performed in oxygenated, bicarbonate buffered Locke's solution supplemented with 10 mm glucose at 35–37°C. The suction pipette contained HEPES-buffered Locke's solution, pH 7.4, connected to a calomel electrode via an agar bridge. An agar bridge also connected the bath to the opposite calomel half-cell.
Brief (10 ms, 500 nm) flashes of calibrated strength were used to elicit electrical responses, which were amplified (AxoPatch 200B, Molecular Devices), filtered at 30 Hz with an 8-pole Bessel (Frequency Devices), and digitized at 200 Hz using custom-written acquisition procedures in IgorPro (WaveMetrics). Responses to a large number (>30) of dim flashes were averaged and used to determine the mean time to peak and integration time (time integral of the response divided by the peak amplitude) if the peak amplitude was <20% of the maximal response amplitude (dark current). Variance to mean analysis of the dim flash responses was used to calculate the average single photon response for each cell, and saturating flash responses were used to calculate the dominant time constant of recovery (τD) as previously described (Krispel et al., 2006; Gross and Burns, 2010). Statistical tests for differences between genotypes were determined for the population average kinetic parameters using a two-tailed t test assuming equal variances. Error bars represent SEM.
Statistical analysis.
Statistical analyses were performed in GraphPad Prism 5 software or Excel 2016. Results are expressed as the mean ± SD, unless noted otherwise. Differences are considered significant if p ≤ 0.05, as determined by a two-tailed homoscedastic t test; n.s. indicates no significant difference (p > 0.05).
Results
Slow photoreceptor degeneration in Mfsd2a−/− mice
We conducted morphometric analysis of Mfsd2a−/− retinas in thin plastic sections at 1 and 6 months of age. The number of photoreceptors was reduced by 12 and 30%, respectively, compared with WT littermates of the same ages (Figs. 1A, 2). The ongoing cellular loss in mutant mice was corroborated by an ∼sevenfold increase in the number of apoptotic nuclei in the outer nuclear layer detected by the TUNEL assay (Fig. 1B). The progressive nature of photoreceptor loss and the presence of apoptotic nuclei argue that the reduction in the photoreceptor cell count in young Mfsd2a−/− mice indeed reflects ongoing photoreceptor degeneration rather than some form of developmental deficiency.
However, the EM analysis of Mfsd2a−/− photoreceptors did not reveal any notable abnormalities in either the overall ultrastructure of outer segments from surviving rods (Fig. 3A) or the spacing of photoreceptor discs (Fig. 3A–C). Also normal were the appearance of the active zone of new disc formation at the outer segment base and the interface between the outer segment tip and the retinal pigment epithelium (Fig. 3B,C). The only distinctions from WT retinas were as follows: (1) an occasional appearance of malformed outer segments, apparently belonging to dying Mfsd2a−/− photoreceptors (Fig. 3A, red asterisk); and (2) the presence of mononuclear phagocytes in the subretinal space (Fig. 3D). The latter are most likely microglia, as shown by Wong et al. (2016). No overt abnormalities in other retinal layers were observed as well. The normal appearance of outer segment ultrastructure was different from (Wong et al., 2016) who observed gross morphological abnormalities in Mfsd2a−/− photoreceptors in both young and aged animals. One potential explanation is that they characterized a Mfsd2a−/− mouse obtained from a different source that contained the rd8 mutation (D. L. Silver, personal communication), whereas our line was back-crossed to the rd8-free C57BL/6J background and specifically tested for the absence of rd8. It is conceivable that rd8 serves as a genetic factor confounding the Mfsd2a−/− phenotype. We also cannot exclude the possibilities that these differences are explained by using different diets or subtle differences in tissue processing protocols.
Consistent with only a minor photoreceptor loss, we did not observe a significant reduction in the expression levels of representative phototransduction proteins, measured in whole retinal lysates of 1-month-old mice (Fig. 4A). We also assessed the localization of rhodopsin in rod photoreceptors of Mfsd2a−/− mice. Considering a previous report that rhodopsin is partially mislocalized from the outer to inner segments of Mfsd2a−/− rods (Wong et al., 2016), we used thin retinal sections, which aids inner segment rhodopsin detection, and optimized image contrast to visualize rhodopsin in this region. Despite these precautions, we observed no difference in inner segment rhodopsin staining between the retinas of Mfsd2a−/− mice and their WT littermates (Fig. 4B), indicating that rhodopsin targeting and trafficking are not affected by the Mfsd2a knock-out.
Mfsd2a expression in the retina
Previous study demonstrated robust expression of Mfsd2a in the RPE and retinal vasculature endothelium (Wong et al., 2016). We investigated whether the Mfsd2a gene could be additionally expressed in other retinal cell types using single-cell transcriptomic analysis. We used a large single-cell transcriptome dataset obtained from an ongoing study of light damaged retinas. Mfsd2a transcript was detected in four clearly separated cell clusters (Fig. 5). As expected, one of them represented RPE (based on expression of the Rpe65 marker) and another represented vascular endothelial cells (based on expression of the Tie1 marker). However, we also detected robust Mfsd2a expression in cell clusters expressing the Muller cell marker, LIF (Joly et al., 2008), and the astrocyte marker, Pdgfra (Mudhar et al., 1993), suggesting that Mfsd2a participates in the lipid transport across the membranes of these glial cells. Expression of Mfsd2a in retinal astrocytes is consistent with single-cell transcriptomics analysis of the brain (Vanlandewijck et al., 2018). Unfortunately, we could not identify antibodies that could corroborate these findings by immunohistochemistry. It is also worth noting that transcriptional analysis performed in a recent study (Benedicto et al., 2017) suggested the presence of Mfsd2a in the choroidal cells, which were absent in our dataset and not investigated here.
Analysis of the retina–blood barrier in Mfsd2a−/− mice
There has been a controversy between the studies by two other laboratories analyzing the brain phenotype of the Mfsd2a−/− mouse. Although Ben-Zvi et al. (2014) claimed that the major consequence of the Mfsd2a knock-out is the blood–brain barrier leakage, Wong et al. (2016) did not observe blood–retina barrier leakage in either the brain or the retina. This prompted us to interrogate whether the lack of Mfsd2a causes a disruption of the blood–retina barrier in our mouse (which was the same strain as used by Ben-Zvi et al., 2014). We used fluorescein angiography, which assesses the integrity of tight junctions in the retinal vasculature endothelium in vivo. Consistent with results reported by Wong et al. (2016), our analysis revealed no evidence of retinal vascular leakage, at least as assessed by this technique (Fig. 6).
Lipid composition of Mfsd2a−/− retinas
We next compared the lipid composition of retinas obtained from 6-week-old Mfsd2a−/− and WT mice. Lipid extracts, pre-fractionated by TLC, were subjected to quantitative top-down shotgun analysis on a high-resolution LTQ Orbitrap mass spectrometer (Schuhmann et al., 2011, 2012). Peaks representing intact molecular ions of lipids were recognized by their accurately determined masses and, therefore, the detection specificity was not biased by the composition of their fatty acid moieties. Notably, alternative analytical approaches (e.g., multiple reaction monitoring) that rely on the detection of carboxylate fragments produced from molecular anions of glycerophospholipids could discriminate against PUFA- containing and, particularly, VLCFA-containing lipid species. Unlike carboxylate fragments of saturated fatty acid moieties, the unsaturated fragments are unstable, undergo facile collisional degradation, and therefore are detected as peaks with much lower amplitudes. Using the shotgun methodology allowed us to identify multiple VLCFA species with >30 carbon atoms and 6 double bonds, whereas none of them were reported by Wong et al. (2016), who relied on targeted mass spectrometry.
Because Mfsd2a is the LPC transporter, we first analyzed the retinal content of this phospholipid and observed a remarkable shift in the saturation profile of its fatty acid moieties: of 17 detected LPC species, 15 were statistically significantly changed (Fig. 7A). These changes could be summarized in four groups: (1) a reduction of LPC with fatty acids containing 22 and more carbons, such as LPC22:6 (containing DHA), LPC24:6, LPC32:6 and LPC34:6; (2) an increase in LPC20:4 (containing AA); (3) an increase in LPC containing MUFAs (LPC16:1, LPC18:1, LPC20:1); and (4) miscellaneous changes in proportions of LPC species containing 0, 2, 3 and 5 double bonds.
Figure 7-1
In a similar pattern, we observed a reduction in the major PC species containing DHA (such as PC38:6, PC40:6, and PC44:12) and those containing 24 and more carbons (VLCFA; PC46:12, PC54:12, PC56:12), as well as an increase in the fraction of PC containing AA (PC36:4, PC38:4) and MUFAs (PC32:1, PC34:1, PC36:1; Fig. 7B). This trend was present for other major phospholipid classes—PE, PS, and PI—although the extent of fatty acid moieties' remodeling was smaller than for PC (Fig. 8).
Consistent with the role of Mfsd2a in LPC transport across the blood-retina barrier, we observed an ∼19% reduction in the total molar fraction of LPC (Fig. 7C). The contents of PE O- and PG species were decreased as well, whereas small changes in other lipid classed did not reach statistical significance.
Major changes in the patterns of membrane lipid composition between WT and Mfsd2a−/− retinas are summarized in Figure 7, D and E. Mfsd2a loss resulted in a 45% reduction of the DHA content in retinal glycerophospholipids and a 57% reduction in the total fraction of VLCFAs, with some individual species reduced by over 10-fold (e.g., PC54:12). On the other hand, the retinal content of AA was increased by 66% and the content of MUFAs by 56%. These changes appear to reflect a compensatory mechanism aiming to maintain a relatively constant overall level of lipid saturation in the retina, perhaps required for maintaining the fluidity of cell membranes.
To confirm that a similar pattern of DHA deficiency takes place specifically in photoreceptor outer segments, we performed lipidomics analysis of rod outer segment membranes isolated from the retinas of both genotypes. Due to a very small amount of available material (particularly from poorly breeding Mfsd2a+/− mice) our quantitative analysis was restricted to four major phospholipid classes (PC, PE, PS, and PI; Fig. 7-1). The changes in PC are illustrated in Figure 9A. Overall, we observed an ∼56% loss of phospholipids containing DHA, and compensatory changes, including an ∼147% increase in AA and ∼84% increase in MUFAs (Fig. 9B). These relative changes are somewhat larger than those observed in whole retinas (for direct comparison, see Fig. 9B), perhaps reflecting the fact that outer segment membranes undergo rapid turnover and are completely replaced in the course of ∼10–12 d in mice (Young, 1967). We also documented a nearly complete loss of major PC species containing two DHA moieties (PC44:12) and VLCFAs containing six double bonds (e.g., PC54:12, PC56:12, PC58:12; Fig. 9A).
This global remodeling of the retinal phospholipid composition could be explained by a combination of altered fatty acid transport, particularly that of DHA, across the blood-retina barrier and an upregulation of desaturase enzymes. The latter was previously documented in the Mfsd2a−/− mouse eyecups analyzed at P13 when retinal development is not yet completed (Wong et al., 2016). To elucidate whether such regulatory mechanisms persist in adult retinas, we used comparative RNA-Seq analysis to evaluate differential gene expression in the retinas of 5- to 6-week-old Mfsd2a−/− mice and their WT littermates (this was the same age at which our lipidome analysis was performed). RNA-Seq data were mapped with DNASTAR, normalized by the RPKM values, and gene expression levels in Mfsd2a−/− retinas were calculated as percentage of mean values in WT retinas. In agreement with previous report (Wong et al., 2016), we observed an upregulation of Stearoyl-CoA desaturases 1 and 2 (Scd1 and Scd2), fatty acid desaturases 1 and 2 (Fads1 and Fads2), as well as sterol regulatory element-binding factor Srebf1, which reported to control their expression in other experimental systems (Seo et al., 2009; Knebel et al., 2012; Angela et al., 2016; Fig. 10 and Fig. 10-1). Also in agreement with Wong et al. (2016), we observed a small but statistically significant ∼15% increase in the level of sterol regulatory element-binding factor Srebf2, which controls the expression of genes regulating cholesterol metabolism. A complementary analysis of the RNA-Seq data using the DeSeq2 normalization corroborated statistically significant upregulation of desaturases (Scd1, Scd2, Fads1, Fads2) and Srebf1, but not Srebf2 (Fig. 10-1).
Figure 10-1
Figure 10-2
Consistent with these findings, the enrichment pathway analysis performed on the data normalized by either RPKM or DeSeq2 (Fig. 10-2) suggested that Mfsd2a−/− retinas have elevated activity of the pathways, which include lipid desaturases (“Oleate Biosynthesis II” and “γ-linolenate Biosynthesis II”) and various lipid pathways which involve the Srebf1 and/or Srebf2 transcription factors (“LXR/RXR Activation”, “FXR/RXR Activation”). We also noted an upregulation of several genes indicative of ongoing photoreceptor stress (Fig. 10-1), including Stat3 (signal transducer and activator of transcription 3; Jiang et al., 2014), which we use as a representative example in Figure 10.
To control that these changes in the expression of genes regulating lipid metabolism do not merely represent a nonspecific transcriptional response of degenerating retinas, we repeated this analysis for retinas of mice bearing the P23H mutation in one copy of the rhodopsin gene (RhoP23H/WT mouse), a model commonly used to study retinal degeneration (Sakami et al., 2011). No changes in transcription levels of any of these genes were found, whereas Stat3 was robustly elevated (Fig. 10). This suggests that upregulation of the genes regulating lipid metabolism represents a specific response to the Mfsd2a knock-out.
Importantly, we did not detect any changes in the mRNA level of the Elovl4 gene (the elongation of very long chain fatty acids protein 4) responsible for VLCFA production. This suggests that the loss of phospholipids containing VLCFA is likely caused by a decreased availability of PUFAs, including DHA, serving among fatty acid precursors in the VLCFA biosynthetic pathway (Yu et al., 2012). Similarly, no changes were detected in the mRNA expression level of the Adipor1 gene (the Adiponectin Receptor 1), which was suggested to play a role in DHA retention in rods (Rice et al., 2015), and the Acsl6 gene coding Acyl-CoA synthetase 6 recently shown to contribute to DHA enrichment in the brain (Fernandez et al., 2018).
Single-cell recordings from Mfsd2a−/− rods
In the last set of experiments, we investigated the functional consequences of DHA loss on the sensitivity and recovery kinetics of rod phototransduction, the latter of which is known to be highly dependent on the diffusional encounter of the RGS9 complex in the disc membrane (Krispel et al., 2006; Gospe et al., 2011). The prevalent, although untested in intact cells, hypothesis is that high DHA content provides exceptional photoreceptor membrane liquidity (Treen et al., 1992; Niu et al., 2004), which supports a rapid rate of protein diffusion on these membranes' surface assuring the rapid diffusion needed for high amplification and timely deactivation. Mfsd2a−/− mice provide a unique opportunity to test this hypothesis directly.
Suction electrode recording from intact rods of Mfsd2a−/− and WT littermate mice revealed light-evoked responses that were highly similar (Fig. 11). Families of responses to flashes that ranged over 3 log units of intensity were qualitatively indistinguishable (Fig. 11A), although on average Mfsd2a−/− rods had smaller maximal dark currents (p = 1.1 × 10−8; Table 1) and captured a smaller fraction of incident photons (“effective collecting area”; p = 4.7 × 10−3; Table 1) than their WT counterparts, both of which may have been influenced by the younger average age of the knock-out mice (34 vs 50 d; Table 1; p = 9.3 × 10−6). Nevertheless, the amplitude and kinetics of the light-evoked responses of Mfsd2a−/− and WT rods were virtually indistinguishable (Table 1), especially when comparing littermates of nearly the exact same ages (Fig. 11B,C).
Discussion
In this study, we analyzed the retinal phenotype of a mouse lacking Mfsd2a, an LPC membrane transporter serving as a major contributor to the DHA delivery across both blood–retina and blood–brain barriers (Nguyen et al., 2014; Wong et al., 2016). Perhaps the most striking observation in our study is how relatively normal the functioning of Mfsd2a−/− retinas is, given that they lack ∼one-half of their normal DHA content and that brain pathology in these mice is much more severe (Nguyen et al., 2014; Chan et al., 2018). The progression of photoreceptor loss in Mfsd2a−/− retinas is relatively slow, particularly compared with many other models of inherited retinal degeneration, and surviving photoreceptors have normal morphology and generate essentially normal light responses.
The retinal phenotype of Mfsd2a−/− mice was previously investigated by (Wong et al., 2016). Whereas our findings are mostly in agreement, we present several unique results and report notable differences that need to be highlighted. Whereas Wong et al. (2016) conducted lipidome analysis using whole eyes, we performed this analysis using isolated retinas and purified outer segments. This allowed us to obtain a more quantitative account of the changes in lipid composition taking place at the retinal side of the blood-retina barrier. Accordingly, we observed an ∼threefold reduction in the overall DHA/AA ratio in Mfsd2a−/− retinas versus an ∼1.9-fold reduction previously measured in whole eyes. Likely for the same reason, we were able to reveal a robust ∼56% increase in the fraction of phospholipids containing MUFAs. In addition, our lipidomics methodology was suited to quantitatively detect lipids containing VLCFA moieties, which revealed an ∼57% reduction in VLCFAs (not analyzed in the previous study). The corresponding changes in outer segment phospholipids were even more pronounced: an ∼5.6-fold reduction in the DHA/AA ratio and ∼84% MUFA increase. Overall, these changes appear to reflect a compensatory response by the retina aiming to maintain a certain saturation level of the fatty acid moieties after losing a normal route of the DHA supply.
The DHA loss in Mfsd2a−/− retinas could be sufficiently explained by impaired LPC transport across the blood–retina barrier, and the reduction in VLCFAs could take place for the same reason because DHA serves as one of the VLCFA biosynthesis precursor (Yu et al., 2012). The other major changes in the retinal lipidome are more complex and could be attributed, at least in part, to a transcriptionally-driven attempt by the retina to rebuild its membrane composition. Indeed, the increased proportion of MUFAs could reflect upregulation of the desaturases, SCD1 and SCD2, which catalyze MUFA formation from saturated fatty acids. Consistently, we observed a small 7% reduction in the saturated fatty acid content, which almost reached statistical significance (p = 0.057). The upregulation of AA and complex changes in profiles of fatty acids with 2, 3 and 5 double bonds could be explained by upregulation of desaturases FADS1 and FADS2.
This type of transcriptional response triggered by deficiency of DHA and other PUFAs has been documented in other experimental systems and is thought to restore membrane fluidity and structure upon global or specific lipid deprivation (for review, see Hagen et al., 2010). Of particular interest is a recent biophysical study directly demonstrating the ability of cells to adjust their lipidome in response to PUFA challenge, resulting in the preservation of membrane packing and fluidity (Levental et al., 2018). This response is also consistent with a feedback mechanism in which LPC-DHA suppresses SREBF1 activity in developing brain, whereas Mfsd2a knock-out upregulates SREBF1 and its downstream targets following LPC-DHA deprivation (Chan et al., 2018).
Another unanswered question relates to the mechanisms accounting for the remaining retinal DHA in Mfsd2a−/− mice. These mechanisms may include developmental DHA accretion before the blood–retina barrier is established with subsequent DHA recycling from the RPE by interphotoreceptor retinoid-binding protein (Chen et al., 1993), or cross-barrier delivery of non-esterified DHA by fatty acid transporters (e.g., CD36 and FATP1; Mitchell and Hatch, 2011; Chouinard-Watkins et al., 2018).
Impaired LPC transport to the retina results in partial photoreceptor loss. However, neither gross morphological abnormalities in Mfsd2a−/− photoreceptors and RPE nor partial rhodopsin mislocalization from rod outer segments of young mice (Wong et al., 2016) were observed in our study. Regardless, both studies documented the invasion of mononuclear phagocytes into the subretinal space, which is a nonspecific hallmark of ongoing photoreceptor pathology (Rathnasamy et al., 2019).
Despite ongoing photoreceptor cell loss, the functional behavior of surviving cells is remarkably normal. This was first suggested by the ERG analysis by Wong et al. (2016). They observed somewhat reduced amplitudes of ERG a- and b-waves but noted that this reduction did not reach statistical significance and could be explained by partial photoreceptor loss. Considering the long history of presumed connection between DHA and normal vision, we were intrigued by this observation and performed single-cell recordings from Mfsd2a−/− rods. Unlike ERG, single-cell recordings are not affected by the fraction of surviving rods and allow precise monitoring of the entire photoresponse time course. The amplitudes and kinetics of light-responses by Mfsd2a−/− rods are strikingly similar to those from WT rods. This result was unexpected, considering that photoreceptor disc membranes harboring the phototransduction machinery are among the most DHA-enriched in our body (Fliesler and Anderson, 1983). Yet, such an unusually high DHA content does not appear to directly serve the needs of phototransduction and may relate to another enigmatic aspect of these cells' function not immediately related to light detection.
So, what causes degeneration of Mfsd2a−/− rods? While obtaining a definitive answer is another subject of future investigation, we can propose several not mutually exclusive mechanisms. A rather trivial explanation is that these cells suffer from a general lipid insufficiency caused by impaired LPC transport. Such insufficiency would particularly affect photoreceptors since they replenish the entire complement of outer segment membranes every 10–12 d (Young, 1967). Consistent with this explanation, the knock-outs of LPCAT1 or LPAAT3, the enzymes regulating the biosynthesis of retinal phospholipids (including but not limited to those containing DHA), cause even more severe photoreceptor degeneration than the Mfsd2a knock-out (Friedman et al., 2010; Shindou et al., 2017). Another potential explanation is that reduced lipid supply to the Mfsd2a−/− retina causes energy deprivation since the retina, and particularly photoreceptors, can oxidize lipid through fatty acid β-oxidation to sustain its energy needs (Joyal et al., 2016). It is also possible that pathology arises from significant changes in fatty acids (MUFAs, DHA, AA and VLCFAs), which act as signaling molecules, either themselves or through their metabolites (Igal, 2016; Rodríguez-Cruz and Serna, 2017; Bazan, 2018). For example, the shift in the DHA/AA ratio may be proinflammatory (SanGiovanni and Chew, 2005). Of particular interest is the observed reduction in VLCFAs, as evolving evidence suggests a complex involvement of VLCFAs in supporting retinal structure and function, particularly in the context of hereditary retinal degenerations, including Stargardt-3 macular degeneration (Hopiavuori et al., 2019). Finally, a recent study reported that Mfsd2a suppresses transcellular transport in endothelial cells of the retinal vasculature (Chow and Gu, 2017). One can imagine that this may be detrimental for the retina, although it is not clear how exactly increased transcytosis in retinal vasculature may lead to photoreceptor degeneration.
Another important conclusion from our study is that the consequences of DHA deficit caused by the Mfsd2a knock-out are different from those caused by dietary DHA restriction. Dietary restriction in mice and rats could lead to a higher degree of DHA deprivation of the retina associated with reductions in ERG responses, outer segment shortening and alteration in the rate of photoreceptor outer segment renewal without notable photoreceptor loss (Benolken et al., 1973; Wheeler et al., 1975; Senapati et al., 2018). However, this diet leads to a major substitution of DHA with docosapentaenoic acid (22:5; Anderson et al., 1974; Niu et al., 2004), which does not happen in Mfsd2a−/− retinas. Therefore, any direct comparison of retinal DHA deprivation through dietary restriction versus Mfsd2a knock-out may be deceiving.
In a similar context, the study by Rice et al. (2015) analyzed the retinal phenotype of the adiponectin receptor 1 (Adipor1) knock-out and found an ∼50% reduction in the DHA content associated with a more severe photoreceptor degeneration than in Mfsd2a knock-out. The authors concluded that the activity of ADIPOR1 is critical for DHA retention in the retina and RPE.
In summary, our results, combined with other observations discussed here, argue that maintaining an efficient phospholipid supply to the retina is more critical for supporting photoreceptor survival than their ability to produce light responses. Furthermore, the degree of photoreceptor pathology does not appear to directly correlate with the degree of retinal DHA deprivation and may be driven by other concurrent changes in the retinal lipidome. Therefore, the challenges of future studies are to pinpoint the precise place of DHA, other unsaturated fatty acids and their metabolites in the complex regulatory networks defining the compositions of retinal and photoreceptor lipidomes and to identify potential targets within these networks for improving photoreceptor health in congenital and age-related diseases of the retina.
Footnotes
This work was supported by the Duke NUS Pilot Award “The role of DHA in supporting retinal function” (V.Y.A.); National Institutes of Health Grants EY022959 (V.Y.A.), EY05722 (V.Y.A.), EY016459 (J.D.A.), EY014047 (M.E.B.), and EY030043(E.S.L.); a Grant from the Knights Templar Eye Foundation (E.S.L.); The University of Florida startup funds (E.S.L.); unrestricted Grants from Research to Prevent Blindness to the Departments of Ophthalmology of University of Florida and Duke University; the German Research Foundation TRR83 Grant project A17 (A.S.); and Federal Ministry of Education and Research LIFS Grant (A.S.). We thank Dr. David L Silver (Duke NUS) for many helpful discussions and critical reading of the paper.
The authors declare no competing financial interests
- Correspondence should be addressed to Ekaterina S. Lobanova at elobanova{at}ufl.edu or Vadim Y. Arshavsky at vadim.arshavsky{at}duke.edu