Impaired glucose-1,6-biphosphate production due to bi-allelic PGM2L1 mutations is associated with a neurodevelopmental disorder

Molecular and Clinical Characterization of a Founder Mutation Causing G6PC3 Deficiency

G6PC3 deficiency is a monogenic immunometabolic disorder that causes severe congenital neutropenia type 4. Patients display heterogeneous extra-hematological manifestations, contributing to delayed diagnosis. Here, we investigated the origin and functional consequence of the G6PC3 c.210delC variant found in patients of Mexican descent. Based on the shared haplotypes amongst mutation carriers, we estimated that this variant originated from a founder effect in a common ancestor. Furthermore, by ancestry analysis, we concluded that it appeared in the indigenous Mexican population. At the protein level, we showed that this frameshift mutation leads to an aberrant protein expression in overexpression and patient-derived Epstein-Barr Virus-immortalized B (EBV-B) cells. The neutropenia observed in G6PC3-deficient patients is driven by the intracellular accumulation of the metabolite 1,5-anhydroglucitol-6-phosphate (1,5-AG6P) that inhibits glycolysis. We characterized how the c.210delC variant impacts glycolysis by performing extracellular flux assays on patient-derived EBV-B cells. When treated with 1,5-anhydroglucitol (1,5-AG), the precursor to 1,5-AG6P, patient cells exhibited markedly reduced engagement of glycolysis. Finally, we compared the clinical presentation of patients with the mutation c.210delC and all other G6PC3-deficient patients reported in the literature, and we found that the c.210delC carriers display all prominent clinical features observed in prior patients. In conclusion, G6PC3 c.210delC is a loss-of-function mutation that arose from a founder effect in the indigenous Mexican population. These findings may facilitate the diagnosis of additional patients in this geographical area. Moreover, the in vitro 1,5-AG-dependent functional assay used in our study could be employed to assess the pathogenicity of additional G6PC3 variants.

Protein kinase N1 deficiency results in upregulation of cerebral energy metabolism and is highly protective in in vivo and in vitro stroke models

BACKGROUND AND AIM

We recently identified protein kinase N1 (PKN1) as a master regulator of brain development. However, its function in the adult brain has not been clearly established. In this study, we assessed the cerebral energetic phenotype of wildtype (WT) and global Pkn1 knockout (Pkn1-/-) animals under physiological and pathophysiological conditions.

METHODS

Cerebral energy metabolism was analyzed by 13C6-glucose tracing in vivo and real time seahorse analysis of extracellular acidification rates as well as mitochondrial oxygen consumption rates (OCR) of brain slice punches in vitro. Isolated WT and Pkn1-/- brain mitochondria were tested for differences in OCR with different substrates. Metabolite levels were determined by mass spectrometric analysis in brain slices under control and energetic stress conditions, induced by oxygen-glucose deprivation and reperfusion, an in vitro model of ischemic stroke. Differences in enzyme activities were assessed by enzymatic assays, western blotting and bulk RNA sequencing. A middle cerebral artery occlusion stroke model was used to analyze lesion volumes and functional recovery in WT and Pkn1-/- mice.

RESULTS

Pkn1 deficiency resulted in a remarkable upregulation of cerebral energy metabolism, in vivo and in vitro. This was due to two separate mechanisms involving an enhanced glycolytic flux and higher pyruvate-induced mitochondrial OCR. Mechanistically we show that Pkn1-/- brain tissue exhibits an increased activity of the glycolysis rate-limiting enzyme phosphofructokinase. Additionally, glucose-1,6-bisphosphate levels, a metabolite that increases mitochondrial pyruvate uptake, were elevated upon Pkn1 deficiency. Consequently, Pkn1-/- brain slices had more ATP and a greater accumulation of ATP degradation metabolites during energetic stress. This translated into increased phosphorylation and activity of adenosine monophosphate (AMP)-activated protein kinase (AMPK) during in vitro stroke. Accordingly, Pkn1-/- brain slices showed a post-ischemic transcriptional upregulation of energy metabolism pathways and Pkn1 deficiency was strongly protective in in vitro and in vivo stroke models. While inhibition of mitochondrial pyruvate uptake only moderately affected the protective phenotype, inhibition of AMPK in Pkn1-/- slices increased post-ischemic cell death in vitro.

CONCLUSION

This is the first study to comprehensively demonstrate an essential and unique role of PKN1 in cerebral energy metabolism, regulating glycolysis and mitochondrial pyruvate-induced respiration. We further uncovered a highly protective phenotype of Pkn1 deficiency in both, in vitro and in vivo stroke models, validating inhibition of PKN1 as a promising new therapeutic target for the development of novel stroke therapies.

The identification of key genes and pathways in polycystic ovary syndrome by bioinformatics analysis of next-generation sequencing data

Background

Polycystic ovary syndrome (PCOS) is a reproductive endocrine disorder. The specific molecular mechanism of PCOS remains unclear. The aim of this study was to apply a bioinformatics approach to reveal related pathways or genes involved in the development of PCOS.

Methods

The next-generation sequencing (NGS) dataset GSE199225 was downloaded from the gene expression omnibus (GEO) database and NGS dataset analyzed is obtained from in vitro culture of PCOS patients’ muscle cells and muscle cells of healthy lean control women. Differentially expressed gene (DEG) analysis was performed using DESeq2. The g:Profiler was utilized to analyze the gene ontology (GO) and REACTOME pathways of the differentially expressed genes. A protein–protein interaction (PPI) network was constructed and module analysis was performed using HiPPIE and cytoscape. The miRNA-hub gene regulatory network and TF-hub gene regulatory network were constructed. The hub genes were validated by using receiver operating characteristic (ROC) curve analysis.

Results

We have identified 957 DEG in total, including 478 upregulated genes and 479 downregulated gene. GO terms and REACTOME pathways illustrated that DEG were significantly enriched in regulation of molecular function, developmental process, interferon signaling and platelet activation, signaling, and aggregation. The top 5 upregulated hub genes including HSPA5, PLK1, RIN3, DBN1, and CCDC85B and top 5 downregulated hub genes including DISC1, AR, MTUS2, LYN, and TCF4 might be associated with PCOS. The hub gens of HSPA5 and KMT2A, together with corresponding predicted miRNAs (e.g., hsa-mir-34b-5p and hsa-mir-378a-5p), and HSPA5 and TCF4 together with corresponding predicted TF (e.g., RCOR3 and TEAD4) were found to be significantly correlated with PCOS.

Conclusions

These study uses of bioinformatics analysis of NGS data to obtain hub genes and key signaling pathways related to PCOS and its associated complications. Also provides novel ideas for finding biomarkers and treatment methods for PCOS and its associated complications.

Absence of the dolichol synthesis gene DHRSX leads to N-glycosylation defects in Lec5 and Lec9 Chinese hamster ovary cells

Glycosylation-deficient Chinese hamster ovary cell lines have been instrumental in the discovery of N-glycosylation machinery. Yet, the molecular causes of the glycosylation defects in the Lec5 and Lec9 mutants have been elusive, even though for both cell lines a defect in dolichol formation from polyprenol was previously established. We recently found that dolichol synthesis from polyprenol occurs in three steps consisting of the conversion of polyprenol to polyprenal by DHRSX, the reduction of polyprenal to dolichal by SRD5A3, and the reduction of dolichal to dolichol, again by DHRSX. This led us to investigate defective dolichol synthesis in Lec5 and Lec9 cells. Both cell lines showed increased levels of polyprenol and its derivatives, concomitant with decreased levels of dolichol and derivatives, but no change in polyprenal levels, suggesting DHRSX deficiency. Accordingly, N-glycan synthesis and changes in polyisoprenoid levels were corrected by complementation with human DHRSX but not with SRD5A3. Furthermore, the typical polyprenol dehydrogenase and dolichal reductase activities of DHRSX were absent in membrane preparations derived from Lec5 and Lec9 cells, while the reduction of polyprenal to dolichal, catalyzed by SRD5A3, was unaffected. Long-read whole genome sequencing of Lec5 and Lec9 cells did not reveal mutations in the ORF of SRD5A3, but the genomic region containing DHRSX was absent. Lastly, we established the sequence of Chinese hamster DHRSX and validated that this protein has similar kinetic properties to the human enzyme. Our work therefore identifies the basis of the dolichol synthesis defect in Chinese hamster ovary Lec5 and Lec9 cells.

Glucose-1,6-bisphosphate: A new gatekeeper of cerebral mitochondrial pyruvate uptake

OBJECTIVE

Glucose-1,6-bisphosphate (G-1,6-BP), a byproduct of glycolysis that is synthesized by phosphoglucomutase 2 like 1 (PGM2L1), is particularly abundant in neurons. G-1,6-BP is sensitive to the glycolytic flux, due to its dependence on 1,3-bisphosphoglycerate as phosphate donor, and the energy state, due to its degradation by inosine monophosphate-activated phosphomannomutase 1. Since the exact role of this metabolite remains unclear, our aim was to elucidate the specific function of G-1,6-BP in the brain.

METHODS

The effect of PGM2L1 on neuronal post-ischemic viability was assessed by siRNA-mediated knockdown of PGM2L1 in primary mouse neurons. Acute mouse brain slices were used to correlate the reduction in G-1,6-BP upon ischemia to changes in carbon metabolism by 13C6-glucose tracing. A drug affinity responsive target stability assay was used to test if G-1,6-BP interacts with the mitochondrial pyruvate carrier (MPC) subunits in mouse brain protein extracts. Human embryonic kidney cells expressing a MPC bioluminescence resonance energy transfer sensor were used to analyze how PGM2L1 overexpression affects MPC activity. The effect of G-1,6-BP on mitochondrial pyruvate uptake and oxygen consumption rates was analyzed in isolated mouse brain mitochondria. PGM2L1 and a predicted upstream kinase were overexpressed in a human neuroblastoma cell line and G-1,6-BP levels were measured.

RESULTS

We found that G-1,6-BP in mouse brain slices was quickly degraded upon ischemia and reperfusion. Knockdown of PGM2L1 in mouse neurons reduced post-ischemic viability, indicating that PGM2L1 plays a neuroprotective role. The reduction in G-1,6-BP upon ischemia was not accompanied by alterations in glycolytic rates but we did see a reduced 13C6-glucose incorporation into citrate, suggesting a potential role in mitochondrial pyruvate uptake or metabolism. Indeed, G-1,6-BP interacted with both MPC subunits and overexpression of PGM2L1 increased MPC activity. G-1,6-BP, at concentrations found in the brain, enhanced mitochondrial pyruvate uptake and pyruvate-induced oxygen consumption rates. Overexpression of a predicted upstream kinase inhibited PGM2L1 activity, showing that besides metabolism, also signaling pathways can regulate G-1,6-BP levels.

CONCLUSIONS

We provide evidence that G-1,6-BP positively regulates mitochondrial pyruvate uptake and post-ischemic neuronal viability. These compelling data reveal a novel mechanism by which neurons can couple glycolysis-derived pyruvate to the tricarboxylic acid cycle. This process is sensitive to the glycolytic flux, the cell's energetic state, and upstream signaling cascades, offering many regulatory means to fine-tune this critical metabolic step.

The N-glycosylation defect in Lec5 and Lec9 CHO cells is caused by absence of the DHRSX gene

Glycosylation-deficient Chinese hamster ovary (CHO) cell lines have been instrumental in the discovery of N-glycosylation machinery. Yet, the molecular causes of the glycosylation defects in the Lec5 and Lec9 mutants have been elusive, even though for both cell lines a defect in dolichol formation from polyprenol was previously established. We recently found that dolichol synthesis from polyprenol occurs in three steps consisting of the conversion of polyprenol to polyprenal by DHRSX, the reduction of polyprenal to dolichal by SRD5A3 and the reduction of dolichal to dolichol, again by DHRSX. This led us to investigate defective dolichol synthesis in Lec5 and Lec9 cells. Both cell lines showed increased levels of polyprenol and its derivatives, concomitant with decreased levels of dolichol and derivatives, but no change in polyprenal levels, suggesting DHRSX deficiency. Accordingly, N-glycan synthesis and changes in polyisoprenoid levels were corrected by complementation with human DHRSX but not with SRD5A3. Furthermore, the typical polyprenol dehydrogenase and dolichal reductase activities of DHRSX were absent in membrane preparations derived from Lec5 and Lec9 cells, while the reduction of polyprenal to dolichal, catalyzed by SRD5A3, was unaffected. Long-read whole genome sequencing of Lec5 and Lec9 cells did not reveal mutations in the ORF of SRD5A3, but the genomic region containing DHRSX was absent. Lastly, we established the sequence of Chinese hamster DHRSX and validated that this protein has similar kinetic properties to the human enzyme. Our work therefore identifies the basis of the dolichol synthesis defect in CHO Lec5 and Lec9 cells.

STW-MD: a novel spatio-temporal weighting and multi-step decision tree method for considering spatial heterogeneity in brain gene expression data

Gene expression during brain development or abnormal development is a biological process that is highly dynamic in spatio and temporal. Previous studies have mainly focused on individual brain regions or a certain developmental stage. Our motivation is to address this gap by incorporating spatio-temporal information to gain a more complete understanding of brain development or abnormal brain development, such as Alzheimer's disease (AD), and to identify potential determinants of response. In this study, we propose a novel two-step framework based on spatial-temporal information weighting and multi-step decision trees. This framework can effectively exploit the spatial similarity and temporal dependence between different stages and different brain regions, and facilitate differential gene analysis in brain regions with high heterogeneity. We focus on two datasets: the AD dataset, which includes gene expression data from early, middle and late stages, and the brain development dataset, spanning fetal development to adulthood. Our findings highlight the advantages of the proposed framework in discovering gene classes and elucidating their impact on brain development and AD progression across diverse brain regions and stages. These findings align with existing studies and provide insights into the processes of normal and abnormal brain development.

Glycan quality control in and out of the endoplasmic reticulum of mammalian cells

The endoplasmic reticulum (ER) is equipped with multiple quality control systems (QCS) that are necessary for shaping the glycoproteome of eukaryotic cells. These systems facilitate the productive folding of glycoproteins, eliminate defective products, and function as effectors to evoke cellular signaling in response to various cellular stresses. These ER functions largely depend on glycans, which contain sugar-based codes that, when needed, function to recruit carbohydrate-binding proteins that determine the fate of glycoproteins. To ensure their functionality, the biosynthesis of such glycans is therefore strictly monitored by a system that selectively degrades structurally defective glycans before adding them to proteins. This system, which is referred to as the glycan QCS, serves as a mechanism to reduce the risk of abnormal glycosylation under conditions where glycan biosynthesis is genetically or metabolically stalled. On the other hand, glycan QCS increases the risk of global hypoglycosylation by limiting glycan availability, which can lead to protein misfolding and the activation of unfolded protein response to maintaining cell viability or to initiate cell death programs. This review summarizes the current state of our knowledge of the mechanisms underlying glycan QCS in mammals and its physiological and pathological roles in embryogenesis, tumor progression, and congenital disorders associated with abnormal glycosylation.

Sleep deprivation and recovery sleep affect healthy male resident’s pain sensitivity and oxidative stress markers: The medial prefrontal cortex may play a role in sleep deprivation model

Sleep is essential for the body’s repair and recovery, including supplementation with antioxidants to maintain the balance of the body’s redox state. Changes in sleep patterns have been reported to alter this repair function, leading to changes in disease susceptibility or behavior. Here, we recruited healthy male physicians and measured the extent of the effect of overnight sleep deprivation (SD) and recovery sleep (RS) on nociceptive thresholds and systemic (plasma-derived) redox metabolism, namely, the major antioxidants glutathione (GSH), catalase (CAT), malondialdehyde (MDA), and superoxide dismutase (SOD). Twenty subjects underwent morning measurements before and after overnight total SD and RS. We found that one night of SD can lead to increased nociceptive hypersensitivity and the pain scores of the Numerical Rating Scale (NRS) and that one night of RS can reverse this change. Pre- and post-SD biochemical assays showed an increase in MDA levels and CAT activity and a decrease in GSH levels and SOD activity after overnight SD. Biochemical assays before and after RS showed a partial recovery of MDA levels and a basic recovery of CAT activity to baseline levels. An animal study showed that SD can cause a significant decrease in the paw withdrawal threshold and paw withdrawal latency in rats, and after 4 days of unrestricted sleep, pain thresholds can be restored to normal. We performed proteomics in the rat medial prefrontal cortex (mPFC) and showed that 37 proteins were significantly altered after 6 days of SD. Current findings showed that SD causes nociceptive hyperalgesia and oxidative stress, and RS can restore pain thresholds and repair oxidative stress damage in the body. However, one night of RS is not enough for repairing oxidative stress damage in the human body.

Glucose-1,6-Bisphosphate, a Key Metabolic Regulator, Is Synthesized by a Distinct Family of α-Phosphohexomutases Widely Distributed in Prokaryotes

The reactions of α-d-phosphohexomutases (αPHM) are ubiquitous, key to primary metabolism, and essential for several processes in all domains of life. The functionality of these enzymes relies on an initial phosphorylation step which requires the presence of α-d-glucose-1,6-bisphosphate (Glc-1,6-BP). While well investigated in vertebrates, the origin of this activator compound in bacteria is unknown. Here we show that the Slr1334 protein from the unicellular cyanobacterium Synechocysitis sp. PCC 6803 is a Glc-1,6-BP-synthase. Biochemical analysis revealed that Slr1334 efficiently converts fructose-1,6-bisphosphate (Frc-1,6-BP) and α-d-glucose-1-phosphate/α-d-glucose-6-phosphate into Glc-1,6-BP and also catalyzes the reverse reaction. As inferred from phylogenetic analysis, the slr1334 product belongs to a primordial subfamily of αPHMs that is present especially in deeply branching bacteria and also includes human commensals and pathogens. Remarkably, the homologue of Slr1334 in the human gut bacterium Bacteroides salyersiae catalyzes the same reaction, suggesting a conserved and essential role for the members of this αPHM subfamily.

Effects of the T337M and G391V disease-related variants on human phosphoglucomutase 1: structural disruptions large and small

Phosphoglucomutase 1 (PGM1) plays a central role in glucose homeostasis in human cells. Missense variants of this enzyme cause an inborn error of metabolism, which is categorized as a congenital disorder of glycosylation. Here, two disease-related variants of PGM1, T337M and G391V, which are both located in domain 3 of the four-domain protein, were characterized via X-ray crystallography and biochemical assays. The studies show multiple impacts resulting from these dysfunctional variants, including both short- and long-range structural perturbations. In the T337M variant these are limited to a small shift in an active-site loop, consistent with reduced enzyme activity. In contrast, the G391V variant produces a cascade of structural perturbations, including displacement of both the catalytic phosphoserine and metal-binding loops. This work reinforces several themes that were found in prior studies of dysfunctional PGM1 variants, including increased structural flexibility and the outsized impacts of mutations affecting interdomain interfaces. The molecular mechanisms of PGM1 variants have implications for newly described inherited disorders of related enzymes.

Understanding Inborn Errors of Metabolism through Metabolomics

Inborn errors of metabolism (IEMs) are rare diseases caused by a defect in a single enzyme, co-factor, or transport protein. For most IEMs, no effective treatment is available and the exact disease mechanism is unknown. The application of metabolomics and, more specifically, tracer metabolomics in IEM research can help to elucidate these disease mechanisms and hence direct novel therapeutic interventions. In this review, we will describe the different approaches to metabolomics in IEM research. We will discuss the strengths and weaknesses of the different sample types that can be used (biofluids, tissues or cells from model organisms; modified cell lines; and patient fibroblasts) and when each of them is appropriate to use.

Clinical implementation of RNA sequencing for Mendelian disease diagnostics

BACKGROUND

Lack of functional evidence hampers variant interpretation, leaving a large proportion of individuals with a suspected Mendelian disorder without genetic diagnosis after whole genome or whole exome sequencing (WES). Research studies advocate to further sequence transcriptomes to directly and systematically probe gene expression defects. However, collection of additional biopsies and establishment of lab workflows, analytical pipelines, and defined concepts in clinical interpretation of aberrant gene expression are still needed for adopting RNA sequencing (RNA-seq) in routine diagnostics.

METHODS

We implemented an automated RNA-seq protocol and a computational workflow with which we analyzed skin fibroblasts of 303 individuals with a suspected mitochondrial disease that previously underwent WES. We also assessed through simulations how aberrant expression and mono-allelic expression tests depend on RNA-seq coverage.

RESULTS

We detected on average 12,500 genes per sample including around 60% of all disease genes-a coverage substantially higher than with whole blood, supporting the use of skin biopsies. We prioritized genes demonstrating aberrant expression, aberrant splicing, or mono-allelic expression. The pipeline required less than 1 week from sample preparation to result reporting and provided a median of eight disease-associated genes per patient for inspection. A genetic diagnosis was established for 16% of the 205 WES-inconclusive cases. Detection of aberrant expression was a major contributor to diagnosis including instances of 50% reduction, which, together with mono-allelic expression, allowed for the diagnosis of dominant disorders caused by haploinsufficiency. Moreover, calling aberrant splicing and variants from RNA-seq data enabled detecting and validating splice-disrupting variants, of which the majority fell outside WES-covered regions.

CONCLUSION

Together, these results show that streamlined experimental and computational processes can accelerate the implementation of RNA-seq in routine diagnostics.