Systemic sclerosis (SSc) is a rheumatic disease characterised by vasculopathy, inflammation and fibrosis. Its aetiopathogenesis is still unknown, and the pathways/mechanisms of the disease are not clarified. This study aimed to perform in silico analysis of the already Mass Spectrometry (MS)-based discovered biomarkers of SSc to extract possible pathways/mechanisms implicated in the disease. We recorded all published candidate MS-based found biomarkers related to SSc. We then selected a number of the candidate biomarkers using specific criteria and performed pathway and cellular component analyses using Enrichr. We used PANTHER and STRING to assess the biological processes and the interactions of the recorded proteins, respectively. Pathway analysis extracted several pathways that are associated with the three different stages of SSc pathogenesis. Some of these pathways are also related to other diseases, including autoimmune diseases. We observe that these biomarkers are located in several cellular components and implicated in many biological processes. STRING analysis showed that some proteins interact, creating significant clusters, while others do not display any evidence of an interaction. All these data highlight the complexity of SSc, and further investigation of the extracted pathways/biological processes and interactions may help study the disease from a different angle.

Distal hereditary motor neuronopathies (dHMNs) form a heterogeneous group of rare disorders characterized by distal weakness and wasting in the limbs with no significant sensory involvement. Harding has classified dHMNs into seven categories based on clinical and genetic criteria. We report a novel form of autosomal recessive dHMN in 7 consanguineous families located in the Jerash region of Jordan. Onset of the disease is between 6 and 10 years of age and is characterized by weakness and atrophy of the lower limbs associated with pyramidal features. Within 2 years, symptoms progress to the upper limbs. Neurophysiological studies typically show normal conduction velocities, reduced compound motor action potential amplitudes, normal sensory nerve action potentials, and chronic neurogenic changes on needle electromyography. No significant abnormalities are seen on sural nerve biopsy. We call this novel form of dHMN Jerash hereditary motor neuronopathy. We studied the families at the molecular genetic level and mapped the Jerash hereditary motor neuronopathy gene to an approximately 0.54-cM region on chromosome 9p21.1-p12, flanked by microsatellite polymorphic marker loci D9S1845 and D9S1791. A maximum LOD score of 19.80 at theta = 0.001 was obtained between the disease and locus D9S1878.

The neuronal ceroid lipofuscinoses (NCLs), also known as Batten disease, are a group of autosomal recessive lysosomal storage disorders that are characterized by neurodegeneration, progressive cognitive decline, motor impairment, ataxia, loss of vision, seizures, and premature death. To date, pathogenic variants in more than 13 genes have been associated with NCLs. CLN6 encodes an endoplasmic reticulum non-glycosylated transmembrane protein, which is involved in lysosomal acidification. Mutations in CLN6 cause late-infantile juvenile NCL (JNCL) adult-onset NCL, and Kufs disease. Members from two available families with JNCL were clinically evaluated, and samples were collected from consenting individuals. The molecular investigation was performed by whole-exome sequencing, Sanger sequencing, and family segregation analysis. Furthermore, in silico prediction analysis and structural modeling of the identified CLN6 variants were performed. We report clinical and genetic findings of three patients from two Greek-Cypriot families (families 915 and 926) with JNCL. All patients were males, and the first symptoms appeared at the age of 6 years. The proband of family 926 presented with loss of motor abilities, ataxia, spasticity, seizure, and epilepsy. The proband of family 915 had ataxia, spasticity, dysarthria, dystonia, and intellectual disability. Both probands did not show initial signs of vision and/or hearing loss. Molecular analysis of family 926 revealed two CLN6 biallelic variants: the novel, de novo p.Tyr295Cys and the known p.Arg136His variants. In family 915, both patients were homozygous for the p.Arg136His CLN6 variant. Prediction analysis of the two CLN6 variants characterized them as probably damaging and disease-causing. Structural modeling of the variants predicted that they probably cause protein structural differentiation. In conclusion, we describe two unrelated Cypriot families with JNCL. Both families had variants in the CLN6 gene; however, they presented with slightly different symptoms, and notably none of the patients has loss of vision. In silico prediction and structural analyses indicate that both variants are most likely pathogenic.

The advancement of scientific and medical research over the past years has generated a wealth of experimental data from multiple technologies, including genomics, transcriptomics, proteomics and other forms of -omics data, which are available for a number of diseases. The integration of such multisource data is a key component towards the success of precision medicine. In this work we are investigating a multisource data integration method developed by our group, regarding its ability to drive to clusters of connected pathways under two different approaches: (a) a disease-centric approach, where we want to integrate data around a disease, and (b) a gene-centric approach, where we want to integrate data around a gene. We have used as a paradigm for the first approach Huntington's disease (HD), a disease with a plethora of available data, whereas for the second approach the GBA2, a gene that is related to Spastic Ataxia (SA), a phenotype with sparse availability of data. Our work shows that valuable information at the level of disease-related pathway clusters can be obtained for both HD and SA. New pathways that classical pathway analysis methods were unable to reveal, are emerging as necessary "connectors" to build connected pathway stories formed as pathway clusters.

Background This study aims to characterize SARS-CoV-2 mutations which are primarily prevalent in the Cypriot population. Moreover, using computational approaches, we assess whether these mutations are associated with changes in viral virulence. Methods We utilize genetic data from 144 sequences of SARS-CoV-2 strains from the Cypriot population obtained between March 2020 and January 2021, as well as all data available from GISAID. We combine this with countries’ regional information, such as deaths and cases per million, as well as COVID-19-related public health austerity measure response times. Initial indications of selective advantage of Cyprus-specific mutations are obtained by mutation tracking analysis. This entails calculating specific mutation frequencies within the Cypriot population and comparing these with their prevalence world-wide throughout the course of the pandemic. We further make use of linear regression models to extrapolate additional information that may be missed through standard statistical analysis. Results We report a single mutation found in the ORF1ab gene (nucleotide position 18,440) that appears to be significantly enriched within the Cypriot population. The amino acid change is denoted as S6059F, which maps to the SARS-CoV-2 NSP14 protein. We further analyse this mutation using regression models to investigate possible associations with increased deaths and cases per million. Moreover, protein structure prediction tools show that the mutation infers a conformational change to the protein that significantly alters its structure when compared to the reference protein. Conclusions Investigating Cyprus-specific mutations for SARS-CoV-2 can lead to a better understanding of viral pathogenicity. Researching these mutations can generate potential links between viral-specific mutations and the unique genomics of the Cypriot population. This can not only lead to important findings from which to battle the pandemic on a national level, but also provide insights into viral virulence worldwide.

The SPG7 gene encodes the paraplegin protein, an inner mitochondrial membrane—localized protease. It was initially linked to pure and complicated hereditary spastic paraplegia with cerebellar atrophy, and now represents a frequent cause of undiagnosed cerebellar ataxia and spastic ataxia. We hereby report the molecular characterization and the clinical features of a large Cypriot family with five affected individuals presenting with spastic ataxia in an autosomal recessive transmission mode, due to a novel SPG7 homozygous missense variant. Detailed clinical histories of the patients were obtained, followed by neurological and neurophysiological examinations. Whole exome sequencing (WES) of the proband, in silico gene panel analysis, variant filtering and family segregation analysis of the candidate variants with Sanger sequencing were performed. RNA and protein expression as well as in vitro protein localization studies and mitochondria morphology evaluation were carried out towards functional characterization of the identified variant. The patients presented with typical spastic ataxia features while some intrafamilial phenotypic variation was noted. WES analysis revealed a novel homozygous missense variant in the SPG7 gene (c.1763C > T, p. Thr588Met), characterized as pathogenic by more than 20 in silico prediction tools. Functional studies showed that the variant does not affect neither the RNA or protein expression, nor the protein localization. However, aberrant mitochondrial morphology has been observed thus indicating mitochondrial dysfunction and further demonstrating the pathogenicity of the identified variant. Our study is the first report of an SPG7 pathogenic variant in the Cypriot population and broadens the spectrum of SPG7 pathogenic variants.

Background Next-generation sequencing (NGS) represents a significant advancement in clinical genetics. However, its use creates several technical, data interpretation and management challenges. It is essential to follow a consistent data analysis pipeline to achieve the highest possible accuracy and avoid false variant calls. Herein, we aimed to compare the performance of twenty-eight combinations of NGS data analysis pipeline compartments, including short-read mapping (BWA-MEM, Bowtie2, Stampy), variant calling (GATK-HaplotypeCaller, GATK-UnifiedGenotyper, SAMtools) and interval padding (null, 50 bp, 100 bp) methods, along with a commercially available pipeline (BWA Enrichment, Illumina®). Fourteen germline DNA samples from breast cancer patients were sequenced using a targeted NGS panel approach and subjected to data analysis. Results We highlight that interval padding is required for the accurate detection of intronic variants including spliceogenic pathogenic variants (PVs). In addition, using nearly default parameters, the BWA Enrichment algorithm, failed to detect these spliceogenic PVs and a missense PV in the TP53 gene. We also recommend the BWA-MEM algorithm for sequence alignment, whereas variant calling should be performed using a combination of variant calling algorithms; GATK-HaplotypeCaller and SAMtools for the accurate detection of insertions/deletions and GATK-UnifiedGenotyper for the efficient detection of single nucleotide variant calls. Conclusions These findings have important implications towards the identification of clinically actionable variants through panel testing in a clinical laboratory setting, when dedicated bioinformatics personnel might not always be available. The results also reveal the necessity of improving the existing tools and/or at the same time developing new pipelines to generate more reliable and more consistent data.

Polyneuropathies are amongst the most common neurological conditions worldwide affecting over 20 million people. However, 40% of patients with primary polyneuropathies have no disease-causing mutation identified.We investigated patients with gene-negative primary polyneuropathies using a combination of whole genome sequencing, homozygosity mapping and segregation analysis. Pathogenicity was confirmed via enzymatic assays and mass spectroscopy on recombinant protein and patient-derived fibroblasts, plasma and erythrocytes. We used circular dichroism to show secondary structure changes and isothermal titration calorimetry to investigate the ATP binding.We report that biallelic mutations in human PDXK are associated with primary axonal polyneuropathy and optic atrophy. Pyridoxal kinase (PDXK) is involved in converting vitamin B6 to its active form, pyridoxal 5’-phosphate (PLP). We show that PDXK mutations lead to disease via decreased plasma PLP concentrations. Our functional studies revealed conformational rearrangement in the mutant enzyme around the kinase ATP-binding pocket with impaired PDXK ability to bind ATP and leading to reduced erythrocyte PDXK activity. We show that both the human clinical picture and biochemical profile in PDXK mutations are rescued by PLP supplementation. Patients regained their ability to walk independently. Furthermore, treatment-led normalisation of plasma PLP levels, correlated with reduction of neurofilament light chain concentrations, a biomarker of axonal breakdown.In conclusion, biallelic mutations in human PDXK are associated with a novel disorder leading to treatable primary axonal polyneuropathy and optic atrophy and identifies PLP as therapeutic target.