Genetic Risk of Reticular Pseudodrusen in Age-Related Macular Degeneration: HTRA1 /lncRNA BX842242.1 dominates, with no evidence for Complement Cascade involvement.

Age-related macular degeneration (AMD) is a multifactorial retinal disease with a large genetic risk contribution. Reticular pseudodrusen (RPD) is a sub-phenotype of AMD with a high risk of progression to late vision threatening AMD. In a genome-wide association study of 2,165 AMD+/RPD+ and 4,181 AMD+/RPD-compared to 7,660 control participants, both chromosomes 1 ( CFH ) and 10 ( ARMS2/HTRA1 ) major AMD risk loci were reidentified. However association was only detected for the chromosome 10 locus when comparing AMD+/RPD+ to AMD+/RPD-cases. The chromosome 1 locus was notably absent. The chromosome 10 RPD risk region contains a long non-coding RNA (ENSG00000285955/BX842242.1) which colocalizes with genetic markers of retinal thickness. BX842242.1 has a strong retinal eQTL signal, pinpointing the parafoveal photoreceptor outer segment layer. Whole genome sequencing of phenotypically extreme RPD cases identified even stronger enrichment for the chromosome 10 risk genotype.

Deep Learning-Based Detection of Reticular Pseudodrusen in Age-Related Macular Degeneration on Optical Coherence Tomography.

Reticular pseudodrusen (RPD) signify a critical phenotype driving vision loss in age-related macular degeneration (AMD). Their detection is paramount in the clinical management of those with AMD, yet they remain challenging to reliably identify. We thus developed a deep learning (DL) model to segment RPD from 9,800 optical coherence tomography B-scans, and this model produced RPD segmentations that had higher agreement with four retinal specialists (Dice similarity coefficient [DSC]=0·76 [95% confidence interval [CI] 0·71-0·81]) than the agreement amongst the specialists (DSC=0·68, 95% CI=0·63-0·73; p <0·001). In five external test datasets consisting of 1,017 eyes from 812 individuals, the DL model detected RPD with a similar level of performance as two retinal specialists (area-under-the-curve of 0·94 [95% CI=0·92-0·97], 0·95 [95% CI=0·92-0·97] and 0·96 [95% CI=0·94-0·98] respectively; p ≥0·32). This DL model enables the automatic detection and quantification of RPD with expert-level performance, which we have made publicly available.

Serine and glycine physiology reversibly modulate retinal and peripheral nerve function.

Metabolic homeostasis is maintained by redundant pathways to ensure adequate nutrient supply during fasting and other stresses. These pathways are regulated locally in tissues and systemically via the liver, kidney, and circulation. Here, we characterize how serine, glycine, and one-carbon (SGOC) metabolism fluxes across the eye, liver, and kidney sustain retinal amino acid levels and function. Individuals with macular telangiectasia (MacTel), an age-related retinal disease with reduced circulating serine and glycine, carrying deleterious alleles in SGOC metabolic enzymes exhibit an exaggerated reduction in circulating serine. A Phgdh+/- mouse model of this haploinsufficiency experiences accelerated retinal defects upon dietary serine/glycine restriction, highlighting how otherwise silent haploinsufficiencies can impact retinal health. We demonstrate that serine-associated retinopathy and peripheral neuropathy are reversible, as both are restored in mice upon serine supplementation. These data provide molecular insights into the genetic and metabolic drivers of neuro-retinal dysfunction while highlighting therapeutic opportunities to ameliorate this pathogenesis.

Decreased Circulating Very Small Low-Density Lipoprotein is Likely Causal for Age-Related Macular Degeneration.

Objective: Abnormal changes in metabolite levels in serum or plasma have been highlighted in several studies in age-related macular degeneration (AMD), the leading cause of irreversible vision loss. Specific changes in lipid profiles are associated with an increased risk of AMD. Metabolites could thus be used to investigate AMD disease mechanisms or incorporated into AMD risk prediction models. However, whether particular metabolites causally affect the disease has yet to be established. Design: A 3-tiered analysis of blood metabolites in the United Kingdom (UK) Biobank cohort to identify metabolites that differ in AMD patients with evidence for a putatively causal role in AMD. Participants: A total of 72 376 donors from the UK Biobank cohort including participants with AMD (N = 1353) and non-AMD controls (N = 71 023). Methods: We analyzed 325 directly measured or derived blood metabolites from the UK Biobank for 72 376 donors to identify AMD-associated metabolites. Genome-wide association studies for 325 metabolites in 98 316 European participants from the UK Biobank were performed. The causal effects of these metabolites in AMD were tested using a 2-sample Mendelian randomization approach. The predictive value of these measurements together with sex and age was assessed by developing a machine learning classifier. Main Outcome Measures: Evaluating metabolic biomarkers associated with AMD susceptibility and investigating their potential causal contribution to the development of the disease. Results: This study noted age to be the prominent risk factor associated with AMD development. While accounting for age and sex, we identified 84 metabolic markers as significantly (false discovery rate-adjusted P value < 0.05) associated with AMD. Lipoprotein subclasses comprised the majority of the AMD-associated metabolites (39%) followed by several lipoprotein to lipid ratios. Nineteen metabolites showed a likely causative role in AMD etiology. Of these, 6 lipoproteins contain very small, very low-density lipoprotein (VLDL), and phospholipids to total lipid ratio in medium VLDL. Based on this we postulate that depletion of circulating very small VLDLs is likely causal for AMD. The risk prediction model constructed from the metabolites, age and sex, identified age as the primary predictive factor with a much smaller contribution by metabolites to AMD risk prediction. Conclusions: This study underscores the pronounced role of lipids in AMD susceptibility and the likely causal contribution of particular subclasses of lipoproteins to AMD. Our study provides valuable insights into the metabopathological mechanisms of AMD disease development and progression.

Mouse microglia express unique miRNA-mRNA networks to facilitate age-specific functions in the developing central nervous system.

Microglia regulate multiple processes in the central nervous system, exhibiting a considerable level of cellular plasticity which is facilitated by an equally dynamic transcriptional environment. While many gene networks that regulate microglial functions have been characterised, the influence of epigenetic regulators such as small non-coding microRNAs (miRNAs) is less well defined. We have sequenced the miRNAome and mRNAome of mouse microglia during brain development and adult homeostasis, identifying unique profiles of known and novel miRNAs. Microglia express both a consistently enriched miRNA signature as well as temporally distinctive subsets of miRNAs. We generated robust miRNA-mRNA networks related to fundamental developmental processes, in addition to networks associated with immune function and dysregulated disease states. There was no apparent influence of sex on miRNA expression. This study reveals a unique developmental trajectory of miRNA expression in microglia during critical stages of CNS development, establishing miRNAs as important modulators of microglial phenotype.

Transcriptomic analysis of albendazole resistance in human diarrheal parasite Giardia duodenalis.

Benzimidazole-2-carbamates (BZ, e.g., albendazole; ALB), which bind ß-tubulin to disrupt microtubule polymerization, are one of two primary compound classes used to treat giardiasis. In most parasitic nematodes and fungi, BZ-resistance is caused by ß-tubulin mutations and its molecular mode of action (MOA) is well studied. In contrast, in Giardia duodenalis BZ MOA or resistance is less well understood, may involve target-specific and broader impacts including cellular damage and oxidative stress, and its underlying cause is not clearly determined. Previously, we identified acquisition of a single nucleotide polymorphism, E198K, in ß-tubulin in ALB-resistant (ALB-R) G. duodenalis WB-1B relative to ALB-sensitive (ALB-S) parental controls. E198K is linked to BZ-resistance in fungi and its allelic frequency correlated with the magnitude of BZ-resistance in G. duodenalis WB-1B. Here, we undertook detailed transcriptomic comparisons of these ALB-S and ALB-R G. duodenalis WB-1B cultures. The primary transcriptional changes with ALB-R in G. duodenalis WB-1B indicated increased protein degradation and turnover, and up-regulation of tubulin, and related genes, associated with the adhesive disc and basal bodies. These findings are consistent with previous observations noting focused disintegration of the disc and associated structures in Giardia duodenalis upon ALB exposure. We also saw transcriptional changes with ALB-R in G. duodenalis WB-1B consistent with prior observations of a shift from glycolysis to arginine metabolism for ATP production and possible changes to aspects of the vesicular trafficking system that require further investigation. Finally, we saw mixed transcriptional changes associated with DNA repair and oxidative stress responses in the G. duodenalis WB-1B line. These changes may be indicative of a role for H2O2 degradation in ALB-R, as has been observed in other G. duodenalis cell cultures. However, they were below the transcriptional fold-change threshold (log2FC > 1) typically employed in transcriptomic analyses and appear to be contradicted in ALB-R G. duodenalis WB-1B by down-regulation of the NAD scavenging and conversion pathways required to support these stress pathways and up-regulation of many highly oxidation sensitive iron-sulphur (FeS) cluster based metabolic enzymes.

iPSC-derived retinal pigmented epithelial cells from patients with macular telangiectasia show decreased mitochondrial function.

Patient-derived induced pluripotent stem cells (iPSCs) provide a powerful tool for identifying cellular and molecular mechanisms of disease. Macular telangiectasia type 2 (MacTel) is a rare, late-onset degenerative retinal disease with an extremely heterogeneous genetic architecture, lending itself to the use of iPSCs. Whole-exome sequencing screens and pedigree analyses have identified rare causative mutations that account for less than 5% of cases. Metabolomic surveys of patient populations and GWAS have linked MacTel to decreased circulating levels of serine and elevated levels of neurotoxic 1-deoxysphingolipids (1-dSLs). However, retina-specific, disease-contributing factors have yet to be identified. Here, we used iPSC-differentiated retinal pigmented epithelial (iRPE) cells derived from donors with or without MacTel to screen for novel cell-intrinsic pathological mechanisms. We show that MacTel iRPE cells mimicked the low serine levels observed in serum from patients with MacTel. Through RNA-Seq and gene set enrichment pathway analysis, we determined that MacTel iRPE cells are enriched in cellular stress pathways and dysregulation of central carbon metabolism. Using respirometry and mitochondrial stress testing, we functionally validated that MacTel iRPE cells had a reduction in mitochondrial function that was independent of defects in serine biosynthesis and 1-dSL accumulation. Thus, we identified phenotypes that may constitute alternative disease mechanisms beyond the known serine/sphingolipid pathway.

Divergent amino acid and sphingolipid metabolism in patients with inherited neuro-retinal disease.

OBJECTIVES: The non-essential amino acids serine, glycine, and alanine, as well as diverse sphingolipid species, are implicated in inherited neuro-retinal disorders and are metabolically linked by serine palmitoyltransferase (SPT), a key enzyme in membrane lipid biogenesis. To gain insight into the pathophysiological mechanisms linking these pathways to neuro-retinal diseases we compared patients diagnosed with two metabolically intertwined diseases: macular telangiectasia type II (MacTel), hereditary sensory autonomic neuropathy type 1 (HSAN1), or both. METHODS: We performed targeted metabolomic analyses of amino acids and broad sphingolipids in sera from a cohort of MacTel (205), HSAN1 (25) and Control (151) participants. RESULTS: MacTel patients exhibited broad alterations of amino acids, including changes in serine, glycine, alanine, glutamate, and branched-chain amino acids reminiscent of diabetes. MacTel patients had elevated 1-deoxysphingolipids but reduced levels of complex sphingolipids in circulation. A mouse model of retinopathy indicates dietary serine and glycine restriction can drive this depletion in complex sphingolipids. HSAN1 patients exhibited elevated serine, lower alanine, and a reduction in canonical ceramides and sphingomyelins compared to controls. Those patients diagnosed with both HSAN1 and MacTel showed the most significant decrease in circulating sphingomyelins. CONCLUSIONS: These results highlight metabolic distinctions between MacTel and HSAN1, emphasize the importance of membrane lipids in the progression of MacTel, and suggest distinct therapeutic approaches for these two neurodegenerative diseases.

Spatial distribution of metabolites in the retina and its relevance to studies of metabolic retinal disorders.

INTRODUCTION: The primate retina has evolved regional specialisations for specific visual functions. The macula is specialised towards high acuity vision and is an area that contains an increased density of cone photoreceptors and signal processing neurons. Different regions in the retina display unique susceptibility to pathology, with many retinal diseases primarily affecting the macula. OBJECTIVES: To better understand the properties of different retinal areas we studied the differential distribution of metabolites across the retina. METHODS: We conducted an untargeted metabolomics analysis on full-thickness punches from three different regions (macula, temporal peri-macula and periphery) of healthy primate retina. RESULTS: Nearly half of all metabolites identified showed differential abundance in at least one comparison between the three regions. Furthermore, mapping metabolomics results from macula-specific eye diseases onto our region-specific metabolite distributions revealed differential abundance defining systemic metabolic dysregulations that were region specific. CONCLUSIONS: The unique metabolic phenotype of different retinal regions is likely due to the differential distribution of different cell types in these regions reflecting the specific metabolic requirements of each cell type. Our results may help to better understand the pathobiology of retinal diseases with region specificity.

VIVID: A Web Application for Variant Interpretation and Visualization in Multi-dimensional Analyses.

Large-scale comparative genomics- and population genetic studies generate enormous amounts of polymorphism data in the form of DNA variants. Ultimately, the goal of many of these studies is to associate genetic variants to phenotypes or fitness. We introduce VIVID, an interactive, user-friendly web application that integrates a wide range of approaches for encoding genotypic to phenotypic information in any organism or disease, from an individual or population, in three-dimensional (3D) space. It allows mutation mapping and annotation, calculation of interactions and conservation scores, prediction of harmful effects, analysis of diversity and selection, and 3D visualization of genotypic information encoded in Variant Call Format on AlphaFold2 protein models. VIVID enables the rapid assessment of genes of interest in the study of adaptive evolution and the genetic load, and it helps prioritizing targets for experimental validation. We demonstrate the utility of VIVID by exploring the evolutionary genetics of the parasitic protist Plasmodium falciparum, revealing geographic variation in the signature of balancing selection in potential targets of functional antibodies.

Neonatal BCG vaccination is associated with a long-term DNA methylation signature in circulating monocytes.

Trained immunity describes the capacity of innate immune cells to develop heterologous memory in response to certain exogenous exposures. This phenomenon mediates, at least in part, the beneficial off-target effects of the BCG vaccine. Using an in vitro model of trained immunity, we show that BCG exposure induces a persistent change in active histone modifications, DNA methylation, transcription, and adenosine-to-inosine RNA modification in human monocytes. By profiling DNA methylation of circulating monocytes from infants in the MIS BAIR clinical trial, we identify a BCG-associated DNA methylation signature that persisted more than 12 months after neonatal BCG vaccination. Genes associated with this epigenetic signature are involved in viral response pathways, consistent with the reported off-target protection against viral infections in neonates, adults, and the elderly. Our findings indicate that the off-target effects of BCG in infants are accompanied by epigenetic remodeling of circulating monocytes that lasts more than 1 year.

Heterozygous PNPT1 Variants Cause Spinocerebellar Ataxia Type 25.

OBJECTIVE: Dominant spinocerebellar ataxias (SCA) are characterized by genetic heterogeneity. Some mapped and named loci remain without a causal gene identified. Here we applied next generation sequencing (NGS) to uncover the genetic etiology of the SCA25 locus. METHODS: Whole-exome and whole-genome sequencing were performed in families linked to SCA25, including the French family in which the SCA25 locus was originally mapped. Whole exome sequence data were interrogated in a cohort of 796 ataxia patients of unknown etiology. RESULTS: The SCA25 phenotype spans a slowly evolving sensory and cerebellar ataxia, in most cases attributed to ganglionopathy. A pathogenic variant causing exon skipping was identified in the gene encoding Polyribonucleotide Nucleotidyltransferase PNPase 1 (PNPT1) located in the SCA25 linkage interval. A second splice variant in PNPT1 was detected in a large Australian family with a dominant ataxia also mapping to SCA25. An additional nonsense variant was detected in an unrelated individual with ataxia. Both nonsense and splice heterozygous variants result in premature stop codons, all located in the S1-domain of PNPase. In addition, an elevated type I interferon response was observed in blood from all affected heterozygous carriers tested. PNPase notably prevents the abnormal accumulation of double-stranded mtRNAs in the mitochondria and leakage into the cytoplasm, associated with triggering a type I interferon response. INTERPRETATION: This study identifies PNPT1 as a new SCA gene, responsible for SCA25, and highlights biological links between alterations of mtRNA trafficking, interferonopathies and ataxia. ANN NEUROL 2022;92:122-137.

Genetics of reticular pseudodrusen in age-related macular degeneration.

Reticular pseudodrusen (RPD) are subretinal deposits that, when observed with age-related macular degeneration (AMD), form a distinct phenotype, often associated with late-stage disease. To date, RPD genetic risk associations overlap six well-established AMD-risk regions. Determining RPD-specific underlying genetic causes by using adequate imaging methods should improve our understanding of the pathophysiology of RPD.

Long-read assembly and comparative evidence-based reanalysis of Cryptosporidium genome sequences reveal expanded transporter repertoire and duplication of entire chromosome ends including subtelomeric regions.

Cryptosporidiosis is a leading cause of waterborne diarrheal disease globally and an important contributor to mortality in infants and the immunosuppressed. Despite its importance, the Cryptosporidium community has only had access to a good, but incomplete, Cryptosporidium parvum IOWA reference genome sequence. Incomplete reference sequences hamper annotation, experimental design, and interpretation. We have generated a new C. parvum IOWA genome assembly supported by Pacific Biosciences (PacBio) and Oxford Nanopore long-read technologies and a new comparative and consistent genome annotation for three closely related species: C. parvum, Cryptosporidium hominis, and Cryptosporidium tyzzeri We made 1926 C. parvum annotation updates based on experimental evidence. They include new transporters, ncRNAs, introns, and altered gene structures. The new assembly and annotation revealed a complete Dnmt2 methylase ortholog. Comparative annotation between C. parvum, C. hominis, and C. tyzzeri revealed that most "missing" orthologs are found, suggesting that the biological differences between the species must result from gene copy number variation, differences in gene regulation, and single-nucleotide variants (SNVs). Using the new assembly and annotation as reference, 190 genes are identified as evolving under positive selection, including many not detected previously. The new C. parvum IOWA reference genome assembly is larger, gap free, and lacks ambiguous bases. This chromosomal assembly recovers all 16 chromosome ends, 13 of which are contiguously assembled. The three remaining chromosome ends are provisionally placed. These ends represent duplication of entire chromosome ends including subtelomeric regions revealing a new level of genome plasticity that will both inform and impact future research.

A survey of RNA editing at single-cell resolution links interneurons to schizophrenia and autism.

Conversion of adenosine to inosine in RNA by ADAR enzymes, termed "RNA editing," is essential for healthy brain development. Editing is dysregulated in neuropsychiatric diseases, but has not yet been investigated at scale at the level of individual neurons. We quantified RNA editing sites in nuclear transcriptomes of 3055 neurons from six cortical regions of a neurotypical female donor, and found 41,930 sites present in at least ten nuclei. Most sites were located within Alu repeats in introns or 3' UTRs, and approximately 80% were cataloged in public RNA editing databases. We identified 9285 putative novel editing sites, 29% of which were also detectable in unrelated donors. Intersection with results from bulk RNA-seq studies provided cell-type and spatial context for 1730 sites that are differentially edited in schizophrenic brain donors, and 910 such sites in autistic donors. Autism-related genes were also enriched with editing sites predicted to modify RNA structure. Inhibitory neurons showed higher overall transcriptome editing than excitatory neurons, and the highest editing rates were observed in the frontal cortex. We used generalized linear models to identify differentially edited sites and genes between cell types. Twenty nine genes were preferentially edited in excitatory neurons, and 43 genes were edited more heavily in inhibitory neurons, including RBFOX1, its target genes, and genes in the autism-associated Prader-Willi locus (15q11). The abundance of SNORD115/116 genes from locus 15q11 was positively associated with editing activity across the transcriptome. We contend that insufficient editing of autism-related genes in inhibitory neurons may contribute to the specific perturbation of those cells in autism.

Multimodal regulation of encystation in Giardia duodenalis revealed by deep proteomics.

Cyst formation in the parasitic protist Giardia duodenalis is critical to its transmission. Existing proteomic data quantifies only 17% of coding genes transcribed during encystation and does not cover the complete process from trophozoite to mature cyst. Using high-resolution mass spectrometry, we have quantified proteomic changes across encystation and compared this with published transcriptomic data. We reproducibly identified 3863 (64.5% of Giardia proteins) and quantified 3382 proteins (56.5% of Giardia proteins) over standard trophozoite growth (TY), during low-bile encystation priming (LB), 16 h into encystation (EC), and at cyst maturation (C). This work provides the first known expanded observation of encystation at the proteomic level and triples the coverage of previous encystation proteomes. One-third (1169 proteins) of the quantified proteome is differentially expressed in the mature cyst relative to the trophozoite, including proteasomal machinery, metabolic pathways, and secretory proteins. Changes in lipid metabolism indicated a shift in lipid species dependency during encystation. Consistent with this, we identified the first, putative lipid transporters in this species, representing the steroidogenic acute regulatory protein-related lipid transfer (StARkin), oxysterol binding protein related protein (ORP/Osh) and glycosphingolipid transfer protein (GLTP) families, and follow their differential expression over cyst formation. Lastly, we undertook correlation analyses of the transcriptome and proteome of trophozoites and cysts, and found evidence of post-transcriptional regulation of key protein classes (RNA binding proteins) and stage-specific genes (encystation markers) implicating translation-repression in encystation. We provide the most extensive proteomic analysis of encystation in Giardia to date and the first known exploration across its complete duration. This work identifies encystation as highly coordinated, involving major changes in proteostasis, metabolism and membrane dynamics, and indicates a potential role for post-transcriptional regulation, mediated through RNA-binding proteins. Together our work provides a valuable resource for Giardia research and the development of transmission-blocking anti-giardials.

Identification of genetic factors influencing metabolic dysregulation and retinal support for MacTel, a retinal disorder.

Macular Telangiectasia Type 2 (MacTel) is a rare degenerative retinal disease with complex genetic architecture. We performed a genome-wide association study on 1,067 MacTel patients and 3,799 controls, which identified eight novel genome-wide significant loci (p < 5 × 10-8), and confirmed all three previously reported loci. Using MAGMA, eQTL and transcriptome-wide association analysis, we prioritised 48 genes implicated in serine-glycine biosynthesis, metabolite transport, and retinal vasculature and thickness. Mendelian randomization indicated a likely causative role of serine (FDR = 3.9 × 10-47) and glycine depletion (FDR = 0.006) as well as alanine abundance (FDR = 0.009). Polygenic risk scoring achieved an accuracy of 0.74 and was associated in UKBiobank with retinal damage (p = 0.009). This represents the largest genetic study on MacTel to date and further highlights genetically-induced systemic and tissue-specific metabolic dysregulation in MacTel patients, which impinges on retinal health.

Genetic disruption of serine biosynthesis is a key driver of macular telangiectasia type 2 aetiology and progression.

BACKGROUND: Macular telangiectasia type 2 (MacTel) is a rare, heritable and largely untreatable retinal disorder, often comorbid with diabetes. Genetic risk loci subtend retinal vascular calibre and glycine/serine/threonine metabolism genes. Serine deficiency may contribute to MacTel via neurotoxic deoxysphingolipid production; however, an independent vascular contribution is also suspected. Here, we use statistical genetics to dissect the causal mechanisms underpinning this complex disease. METHODS: We integrated genetic markers for MacTel, vascular and metabolic traits, and applied Mendelian randomisation and conditional and interaction genome-wide association analyses to discover the causal contributors to both disease and spatial retinal imaging sub-phenotypes. RESULTS: Genetically induced serine deficiency is the primary causal metabolic driver of disease occurrence and progression, with a lesser, but significant, causal contribution of type 2 diabetes genetic risk. Conversely, glycine, threonine and retinal vascular traits are unlikely to be causal for MacTel. Conditional regression analysis identified three novel disease loci independent of endogenous serine biosynthetic capacity. By aggregating spatial retinal phenotypes into endophenotypes, we demonstrate that SNPs constituting independent risk loci act via related endophenotypes. CONCLUSIONS: Follow-up studies after GWAS integrating publicly available data with deep phenotyping are still rare. Here, we describe such analysis, where we integrated retinal imaging data with MacTel and other traits genomics data to identify biochemical mechanisms likely causing this disorder. Our findings will aid in early diagnosis and accurate prognosis of MacTel and improve prospects for effective therapeutic intervention. Our integrative genetics approach also serves as a useful template for post-GWAS analyses in other disorders.

Systemic lipid dysregulation is a risk factor for macular neurodegenerative disease.

Macular Telangiectasia type 2 (MacTel) is an uncommon bilateral retinal disease, in which glial cell and photoreceptor degeneration leads to central vision loss. The causative disease mechanism is largely unknown, and no treatment is currently available. A previous study found variants in genes associated with glycine-serine metabolism (PSPH, PHGDH and CPS1) to be associated with MacTel, and showed low levels of glycine and serine in the serum of MacTel patients. Recently, a causative role of deoxysphingolipids in MacTel disease has been established. However, little is known about possible other metabolic dysregulation. Here we used a global metabolomics platform in a case-control study to comprehensively profile serum from 60 MacTel patients and 58 controls. Analysis of the data, using innovative computational approaches, revealed a detailed, disease-associated metabolic profile with broad changes in multiple metabolic pathways. This included alterations in the levels of several metabolites that are directly or indirectly linked to glycine-serine metabolism, further validating our previous genetic findings. We also found changes unrelated to PSPH, PHGDH and CPS1 activity. Most pronounced, levels of several lipid groups were altered, with increased phosphatidylethanolamines being the most affected lipid group. Assessing correlations between different metabolites across our samples revealed putative functional connections. Correlations between phosphatidylethanolamines and sphingomyelin, and glycine-serine and sphingomyelin, observed in controls, were reduced in MacTel patients, suggesting metabolic re-wiring of sphingomyelin metabolism in MacTel patients. Our findings provide novel insights into metabolic changes associated with MacTel and implicate altered lipid metabolism as a contributor to this retinal neurodegenerative disease.