- 1Zentrum für Sozial-und Neuropädiatrie (DBZ), Vivantes Klinikum Neukölln, Berlin, Germany
- 2Institute of Cell Biology and Neurobiology, Charité - Universitätsmedizin Berlin, corporate member of Freie Universität Berlin und Humboldt-Universität zu Berlin, Berlin, Germany
- 3Praxis für Humangenetik, Berlin, Germany
- 4Diagenom GmbH, Rostock, Germany
- 5Zentrum für Humangenetik, SYNLAB MVZ Humangenetik Freiburg, Tübingen, Germany
- 6Zentrum für Humangenetik Tübingen, Tübingen, Germany
- 7CeGaT, Tübingen, Germany
- 8Experimental and Clinical Research Center (ECRC), a cooperation been the Max Delbrück Center for Molecular Medicine in the Helmholtz Association and Charité—Universitätsmedizin Berlin, Berlin, Germany
- 9Institute of Physiology, Brandenburg Medical School (MHB) Theodor Fontane, Brandenburg an der Havel, Germany
- 10Institute for Medical Genetics and Human Genetics, Charité—Universitätsmedizin Berlin, corporate member of Freie Universität Berlin und Humboldt Universität zu Berlin, Berlin, Germany
Introduction: Cohen syndrome (CS) is an early-onset pediatric neurodevelopmental disorder characterized by postnatal microcephaly and intellectual disability. An accurate diagnosis for individuals with CS is crucial, particularly for their caretakers and future prospects. CS is predominantly caused by rare homozygous or compound heterozygous pathogenic variants in the vacuolar protein sorting-associated 13B (VPS13B) gene, which disrupt protein translation and lead to a loss of function (LoF) of the encoded VPS13B protein.
Methods: The widespread incorporation of next-generation sequencing approaches in genetic diagnostics increases the number of individuals carrying VPS13B mutant alleles. At the same time, it increases the detection of variants of unknown clinical significance, necessitating further functional pathogenicity validation.
Results: In this study, we present a family with two CS patients. Within this family, four rare VPS13B variants were detected: c.710G > C, p.Arg237Pro; c.6804delT, p.Phe2268Leufs*24; c.7304C > T, p.Ala2435Val; and c.10302T > A, p.Tyr3434*. These variants challenge the interpretation of their disease-causing role. Specifically, the variants c.6804delT, p.Phe2268Leufs*24 and c.710G > C, p.Arg237Pro were detected in trans configuration and are considered to be causing CS genetically. The functional characterization of the missense variant c.710G > C, p.Arg237Pro shows diminished localization at the Golgi complex, highlighting its clinical relevance and supporting its classification by the American College of Medical Genetics and Genomics (ACMG) as likely pathogenic, class 4.
Discussion: Overall, we emphasize the need for combining genetic and functional testing of VPS13B missense variants to ensure accurate molecular diagnosis and personalized medical care for CS patients.
Introduction
In recent decades, human genetics has enabled the rapid and comprehensive identification of genetic variants and their association with neurological diseases (Satam et al., 2023). However, we emphasize the importance of systematic phenotyping and functional studies, especially when the interpretation of variants of unknown significance (VUS) leaves uncertainty regarding their pathogenic character.
In this study, we present two brothers with the phenotype of Cohen syndrome (CS, MIM: #216550). CS is a rare autosomal recessive genetic condition that is primarily characterized by developmental delay, postnatal microcephaly, and distinctive facial gestalt with wave-shaped eyelids and a short philtrum (Kolehmainen et al., 2003; El Chehadeh-Djebbar et al., 2013). The clinical diagnosis of CS is facilitated by the presence of progressive retinal dystrophy and/or neutropenia. Additional facultative clinical signs include myopia, childhood hypotonia, joint laxity, a cheerful disposition, and autism spectrum disorder (ASD)-like behavior (Katzaki et al., 2007; Gueneau et al., 2014). Genetically, CS is predominantly caused by splice site mutations, frameshift indels, or nonsense variants, which typically results in a complete loss of function (LoF) of the encoded vacuolar protein sorting-associated 13B (VPS13B) (Seifert et al., 2009). VPS13B belongs to the VPS13 family of bridge-like lipid transfer proteins, suggesting that it facilitates lipid transport between adjacent membranes (Vacca et al., 2024; Mochida et al., 2004; Hanna et al., 2023). At the cellular level, VPS13B is primarily localized to the Golgi complex. It plays a crucial role in preserving the morphology and function of the Golgi, potentially influencing membrane organization and intracellular membrane transport (Vacca et al., 2024; Seifert et al., 2011). Recent studies using neuronal models, including iPSC-derived neurons and animal models of CS, have indicated that VPS13B depletion is linked to neurodevelopmental abnormalities, such as axonal elongation defects, defective synaptogenesis, and increased autophagic flux (Seifert et al., 2015; Lee et al., 2020; Lee et al., 2020). However, a clear genotype–phenotype correlation or functional correlation regarding the differences in the clinical expression of CS has not yet been established.
In contrast, the interpretation of VPS13B missense variants is challenging due to the lack of biochemical and pathomechanistic characterization of most VPS13B protein regions. Although in silico analysis supports predictions of pathogenicity, functional evidence provides much stronger validation for the pathogenic character of a genetic variant. To overcome this limitation, we established a cellular detection assay that measures the subcellular distribution of the encoded VPS13B protein in the Golgi apparatus and the cytosol after immunostaining (Zorn et al., 2022). In this study, we employed a quantitative subcellular analysis of VPS13B mutant protein distribution to classify the familial VPS13B missense VUS, thereby validating the genetic diagnosis for CS.
Methods
Family consent
Written informed consent for genetic testing and the use of data for publication was obtained from each participant. Parental consent was secured for both siblings under 18 years of age.
Genetic testing
DNA was isolated from peripheral EDTA blood. Both patients underwent initial Sanger sequencing of the entire coding region of the VPS13B gene (exon 2–62, including 20 bp of the 3′ and 5′ intron regions) using standard procedures with a 3130XL Genetic Analyzer (Applied Biosystems). All details can be obtained upon request. Quadro analysis of the VPS13B gene using NGS-based whole exome sequencing was performed on both patients and their parents to reevaluate the variants, technically confirm allelic distribution, and ascertain the de novo status. Enrichment of the coding regions, adjacent intronic regions, and other non-coding, disease-relevant regions of the VPS13B candidate region was performed using in-solution hybridization. For exome capture, we used enrichment kits from Twist Bioscience. High-throughput sequencing was conducted using the Illumina NovaSeq 6,000/NovaSeq X Plus System. Adapter sequences were removed using Skewer, a tool designed for trimming high-throughput sequencing reads. The resulting sequences were aligned to the human reference genome (hg19) using the Burrows–Wheeler Aligner sequences. The sequences that were not uniquely aligned or were identified as duplicates, most likely due to PCR amplification, were removed to ensure data accuracy. Sequencing data were processed using Illumina bcl2fastq2 software, which converts raw BCL files generated by the sequencer into FASTQ files for further analysis. The remaining high-quality sequences were annotated for sequence variants using various internal and external databases, including the Illumina database and the VPS13B MANE PLUS CLINICAL transcript ENST00000358544, NM_017890.5, comprising 12,069 bp. However, in this study, all variants were aligned to the ubiquitously expressed and conserved variant VPS13B ENST00000357162.7, NM_152564.5, comprising 11,994 bp and encoding 3,997 amino acids (aa) as the reference.
Reagents and antibodies
All reagents were obtained from Sigma-Aldrich, Roth or Merck, unless stated otherwise. The following commercial antibody was used: Mouse anti-GM130 (BD Transductions Laboratories, Cat. No. 610822). Rabbit anti-VPS13B (442) was described earlier (Seifert et al., 2011). Secondary antibodies for immunofluorescence included donkey anti-mouse-Cy3, anti-rabbit-Cy3, and anti-rabbit-Cy2 (all from Dianova GmbH). 6-Diamidino-2-phenylindole (DAPI) (Invitrogen) was used for nuclear DNA staining.
Cloning of the VPS13B variant expression vectors by site-directed mutagenesis
Missense variants were annotated using the VPS13B transcript [Ensembl, ENST00000357162.6 (VPS13B-202); RefSeq, NM_152564.5]. For transient expression, the full-length human VPS13B construct encoding the VPS13B missense variant p.Arg237Pro (pcDNA3.1_VPS13B-mutR237P) was cloned using site-directed mutagenesis PCR (Reikofski and Tao, 1992). Briefly, the full-length wild-type human construct pcDNA3.1_VPS13B-wt (Seifert et al., 2011) was amplified using primer pairs containing the corresponding missense variants in VPS13B using PrimeSTAR GXL DNA Polymerase (NEB Inc.) according to the manufacturer’s instructions. The PCR products were subsequently treated with the KLD enzyme mix (NEB Inc.), which eliminated the methylated pcDNA3.1_VPS13B-wt vectors and prepared them for the DNA ligation reaction. The integrity and correct mutagenesis of all cloned vectors were confirmed through Sanger sequencing.
Cell culture and transfection
HeLa-cells (ATCC) were cultured at 37°C in a 5% CO2 atmosphere using DMEM culture medium supplemented with 5% FCS and 1% glutamine. Transient transfection of the expression vectors was performed at approximately 50% confluence using Lipofectamine (Invitrogen, Thermo Fisher) according to the manufacturer’s instructions. The variant of interest was pcDNA3.1_VPS13B-mutArg237Pro, while pcDNA3.1_VPS13B-mutAla590Thr and pcDNA3.1_hVPS13B-mutGlu2704Arg served as negative and positive controls, respectively. Cell harvesting for further analysis was performed 24–48 h post-transfection.
Immunofluorescence
For visualization of intracellular protein expression and localization, the cells were cultured on coverslips and, if required, transfected (see above). The cells on the coverslips were fixed using 3% (w/v) paraformaldehyde in 1x PBS at 4°C and permeabilized in 0.5% (v/v) Triton X-100 and 1% (w/v) BSA in 1x PBS. Primary antibodies were incubated in 1% (w/v) BSA diluted in 1x PBS for at least 8 h at 4°C. After rinsing with PBS, secondary antibodies were incubated in 1% (w/v) BSA diluted in 1x PBS for 2 h at 4°C. Subsequently, the coverslips were attached to slides using ImmuMount (ThermoFisher). The visualization was performed using a spinning disk microscope (Zeiss), with a comparable cellular imaging level and identical imaging conditions.
ImageJ and statistical analysis
Image analysis and Golgi enrichment were performed using ImageJ. Two different regions of interest (ROIs) were defined for each cell separately, as previously described (Zorn et al., 2022). Briefly, one ROI outlining the cell border measured the total VPS13B immunofluorescence of the cell (total cell ROI), while the second ROI outlining the GM130-positive Golgi structure (Golgi ROI) measured the Golgi-associated immunofluorescence of VPS13B. The percentage of Golgi-associated VPS13B fluorescence compared to the total VPS13 cell fluorescence intensity was calculated. The graphical representation and statistical analysis were carried out using GraphPad Prism software. The VPS13B variants pcDNA3.1_VPS13B-mutAla590Thr and pcDNA3.1_hVPS13B-mutGlu2704Arg were previously described (Zorn et al., 2022) and included in the study as Golgi enrichment-affecting negative and positive controls for the VPS13B missense variant, respectively.
In silico prediction and variant classification for the evaluation of pathogenicity
All variant pathogenicity prediction algorithms and the ACMG classification (ACMG/ACGS 2020 guidelines, version 4.01) were conducted as previously described (Zorn et al., 2022). Truncating VPS13B variants were classified as decribed (Abou Tayoun et al., 2018).
Results
Patients and family
Here, we report the phenotype of two brothers with clinical manifestations indicating CS.
Patient 1 was born at term (41 weeks) via spontaneous delivery, with a relatively low birth weight of 3,030 g (P16, −1.01 SD). He exhibited respiratory adaptation disorder and muscular hypotonia. A cranial MRI scan at birth revealed a structurally unremarkable brain with slightly enlarged cerebrospinal fluid spaces. After birth, the patient developed postnatal microcephaly, with a head circumference of 41 cm at the age of 1 (<P1, −5.19 SD), and a psychomotor developmental delay. All developmental milestones were met at delayed time points. The patient was able to sit independently at the age of 1.5 years, walked independently by the age of 4, and has remained non-verbal to date. Myopia was diagnosed at the age of 4. He is usually in a friendly mood and shows limited understanding of simple content. Moreover, within the first year, the patient was noted to have recurrent bronchial infections and neutropenia. Other symptoms included feeding difficulties, small genitalia, and bilateral undescended testis.
Patient 2, the younger brother, was born via spontaneous delivery at term after a pregnancy complicated by maternal diabetes requiring insulin therapy. At 6 months of age, developmental delay, muscular hypotonia, and feeding difficulties became evident. He developed distinct postnatal microcephaly and, despite intensive supportive therapy, showed slow developmental progress. He was able to sit unaided at 21 months and has been walking with support since the age of 3 years. Sleeping disorders and feeding difficulties have been ongoing to this day. The boy frequently shows high irritability and exhibits autistic features, including reduced eye contact, stereotypical movements, and social withdrawal. His mother reports that he often appears to be in his own world and lacks a recognizable need for social contact, which is in line with previous reports of autistic features in CS (Table 1).
Sequencing results
Chromosomal analysis and microarray comparative genomic hybridization (CGH) were unremarkable, indicating no major chromosomal abnormalities. Initial Sanger sequencing of the VPS13B gene was performed based on the clinically distinct phenotype of CS in Patient 1. Further molecular testing of both patients and the parents was performed using candidate region in-solution hybridization enrichment and high-throughput sequence analysis using the Illumina NovaSeq 6000/NovaSeq X Plus System. Next-generation sequencing-based copy number variations were calculated based on the sequences that could be assigned to a genomic position using an internally developed method based on sequencing depth. No genomic copy number imbalances were detected that could explain the patients’ phenotype. Further annotation of the high-throughput sequencing data using various internal and external databases for the sequence variants in the candidate region of VPS13B confirmed two pathogenic heterozygous variants and two heterozygous VUS in the older patient (Patient 1): c.6804delT, p.Phe2268fs*24 (Rafiq et al., 2015); c.10302T > A, p.Tyr3434*; c.710G > C, p. Arg237Pro; and c.7304C > T, p.Ala2435Val. In the younger patient (Patient 2), sequencing of the candidate region of VPS13B identified one pathogenic heterozygous variant and two heterozygous VUS: c.6804delT, p.Phe2268fs*24 (Rafiq et al., 2015); c.710G > C, p. Arg237Pro; and c.7304C > T, p.Ala2435Val. Subsequent analysis of the parents suggested that the mother carried the pathogenic variant c.6804delT, p.Phe2268fs*24 and the VUS c.7304C > T, p.Ala2435Val, which were inherited by both brothers. The second pathogenic variant, c.10302T > A, p.Tyr3434*, in Patient 1, was not inherited from the mother or the father, suggesting a de novo agent, germ cell mosaicism, or a somatic event. The father carried the heterozygous VPS13B VUS c.710G > C, p. Arg237Pro, which was inherited by both brothers (Figure 1; Table 2).
Figure 1. Pedigree and genetic findings in the family with CS. (A) Squares and circles symbolize males and females, respectively. Filled symbols indicate affected individuals, while unfilled symbols represent healthy individuals. The VPS13B genotypes are depicted below each symbol (n.a., no DNA available). (B) A partial screenshot of the VPS13B gene region on chr8:100,123,423–100,123,486 from the UCSC genome browser shows the conservation plot of the VPS13B p.Arg237 variant (highlighted in a red box), with conservation extending down to Xenopus tropicalis.
Table 2. Variant pathogenicity interpretation was conducted according to the ACMG standards and guidelines (Richards et al., 2015).
Variant pathogenicity prediction algorithms
To understand the complex genetic and clinical data, we interpreted the VPS13B variants according to the ACMG standards and guidelines (Deignan et al., 2019). The Phe2268Leufs*24 variant has been previously reported to cause CS (Table 2, ClinVarID: 817631) (Rafiq et al., 2015; Zare Ashrafi et al., 2023). The p.Tyr3434* variant was detected in Patient 1, and it is classified as a pathogenic variant due to its deleterious nature. However, this variant did not co-segregate with the CS phenotype within the family, likely indicating a de novo mutation, germ cell mosaicism, or a somatic event. Paternity and maternity for both patients were confirmed. The missense variants were further processed using standard computational and predictive in silico algorithms. In this study, the major difference between both missense variants was predicted using SNPs&GO analyses (Calabrese et al., 2009), which indicated a disease-causing prediction for p.Arg237Pro and a neutral character for p.Ala2435Val (Table 2). The disease-causing prediction of pArg237Pro was further supported by its conservation down to Xenopus tropicalis (Figure 1B), while p.Ala2435Val was excluded as a cause of the CS phenotype due to its co-segregation as a complex maternal allele with p.Phe2268Leufs*24. Overall, the segregation of p.Arg237Pro with the CS phenotype suggested that this missense variant caused autosomal recessive CS in both affected brothers, along with the p.Phe2268Leufs*24 variant.
Golgi localization
To date, the analysis of cellular VPS13B distribution is the only simple and reliable method to investigate the effects of VPS13B missense variants (Zorn et al., 2022). This method involves cloning the respective missense variant into a mammalian expression vector (pcDNA3.1_VPS13B) and subsequent overexpression in standard cell cultures, such as HeLa cells (Zorn et al., 2022). Using this approach, we showed that the wild-type VPS13B protein localizes to the Golgi complex, as indicated by its co-localization with the Golgi matrix protein GM130 (Figure 2A). As controls, we analyzed two previously reported VPS13B missense variants regarding their subcellular localization (Zorn et al., 2022). The p.Ala590Thr variant exhibited a subcellular distribution similar to the wild-type VPS13B protein and was therefore used as a negative control. In contrast, the p.Gly2704Arg variant showed reduced localization to the Golgi complex and was used as a positive control, indicating a disruption in subcellular targeting. The analysis of the p.Arg237Pro variant revealed a significant reduction in Golgi enrichment of the mutated VPS13B protein (Figures 2A,B). It is important to note that the overall area covered by the Golgi complex was not affected; therefore, all variants are comparable to the wild-type VPS13B protein in this respect (Figure 2C).
Figure 2. Golgi localization. (A) HeLa cells were transfected with different pcDNA3.1_VPS13B constructs, encoding either wild-type (wt) or mutated VPS13B (missense variants p.Arg237Pro, p.Ala590Thr, p.Glu2704Arg). Twenty-four hours post-transfection, all cells were processed for immunostaining of VPS13B (green) and GM130 (red). Imaging was performed using a Zeiss spinning disk microscope. Nuclei were stained with DAPI (blue). The scale bar represents 20 μm. (B) Quantitative analysis of the subcellular Golgi enrichment of VPS13B was conducted using ImageJ. Two ROIs were defined for each analyzed cell: one total cell ROI, drawn around the cell area, and one Golgi ROI, drawn around the GM130-positive Golgi area. Statistical analysis was conducted on the Golgi enrichment by calculating the percentage of the VPS13B fluorescence intensity at the Golgi compared to the total cellular fluorescence intensity. (C) Analysis of the Golgi area of VPS13B was conducted using ImageJ. Both ROIs were used to measure the area. Statistical analysis was performed by calculating the percentage of the occupied Golgi area compared to the total cellular area.
Discussion
CS is a rare neurodevelopmental autosomal recessive condition mainly characterized by developmental delay and postnatal microcephaly. In this study, we present a family with two affected brothers who exhibit the obligate clinical CS signs, such as postnatal microcephaly, developmental delay, and distinctive facial gestalt. When we compared the two brothers, we found divergent clinical signs: for example, certain autistic features were exhibited by Patient 2 whereas myopia and neutropenia were detected in Patient 1 only. Genetic testing of VPS13B in the family revealed two LoF variants, p.Phe2268fs*and p.Tyr3434*, as well as two missense variants, p.Arg237Pro and p.Ala2435Val. The heterozygous variant p.Phe2268fs*24 was found to be present in both Patient 1 and Patient 2, which was inherited from their phenotypically unaffected mother. This variant was recently associated with CS and classified as pathogenic (class 5) according to the ACMG criteria (Deignan et al., 2019). The heterozygous variant p.Ala2435Val was also found to be present in both Patient 1 and Patient 2, which was inherited from their phenotypically unaffected mother. The co-segregation of p.Ala2435Val with the heterozygous p.Phe2268fs*24 in all three individuals suggested the presence of a complex VPS13B allele. Based on this, p.Ala2435Val was classified as a VUS (class 3) according to the ACMG guidelines, and current evidence suggests that it may not have a pathogenic impact. Patient 1 had the heterozygous LoF variant p.Tyr3434*, which had not been previously described. This variant appeared de novo in Patient 1 and was classified as pathogenic (class 5) due to the presence of a premature stop codon. The rare heterozygous variant p.Arg237Pro was found to be present in both Patient 1 and Patient 2 and was inherited from their phenotypically unaffected father. So far, this variant has been annotated in ClinVar as a VUS. Our analysis activated the ACMG criteria PM2, PP3, and PM3 for this variant. Importantly, diminished Golgi enrichment in the subcellular analysis of VPS13B p.Arg237Pro allowed for the additional activation of the PS3 term, resulting in the ACMG classification as likely pathogenic, class 4. In addition, a comparable substitution at the same amino acid position (p.Arg237His) has previously been described as a variant of uncertain significance in a compound heterozygous context, based on a trio exome analysis of a patient with CS (Liu et al., 2021). The combination of genetic and functional analyses helps to elucidate how different variants affect protein function, which will be crucial in the future for explaining phenotypic variability. Overall, genetic analysis suggested that CS in both Patients 1 and 2 was caused by the compound heterozygous VPS13B variants p.Phe2268fs*24 and p.Arg237Pro. We speculated that the de novo event p.Tyr3434* dominated the missense variant p.Arg237Pro in Patient 1, resulting in a complete VPS13B LoF and, consequently, a CS phenotype that includes neutropenia. Furthermore, we hypothesized that an incomplete LoF of VPS13B, due to exon-skipping splice variants, larger in-frame indels, and missense variants of uncertain significance (VUS), likely contributes to the phenotypic complexity and heterogeneity of both CS and VPS13B-associated ASD-like behavior (Zorn et al., 2022; Afridi et al., 2024; AbdelAleem et al., 2023). This hypothesis was supported by the altered CS phenotype observed in Patient 2 with ASD-like behavior. Together, this also suggests that establishing a phenotype–genotype correlation will improve patient management by anticipating specific clinical manifestations.
VPS13B missense variants have been rarely linked to CS due to low detection rates and limited biological evidence. A total of 4,973 coding missense variants have been reported to date, of which 109 have been found in homozygosity. The clinical significance of these variants in ClinVar is categorized as follows: 27 classified as benign/likely benign, 42 with conflicting interpretations, 28 as variants of uncertain significance (VUS), and 12 without sufficient documentation. According to data from the patient cohort, literature, and ClinVar, 42 missense variants of VPS13B have been associated with CS so far (Zorn et al., 2022). VPS13B missense variants have also been linked to autism spectrum disorders (ASDs), as noted in a previous study (Yu et al., 2013). The connection between VPS13B and ASD, particularly in patients with CS, raises questions about whether some variants may contribute to ASD-like behavior (Rafiq et al., 2015; Ionita-Laza et al., 2014; Howlin, 2001). Further exploration of this possibility could help identify more specific clinical subtypes of CS, potentially allowing for personalized care and management strategies. Notably, the majority of established CS-causing VSP13B missense variants cluster within or near the Vps13 adaptor binding (VAB) domain, which facilitates the membrane recruitment of the VPS13B C-terminus to the Golgi complex (Zorn et al., 2022; Levine, 2022). The cell biological contribution of the previously described N-terminal missense variant VPS13B p.Arg237Pro, which localizes to the highly conserved N-terminal VPS13 region (VPS13_N, also known as RBG2), to CS and the ADS-like behavior is currently unclear (Levine, 2022; Velayos-Baeza et al., 2004). The VPS13B region around the position p.Arg237 is highly conserved in vertebrates. The human VPS13_N domain has not been characterized; however, a study utilizing cryo-electron microscopy analyzed the VPS13 N-terminal regions of the fungi Chaetomium thermophilum and yeast Saccharomyces cerevisiae (Li et al., 2020). The VPS13 N-terminal region was postulated to mediate the lipid membrane attachment toward the donor lipid membrane. Mutation of critical hydrophobic amino acids into charged amino acids in the Vps13p N-terminal region of S. cerevisiae reduced yeast sporulation, indicating defective lipid transport (Li et al., 2020). Overall, genetic analysis of VPS13B p.Arg237Pro suggested that the VPS13_N domain is critical for VPS13B function, impacts VPS13B Golgi enrichment, and may affect the stability of the VPS13B lipid membrane complex.
In conclusion, while phenotyping and genetic testing are powerful tools for the initial identification of genetic variants, functional studies are indispensable for the accurate interpretation of missense variants. By combining these strengths, especially in the era of NGS, rigorous functional analysis will enhance our ability to diagnose accurately and consequently improve our understanding of molecular pathomechanisms. Furthermore, this study highlights the importance of integrating accurate genetic testing, comprehensive evaluation, and functional analyses to elucidate the complex phenotype–genotype correlation and the pathogenic mechanisms underlying CS. Such integration not only improves accurate diagnosis but also helps identify new therapeutic targets and improve personalized patient care. To increase understanding of CS and related neurodevelopmental disorders, future research should focus on larger cohort studies to better assess the prevalence and spectrum of VPS13B variants. Functional assays, particularly in vitro and in vivo models, are essential for validating the effects of newly identified missense variants on VSP13B function and may help uncover the molecular mechanisms driving phenotypic variability. In addition, longitudinal studies, including standardized clinical assessments and precise phenotyping of patients with CS, may assist in identifying potential genetic modifiers.
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
Ethical approval was not required for the study involving human samples in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the each adult and participants’ legal guardians.
Author contributions
GS: Investigation, Writing – original draft. CM: Formal analysis, Writing – original draft, Data curation. JKo: Data curation, Formal analysis, Writing – original draft. ES: Investigation, Writing – original draft. ZM: Data curation, Methodology, Writing – original draft. JK: Writing – review & editing. NH: Data curation, Investigation, Methodology, Writing – original draft. JKü: Formal analysis, Validation, Writing – original draft, Writing – review & editing. C-EO: Investigation, Writing – original draft. WS: Conceptualization, Formal analysis, Funding acquisition, Investigation, Project administration, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. CM received an ERASMUS scholarship for the duration of the project. For cell biological research on VPS13B, WS was funded by a DFG grant (SE2266_1-1).
Acknowledgments
We would like to thank Magdalena Feldhahn, CeGaT, Tübingen, Germany, for bioinformatical support.
Conflict of interest
ES was employed by company Praxis für Humangenetik. ZM was employed by company Diagenom GmbH. JK was employed by company Zentrum für Humangenetik, SYNLAB MVZ Humangenetik Freiburg. NH was employed by companies Zentrum für Humangenetik Tübingen and CeGaT.
The remaining 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.
References
AbdelAleem, A., Haddad, N., Al-Ettribi, G., Crunk, A., and Elsotouhy, A. (2023). Cohen syndrome and early-onset epileptic encephalopathy in male triplets: two disease-causing mutations in VPS13B and NAPB. Neurogenetics 24, 103–112. doi: 10.1007/s10048-023-00710-2
Abou Tayoun, A. N., Pesaran, T., DiStefano, M. T., Oza, A., Rehm, H. L., Biesecker, L. G., et al. (2018). Recommendations for interpreting the loss of function PVS1 ACMG/AMP variant criterion. Hum. Mutat. 39, 1517–1524. doi: 10.1002/humu.23626
Adzhubei, I., Jordan, D. M., and Sunyaev, S. R. (2013). Predicting functional effect of human missense mutations using PolyPhen-2. Curr. Prot. 7:Unit7.20. doi: 10.1002/0471142905.hg0720s76
Afridi, T. U. K., Fatima, A., Satti, H. S., Akram, Z., Yousafzai, I. K., Naeem, W. B., et al. (2024). Exome sequencing in four families with neurodevelopmental disorders: genotype-phenotype correlation and identification of novel disease-causing variants in VPS13B and RELN. Mol. Gen. Genomics. 299:55. doi: 10.1007/s00438-024-02149-y
Calabrese, R., Capriotti, E., Fariselli, P., Martelli, P. L., and Casadio, R. (2009). Functional annotations improve the predictive score of human disease-related mutations in proteins. Hum. Mutat. 30, 1237–1244. doi: 10.1002/humu.21047
Chen, S., Francioli, L. C., Goodrich, J. K., Collins, R. L., Kanai, M., Wang, Q., et al. (2024). A genomic mutational constraint map using variation in 76,156 human genomes. Nature 625, 92–100. doi: 10.1038/s41586-023-06045-0
Deignan, J. L., Chung, W. K., Kearney, H. M., Monaghan, K. G., Rehder, C. W., Chao, E. C., et al. (2019). Points to consider in the reevaluation and reanalysis of genomic test results: a statement of the American College of Medical Genetics and Genomics (ACMG). Genet. Med. 21, 1267–1270. doi: 10.1038/s41436-019-0478-1
El Chehadeh-Djebbar, S., Blair, E., Holder-Espinasse, M., Moncla, A., Frances, A. M., Rio, M., et al. (2013). Changing facial phenotype in Cohen syndrome: towards clues for an earlier diagnosis. Eur. J. Hum. Genet. 21, 736–742. doi: 10.1038/ejhg.2012.251
Grantham, R. (1974). Amino acid difference formula to help explain protein evolution. Science 185, 862–864. doi: 10.1126/science.185.4154.862
Gueneau, L., Duplomb, L., Sarda, P., Hamel, C., Aral, B., Chehadeh, S. E., et al. (2014). Congenital neutropenia with retinopathy, a new phenotype without intellectual deficiency or obesity secondary to VPS13B mutations. Am. J. Med. Genet. A 164A, 522–527. doi: 10.1002/ajmg.a.36300
Hanna, M., Guillén-Samander, A., and De Camilli, P. (2023). RBG motif bridge-like lipid transport proteins: structure, functions, and open questions. Annu. Rev. Cell Dev. Biol. 39, 409–434. doi: 10.1146/annurev-cellbio-120420-014634
Howlin, P. (2001). Autistic features in Cohen syndrome: a preliminary report. Dev. Med. Child Neurol. 43, 692–696. doi: 10.1111/j.1469-8749.2001.tb00143.x
Ionita-Laza, I., Capanu, M., De Rubeis, S., McCallum, K., and Buxbaum, J. D. (2014). Identification of rare causal variants in sequence-based studies: methods and applications to VPS13B, a gene involved in Cohen syndrome and autism. PLoS Genet. 10:e1004729. doi: 10.1371/journal.pgen.1004729
Katzaki, E., Pescucci, C., Uliana, V., Papa, F. T., Ariani, F., Meloni, I., et al. (2007). Clinical and molecular characterization of Italian patients affected by Cohen syndrome. J. Hum. Genet. 52, 1011–1017. doi: 10.1007/s10038-007-0208-4
Kolehmainen, J., Black, G. C., Saarinen, A., Chandler, K., Clayton-Smith, J., Träskelin, A. L., et al. (2003). Cohen syndrome is caused by mutations in a novel gene, COH1, encoding a transmembrane protein with a presumed role in vesicle-mediated sorting and intracellular protein transport. Am. J. Hum. Genet. 72, 1359–1369. doi: 10.1086/375454
Kumar, P., Henikoff, S., and Ng, P. C. (2009). Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat. Protoc. 4, 1073–1081. doi: 10.1038/nprot.2009.86
Lee, Y. K., Hwang, S. K., Lee, S. K., Yang, J. E., Kwak, J. H., Seo, H., et al. (2020). Cohen syndrome patient iPSC-derived Neurospheres and forebrain-like glutamatergic neurons reveal reduced proliferation of neural progenitor cells and altered expression of synapse genes. J. Clin. Med. 9:1886. doi: 10.3390/jcm9061886
Lee, Y. K., Lee, S. K., Choi, S., Huh, Y. H., Kwak, J. H., Lee, Y. S., et al. (2020). Autophagy pathway upregulation in a human iPSC-derived neuronal model of Cohen syndrome with VPS13B missense mutations. Mol. Brain 13:69. doi: 10.1186/s13041-020-00611-7
Levine, T. P. (2022). Sequence analysis and structural predictions of lipid transfer bridges in the repeating Beta groove (RBG) superfamily reveal past and present domain variations affecting form, function and interactions of VPS13, ATG2, SHIP164, hobbit and Tweek. Contact (Thousand Oaks) 5:251525642211343. doi: 10.1177/25152564221134328
Li, P., Lees, J. A., Lusk, C. P., and Reinisch, K. M. (2020). Cryo-EM reconstruction of a VPS13 fragment reveals a long groove to channel lipids between membranes. J. Cell Biol. 219:1161. doi: 10.1083/jcb.202001161
Liu, Y., Liu, X., Qin, D., Zhao, Y., Cao, X., Deng, X., et al. (2021). Clinical utility of next-generation sequencing for developmental disorders in the rehabilitation department: experiences from a single Chinese center. J. Mol. Neurosci. 71, 845–853. doi: 10.1007/s12031-020-01707-4
Mochida, G. H., Rajab, A., Eyaid, W., Lu, A., Al-Nouri, D., Kosaki, K., et al. (2004). Broader geographical spectrum of Cohen syndrome due to COH1 mutations. J. Med. Genet. 41:e87. doi: 10.1136/jmg.2003.014779
Pejaver, V., Urresti, J., Lugo-Martinez, J., Pagel, K. A., Lin, G. N., Nam, H. J., et al. (2020). Inferring the molecular and phenotypic impact of amino acid variants with MutPred2. Nat. Commun. 11:5918. doi: 10.1038/s41467-020-19669-x
Rafiq, M. A., Leblond, C. S., Saqib, M. A., Vincent, A. K., Ambalavanan, A., Khan, F. S., et al. (2015). Novel VPS13B mutations in three large Pakistani Cohen syndrome families suggests a Baloch variant with autistic-like features. BMC Med. Genet. 16:41. doi: 10.1186/s12881-015-0183-0
Reikofski, J., and Tao, B. Y. (1992). Polymerase chain reaction (PCR) techniques for site-directed mutagenesis. Biotechnol. Adv. 10, 535–547. doi: 10.1016/0734-9750(92)91451-J
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, 405–424. doi: 10.1038/gim.2015.30
Satam, H., Joshi, K., Mangrolia, U., Waghoo, S., Zaidi, G., Rawool, S., et al. (2023). Next-generation sequencing technology: current trends and advancements. Biology (Basel) 12:997. doi: 10.3390/biology12070997
Schwarz, J. M., Cooper, D. N., Schuelke, M., and Seelow, D. (2014). MutationTaster2: mutation prediction for the deep-sequencing age. Nat. Methods 11, 361–362. doi: 10.1038/nmeth.2890
Seifert, W., Holder-Espinasse, M., Kühnisch, J., Kahrizi, K., Tzschach, A., Garshasbi, M., et al. (2009). Expanded mutational spectrum in Cohen syndrome, tissue expression, and transcript variants of COH1. Hum. Mutat. 30, E404–E420. doi: 10.1002/humu.20886
Seifert, W., Kuhnisch, J., Maritzen, T., Horn, D., Haucke, V., and Hennies, H. C. (2011). Cohen syndrome-associated protein, COH1, is a novel, giant Golgi matrix protein required for Golgi integrity. J. Biol. Chem. 286, 37665–37675. doi: 10.1074/jbc.M111.267971
Seifert, W., Kuhnisch, J., Maritzen, T., Lommatzsch, S., Hennies, H. C., Bachmann, S., et al. (2015). Cohen syndrome-associated protein COH1 physically and functionally interacts with the small GTPase RAB6 at the Golgi complex and directs neurite outgrowth. J. Biol. Chem. 290, 3349–3358. doi: 10.1074/jbc.M114.608174
Vacca, F., Yalcin, B., and Ansar, M. (2024). Exploring the pathological mechanisms underlying Cohen syndrome. Front. Neurosci. 18:1431400. doi: 10.3389/fnins.2024.1431400
Velayos-Baeza, A., Vettori, A., Copley, R. R., Dobson-Stone, C., and Monaco, A. P. (2004). Analysis of the human VPS13 gene family. Genomics 84, 536–549. doi: 10.1016/j.ygeno.2004.04.012
Yu, T. W., Chahrour, M. H., Coulter, M. E., Jiralerspong, S., Okamura-Ikeda, K., Ataman, B., et al. (2013). Using whole-exome sequencing to identify inherited causes of autism. Neuron 77, 259–273. doi: 10.1016/j.neuron.2012.11.002
Zare Ashrafi, F., Akhtarkhavari, T., Fattahi, Z., Asadnezhad, M., Beheshtian, M., Arzhangi, S., et al. (2023). Emerging epidemiological data on rare intellectual disability syndromes from analyzing the data of a large Iranian cohort. Arch. Iran. Med. 26, 186–197. doi: 10.34172/aim.2023.29
Keywords: VPS13B, Cohen syndrome, missense variant, functional testing, Golgi complex
Citation: Schottmann G, Martínez Almudéver C, Knop JCM, Suk EK, Meyer Z, Kohlhase J, Himmelreich N, Kühnisch J, Ott C-E and Seifert W (2024) Impact of genetic test interpretation on a VPS13B missense variant in Cohen syndrome. Front. Neurosci. 18:1488133. doi: 10.3389/fnins.2024.1488133
Edited by:
Kevin Peikert, University Hospital Rostock, GermanyReviewed by:
Saurabh Srivastav, Texas Children’s Hospital, United StatesMuhammad Ansar, Fondation Asile des Aveugles, Switzerland
Copyright © 2024 Schottmann, Martínez Almudéver, Knop, Suk, Meyer, Kohlhase, Himmelreich, Kühnisch, Ott and Seifert. 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: Wenke Seifert, d2Vua2Uuc2VpZmVydEBjaGFyaXRlLmRl