Major depressive disorder (MDD) is a leading cause of global disability, and a substantial proportion of patients develop treatment-resistant depression (TRD), highlighting the urgent need for mechanistically novel and fast-acting antidepressants. Ketamine has emerged as a paradigm-shifting therapy due to its rapid and robust antidepressant effects, yet its molecular mechanisms in human neurons remain incompletely understood. Here, we investigated the cellular and transcriptomic effects of ketamine in human induced pluripotent stem cell (hiPSC)-derived forebrain neurons. Neuronal cultures were treated with ketamine (1 µM) for 72 hours and analyzed using immunocytochemistry, western blotting, and RNA sequencing. Ketamine treatment did not alter gross neuronal morphology or presynaptic organization but significantly increased synaptic localization of the postsynaptic scaffolding protein PSD-95, suggesting rapid functional strengthening of existing excitatory synapses. At the signaling level, ketamine selectively modulated the Akt pathway without inducing changes in BDNF protein levels or MAPK activation, and without broadly altering excitatory or inhibitory synaptic protein expression. Transcriptomic profiling revealed 604 differentially expressed genes, with significant enrichment of pathways related to serotonergic and dopaminergic neurotransmission. Notably, ketamine upregulated genes involved in monoamine synthesis and signaling, including DDC and DRD1 , while downregulating genes associated with GPCR ligand binding and neurotransmitter transport. SynGO analysis further identified selective regulation of synapse-associated genes, predominantly affecting presynaptic and neurotransmitter uptake mechanisms. Together, these findings indicate that ketamine exerts rapid antidepressant-relevant effects in human neurons by enhancing excitatory synaptic organization and engaging downstream transcriptional programs linked to monoaminergic signaling. This hiPSC-based model provides a mechanistic framework for understanding ketamine’s rapid action and a foundation for future studies on differential responses in TRD patients.

Chromosomal triplications are rare structural variations often associated with complex phenotypes. We report the molecular characterization of a novel intrachromosomal triplication at 18q12.1q21.2 identified in a fetus with ultrasound abnormalities. Conventional karyotyping and array-CGH revealed a partial tetrasomy and a 26 Mb region of loss of homozygosity (LOH), extending from the triplication to the telomere. Long-read sequencing (LRS) identified breakpoint junctions revealing a direct-inverted-direct triplication structure. Breakpoint analysis suggested that this rearrangement arose through a U-type exchange between sister chromatids, likely mediated by microhomology-based mechanisms. This process likely generated a transient dicentric chromosome that subsequently broke during mitosis. The resulting duplicated chromosome may have been stabilized by telomere capture, consistent with the triplicated 18q12.1q21.2 region followed by the 18q21.2q23 LOH. Nine genes within the triplicated region, including SMAD2 and SMAD4, showed high predicted sensitivity to increased dosage, possibly contributing to the clinical phenotype. This study highlights the utility of LRS in defining complex chromosomal rearrangements and emphasizes the importance of molecular breakpoint analysis for understanding pathogenic mechanisms and improving genetic prenatal diagnosis.

NF1 microdeletion syndrome, accounting for 5-11% of NF1 patients, is caused by a deletion in the NF1 region and it is generally characterized by a severe phenotype. Although 70% of NF1 microdeletion patients presents the same 1.4 Mb type-I deletion, some patients may show additional clinical features. Therefore, the contribution of several pathogenic mechanisms, besides haploinsufficiency of some genes within the deletion interval, is expected and needs to be defined. We investigated an altered expression of deletion flanking genes by qPCR in patients with type-1 NF1 deletion, compared to healthy donors, possibly contributing to the clinical traits of NF1 microdeletion syndrome. In addition, the 1.4-Mb deletion leads to changes in the 3D chromatin structure in the 17q11.2 region. Specifically, this deletion alters DNA-DNA interactions in the regions flanking the breakpoints, as demonstrated by our 4C-seq analysis. This alteration likely causes position effect on the expression of deletion flanking genes.Interestingly, 4C-seq analysis revealed that in microdeletion patients, an interaction was established between the RHOT1 promoter and the SLC6A4 gene, which showed increased expression. We performed NGS on putative modifier genes, and identified two "likely pathogenic" rare variants in RAS pathway, possibly contributing to incidental phenotypic features.This study provides new insights into understanding the pathogenesis of NF1 microdeletion syndrome and suggests a novel pathomechanism that contributes to the expression phenotype in addition to haploinsufficiency of genes located within the deletion.This is a pivotal approach that can be applied to unravel microdeletion syndromes, improving precision medicine, prognosis and patients' follow-up.

Spinal neurofibromatosis (SNF), a phenotypic subclass of neurofibromatosis 1 (NF1), is characterized by bilateral neurofibromas involving all spinal roots. In order to deepen the understanding of SNF’s clinical and genetic features, we identified 81 patients with SNF, 55 from unrelated families, and 26 belonging to 19 families with at least 1 member affected by SNF, and 106 NF1 patients aged >30 years without spinal tumors. A comprehensive NF1 mutation screening was performed using NGS panels, including NF1 and several RAS pathway genes. The main features of the SNF subjects were a higher number of internal neurofibromas (p < 0.001), nerve root swelling (p < 0.001), and subcutaneous neurofibromas (p = 0.03), while hyperpigmentation signs were significantly less frequent compared with the classical NF1-affected cohorts (p = 0.012). Fifteen patients underwent neurosurgical intervention. The histological findings revealed neurofibromas in 13 patients and ganglioneuromas in 2 patients. Phenotypic variability within SNF families was observed. The proportion of missense mutations was higher in the SNF cases than in the classical NF1 group (21.40% vs. 7.5%, p = 0.007), conferring an odds ratio (OR) of 3.34 (CI = 1.33–10.78). Two unrelated familial SNF cases harbored in trans double NF1 mutations that seemed to have a subclinical worsening effect on the clinical phenotype. Our study, with the largest series of SNF patients reported to date, better defines the clinical and genetic features of SNF, which could improve the management and genetic counseling of NF1.

Spinal neurofibromatosis (SNF) is a form of neurofibromatosis type 1 (NF1) characterized by bilateral neurofibromas involving all spinal roots. The pathogenic mechanisms determining the SNF form are currently unknown. To verify the presence of genetic variants possibly related to SNF or classic NF1, we studied 106 sporadic NF1 and 75 SNF patients using an NGS panel of 286 genes encoding RAS pathway effectors and neurofibromin interactors and evaluated the expression of syndecans (SDC1, SDC2, SDC3, SDC4), the NF1 3' tertile interactors, by quantitative real-time PCR. We previously identified 75 and 106 NF1 variants in SNF and NF1 cohorts, respectively. The analysis of the distribution of pathogenic NF1 variants in the three NF1 tertiles showed a significantly higher prevalence of NF1 3' tertile mutations in SNF than in the NF1 cohort. We hypothesized a potential pathogenic significance of the 3' tertile NF1 variants in SNF. The analysis of syndecan expression on PBMCs RNAs from 16 SNF, 16 classic NF1 patients and 16 healthy controls showed that the expression levels of SDC2 and SDC3 were higher in SNF and NF1 patients than in controls; moreover, SDC2, SDC3 and SDC4 were significantly over expressed in patients mutated in the 3' tertile compared to controls. Two different mutational NF1 spectra seem to characterize SNF and classic NF1, suggesting a pathogenic role of NF1 3' tertile and its interactors, syndecans, in SNF. Our study, providing new insights on a possible role of neurofibromin C-terminal in SNF, could address effective personalized patient management and treatments.

Over 200 genetic loci have been associated with multiple sclerosis (MS) explaining ~ 50% of its heritability, suggesting that additional mechanisms may account for the “missing heritability” phenomenon. To analyze a large cohort of Italian individuals to identify markers associated with MS with potential functional impact in the disease. We studied 2571 MS and 3234 healthy controls (HC) of continental Italian origin. Discovery phase included a genome wide association study (1727 MS, 2258 HC), with SNPs selected according to their association in the Italian cohort only or in a meta-analysis of signals with a cohort of European ancestry (4088 MS, 7144 HC). Top associated loci were then tested in two Italian cohorts through array-based genotyping (903 MS, 884 HC) and pool-based target sequencing (588 MS, 408 HC). Finally, functional prioritization through conditional eQTL and mQTL has been performed. Top associated signals overlap with already known MS loci on chromosomes 3 and 17. Three SNPs (rs4267364, rs8070463, rs67919208), all involved in the regulation of TBKBP1, were prioritized to be functionally relevant. No evidence of novel signal of association with MS specific for the Italian continental population has been found; nevertheless, two MS loci seems to play a relevant role, raising the interest to further investigations for TBKBP1 gene.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Genome-wide association studies identified over 200 risk loci for multiple sclerosis (MS) focusing on common variants, which account for about 50% of disease heritability. The goal of this study was to investigate whether low-frequency and rare functional variants, located in MS-established associated loci, may contribute to disease risk in a relatively homogeneous population, testing their cumulative effect (burden) with gene-wise tests. We sequenced 98 genes in 588 Italian patients with MS and 408 matched healthy controls (HCs). Variants were selected using different filtering criteria based on allelic frequency and in silico functional impacts. Genes showing a significant burden (n = 17) were sequenced in an independent cohort of 504 MS and 504 HC. The highest signal in both cohorts was observed for the disruptive variants (stop-gain, stop-loss, or splicing variants) located in EFCAB13, a gene coding for a protein of an unknown function (p < 10-4). Among these variants, the minor allele of a stop-gain variant showed a significantly higher frequency in MS versus HC in both sequenced cohorts (p = 0.0093 and p = 0.025), confirmed by a meta-analysis on a third independent cohort of 1298 MS and 1430 HC (p = 0.001) assayed with an SNP array. Real-time PCR on 14 heterozygous individuals for this variant did not evidence the presence of the stop-gain allele, suggesting a transcript degradation by non-sense mediated decay, supported by the evidence that the carriers of the stop-gain variant had a lower expression of this gene (p = 0.0184). In conclusion, we identified a novel low-frequency functional variant associated with MS susceptibility, suggesting the possible role of rare/low-frequency variants in MS as reported for other complex diseases.

Among multiple sclerosis (MS) susceptibility genes, the strongest non-human leukocyte antigen (HLA) signal in the Italian population maps to the TNFSF14 gene encoding LIGHT, a glycoprotein involved in dendritic cell (DC) maturation. Through fine-mapping in a large Italian dataset (4,198 patients with MS and 3,903 controls), we show that the TNFSF14 intronic SNP rs1077667 is the primarily MS-associated variant in the region. Expression quantitative trait locus (eQTL) analysis indicates that the MS risk allele is significantly associated with reduced TNFSF14 messenger RNA levels in blood cells, which is consistent with the allelic imbalance in RNA-Seq reads (P < 0.0001). The MS risk allele is associated with reduced levels of TNFSF14 gene expression (P < 0.01) in blood cells from 84 Italian patients with MS and 80 healthy controls (HCs). Interestingly, patients with MS are lower expressors of TNFSF14 compared to HC (P < 0.007). Individuals homozygous for the MS risk allele display an increased percentage of LIGHT-positive peripheral blood myeloid DCs (CD11c+, P = 0.035) in 37 HCs, as well as in in vitro monocyte-derived DCs from 22 HCs (P = 0.04). Our findings suggest that the intronic variant rs1077667 alters the expression of TNFSF14 in immune cells, which may play a role in MS pathogenesis.

BACKGROUND Triadin is a protein expressed in cardiac and skeletal muscle that has an essential role in the structure and functional regulation of calcium release units and excitation-contraction coupling. Mutations in the triadin gene (TRDN) have been described in different forms of human arrhythmia syndromes with early onset and severe arrhythmogenic phenotype, including triadin knockout syndrome. OBJECTIVE The purpose of this study was to characterize the pathogenetic mechanism underlying a case of severe pediatric malignant arrhythmia associated with a defect in the TRDN gene. METHODS We used a trio whole exome sequencing approach to identify the genetic defect in a 2-year-old boy who had been resuscitated from sudden cardiac arrest and had frequent episodes of ventricular fibrillation and a family history positive for sudden death. We then performed in vitro functional analysis to investigate possible pathogenic mechanisms underlying this severe phenotype. RESULTS We identified a novel homozygous missense variant (p.L56P) in the TRDN gene in the proband that was inherited from the heterozygous unaffected parents. Expression of a green fluorescent protein (GFP)-tagged mutant human cardiac triadin isoform (TRISK32-L56P-GFP) in heterologous systems revealed that the mutation alters protein dynamics. Furthermore, when co-expressed with the type 2 ryanodine receptor, caffeine-induced calcium release from TRISK32-L56P-GFP was relatively lower compared to that observed with the wild-type construct. CONCLUSION The results of this study allowed us to hypothesize a pathogenic mechanism underlying this rare arrhythmogenic recessive form, suggesting that the mutant protein potentially can trigger arrhythmias by altering calcium homeostasis.

Transcriptional changes normally occur during development but also underlie differences between healthy and pathological conditions. Transcription factors or chromatin modifiers are involved in orchestrating gene activity, such as the cohesin genes and their regulator NIPBL. In our previous studies, using a zebrafish model for nipblb knockdown, we described the effect of nipblb loss-of-function in specific contexts, such as central nervous system development and hematopoiesis. However, the genome-wide transcriptional impact of nipblb loss-of-function in zebrafish embryos at diverse developmental stages remains under investigation. By RNA-seq analyses in zebrafish embryos at 24 h post-fertilization, we examined genome-wide effects of nipblb knockdown on transcriptional programs. Differential gene expression analysis revealed that nipblb loss-of-function has an impact on gene expression at 24 h post fertilization, mainly resulting in gene inactivation. A similar transcriptional effect has also been reported in other organisms, supporting the use of zebrafish as a model to understand the role of Nipbl in gene regulation during early vertebrate development. Moreover, we unraveled a connection between nipblb-dependent differential expression and gene expression patterns of hematological cell populations and AML subtypes, enforcing our previous evidence on the involvement of NIPBL-related transcriptional dysregulation in hematological malignancies.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

BACKGROUND:The application of single-cell RNA sequencing (scRNASeq) represents a unique approach to identify hundreds to millions of cells in mammalian cortical multilayers at different stages of embryogenesis. ScRNASeq technology applied to neurological studies requires the use of fresh starting materials because standard cryopreservation methods do not guarantee high viability of cortical primary cells derived from dissected brain areas.NEW METHOD:Here we set up and validate an innovative strategy to perform scRNASeq studies in cryopreserved primary cortical cells isolated from E15.5 mouse embryo. In order to freeze cortical primary cells, we have employed Neurostore, a medium able to guarantee high viability and cell composition of embryonic cortex after thawing.COMPARISON WITH EXISTING METHODS:We showed for the first time the possibility to run scRNASeq experiments on primary cortical cells in an off-line set-up, ensuring cellular integrity and diversity.RESULTS:By trypan blue assay and flow cytometry analysis, we found that Neurostore-cryopreserved cortical cells showed approximately 95 % of viability. Satisfactory RNA recovery and cDNA libraries were achieved. Transcriptome sequencing of 35,763 cryoconserved single cells yielded a robust data-set, identifying 25 cell clusters in three biological samples. Prevalence of peculiar neural populations before and after the cryopreservation-resuscitation procedure was verified by marker gene expression and immunofluorescence analysis.CONCLUSIONS:Our findings support the evidence that frozen primary cortical cells can be successfully employed in scRNASeq experiments allowing an unprecedented flexibility in experimental procedures, such as sample preparation and subsequent processing steps performed in different locations.

Spinal Neurofibromatosis (SNF), a distinct clinical entity of NF1, characterized by bilateral neurofibromas involving all spinal roots and a few, if any, cutaneous manifestations, entails greater morbidity than the classical form of disease. Nevertheless, there are no reliable patterns to sort out patients at risk for developing SNF. We investigated 19 NF1 families with at least one SNF member, 37 sporadic SNF patient and 100 NF1 patients with classical form of disease. We applied Targeted NGS using a panel consisting of 139 genes encoding RAS pathway effectors, neurofibromin interactors and genes mapping at 17q11.2 region. In SNF patients we found a higher percentage of missense (21% versus 8%, p=0 0.016, OR 3.13 (95% CI 01.1 -8.2) and a lower percentage of nonsense NF1 mutations (12.5% versus 28%,, p= 0.026, OR 0.36 (95% CI 0.14–0.9) than in classical NF1 cases. Furthermore, we evaluated rare variants with damaging potential predictors in genes of the RAS pathway and in neurofibromin interactors. In more than one sporadic case possible pathogenic variants were found in LIMK2 (neurofibromin interactor), RASAL1, RASAL3, SOS1, A2ML1, MAP3K1 (RAS pathway effectors), while in more than one SNF family were detected RASAL1, RASAL3, MAP3K1 genes variations. Our results confirm the correlation between NF1 genotype and SNF phenotype as previously reported (Ruggieri, 2015), suggesting that neurofibromin gain-of-functions mutations are associated to SNF. In some patients, the co-occurrence of potential pathogenic variants in NF1 related genes with severe phenotypes was detected supporting their role as modifier genes and promising therapeutic targets.

Noonan syndrome (NS) is an autosomal-dominant disorder with variable expressivity and locus heterogeneity. Despite several RAS pathway genes were implicated in NS, 20-30% of patients remain without molecular diagnosis, suggesting the involvement of further genes or multiple mechanisms. Eight patients out of 60, negative for conventional NS mutation analysis, with heterogeneous NS phenotype were investigated by means of target resequencing of 26 RAS/MAPK pathway genes. A trio was further characterized by means of whole-exome sequencing. Protein modeling and in silico prediction of protein stability allowed to identify possible pathogenic RAS pathway variants in four NS patients. A new c.355T>C variant in LZTR1 was found in patient 43. Two patients co-inherited variants in LRP1 and LZTR1 (patient 53), or LRP1 and SOS1 genes (patient 67). The forth patient (56) carried a compound heterozygote of RASAL3 gene variants and also an A2ML1 variant. While these subclinical variants are singularly present in healthy parents, they co-segregate in patients, suggesting their addictive effect and supporting a digenic inheritance, as alternative model to a more common monogenic transmission. The ERK1/2 and SAPK/JNK activation state, assessed on immortalized lymphocytes from patients 53 and 67 showed highest phosphorylation levels compared to their asymptomatic parents. These findings together with the lack of their co-occurrence in the 1000Genomes database strengthen the hypothesis of digenic inheritance in a subset of NS patients. This study suggests caution in the exclusion of subclinical variants that might play a pathogenic role providing new insights for alternative hereditary mechanisms.