In Silico Modeling of Fabry Disease Pathophysiology for the Identification of Early Cellular Damage Biomarker Candidates.

Fabry disease (FD) is an X-linked lysosomal disease whose ultimate consequences are the accumulation of sphingolipids and subsequent inflammatory events, mainly at the endothelial level. The outcomes include different nervous system manifestations as well as multiple organ damage. Despite the availability of known biomarkers, early detection of FD remains a medical need. This study aimed to develop an in silico model based on machine learning to identify candidate vascular and nervous system proteins for early FD damage detection at the cellular level. A combined systems biology and machine learning approach was carried out considering molecular characteristics of FD to create a computational model of vascular and nervous system disease. A data science strategy was applied to identify risk classifiers by using 10 K-fold cross-validation. Further biological and clinical criteria were used to prioritize the most promising candidates, resulting in the identification of 36 biomarker candidates with classifier abilities, which are easily measurable in body fluids. Among them, we propose four candidates, CAMK2A, ILK, LMNA, and KHSRP, which have high classification capabilities according to our models (cross-validated accuracy ≥ 90%) and are related to the vascular and nervous systems. These biomarkers show promise as high-risk cellular and tissue damage indicators that are potentially applicable in clinical settings, although in vivo validation is still needed.

Machine Learning-Driven Biomarker Discovery for Skeletal Complications in Type 1 Gaucher Disease Patients.

Type 1 Gaucher disease (GD1) is a rare, autosomal recessive disorder caused by glucocerebrosidase deficiency. Skeletal manifestations represent one of the most debilitating and potentially irreversible complications of GD1. Although imaging studies are the gold standard, early diagnostic/prognostic tools, such as molecular biomarkers, are needed for the rapid management of skeletal complications. This study aimed to identify potential protein biomarkers capable of predicting the early diagnosis of bone skeletal complications in GD1 patients using artificial intelligence. An in silico study was performed using the novel Therapeutic Performance Mapping System methodology to construct mathematical models of GD1-associated complications at the protein level. Pathophysiological characterization was performed before modeling, and a data science strategy was applied to the predicted protein activity for each protein in the models to identify classifiers. Statistical criteria were used to prioritize the most promising candidates, and 18 candidates were identified. Among them, PDGFB, IL1R2, PTH and CCL3 (MIP-1α) were highlighted due to their ease of measurement in blood. This study proposes a validated novel tool to discover new protein biomarkers to support clinician decision-making in an area where medical needs have not yet been met. However, confirming the results using in vitro and/or in vivo studies is necessary.

Frequency of rs1051338 and rs116928232 Variants in Individuals from Northwest Mexico.

BACKGROUND: LIPA, situated on chromosome 10q23.2-q23.3, encodes the enzyme lysosomal acid lipase (LAL) (EC 3.1.1.13). Genetic alterations in LIPA lead to lysosomal acid lipase deficiency (LALD), an inborn error causing lipid metabolism anomalies and impairing cholesterol and triacylglyceride degradation. Over 40 LIPA variants have been documented, yet this study focuses on just two. The rs1051338 variant (NM_000235:c.46A>C) affects the signal peptide in Exon 2, whereas rs116928232, located in Exon 8, alters the splice site (NM_000235:c.894G>A), impacting lysosomal acid lipase activity. Considering the diverse clinical manifestations of LALD and the rising hepatic steatosis prevalence in Mexican population, mainly due to diet, these variants were investigated within this demographic to uncover potential contributing factors. This study aimed to reveal the frequency of rs1051338 and rs116928232 among healthy mestizo individuals in Northwest Mexico, marking a significant genetic exploration in this demographic. METHODS: Three hundred ten healthy mestizo individuals underwent PCR-RFLP analysis for both variants, and Sanger sequencing was performed for variant rs116928232. Bioinformatic analysis was also performed to predict protein changes. RESULTS: Allele frequencies for rs1051338 (FA = 0.39, p value = 0.15) and rs116928232 (FA = 0.0016, p value = 0.49) aligned with reported data, while bioinformatic analysis allowed us to identify the protein alteration observed in both variants; finally, the variants showed no linkage between them (normalized D' = 1.03, p value = 0.56). CONCLUSIONS: Allelic frequencies closely matched reported data, and protein structure analysis confirmed variant impacts on LAL enzyme function. Notably, this study marks the first analysis of rs1051338 and rs116928232 in a healthy Mexican mestizo population.

Identification of risk features for complication in Gaucher's disease patients: a machine learning analysis of the Spanish registry of Gaucher disease.

BACKGROUND: Since enzyme replacement therapy for Gaucher disease (MIM#230800) has become available, both awareness of and the natural history of the disease have changed. However, there remain unmet needs such as the identification of patients at risk of developing bone crisis during therapy and late complications such as cancer or parkinsonism. The Spanish Gaucher Disease Registry has worked since 1993 to compile demographic, clinical, genetic, analytical, imaging and follow-up data from more than 400 patients. The aims of this study were to discover correlations between patients' characteristics at diagnosis and to identify risk features for the development of late complications; for this a machine learning approach involving correlation networks and decision trees analyses was applied. RESULTS: A total of 358 patients, 340 type 1 Gaucher disease and 18 type 3 cases were selected. 18% were splenectomyzed and 39% had advanced bone disease. 81% of cases carried heterozygous genotype. 47% of them were diagnosed before the year 2000. Mean age at diagnosis and therapy were 28 and 31.5 years old (y.o.) respectively. 4% developed monoclonal gammopathy undetermined significance or Parkinson Disease, 6% cancer, and 10% died before this study. Previous splenectomy correlates with the development of skeletal complications and severe bone disease (p = 0.005); serum levels of IgA, delayed age at start therapy (> 9.5 y.o. since diagnosis) also correlates with severe bone disease at diagnosis and with the incidence of bone crisis during therapy. High IgG (> 1750 mg/dL) levels and age over 60 y.o. at diagnosis were found to be related with the development of cancer. When modelling the decision tree, patients with a delayed diagnosis and therapy were the most severe and with higher risk of complications. CONCLUSIONS: Our work confirms previous observations, highlights the importance of early diagnosis and therapy and identifies new risk features such as high IgA and IgG levels for long-term complications.

Serum protein profile analysis in lysosomal storage disorders patients.

INTRODUCTION: Serum protein electrophoresis (SPE) is a well-established technique to identify alterations in plasma protein profiles, caused by diseases as multiple myeloma (MM). In addition, it could be a cost-effective technique to discover new plasma biomarkers. Relation between MM and lysosomal storage diseases (LSDs) as Gaucher disease has been set out but, it has not been evaluated on other LSDs nor the utility of the SPE as first step on LSDs biomarkers discovery projects. MATERIALS AND METHODS: Stored plasma samples at diagnosis from several LSDs patients underwent analysis. Quality control was checked prior to the SPE was analyzed by capillary electrophoresis. The analysis for monoclonal spikes and the differences between each fraction on patients' samples vs the control data previously published, were evaluated. Furthermore, immunoprotein quantification and free light chains ratio were done by nephelometry and turbidimetry. RESULTS: Seventy-five samples of LSD patients at diagnosis, were assessed. The frequency of the MGUS on LSDs patients was not higher than in general population whereas one lysosomal acid lipase deficiency infant showed increased IgA and kappa deviation. Regarding to the usefulness of SPE in biomarkers discovery, statistically significant differences were observed on SPE fractions between LSDs and healthy population. DISCUSSION: The evaluation of SPE fractions can be a useful tool to understand pathophysiologic aspects in LDSs and, to simplify new marker discovery projects. In some of them, the MGUS appearance is a risk factor for the MM development despite its frequency is not increased on the studied LSDs at diagnosis.

LC-MS/MS analysis of plasma glucosylsphingosine as a biomarker for diagnosis and follow-up monitoring in Gaucher disease in the Spanish population.

Background Gaucher disease (GD), caused by a deficiency in acid ß-glucosidase, leads to the accumulation of glucosylsphingosine (GluSph), which has been used as a powerful biomarker for the diagnosis and follow-up of GD. Our aim was to perform the first retrospective study of GluSph in Spanish patients, analyzing its relationship with classical biomarkers and other parameters of disease and its utility regarding treatment monitoring. Methods Classical biomarkers were evaluated retrospectively by standard methods in a total of 145 subjects, including 47 GD patients, carriers, healthy controls and patients suffering from other lysosomal lipidoses. GluSph was also measured using a liquid chromatography-tandem mass spectrometry (LC-MS/MS) method developed as part of the present study. Results The optimized method presented intra- and inter-assay variations of 3.1 and 11.5%, respectively, overall recovery higher than 96% and linearity up to plasma concentrations of 1000 ng/mL with 100% specificity and sensitivity. Only GD patients displayed GluSph levels above 5.4 ng/mL at diagnosis and this was significantly correlated with the classical biomarkers chitotriosidase (r = 0.560) and the chemokine CCL18/PARC (CCL18/PARC) (ρ = 0.515), as well as with the Spanish magnetic resonance imaging index (S-MRI, r = 0.364), whereas chitotriosidase correlated with liver volume (r = 0.372) and CCL18/PARC increased in patients with bone manifestations (p = 0.005). GluSph levels decreased with treatment in naïve patients. Conclusions Plasma GluSph is the most disease-specific biomarker for GD with demonstrated diagnostic value and responsiveness to therapy. GluSph in the present series of patients failed to demonstrate better correlations with clinical characteristics at onset than classical biomarkers.

New variants in Spanish Niemann-Pick type c disease patients.

Niemann-Pick type C (NPC) disease is a rare inherited disease, with progressive neurodegeneration as the main symptom. It is a lysosomal storage disorder caused by mutations in NPC1 or NPC2 genes, leading to a lysosomal cholesterol trafficking impairment. Disease indicators are the clinical suspicion and biomarker levels. However, a genetic study is mandatory for the diagnosis, which is complicated due to the different variants with unknown significance. The aim of this work was to identify the variants responsible for NPC in our pediatric population. Twenty-two samples from non-related infants believed to have NPC disease were analyzed during the last 3 years. Surrogate biomarkers of the disease were evaluated whenever possible. Sanger sequencing for both genes is reported for all samples. Complementary genetic studies were performed when necessary. NPC disease was confirmed in 31.8% of subjects due to homozygous or compound heterozygous genetic variants in NPC1. The following four novel variants were identified: a gross deletion variant composed of the gene promoter and the first exon, NM_000271.3:c.385delT, NM_000271.3:c.1553+1342_1655-291del, and NM_000271.3:c.1757delA. None had functional activity and all resulted in important structural changes in the protein.

Muscle-tendon weakness contributes to chronic fatigue syndrome in Gaucher's disease.

BACKGROUND: Chronic fatigue (CFg) is a prevalent symptom in Gaucher disease (GD) at diagnosis (79%) and remains in a quarter of patients after years of therapy. Bone abnormalities are present in over 70% and peripheral neuropathy in about 11% of the patients, which contributes to the disabling and debilitating complications. Our hypothesis is that other factors such as muscle-tendinous weakness could have influence in the development of CFg. METHODS: We have evaluated the fiber structure and elasticity of muscle-tendinous unit by strain-elastography (S-ELA) and analyzed their influence in the CFg. S-ELA study was performed in Achilles tendon in 25 type 1 and two type 3 GD patients, all of them with fatigue and were on enzymatic replacement therapy for mean 13 years; simultaneously, bone marrow burden by MRI and calcaneus ultrasound densitometry were evaluated. Blood cell counts, plasma biomarkers, GBA1 genotyping, and SF36 quality of life scale (QoL) were also performed. STATISTICAL ANALYSIS: descriptive and comparative test. RESULTS: All patients showed a normal Achilles tendinous structure. Abnormal stiff grade 2-3 was found in 17/27 (62.9%); in 11/27 (40.7%) of patients, the alteration was bilateral. There were no correlations between the S-ELA results to other variables; nevertheless, a significant correlation between the degree of tendon hardness and the low score on the QoL scales (p = 0.0035) was found. The S-ELA is a sensitive painless, fast, and low cost method to detect muscle-tendinous subclinical dysfunction that could contribute to CFg in GD. The identification of subclinical tendon alteration would be a sign of alarm, focused on the risk of development of bone complications. CONCLUSION: Intratendinous alteration in strain-elastography is an independent variable in GD patients with persistent fatigue.

Web-Based Bioinformatics Predictors: Recommendations to Assess Lysosomal Cholesterol Trafficking Diseases-Related Genes.

INTRODUCTION: The growing number of genetic variants of unknown significance (VUS) and availability of several in silico prediction tools make the evaluation of potentially deleterious gene variants challenging. MATERIALS AND METHODS: We evaluated several programs and software to determine the one that can predict the impact of genetic variants found in lysosomal storage disorders (LSDs) caused by defects in cholesterol trafficking best. We evaluated the sensitivity, specificity, accuracy, precision, and Matthew's correlation coefficient of the most common software. RESULTS: Our findings showed that for exonic variants, only MutPred1 reached 100% accuracy and generated the best predictions (sensitivity and accuracy = 1.00), whereas intronic variants, SROOGLE or Human Splicing Finder (HSF) generated the best predictions (sensitivity = 1.00, and accuracy = 1.00). DISCUSSION: Next-generation sequencing substantially increased the number of detected genetic variants, most of which were considered to be VUS, creating a need for accurate pathogenicity prediction. The focus of the present study is the importance of accurately predicting LSDs, with majority of previously unreported specific mutations. CONCLUSION: We found that the best prediction tool for the NPC1, NPC2, and LIPA variants was MutPred1 for exonic regions and HSF and SROOGLE for intronic regions and splice sites.

Mutations identified in a cohort of Mexican patients with lysosomal acid lipase deficiency.

INTRODUCTION AND OBJECTIVES: Lysosomal acid lipase deficiency (LAL-D) is an autosomal recessive disease caused by mutations in the LIPA gene, located on the long arm of chromosome 10 (10q23.31). Up until now, more than 59 mutations have been described and which are the cause of a very wide clinical spectrum. The goal of this study was to identify the mutations present in Mexican pediatric patients with a diagnosis of LAL-D. MATERIALS AND METHODS: A cross-sectional study was carried out which included all the pediatric patients with LAL-D treated in a tertiary hospital in Mexico from January 2000 to June 2017. RESULTS: Sixteen patients with LAL-D were identified with a disease phenotype marked by the accumulation of cholesteryl esters. Eight distinct variants in the LIPA gene sequence were found, four pathogenic variants and four probably pathogenic. In six individuals, the variants were found in the homozygous state and ten were compound heterozygous. The eight variants were inverted, with five found on exon 4 and the others on exons 2, 8 and 10. The variant c.386A>G;p.His129Arg was the most common, being found in six of the 16 individuals (37.5%), making it much more frequent than what had previously been reported in the literature in proportion to the rest of the variants. The mutation known as E8SJM, which has been the mostly frequently found at the international level, was not the most common among this group of Mexican patients. In conclusion, Mexican patients present a different frequency of mutations associated with LAL-D in comparison to European populations.

Evaluation of two approaches to lysosomal acid lipase deficiency patient identification: An observational retrospective study.

BACKGROUND AND AIMS: Lysosomal acid lipase deficiency (LALD) leads to the accumulation of cholesteryl esters and/or triglycerides (TG) in lysosomes due to the lack of the enzyme codified by the LIPA gene. The most common symptoms are dyslipidaemia and hypertransaminasemia, together with manifestations common to other lysosomal storage disorders (LSDs), including visceromegalies and elevated plasma biomarkers. Alteration of the lipid-liver profile (LLP) has been widely applied as a criterion for LALD screening, but the usefulness of biomarkers has not yet been explored. Our purpose was to explore the utility of plasma chitotriosidase activity (ChT) and CCL18/PARC concentration in addition to LLP to identify LALD patients in an observational retrospective study of two different sample collections. METHODS: Biological samples refining: Collection 1 (primary hypercholesterolemia suspected) included unrelated individuals with hyperlipidaemia and without LDLR, APOB and PCSK9 gene mutations (Set 1), and Collection 2 (LSD suspected) included individuals without definitive LSD diagnosis (Set 2). We assessed plasma LLP (total cholesterol and its fractions, TG concentration and transaminases activities), as well as plasma ChT and CCL18/PARC. All subjects with anomalous LLP and/or biomarker levels were LIPA sequenced. RESULTS: Twenty-four subjects showed altered LLP and/or biomarkers. We identified two LALD patients (one homozygous and one compound heterozygous) and one carrier of a novel LIPA variant. CONCLUSIONS: The measurement of plasma ChT and CCL18/PARC combined with LLP will be a useful approach to identifying LALD patients in retrospective LALD patient studies.

The erythrocyte osmotic resistance test as screening tool for cholesterol-related lysosomal storage diseases.

BACKGROUND: Erythrocyte volume regulation and membrane elasticity are essential for adaptation to osmotic and mechanical stress, and life span. Here, we evaluated whether defective cholesterol trafficking caused by the rare lysosomal storages diseases (LSDs), Niemann-Pick type C (NPC) and Lysosomal acid lipase (LAL) deficiency (LALD) impairs these properties. Moreover, we tested whether measurements of cholesterol membrane content and osmotic resistance serve as a screening test for these LSDs. METHODS: Patients were genotyped for mutations in NPC1, NPC2, or LIPA genes. We measured LSD plasma biomarkers and LAL activity. Red blood cells (RBC) membrane cholesterol content was evaluated in 73 subjects. Osmotic resistance tests (ORT) were conducted in 121 blood samples from LSD suspected patients and controls. RESULTS: We did not find statistically significant differences between RBC cholesterol content between subjects and controls. However, the ORT, particularly at 0.49% (w/v) hypotonic sodium chloride solution, revealed a significant higher osmotic resistance in LSDs patients than in controls. We established a cut-off value of ≤51% of haemolysis with sensibility and specificity values of 80% and 70%, respectively. CONCLUSIONS: NPC and LALD do not alter cholesterol content in the RBC membrane but increase osmotic resistance. Therefore, ORT serves as screening test for the studied LSDs.

Update on lysosomal acid lipase deficiency: Diagnosis, treatment and patient management. / Actualización en deficiencia de lipasa ácida lisosomal: diagnóstico, tratamiento y seguimiento de los pacientes.

Lysosomal acid lipase deficiency (LALD) is an ultra-rare disease caused by a congenital disorder of the lipid metabolism, characterized by the deposition of cholesterol esters and triglycerides in the organism. In patients with no enzyme function, the disease develops during the perinatal period and is invariably associated with death during the first year of life. In all other cases, the phenotype is heterogeneous, although most patients develop chronic liver diseases and may also develop an early cardiovascular disease. Treatment for LALD has classically included the use of supportive measures that do not prevent the progression of the disease. In 2015, regulatory agencies approved the use of a human recombinant LAL for the treatment of LALD. This long-term enzyme replacement therapy has been associated with significant improvements in the hepatic and lipid profiles of patients with LALD, increasing survival rates in infants with a rapidly progressive disease. Both the severity of LALD and the availability of a specific treatment highlight the need to identify these patients in clinical settings, although its low prevalence and the existing clinical overlap with other more frequent pathologies limit its diagnosis. In this paper we set out practical recommendations to identify and monitor patients with LALD, including a diagnostic algorithm, along with an updated treatment.

Inhibition of Intermediate-Conductance Calcium-Activated K Channel (KCa3.1) and Fibroblast Mitogenesis by α-Linolenic Acid and Alterations of Channel Expression in the Lysosomal Storage Disorders, Fabry Disease, and Niemann Pick C.

The calcium/calmodulin-gated KCa3.1 channel regulates normal and abnormal mitogenesis by controlling K+-efflux, cell volume, and membrane hyperpolarization-driven calcium-entry. Recent studies suggest modulation of KCa3.1 by omega-3 fatty acids as negative modulators and impaired KCa3.1 functions in the inherited lysosomal storage disorder (LSD), Fabry disease (FD). In the first part of present study, we characterize KCa3.1 in murine and human fibroblasts and test the impact of omega-3 fatty acids on fibroblast proliferation. In the second, we study whether KCa3.1 is altered in the LSDs, FD, and Niemann-Pick disease type C (NPC). Our patch-clamp and mRNA-expression studies on murine and human fibroblasts show functional expression of KCa3.1. KCa currents display the typical pharmacological fingerprint of KCa3.1: Ca2+-activation, potentiation by the positive-gating modulators, SKA-31 and SKA-121, and inhibition by TRAM-34, Senicapoc (ICA-17043), and the negative-gating modulator, 13b. Considering modulation by omega-3 fatty acids we found that α-linolenic acid (α-LA) and docosahexanenoic acid (DHA) inhibit KCa3.1 currents and strongly reduce fibroblast growth. The α-LA-rich linseed oil and γ-LA-rich borage oil at 0.5% produce channel inhibition while α-LA/γ-LA-low oils has no anti-proliferative effect. Concerning KCa3.1 in LSD, mRNA expression studies, and patch-clamp on primary fibroblasts from FD and NPC patients reveal lower KCa3.1-gene expression and membrane expression than in control fibroblasts. In conclusion, the omega-3 fatty acid, α-LA, and α-LA/γ-LA-rich plant oils, inhibit fibroblast KCa3.1 channels and mitogenesis. Reduced fibroblast KCa3.1 functions are a feature and possible biomarker of cell dysfunction in FD and NPC and supports the concept that biased lipid metabolism is capable of negatively modulating KCa3.1 expression.

Increased glycolipid storage produced by the inheritance of a complex intronic haplotype in the α-galactosidase A (GLA) gene.

BACKGROUND: Accumulation of galactosphingolipids is a general characteristic of Fabry disease, a lysosomal storage disorder caused by the deficient activity of α-galactosidase A encoded by the GLA gene. Although many polymorphic GLA haplotypes have been described, it is still unclear whether some of these variants are causative of disease symptoms. We report the study of an inheritance of a complex intronic haplotype (CIH) (c.-10C > T, c.369 + 990C > A, c.370-81_370-77delCAGCC, c.640-16A > G, c.1000-22C > T) within the GLA gene associated with Fabry-like symptoms and galactosphingolipid accumulation. We analysed α-Gal A activity in plasma, leukocytes and skin fibroblasts in patients, and measured accumulation of galactosphingolipids by enzymatic methods and immunofluorescence techniques. Additionally, we evaluated GLA expression using quantitative PCR, EMSA, and cDNA cloning. RESULTS: CIH carriers had an altered GLA expression pattern, although most of the carriers had high residual enzyme activity in plasma, leukocytes and in skin fibroblasts. Nonetheless, CIH carriers had significant galactosphingolipid accumulation in fibroblasts in comparison with controls, and also glycolipid deposits in renal tubules and glomeruli. EMSA assays indicated that the c.-10C > T variant in the promoter affected a nuclear protein binding site. CONCLUSIONS: Thus, inheritance of the CIH caused an mRNA deregulation altering the GLA expression pattern, producing a tissue glycolipid storage.