- 1Shanghai National Clinical Research Center for Metabolic Diseases, Key Laboratory for Endocrine and Metabolic Diseases of the National Health Commission of the PR China, Shanghai National Center for Translational Medicine, Shanghai, China
- 2Department of Endocrine and Metabolic Diseases, Shanghai Institute of Endocrine and Metabolic Diseases, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China
Background: The WFS1 gene encodes the protein wolframin, which is crucial for maintaining endoplasmic reticulum homeostasis. Variants in this gene are predominantly associated with Wolfram syndrome and have been implicated in other disorders such as diabetes mellitus and psychiatric diseases, which increases the rate of clinical misdiagnosis.
Methods: Patients were diagnosed with early-onset unclassified diabetes according to their clinical and laboratory data. We performed whole-exome sequencing (WES) in 165 patients, interpreting variants according to the American College of Medical Genetics/Association for Molecular Pathology (ACMG/AMP) 2015 guidelines. Variant verification was done by Sanger sequencing. In vitro experiments were conducted to evaluate the effects of WFS1 compound heterozygous variants.
Results: We identified WFS1 compound heterozygous variants (p.A214fs*74/p.F329I and p.I427S/p.I304T) in two patients with Wolfram Syndrome-Like disorders (WSLD). Both WFS1 compound heterozygous variants were associated with increased ER stress, reduced cell viability, and decreased SERCA2b mRNA levels. Additionally, pathogenic or likely pathogenic WFS1 heterozygous variants were identified in the other three patients.
Conclusion: Our results underscore the importance of early genetic testing for diagnosing young-onset diabetes and highlight the clinical relevance of WFS1 variants in increasing ER stress and reducing cell viability. Incorporating these genetic insights into clinical practice can reduce misdiagnoses and improve treatment strategies for related disorders.
Introduction
Wolfram syndrome (WS) is a rare autosomal recessive disorder, caused by variants in WFS1, which encodes the wolframin protein (Rigoli et al., 2011). The primary phenotypes of WS are diabetes insipidus (DI), diabetes mellitus (DM), optic atrophy (OA), and deafness (D), also known as DIDMOAD syndrome (Barrett et al., 1995). Over 98% of these patients develop diabetes, which is often the first symptom in nearly 80% of cases (Lusk et al., 2020). The clinical diagnosis of WS typically requires the coexistence of diabetes and bilateral optic atrophy (Cano et al., 2007). However, only about half of the cases have the complete DIDMOAD phenotype (de Heredia et al., 2013). Some reports are suggesting that WFS1 can also cause diabetes in an autosomal dominant mode (Bonnycastle et al., 2013), that WFS1 heterozygous variants can also cause WS (Gong et al., 2021; Morikawa et al., 2017), and that WFS1 compound heterozygous variants can also cause diabetes only (Li M. et al., 2020), complicating genotype-phenotype correlations and risking misdiagnosis. Therefore, these heterozygous variants often present with a diverse range of clinical manifestations, from traditional WS symptoms to subtler diabetic features, complicating their diagnosis and subsequent management.
The WFS1 gene is located on chromosome 4p16.1 and spans a length of 33.4 kb, encoding wolframin protein comprising of 890 amino acids. It consists 8 exons, of which exon 1 is non-coding, and exons 2 to 7 are responsible for encoding the protein, with exon 8 mainly encoding the transmembrane region and the carboxyl terminus (Li L. et al., 2020). Islet β-cells are the primary expression site for WFS1 (Morikawa et al., 2017). Wolframin is a resident protein in the endoplasmic reticulum (ER) that maintains Ca2+ homeostasis and regulates ER stress (Fonseca et al., 2010). The development of diabetes in WS is associated with elevated ER stress in β-cells. The ER is responsible for proper protein folding and degradation of misfolded proteins.
In our study of 165 patients with early-onset unclassified diabetes, we identified two individuals with compound heterozygous WFS1 variants. We conducted in vitro experiments to evaluate the functional impact of compound heterozygous WFS1 variants. We also found three patients with WFS1 heterozygous variants considered pathogenic or likely pathogenic (P/LP). This research broadens our understanding of the WFS1 gene and its phenotypes.
Materials and methods
Patients
Between January 2021 and December 2022, we enrolled 225 patients with young-onset diabetes at Ruijin Hospital, affiliated with Shanghai Jiao Tong University School of Medicine. Inclusion criteria were an onset age of ≤30 years and the absence of pancreatic exocrine diseases or other conditions inducing diabetes. Patients were diagnosed using the 2019 WHO guidelines (World Health Organization, 2019). 225 patients were aged 20 [15, 16] and the youngest age of onset was 8 years. All patients were tested for GADA and GADA <7 IU/mL was considered negative. Subsequently, 39 patients with T1DM and 21 with T2DM were excluded. Among the remaining 165 patients with unclassified diabetes, all tested negative for GADA and exhibited no typical T1DM or T2DM phenotypes. Clinical data were collected from the medical records. HbA1c and glucose levels were measured in the central laboratory of Ruijin Hospital using high-performance liquid chromatography and an autoanalyzer, respectively. The study received approval from the Ethics Committee of Ruijin Hospital, and informed consent was obtained from all participants and their relatives.
Genetic analysis
We used the QIAamp DNA Blood Mini Kit (Qiagen, Germany) to extract genomic DNA according to the manufacturer’s protocol. The integrity and concentration of the extracted DNA were assessed using a NanoDrop spectrophotometer and agarose gel electrophoresis. Library capture was performed using custom probes from Integrated DNA Technologies, Inc. (Integrated DNA Technologies, United States), which were biotinylated to allow for sequence enrichment using the xGenTM Hybridization and Wash Kit (Integrated DNA Technologies, United States). We used the Illumina NextSeq 500 system (Illumina, United States) to sequence captured libraries, and 150 bp paired-end reads were generated, aiming for a minimum average coverage depth of 100x to ensure adequate coverage of exonic regions. We only focused on 52 genes (Supplementary Table S1) and interpreted those genes (Caruso et al., 2014; Hosoe et al., 2017; Guillín-Amarelle et al., 2018; Rutkowska et al., 2022; Bonnefond et al., 2023). Uncommon coding or splicing variants (MAF<1%) in those genes were analyzed according to the ACMG/AMP 2015 guidelines (Richards et al., 2015). All patients and their family members were verified by Sanger sequencing.
Cell culture
The medium for HEK-293T cells was Dulbecco’s Modified Eagle’s Medium (DMEM) (Meilunbio, China) supplemented with 10% fetal bovine serum (FBS) (Gibco, Canada). HEK-293T cells were cultured at 37°C, under 95% air and 5% CO2. We used Lipofectamine 2000 (Invitrogen, United States) to transfect cells with indicated plasmids according to the manufacturer’s instructions.
Plasmids and luciferase reporter assay
Wild type (WT), p.A214fs*74, p.F329I, p.I427S, p.I304T, p.W690fs*706 and p.E385K were constructed and cloned into pCMV vector (Promega, United States), respectively. Variants of p.W690fs*706 and p.E385K were from two unrelated individuals with WS (Gong et al., 2021; Hu et al., 2022). We used QuickChange Site-Directed Mutagenesis Kit (Stratagene, United States) to generate all WFS1 variants and Sanger sequencing to verify the full-length coding sequences of all plasmids. HEK-293T cells were seeded in 48-well plates, and 100 ng of ERSE luciferase plasmid and 20 ng of PRL-SV40 plasmid (expressing Renilla luciferase for normalization) were transfected. HEK-293T cells were also co-transfected with 0.2 μg each of WT, p.A214fs*74, p.F329I, p.I427S, p.I304T, p.W690fs*706, p.E385K, and 0.1 μg each of p.A214fs*74, p.F329I, p.I427S, p.I304T. Twenty-4 hours after transfection, cells were stimulated with thapsigargin (TG, ER stress inducer) (10 nM) for 6 h and then harvested. We used the dual-luciferase reporter assay system (Promega, United States) to measure luciferase activity according to the manufacturer’s protocol.
RNA extraction and real-time PCR
HEK-293T cells were seeded in 48-well plates, transfected with 350 ng of the indicated plasmids (WT, p.A214fs*74, p.F329I, p.I427S, p.I304T, p.W690fs*706, p.E385K) and 175 ng each of p.A214fs*74, p.F329I, p.I427S, p.I304T. RNA was extracted using the EZ-press RNA Purification Kit (EZBioscience, United States) according to the manufacturer’s protocol, and was transcribed into cDNA using the Reverse Transcription System (Promega, United States). An ABI system (Life Technology, United States) and SYBR Green Supermix (Takara, Japan) were used to perform real-time PCR. Primers are shown in the table (Supplementary Table S2).
Protein extraction and Western blotting
HEK-293T cells were seeded in 24-well plates, transfected with 500 ng each of the indicated plasmids (WT, p.A214fs*74, p.F329I, p.I427S, p.I304T, p.W690fs*706, p.E385K) and 250 ng each of p.A214fs*74, p.F329I, p.I427S, p.I304T. After 48 h, protein samples were isolated with radioimmuno-precipitation assay (RIPA) buffer plus a protease inhibitor cocktail (Thermo Scientific, United States). The primary antibodies were incubated overnight at 4°C and the secondary antibody was incubated for 1 h at room temperature. We used eBlot Touch Imager (eBlot, United States) to visualize bands. The antibodies used are shown in the table (Supplementary Table S3).
Cell proliferation analysis
HEK-293T cells were seeded in 96-well plates, transfected with 50 ng each of the indicated plasmids (WT, p.A214fs*74, p.F329I, p.I427S, p.I304T, p.W690fs*706, p.E385K) and 25 ng each of p.A214fs*74, p.F329I, p.I427S, p.I304T. After 24 h, we used Cell Counting Kit-8 (CCK-8) (MCE, United States) to measure cell viability. A microplate reader (Bio-Rad, United States) was used to measure OD values at 450 nm wavelength.
Flow cytometric
HEK-293T cells were seeded in 12-well plates, transfected with 1 μg each of the indicated plasmids (WT, p.A214fs*74, p.F329I, p.I427S, p.I304T, p.W690fs*706, p.E385K) and 500 ng each of p.A214fs*74, p.F329I, p.I427S, p.I304T. We used Annexin V-FITC/PI apoptosis detection kit (YEASEN, China) to stain the cells according to the manufacturer’s protocols. A FACSCalibur™ flow cytometry (BD Biosciences, United States) was used to measure the fluorescent intensity and data were analyzed using FlowJo X 10.0.7 software.
Bisulfite treatment of DNA and PCR analysis
We used the QIAamp DNA Blood Mini Kit (Qiagen, Germany) to extract genomic DNA according to the manufacturer’s protocol. Isolated genomic DNA was treated with an EZ DNA methylation kit (Zymo Research, United States), and then analyzed by ddPCR. We used ddPCR Supermix for Probes (Bio-Rad, United States) to perform ddPCR, cycling conditions were as follows: 95°C for 10 min, 94°C for 30 s, and 57.5°C for 1 min for 40 cycles. The PCR-amplified 96-well plate was placed into the microtiter analyzer of the QX200 Droplet Reader (Bio-Rad, United States) to detect the fluorescence signals of FAM and VIC (HEX). Then, QuantaSoft Software (Bio-Rad, United States) automatically processed the data to obtain a concentration of methylated and unmethylated INS DNA (copies/μL). Then we calculated the ratio of unmethylated to methylated INS DNA.
Statistical analysis
Data are expressed as mean ± SEM. Comparisons between two groups were made by Student’s t -test and one-way ANOVA was used to compare three or more groups. Statistics were analyzed by Student’s t-test. We used GraphPad Prism V9.0 to analyze all statistics. P-value < 0.05 was statistically different.
Results
In our study, we performed WES on 165 patients with young-onset unclassified diabetes (Figure 1). We identified pathogenic/likely pathogenic (P/LP) WFS1 variants in five patients (Table 1) and other P/LP genes linked to MODY, insulin resistance, or lipoatrophic diabetes in seven patients (Supplementary Table S4). We then focused on WFS1. Among these, three patients (2%) exhibited heterozygous variants, while the remaining two patients (1%) had compound heterozygous variants in WFS1. Notably, five variants were localized to exon 8, with the remaining two found in exon 6 (Figure 2A). Our analysis classified five variants as P/LP, and two were classified as variants of uncertain significance (VUS) (Table 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1433060/fgene-15-1433060-HTML/image_m/fgene-15-1433060-t001.jpg)
Table 1. A summary of WFS1 compound heterozygous and P/LP WFS1 heterozygous variants identified in this study.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1433060/fgene-15-1433060-HTML/image_m/fgene-15-1433060-g002.jpg)
Figure 2. Pedigrees and Sanger sequencing of the two patients (P1 and P2) and locations of variants. (A) Locations of variants. (B) Pedigree of the S1 family. (C) Pedigree of the S2 family. (D) Sanger sequencing peak map of the patient 1. (E) Sanger sequencing peak map of the patient two.
Clinical characteristics of patients with WFS1 compound heterozygous variants
Diagnosed with diabetes at 25, patient one has never experienced ketoacidosis since onset, with GADA negative. Managed effectively with oral medication, she maintains satisfactory HbA1c and C-peptide levels (Table 2). Notably, she exhibits none of the symptoms typically associated with WS, such as DI, OA, or D, and has had no diabetes-related complications to date. Her family history reveals a prominent pattern of diabetes mellitus, with all members developing diabetes at age 35 or older (Figure 2B). Genetic testing showed that this patient carries the WFS1 compound heterozygous variants (Figure 2D). Sanger sequencing identified heterozygous variants p.A214fs*74 in her father and uncle, and p.F329I in her mother, while her grandmother, aunt, and maternal aunt showed no variants (Supplementary Figure S1). Patient two was first diagnosed with diabetes at 13, presenting with ketoacidosis. At 14, a hospital evaluation indicated average islet function and negative insulin autoantibodies (Table 2). Her clinical profile includes normal urine output and no neurological, psychiatric, or high-frequency hearing issues, or optic nerve atrophy. However, she has a refractive error in her left eye and difficulty distinguishing blue-violet colors. Her medical history includes hydronephrosis, revealed by abdominal ultrasound as bilateral ureteral dilatation, but pancreatic MR scans showed no occupying lesions. Her maternal grandfather, diagnosed with diabetes at age 40, was treated with insulin and medication (Figure 2C). Genetic analysis identified WFS1 compound heterozygous variants in the patient (Figure 2E). Sanger sequencing in her family revealed heterozygous variants: p.I427S in her father and p.I304T in her mother and brother (Supplementary Figure S2). Other family members were unavailable for genetic testing. According to the guidelines from EURO-WABB (http://euro-wabb.org/guidelines/guidelines/), we can diagnose these patients as WSLD at this moment. Recent findings indicate that certain adolescents exhibiting WSLD present solely with diabetes, leading to frequent misdiagnoses as either T1DM or T2DM (Zalloua et al., 2008; Bonnycastle et al., 2013; Li M. et al., 2020).
Functional analysis of WFS1 compound heterozygous variants
Protein expression of wolframin
We evaluated wolframin protein expression in HEK cells transfected with WFS1 variant-carrying plasmids. All four variants, particularly the frameshift variant p.A214fs*74 and the missense variant p.I427S, significantly reduced wolframin levels, with p.A214fs*74 causing complete absence and p.I427S showing near absence of the protein. Combinations of p.A214fs*74 with p.F329I and p.I304T with p.I427S also notably decreased wolframin expression compared to WT (Figure 3). These results suggest that these variants disrupt wolframin protein production, potentially contributing to the development of the associated disease phenotype.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1433060/fgene-15-1433060-HTML/image_m/fgene-15-1433060-g003.jpg)
Figure 3. Protein expression levels of Wolframin. (MOCK: transfected with pCMV vector; WT: wild type).
Effects on ER stress and UPR
Wolframin protein is mainly located in the ER, and variants in the WFS1 gene cause an imbalance in ER homeostasis, which activates the unfolded protein response (UPR) (Wu and Kaufman, 2006). We performed a dual luciferase reporter gene assay to determine whether WFS1 variants increase ER stress. The p.I427S and the compound heterozygous p.I427S/p.I304T variants significantly activated the ERSE reporter, both with and without TG stimulation, indicating increased ER stress (Figure 4A). However, the p.I304T variant did not change ERSE activity from the WT. The p.A214fs74 variant significantly increased ER stress enhancer (ERSE) activity under thapsigargin (TG) stimulation, in contrast to the compound heterozygous p.A214fs74/p.F329I and the p.F329I variants (Figure 4A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1433060/fgene-15-1433060-HTML/image_m/fgene-15-1433060-g004.jpg)
Figure 4. Effects of WFS1 variants on ER stress and UPR. (A) Luciferase reporter assay in HEK-293T cells (n = 4). (B) HEK-293T cells transfected with mutant plasmids upregulated ER stress-related molecules. (C) Protein expression of CHOP. (D) Quantitative PCR for mRNA levels of XBP-1, ATF4α, BiP/GRP78 and CHOP in HEK-293T cells (n = 4). (MOCK/Empty: transfected with pCMV vector; WT: wild type; TG: thapsigargin, ER stress inducer). All statistics were represented as mean ± SEM. *P < 0.05, **P < 0.01 and ***P < 0.001.
Further, protein levels of ER stress markers BiP/GRP78 and XBP-1 were markedly increased in cells with compound heterozygous variants (Figure 4B). Specifically, the p.I427S/p.I304T variant also slightly raised CHOP protein levels (Figure 4C). The single p.I427S and p.I304T variants each upregulated BiP/GRP78 and XBP-1 proteins, with p.I427S additionally boosting CHOP levels (Figures 4B, C). Correspondingly, mRNA expressions of CHOP, BiP/GRP78, and XBP-1 were higher in cells with the compound heterozygous variants (Figure 4D). These results suggest that WFS1 compound heterozygous variants increase ER stress and activate the UPR pathway.
Evaluation of calcium homeostasis and cell apoptosis
ER stress is triggered by increased Ca2+ efflux from the ER, which leads to an increase in the cytoplasmic Ca2+ concentration and consequent activation of calpain-2, resulting in β-cell death (Hara et al., 2014). Wolframin protein is mainly located in the ER and is responsible for maintaining intracellular Ca2+ homeostasis by repressing the expression of sarcoendoplasmic reticulum Ca2+-ATPase 2b (SERCA2b) (Takei et al., 2006). We examined the impact of WFS1 compound heterozygous variants (p.A214fs*74/p.F329I and p.I304T/p.I427S) on ER calcium homeostasis. WFS1 compound heterozygous variants significantly lowered SERCA2b mRNA levels (Figure 5A), disrupting ER calcium balance. Cell viability, measured by CCK-8 assays, was reduced in cells with WFS1 compound heterozygous variants compared to WT (Figure 5B). Flow cytometry showed no significant increase in apoptosis in cells with WFS1 compound heterozygous variants compared to WT (Figure 5C). Previous studies have shown that the extent of β-cell death in pancreatic islet cells can be analyzed indirectly using droplet digital PCR (ddPCR). Non-methylated INS DNA/methylated INS DNA is significantly higher in the diabetic population compared to the normal population (Usmani-Brown et al., 2014). We then analyzed β cell death by ddPCR, it showed no significant increase in β cell death (Figure 5D). This suggests that while WFS1 compound heterozygous variants impair ER calcium stability and cell viability, they do not significantly heighten β-cell death, providing insights into the cellular effects of WFS1 variants in diabetes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1433060/fgene-15-1433060-HTML/image_m/fgene-15-1433060-g005.jpg)
Figure 5. Evaluation of calcium flux and cell proliferation. (A) Quantitative PCR analysis of SERCA2b (n = 4). (B) Cell proliferation was measured by CCK-8 assay (n = 5). (C) Flow cytometry to assess cell apoptosis. (D) The extent of pancreatic β cell damage was measured through ddPCR. (Vector/Empty: transfected with pCMV vector. WT: wild type). All data were expressed as mean ± SEM. *P < 0.05, **P < 0.01 and ***P < 0.001.
Clinical characteristics of patients with heterozygous WFS1 variants
In addition to the compound heterozygous variants, we identified three different heterozygous WFS1 variants in three patients, each classified as P/LP (Table 1). Bonnefond et al. (2023) included the WFS1 gene in the MODY. Notably, these patients had a significant family history of diabetes, all were diagnosed with diabetes at 25 years old (Supplementary Figure S3; Table 2). They were treated with oral antidiabetic drugs and showed no symptoms of WS. Therefore, these patients were clinically suspected of MODY. The implications of these variants warrant further genetic screening and detailed evaluation of additional family members.
Discussion
Approximately half of WFS1 variants are homozygous, but autosomal dominant pathogenic WFS1 variants have also been identified (Astuti et al., 2017), recessive WFS1 variants can also lead to syndromic and nonsyndromic diabetes. In patients with WS, homozygous loss-of-function variants are the most common, with fewer missense homozygous variants observed (Bansal et al., 2018). This variation in variant types can lead to clinical misdiagnosis or underdiagnosis due to their varying features.
WS is a rare autosomal recessive neurodegenerative disorder that presents with diabetes insipidus, insulin-deficient diabetes mellitus, optic atrophy, and deafness (Khanim et al., 2001). The most typical manifestation in WS patients is diabetes, but neurological diseases and urinary complications are major contributors to morbidity and mortality in these individuals (Kinsley et al., 1995). Diabetes is commonly diagnosed in the first decade of life, while optic atrophy often becomes evident in the second decade. Urinary and neurological abnormalities generally emerge between the ages of 10 and 30 years (Sherif et al., 2021). Traditionally, the concurrent presence of DM and OA has been considered essential for WS diagnosis, some cases lack early optic atrophy, leading to misdiagnosis as T1DM or T2DM (Zmyslowska et al., 2014; Li et al., 2023). In this research, we detected two patients with compound heterozygous variants in WFS1 (P1: p.A214fs*74/p.F329I, P2: p.I427S/p.I304T). The p.A214fs*74 and p.I427S variants were classified as LP and P respectively, the other two were classified as VUS, and computational evidences predicted deleterious impacts from these variants. Both patients presented with diabetes mellitus but without optic atrophy. Considering their family history and patient two showed bilateral ureteral dilatation with pelvic separation, it underscored the necessity for comprehensive clinical evaluations in WS and functional analysis to verify the impacts of the two compound heterozygous variants.
Variants in the WFS1 gene are known to disrupt ER function, linked to β-cell death in permanent neonatal diabetes (Fonseca et al., 2005; Edghill et al., 2008; Bonfanti et al., 2009; Fonseca et al., 2010). The UPR is activated when ER stress occurs, which involves the upregulation of CHOP and transactivation of ATF6α to maintain homeostasis within the ER. As a critical transcription factor, ATF6α activates target genes of the UPR and facilitates the expression of protein-coding genes such as Bip/GRP78, thereby restoring protein folding within the ER lumen (Ye et al., 2000; Wu and Kaufman, 2006). However, WFS1, a component of the UPR, negatively regulates the ER stress signaling network by recruiting ATF6α to HRD1, which leads to the degradation of ATF6α and suppression of the UPR. This ultimately leads to a reduction in CHOP, BiP/GRP78, and XBP-1 (Fonseca et al., 2010; Guo et al., 2011). CHOP, a transcriptional repressor, is typically low under homeostatic conditions and is triggered by ER stress (Mv et al., 1994; McCullough et al., 2001; Kwok and Daskal, 2008). Loss-of-function of the WFS1 gene may result in Ca2+ deletion in the ER (Hara et al., 2014), Ca2+ homeostasis is essential for maintaining cellular function. SERCA2b is an important gene for maintaining Ca2+ balance in the ER, studies have shown that exogenous supplementation of SERCA2b reversed ER Ca2+ efflux and prevented cell death (Hara et al., 2014). Thus, SERCA2b downregulation may be detrimental to maintaining ER Ca2+ homeostasis. Subsequently, elevated cytosolic Ca2+ induced cell death. Our in vitro studies showed that protein expression of WFS1 compound heterozygous variants (p.A214fs*74/p.F329I and p.I304T/p.I427S) significantly decreased compared to WT. The apparent decrease in wolframin protein in the double p.I304T/p.I427S transfection versus the single p.I304T transfection may include a dominant negative effect of the p.I427S variant when co-expressed with p.I304T, potentially through misfolding and subsequent degradation of the protein, or by interference with the translation machinery. WFS1 compound heterozygous variants (p.A214fs*74/p.F329I and p.I304T/p.I427S) upregulated the UPR pathway and decreased SERCA2b mRNA levels, indicating increased calcium efflux. This led to reduced cell viability, while flow cytometry and INS DNA methylation analyses showed no significant differences in apoptosis or β-cell death. Reduced cell viability, in the absence of a noticeable increase in apoptosis or beta-cell death, may be attributed to mitochondrial dysfunction. This dysfunction can diminish cellular viability by impairing ATP production and elevating oxidative stress. Additionally, alterations in the cell cycle can induce a quiescent state that inhibits cell proliferation without leading to cell death. In addition, p.I304T/p.I427S resulted in more severe functional impairment than p.A214fs*74/p.F329I, which may explain why patient two showed a more severe phenotype than patient one.
The pathogenicity of WFS1 variants correlates with their location in the gene. Variants in the luminal domain often lead to early-onset diabetes and typical WS symptoms, while those in the transmembrane and cytoplasmic domains are associated with milder, later-onset forms (Qian et al., 2015). Patient one, with the p.A214fs*74 variant in the cytoplasmic domain and the p.F329I variant in the transmembrane domain (Supplementary Figure S4), was diagnosed with diabetes at 25 but showed no WS symptoms for 7 years. Patient two, carried the p.I304T variant in the cytoplasmic domain and the p.I427S in the luminal domain (Supplementary Figure S4), which may explain why she had an earlier onset and a heavier phenotype than patient one. One report indicated delayed onset of optic atrophy at age 53 in a patient first diagnosed with diabetes at age 33 (Lieber et al., 2012), suggesting the potential for later emergence of WS symptoms. Another study demonstrated that an individual diagnosed with early-onset diabetes was found to be homozygous for the WFS1 gene, yet did not exhibit the additional clinical features typically associated with WS (Bansal et al., 2017). Therefore, we verified the two WFS1 compound heterozygous variants indeed lead to impaired function. We still cannot diagnose them as WS, but now as WSLD. Until now, no studies have shown that isolated diabetes progresses to WS, so this needs to be followed up closely with our patients later.
We also identified three patients with P/LP heterozygous variants in WFS1. Notably, these patients had a strong family history of diabetes and were managed with oral antidiabetic drugs. We, therefore, suspected them with MODY, follow-up is needed as well as family line validation. The demographic and clinical analysis of the five patients with P/LP WFS1 variants in our study has uncovered significant patterns that contribute to the understanding of diabetes phenotypes associated with WFS1 variants. Heterozygous variants in the WFS1 gene present considerable diagnostic challenges in WS, mainly because their clinical manifestations are often milder and less distinct than those seen in homozygous cases. This variability complicates clinical identification and timely diagnosis of WS, potentially delaying appropriate management and counseling, which could impact patient outcomes. Therefore, a deeper understanding of the impact of these heterozygous variants is essential for refining diagnostic criteria and developing targeted therapeutic strategies.
Wolfram syndrome is also associated with optic atrophy, hearing loss, and neurodegeneration, among other issues. Early genetic diagnosis could facilitate the monitoring of these comorbidities and the initiation of supportive therapies. Several classes of drugs have been reported for the treatment of WS, including dantrolene sodium (Nguyen et al., 2020), 4-phenylbutyric acid (4-PBA) and tauroursodeoxycholic acid (TUDCA) (Kitamura et al., 2022), glucagon-like peptide (GLP)-1 receptor agonists (Sedman et al., 2016; Kondo et al., 2018), and valproic acid (Kakiuchi et al., 2009). We are going to closely follow up on the two patients identified as harboring WFS1 compound heterozygous variants with targeted examinations as well as targeted medications.
However, the study’s limitations include a small sample size and a lack of ethnic diversity, which may affect the generalizability of the findings. Future studies involving larger and more diverse cohorts are essential to further validate and enhance our understanding of the role of WFS1 variants in diabetes.
In conclusion, we identified two WFS1 compound heterozygous variants in patients with WSLD. The functional analysis verified the impaired function of the WFS1 compound heterozygous variants. The WFS1 gene is associated with diverse phenotypes, ranging from nonsyndromic diabetes to syndromic conditions. Early genetic screening is crucial for optimal patient management and timely intervention as the disease progresses.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving humans were approved by the ethics committee of Ruijin Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’; legal guardians/next of kin. Written informed consent was obtained from the individual(s), and minor(s)’; legal guardian/next of kin, for the publication of any potentially identifiable images or data included in this article.
Author contributions
LW: Conceptualization, Data curation, Investigation, Methodology, Software, Writing–original draft, Writing–review and editing. JZ: Conceptualization, Investigation, Methodology, Supervision, Writing–review and editing. DL: Data curation, Methodology, Software, Supervision, Writing–review and editing. ZZ: Methodology, Supervision, Writing–review and editing. QN: Conceptualization, Investigation, Methodology, Writing–review and editing. RH: Conceptualization, Methodology, Software, Writing–review and editing. LY: Investigation, Supervision, Writing–review and editing. YZ: Supervision, Writing–review and editing. JH: Supervision, Investigation, Writing–review and editing. WW: Supervision, Writing–review and editing. GN: Supervision, Writing–review and editing. WG: Supervision, Writing–review and editing, Data curation, Funding acquisition, Investigation, Methodology.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was supported by National Natural Science Foundation of China (82070864), Shanghai Special Project "Science and Technology Innovation Action Plan" Medical Innovation Research ( 22Y11904800), and National Key R and D Program of China (2023YFC2506000, 2023YFC2506001).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2024.1433060/full#supplementary-material
References
Astuti, D., Sabir, A., Fulton, P., Zatyka, M., Williams, D., Hardy, C., et al. (2017). Monogenic diabetes syndromes: locus-specific databases for Alström, Wolfram, and Thiamine-responsive megaloblastic anemia. Hum. Mutat. 38 (7), 764–777. doi:10.1002/humu.23233
Bansal, V., Boehm, B. O., and Darvasi, A. (2018). Identification of a missense variant in the WFS1 gene that causes a mild form of Wolfram syndrome and is associated with risk for type 2 diabetes in Ashkenazi Jewish individuals. Diabetologia 61 (10), 2180–2188. doi:10.1007/s00125-018-4690-3
Bansal, V., Gassenhuber, J., Phillips, T., Oliveira, G., Harbaugh, R., Villarasa, N., et al. (2017). Spectrum of mutations in monogenic diabetes genes identified from high-throughput DNA sequencing of 6888 individuals. BMC Med. 15 (1), 213. doi:10.1186/s12916-017-0977-3
Barrett, T. G., Bundey, S. E., and Macleod, A. F. (1995). Neurodegeneration and diabetes: UK nationwide study of Wolfram (DIDMOAD) syndrome. Lancet 346 (8988), 1458–1463. doi:10.1016/s0140-6736(95)92473-6
Bonfanti, R., Colombo, C., Nocerino, V., Massa, O., Lampasona, V., Iafusco, D., et al. (2009). Insulin gene mutations as cause of diabetes in children negative for five type 1 diabetes autoantibodies. Diabetes Care 32 (1), 123–125. doi:10.2337/dc08-0783
Bonnefond, A., Unnikrishnan, R., Doria, A., Vaxillaire, M., Kulkarni, R. N., Mohan, V., et al. (2023). Monogenic diabetes. Nat. Rev. Dis. Prim. 9 (1), 12. doi:10.1038/s41572-023-00421-w
Bonnycastle, L. L., Chines, P. S., Hara, T., Huyghe, J. R., Swift, A. J., Heikinheimo, P., et al. (2013). Autosomal dominant diabetes arising from a Wolfram syndrome 1 mutation. Diabetes 62 (11), 3943–3950. doi:10.2337/db13-0571
Cano, A., Rouzier, C., Monnot, S., Chabrol, B., Conrath, J., Lecomte, P., et al. (2007). Identification of novel mutations inWFS1 and genotype–phenotype correlation in Wolfram syndrome. Am. J. Med. Genet. 143A (14), 1605–1612. doi:10.1002/ajmg.a.31809
Caruso, M., Ma, D., Msallaty, Z., Lewis, M., Seyoum, B., Al-janabi, W., et al. (2014). Increased interaction with insulin receptor substrate 1, a novel abnormality in insulin resistance and type 2 diabetes. Diabetes 63 (6), 1933–1947. doi:10.2337/db13-1872
de Heredia, M. L., Clèries, R., and Nunes, V. (2013). Genotypic classification of patients with Wolfram syndrome: insights into the natural history of the disease and correlation with phenotype. Genet. Med. 15 (7), 497–506. doi:10.1038/gim.2012.180
Edghill, E. L., Flanagan, S. E., Patch, A.-M., Boustred, C., Parrish, A., Shields, B., et al. (2008). Insulin mutation screening in 1044 patients with diabetes: mutations in the INS gene are a common cause of neonatal diabetes but a rare cause of diabetes diagnosed in childhood or adulthood. Diabetes 57 (4), 1034–1042. doi:10.2337/db07-1405
Fonseca, S. G., Fukuma, M., Lipson, K. L., Nguyen, L. X., Allen, J. R., Oka, Y., et al. (2005). WFS1 is a novel component of the unfolded protein response and maintains homeostasis of the endoplasmic reticulum in pancreatic beta-cells. J. Biol. Chem. 280 (47), 39609–39615. doi:10.1074/jbc.M507426200
Fonseca, S. G., Ishigaki, S., Oslowski, C. M., Lu, S., Lipson, K. L., Ghosh, R., et al. (2010). Wolfram syndrome 1 gene negatively regulates ER stress signaling in rodent and human cells. J. Clin. Invest 120 (3), 744–755. doi:10.1172/JCI39678
Gong, Y., Xiong, L., Li, X., Su, L., and Xiao, H. (2021). A novel mutation of WFS1 gene leading to increase ER stress and cell apoptosis is associated an autosomal dominant form of Wolfram syndrome type 1. BMC Endocr. Disord. 21, 76. doi:10.1186/s12902-021-00748-z
Guillín-Amarelle, C., Fernández-Pombo, A., Sánchez-Iglesias, S., and Araújo-Vilar, D. (2018). Lipodystrophic laminopathies: diagnostic clues. Nucleus 9 (1), 249–260. doi:10.1080/19491034.2018.1454167
Guo, X., Shen, S., Song, S., He, S., Cui, Y., Xing, G., et al. (2011). The E3 ligase Smurf1 regulates Wolfram syndrome protein stability at the endoplasmic reticulum. J. Biol. Chem. 286 (20), 18037–18047. doi:10.1074/jbc.M111.225615
Hara, T., Mahadevan, J., Kanekura, K., Hara, M., Lu, S., and Urano, F. (2014). Calcium efflux from the endoplasmic reticulum leads to β-cell death. Endocrinology 155 (3), 758–768. doi:10.1210/en.2013-1519
Hosoe, J., Kadowaki, H., Miya, F., Aizu, K., Kawamura, T., Miyata, I., et al. (2017). Structural basis and genotype–phenotype correlations of INSR mutations causing severe insulin resistance. Diabetes 66 (10), 2713–2723. doi:10.2337/db17-0301
Hu, K., Zatyka, M., Astuti, D., Beer, N., Dias, R. P., Kulkarni, A., et al. (2022). WFS1 protein expression correlates with clinical progression of optic atrophy in patients with Wolfram syndrome. J. Med. Genet. 59 (1), 65–74. doi:10.1136/jmedgenet-2020-107257
Kakiuchi, C., Ishigaki, S., Oslowski, C. M., Fonseca, S. G., Kato, T., and Urano, F. (2009). Valproate, a mood stabilizer, induces WFS1 expression and modulates its interaction with ER stress protein GRP94. PLoS One 4 (1), e4134. doi:10.1371/journal.pone.0004134
Khanim, F., Kirk, J., Latif, F., and Barrett, T. G. (2001). WFS1/wolframin mutations, Wolfram syndrome, and associated diseases. Hum. Mutat. 17 (5), 357–367. doi:10.1002/humu.1110
Kinsley, B. T., Swift, M., Dumont, R. H., and Swift, R. G. (1995). Morbidity and mortality in the Wolfram syndrome. Diabetes Care 18 (12), 1566–1570. doi:10.2337/diacare.18.12.1566
Kitamura, R. A., Maxwell, K. G., Ye, W., Kries, K., Brown, C. M., Augsornworawat, P., et al. (2022). Multidimensional analysis and therapeutic development using patient iPSC-derived disease models of Wolfram syndrome. JCI Insight 7 (18), e156549. doi:10.1172/jci.insight.156549
Kondo, M., Tanabe, K., Amo-Shiinoki, K., Hatanaka, M., Morii, T., Takahashi, H., et al. (2018). Activation of GLP-1 receptor signalling alleviates cellular stresses and improves beta cell function in a mouse model of Wolfram syndrome. Diabetologia 61 (10), 2189–2201. doi:10.1007/s00125-018-4679-y
Kwok, S. C. M., and Daskal, I. (2008). Brefeldin A activates CHOP promoter at the AARE, ERSE and AP-1 elements. Mol. Cell Biochem. 319 (1), 203–208. doi:10.1007/s11010-008-9893-3
Li, L., Venkataraman, L., Chen, S., and Fu, H. (2020b). Function of WFS1 and WFS2 in the central nervous system: implications for Wolfram syndrome and alzheimer’s disease. Neurosci. Biobehav Rev. 118, 775–783. doi:10.1016/j.neubiorev.2020.09.011
Li, M., Wang, S., Xu, K., Chen, Y., Fu, Q., Gu, Y., et al. (2020a). High prevalence of a monogenic cause in han Chinese diagnosed with type 1 diabetes, partly driven by nonsyndromic recessive WFS1 mutations. Diabetes 69 (1), 121–126. doi:10.2337/db19-0510
Li, Y., Gong, S., Li, M., Cai, X., Liu, W., Zhang, S., et al. (2023). The genetic and clinical characteristics of WFS1 related diabetes in Chinese early onset type 2 diabetes. Sci. Rep. 13 (1), 9127. doi:10.1038/s41598-023-36334-7
Lieber, D. S., Vafai, S. B., Horton, L. C., Slate, N. G., Liu, S., Borowsky, M. L., et al. (2012). Atypical case of Wolfram syndrome revealed through targeted exome sequencing in a patient with suspected mitochondrial disease. BMC Med. Genet. 13, 3. doi:10.1186/1471-2350-13-3
Lusk, L., Black, E., and Vengoechea, J. (2020). Segregation of two variants suggests the presence of autosomal dominant and recessive forms of WFS1 -related disease within the same family: expanding the phenotypic spectrum of Wolfram Syndrome. J. Med. Genet. 57 (2), 121–123. doi:10.1136/jmedgenet-2018-105782
Matsunaga, K., Tanabe, K., Inoue, H., Okuya, S., Ohta, Y., Akiyama, M., et al. (2014). Wolfram syndrome in the Japanese population; molecular analysis of WFS1 gene and characterization of clinical features. PLoS One 9 (9), e106906. doi:10.1371/journal.pone.0106906
McCullough, K. D., Martindale, J. L., Klotz, L.-O., Aw, T.-Y., and Holbrook, N. J. (2001). Gadd153 sensitizes cells to endoplasmic reticulum stress by down-regulating Bcl2 and perturbing the cellular redox state. Mol. Cell Biol. 21 (4), 1249–1259. doi:10.1128/MCB.21.4.1249-1259.2001
Morikawa, S., Tajima, T., Nakamura, A., Ishizu, K., and Ariga, T. (2017). A novel heterozygous mutation of the WFS1 gene leading to constitutive endoplasmic reticulum stress is the cause of Wolfram syndrome: MORIKAWA et al. Pediatr. Diabetes 18 (8), 934–941. doi:10.1111/pedi.12513
Mv, B., A, C., A, T., L, P., and D, R. (1994). CHOP (GADD153) and its oncogenic variant, TLS-CHOP, have opposing effects on the induction of G1/S arrest. Genes and Dev. 8 (4), 453–464. doi:10.1101/gad.8.4.453
Nguyen, L. D., Fischer, T. T., Abreu, D., Arroyo, A., Urano, F., and Ehrlich, B. E. (2020). Calpain inhibitor and ibudilast rescue β cell functions in a cellular model of Wolfram syndrome. Proc. Natl. Acad. Sci. U. S. A. 117 (29), 17389–17398. doi:10.1073/pnas.2007136117
Qian, X., Qin, L., Xing, G., and Cao, X. (2015). Phenotype prediction of pathogenic nonsynonymous single nucleotide polymorphisms in WFS1. Sci. Rep. 5, 14731. doi:10.1038/srep14731
Richards, S., Aziz, N., Bale, S., Bick, D., Das, S., Gastier-Foster, J., et al. (2015). Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of medical genetics and genomics and the association for molecular Pathology. Genet. Med. 17 (5), 405–424. doi:10.1038/gim.2015.30
Rigoli, L., Lombardo, F., and Di Bella, C. (2011). Wolfram syndrome and WFS1 gene. Clin. Genet. 79 (2), 103–117. doi:10.1111/j.1399-0004.2010.01522.x
Rohayem, J., Ehlers, C., Wiedemann, B., Holl, R., Oexle, K., Kordonouri, O., et al. (2011). Diabetes and neurodegeneration in Wolfram syndrome: a multicenter study of phenotype and genotype. Diabetes Care 34 (7), 1503–1510. doi:10.2337/dc10-1937
Rutkowska, L., Salachna, D., Lewandowski, K., Lewiński, A., and Gach, A. (2022). Familial partial lipodystrophy—literature review and report of a novel variant in PPARG expanding the spectrum of disease-causing alterations in FPLD3. Diagn. (Basel) 12 (5), 1122. doi:10.3390/diagnostics12051122
Sedman, T., Rünkorg, K., Krass, M., Luuk, H., Plaas, M., Vasar, E., et al. (2016). Exenatide is an effective antihyperglycaemic agent in a mouse model of Wolfram syndrome 1. J. Diabetes Res. 2016, 9239530. doi:10.1155/2016/9239530
Sherif, M., Demirbilek, H., Çayır, A., Tahir, S., Çavdarlı, B., Demiral, M., et al. (2021). Identification of three novel and one known mutation in the WFS1 gene in four unrelated Turkish families: the role of homozygosity mapping in the early diagnosis. J. Clin. Res. Pediatr. Endocrinol. 13 (1), 34–43. doi:10.4274/jcrpe.galenos.2020.2020.0152
Strom, T., Hörtnagel, K., Hofmann, S., Gekeler, F., Scharfe, C., Rabl, W., et al. (1998). Diabetes insipidus, diabetes mellitus, optic atrophy and deafness (DIDMOAD) caused by mutations in a novel gene (wolframin) coding for a predicted transmembrane protein. Hum. Mol. Genet. 7 (13), 2021–2028. doi:10.1093/hmg/7.13.2021
Takei, D., Ishihara, H., Yamaguchi, S., Yamada, T., Tamura, A., Katagiri, H., et al. (2006). WFS1 protein modulates the free Ca2+ concentration in the endoplasmic reticulum. FEBS Lett. 580 (24), 5635–5640. doi:10.1016/j.febslet.2006.09.007
Usmani-Brown, S., Lebastchi, J., Steck, A. K., Beam, C., Herold, K. C., and Ledizet, M. (2014). Analysis of β-cell death in type 1 diabetes by droplet digital PCR. Endocrinology 155 (9), 3694–3698. doi:10.1210/en.2014-1150
World Health Organization (2019). Classification of diabetes mellitus. Geneva: World Health Organization. Available at: https://apps.who.int/iris/handle/10665/325182 (Accessed June 7, 2023).
Wu, J., and Kaufman, R. J. (2006). From acute ER stress to physiological roles of the Unfolded Protein Response. Cell Death Differ. 13 (3), 374–384. doi:10.1038/sj.cdd.4401840
Ye, J., Rawson, R. B., Komuro, R., Chen, X., Davé, U. P., Prywes, R., et al. (2000). ER stress induces cleavage of membrane-bound ATF6 by the same proteases that process SREBPs. Mol. Cell 6 (6), 1355–1364. doi:10.1016/s1097-2765(00)00133-7
Zalloua, P. A., Azar, S. T., Delépine, M., Makhoul, N. J., Blanc, H., Sanyoura, M., et al. (2008). WFS1 mutations are frequent monogenic causes of juvenile-onset diabetes mellitus in Lebanon. Hum. Mol. Genet. 17 (24), 4012–4021. doi:10.1093/hmg/ddn304
Keywords: WFS1 gene, Wolfram syndrome-like disorders (WSLD), ACMG classification, functional analysis, whole-exome sequencing (WES)
Citation: Wu L, Zhang J, Li D, Zhang Z, Ni Q, Han R, Ye L, Zhang Y, Hong J, Wang W, Ning G and Gu W (2024) Novel WFS1 variants are associated with different diabetes phenotypes. Front. Genet. 15:1433060. doi: 10.3389/fgene.2024.1433060
Received: 15 May 2024; Accepted: 31 July 2024;
Published: 16 August 2024.
Edited by:
Maurizio Delvecchio, University of L’Aquila, ItalyReviewed by:
Fabrizio Barbetti, University of Rome Tor Vergata, ItalyLuciana Rigoli, University of Messina, Italy
Copyright © 2024 Wu, Zhang, Li, Zhang, Ni, Han, Ye, Zhang, Hong, Wang, Ning and Gu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Weiqiong Gu, weiqionggu@163.com
†These authors have contributed equally to this work and share first authorship