Mutations in the coding region of the WFS1 gene cause Wolfram syndrome, a rare multisystem neurodegenerative disorder of autosomal recessive inheritance. Patients with Wolfram syndrome display considerable clinical pleiomorphism, and symptoms such as neurological complications and psychiatric disorders are common. In the present study we have characterized Wfs1 expression pattern in the mouse central nervous system by using a combination of immunohistochemistry on wild-type mice and X-Gal staining of Wfs1 knockout mice with targeted insertion of the lacZ reporter. We identified a robust enrichment of Wfs1 protein in the central extended amygdala and ventral striatum. Prominent Wfs1 expression was seen in the hippocampal CA1 region, parasubiculum, superficial part of the second and third layers of the prefrontal cortex and proisocortical areas, hypothalamic magnocellular neurosecretory system, and central auditory pathway. Wfs1 expression was also detected in numerous brainstem nuclei and in laminae VIII and IX of the spinal cord. Wfs1-positive nerve fibers were found in the medial forebrain bundle, reticular part of the substantia nigra, globus pallidus, posterior caudate putamen, lateral lemniscus, alveus, fimbria, dorsal hippocampal commissure, subiculum, and to a lesser extent in the central sublenticular extended amygdala, compact part of substantia nigra, and ventral tegmental area. The neuroanatomical findings suggest that the lack of Wfs1 protein function can be related to several neurological and psychiatric symptoms found in Wolfram syndrome. Enrichment of Wfs1 protein in the central extended amygdala suggests a role in the modulation of anxiety and fear.

Mitochondrial dysfunction has been implicated in the etiology of monogenic Parkinson's disease (PD). Yet the role that mitochondrial processes play in the most common form of the disease; sporadic PD, is yet to be fully established. Here, we comprehensively assessed the role of mitochondrial function-associated genes in sporadic PD by leveraging improvements in the scale and analysis of PD GWAS data with recent advances in our understanding of the genetics of mitochondrial disease. We calculated a mitochondrial-specific polygenic risk score (PRS) and showed that cumulative small effect variants within both our primary and secondary gene lists are significantly associated with increased PD risk. We further reported that the PRS of the secondary mitochondrial gene list was significantly associated with later age at onset. Finally, to identify possible functional genomic associations we implemented Mendelian randomization, which showed that 14 of these mitochondrial function-associated genes showed functional consequence associated with PD risk. Further analysis suggested that the 14 identified genes are not only involved in mitophagy, but implicate new mitochondrial processes. Our data suggests that therapeutics targeting mitochondrial bioenergetics and proteostasis pathways distinct from mitophagy could be beneficial to treating the early stage of PD.

Failure of conventional anti-inflammatory therapies in osteoarthritis (OA) underlines the insufficient knowledge about inflammatory mechanisms, patterns and their relationship with cartilage degradation. Considering non-linear nature of cartilage loss in OA, a better understanding of inflammatory milieu and MMP status at different stages of OA is required to design early-stage therapies or personalized disease management. For this, an investigation based on a synovium-synovial fluid (SF) axis was planned to study OA associated changes in synovium and SF along the progressive grades of OA. Gene expressions in synovial-biopsies from different grades OA patients (N = 26) revealed a peak of IL-1β, IL-15, PGE2 and NGF in early OA (Kellgren-Lawrence (KL) grade-I and II); the highest MMP levels were found in advanced stages (KL grade-III and IV). MMPs (MMP-1, 13, 2 and 9) abundance and FALGPA activity estimated in forty SFs of progressive grades showed the maximum protein levels and activity in KL grade-II and III. In an SF challenge test, SW982 and THP1 cells were treated with progressive grade SFs to study the dynamics of MMPs modulation in inflammatory microenvironment; the test yielded a result pattern, which matched with FALGPA and the protein-levels estimation. Inflammatory mediators in SFs served as steering factor for MMP up-regulation. A correlation-matrix of IL-1β and MMPs revealed expressional negative correlation.

Gene Expression Omnibus (GEO) is a database repository hosting a substantial proportion of publicly available high throughput gene expression data. Gene expression analysis is a powerful tool to gain insight into the mechanisms and processes underlying the biological and phenotypic differences between sample groups. Despite the wide availability of gene expression datasets, their access, analysis, and integration are not trivial and require specific expertise and programming proficiency. We developed the GEOexplorer webserver to allow scientists to access, integrate and analyse gene expression datasets without requiring programming proficiency. Via its user-friendly graphic interface, users can easily apply GEOexplorer to perform interactive and reproducible gene expression analysis of microarray and RNA-seq datasets, while producing a wealth of interactive visualisations to facilitate data exploration and interpretation, and generating a range of publication ready figures. The webserver allows users to search and retrieve datasets from GEO as well as to upload user-generated data and combine and harmonise two datasets to perform joint analyses. GEOexplorer, available at https://geoexplorer.rosalind.kcl.ac.uk, provides a solution for performing interactive and reproducible analyses of microarray and RNA-seq gene expression data, empowering life scientists to perform exploratory data analysis and differential gene expression analysis on-the-fly without informatics proficiency.

Considerable heterogeneity exists in the literature concerning genetic determinants of the age at onset (AAO) of Parkinson disease (PD), which could be attributed to a lack of well-powered replication cohorts. The previous largest genome-wide association studies (GWAS) identified SNCA and TMEM175 loci on chromosome (Chr) 4 with a significant influence on the AAO of PD; these have not been independently replicated. This study aims to conduct a meta-analysis of GWAS of PD AAO and validate previously observed findings in worldwide populations.

Psoriasis is a chronic autoimmune disease with complex genetic architecture. Previous genome-wide association studies (GWAS) and a recent meta-analysis using Immunochip data have uncovered 36 susceptibility loci. Here, we extend our previous meta-analysis of European ancestry by refined genotype calling and imputation and by the addition of 5,033 cases and 5,707 controls. The combined analysis, consisting of over 15,000 cases and 27,000 controls, identifies five new psoriasis susceptibility loci at genome-wide significance (P<5 × 10(-8)). The newly identified signals include two that reside in intergenic regions (1q31.1 and 5p13.1) and three residing near PLCL2 (3p24.3), NFKBIZ (3q12.3) and CAMK2G (10q22.2). We further demonstrate that NFKBIZ is a TRAF3IP2-dependent target of IL-17 signalling in human skin keratinocytes, thereby functionally linking two strong candidate genes. These results further integrate the genetics and immunology of psoriasis, suggesting new avenues for functional analysis and improved therapies.

The genetics of an individual is a crucial factor in understanding the risk of developing the neurodegenerative disease amyotrophic lateral sclerosis (ALS). There is still a large proportion of the heritability of ALS, particularly in sporadic cases, to be understood. Among others, active transposable elements drive inter-individual variability, and in humans long interspersed element 1 (LINE1, L1), Alu and SINE-VNTR-Alu (SVA) retrotransposons are a source of polymorphic insertions in the population. We undertook a pilot study to characterise the landscape of non-reference retrotransposon insertion polymorphisms (non-ref RIPs) in 15 control and 15 ALS individuals' whole genomes from Project MinE, an international project to identify potential genetic causes of ALS. The combination of two bioinformatics tools (mobile element locator tool (MELT) and TEBreak) identified on average 1250 Alu, 232 L1 and 77 SVA non-ref RIPs per genome across the 30 analysed. Further PCR validation of individual polymorphic retrotransposon insertions showed a similar level of accuracy for MELT and TEBreak. Our preliminary study did not identify a specific RIP or a significant difference in the total number of non-ref RIPs in ALS compared to control genomes. The use of multiple bioinformatic tools improved the accuracy of non-ref RIP detection and our study highlights the potential importance of studying these elements further in ALS.

We describe the characterisation of a variable number tandem repeat (VNTR) domain within intron 1 of the amyotrophic lateral sclerosis (ALS) risk gene CFAP410 (Cilia and flagella associated protein 410) (previously known as C21orf2), providing insight into how this domain could support differential gene expression and thus be a modulator of ALS progression or risk. We demonstrated the VNTR was functional in a reporter gene assay in the HEK293 cell line, exhibiting both the properties of an activator domain and a transcriptional start site, and that the differential expression was directed by distinct repeat number in the VNTR. These properties embedded in the VNTR demonstrated the potential for this VNTR to modulate CFAP410 expression. We extrapolated these findings in silico by utilisation of tagging SNPs for the two most common VNTR alleles to establish a correlation with endogenous gene expression. Consistent with in vitro data, CFAP410 isoform expression was found to be variable in the brain. Furthermore, although the number of matched controls was low, there was evidence for one specific isoform being correlated with lower expression in those with ALS. To address if the genotype of the VNTR was associated with ALS risk, we characterised the variation of the CFAP410 VNTR in ALS cases and matched controls by PCR analysis of the VNTR length, defining eight alleles of the VNTR. No significant difference was observed between cases and controls, we noted, however, the cohort was unlikely to contain sufficient power to enable any firm conclusion to be drawn from this analysis. This data demonstrated that the VNTR domain has the potential to modulate CFAP410 expression as a regulatory element that could play a role in its tissue-specific and stimulus-inducible regulation that could impact the mechanism by which CFAP410 is involved in ALS.

Transposable elements (TEs) are repetitive elements that belong to a variety of functional classes and have an important role in shaping genome evolution. Around 50% of the human genome contains TEs, and they have been termed the "dark matter" of the genome because relatively little is known about their function. While TEs have been shown to participate in aberrant gene regulation and the pathogenesis of diseases, only a few studies have explored the systemic effect of TEs on gene expression. In the present study, we analysed whole genome sequences and blood whole transcriptome data from 570 individuals within the Parkinson's Progressive Markers Initiative (PPMI) cohort to identify expression quantitative trait loci (eQTL) regulating genome-wide gene expression associated with TEs. We identified 2132 reference TEs that were polymorphic for their presence or absence in our study cohort. The presence or absence of the TE element could change the expression of the gene or gene clusters from zero to tens of thousands of copies of RNA. The main finding is that many TEs possess very strong regulatory effects, and they have the potential to modulate large genetic networks with hundreds of target genes over the genome. We illustrate the plethora of regulatory mechanisms using examples of their action at the HLA gene cluster and data showing different TEs' convergence to modulate WFS1 gene expression. In conclusion, the presence or absence of polymorphisms of TEs has an eminent genome-wide regulatory function with large effect size at the level of the whole transcriptome. The role of TEs in explaining, in part, the missing heritability for complex traits is convincing and should be considered.

Macrophage polarization is a steering factor of osteoarthritis (OA) progression. Synovial fluid (SF) obtained from OA patients with different Kellgren-Lawrence grades (KL grades) holds several proinflammatory factors and was hypothesized to induce macrophage differentiation and polarization by providing the needed microenvironment. U937 cells and peripheral-blood-mononuclear-cell-derived monocytes (PBMC-derived CD14+ cells) were induced with SFs of progressive KL grades for 48 h, and the status of the differentiated cells was evaluated by cell surface markers representing M1 and M2 macrophage phenotypes. Functional viability assessment of the differentiated cells was performed by cytokine estimation. The fraction of macrophages and their phenotypes were estimated by immunophenotyping of SF-isolated cells of different KL grades. A grade-wise proteome analysis of SFs was performed in search of the factors which are influential in macrophage differentiation and polarization. In the assay on U937 cells, induction with SF of KL grade III and IV showed a significant increase in M1 type (CD86+). The percentage of M2 phenotype (CD163+) was significantly higher after the induction with SF of KL grade II. A Significantly higher M1/M2 ratio was estimated in the cells induced with KL grade III and IV. The cell differentiation pattern in the assay on PBMC-derived CD14+ cells showed a grade-wise decline in both M1 (CD11C+, CD86+) and M2 phenotype (CD163+). Cytokine estimation specific to M1 (TNF-α, IL-6, IL-1β, IFN-γ) and M2 (IL-4 and IL-10) macrophages corelated with the differentiation pattern in the U937 cell assay, while it did not reveal any significant changes in the PBMC-derived CD14+ cells assay. SF cells' immunophenotyping showed the highest percentage of CD14+ macrophages in KL grade II; CD86+ and CD163+ cells were minimal in all KL grades' SFs. The proteome analysis revealed significantly expressed MIF, CAPG/MCP, osteopontin, and RAS-related RAB proteins in KL grade III and IV samples, which are linked with macrophages' movement, polarization, and migration-behavior. In conclusion, this study demonstrated that SF in OA joints acts as a niche and facilitates M1 phenotype polarization by providing a proinflammatory microenvironment.

Long interspersed element-1 (LINE-1/L1s) contributes 17% of the human genome with more than 1 million elements present; however, fewer than 100 of these have evidence for being retrotransposition competent (RC). In addition to those RC-L1s present in the reference genome, there are a small number of known non-reference L1 insertions that are also retrotransposition competent. L1 activity, whether through the potentially detrimental effects of their mRNA or protein expression or somatic retrotransposition events, has been linked to several neurological conditions. The polymorphic nature of both reference and non-reference RC-L1s in terms of their presence or absence will result in individuals harboring a different combination of these elements and it is currently unknown if this type of germline variation contributes to the risk of neurological disease. Here, we utilized whole-genome sequencing data from 178 healthy controls and 372 Parkinson's disease (PD) subjects from the Parkinson's Progression Markers Initiative (PPMI) to investigate the role of RC-L1s in PD. In the PPMI cohort, we identified 22 reference and 50 non-reference polymorphic RC-L1 loci. Focusing on 16 highly active RC-L1 loci, an increased burden of these elements (≥9) was associated with PD (OR 1.25, 95% CI 1.03-1.51, p = 0.02). In addition, we identified significant associations of progression markers of PD and the burden of highly active RC-L1s. This study has identified a novel type of genetic element associated with PD risk and disease progression.

The development of Parkinson's disease (PD) involves a complex interaction of genetic and environmental factors. Genome-wide association studies using extensive single nucleotide polymorphism datasets have identified many loci involved in disease. However much of the heritability of Parkinson's disease is still to be identified and the functional elements associated with the risk to be determined and understood. To investigate the component of PD that may involve complex genetic variants we characterised the hominid specific retrotransposon SINE-VNTR-Alus (SVAs) in the Parkinson's Progression Markers Initiative cohort utilising whole genome sequencing. We identified 81 reference SVAs polymorphic for their presence/absence, seven of which were associated with the progression of the disease and with differential gene expression in whole blood RNA sequencing data. This study highlights the importance of addressing SVA variants and potentially other types of retrotransposons in PD genetics, furthermore, these SVA elements should be considered as regulatory domains that could play a role in disease progression.

Levodopa-induced dyskinesia (LID) is a common adverse effect of levodopa, one of the main therapeutics used to treat the motor symptoms of Parkinson's disease (PD). Previous evidence suggests a connection between LID and a disruption of the dopaminergic system as well as genes implicated in PD, including GBA1 and LRRK2.

Psoriatic arthritis (PsA) is a chronic inflammatory musculoskeletal disease affecting up to 30% of psoriasis vulgaris (PsV) cases and approximately 0.25 to 1% of the general population. To identify common susceptibility loci, we performed a meta-analysis of three imputed genome-wide association studies (GWAS) on psoriasis, stratified for PsA. A total of 1,160,703 single-nucleotide polymorphisms (SNPs) were analyzed in the discovery set consisting of 535 PsA cases and 3,432 controls from Germany, the United States, and Canada. We followed up two SNPs in 1,931 PsA cases and 6,785 controls comprising six independent replication panels from Germany, Estonia, the United States, and Canada. In the combined analysis, a genome-wide significant association was detected at 2p16 near the REL locus encoding c-Rel (rs13017599, P=1.18 × 10(-8), odds ratio (OR)=1.27, 95% confidence interval (CI)=1.18-1.35). The rs13017599 polymorphism is known to associate with rheumatoid arthritis (RA), and another SNP near REL (rs702873) was recently implicated in PsV susceptibility. However, conditional analysis indicated that rs13017599, rather than rs702873, accounts for the PsA association at REL. We hypothesize that c-Rel, as a member of the Rel/NF-κB family, is associated with PsA in the context of disease pathways that involve other identified PsA and PsV susceptibility genes including TNIP1, TNFAIP3, and NFκBIA.

To gain further insight into the genetic architecture of psoriasis, we conducted a meta-analysis of 3 genome-wide association studies (GWAS) and 2 independent data sets genotyped on the Immunochip, including 10,588 cases and 22,806 controls. We identified 15 new susceptibility loci, increasing to 36 the number associated with psoriasis in European individuals. We also identified, using conditional analyses, five independent signals within previously known loci. The newly identified loci shared with other autoimmune diseases include candidate genes with roles in regulating T-cell function (such as RUNX3, TAGAP and STAT3). Notably, they included candidate genes whose products are involved in innate host defense, including interferon-mediated antiviral responses (DDX58), macrophage activation (ZC3H12C) and nuclear factor (NF)-κB signaling (CARD14 and CARM1). These results portend a better understanding of shared and distinctive genetic determinants of immune-mediated inflammatory disorders and emphasize the importance of the skin in innate and acquired host defense.

Long interspersed element-1 (LINE-1/L1) is the only autonomous transposable element in the human genome that currently mobilises in both germline and somatic tissues. Recent studies have identified correlations between altered retrotransposon expression and the fatal neurodegenerative disease amyotrophic lateral sclerosis (ALS) in a subset of patients. The risk of an individual developing ALS is dependent on an interaction of genetic variants and subsequent modifiers during life. These modifiers could include environmental factors, which can lead to epigenetic and genomic changes, such as somatic mutations, occurring in the neuronal cells that degenerate as the disease develops. There are more than 1 million L1 copies in the human genome today, but only 80-100 L1 loci in the reference genome are considered to be retrotransposition-competent (RC) and an even smaller number of these RC-L1s loci are highly active. We hypothesise that RC-L1s could affect normal cellular function through their mutagenic potential conferred by their ability to retrotranspose in neuronal cells and through DNA damage caused by the endonuclease activity of the L1-encoded ORF2 protein. To investigate whether either an increase in the genomic burden of RC-L1s or epigenetic changes to RC-L1s altering their expression, could play a role in disease development, we chose a set of seven well characterised genomic RC-L1 loci that were reported earlier to be highly active in a cellular L1 retrotransposition reporter assay or serve as major source elements for germline and/or somatic retrotransposition events. Analysis of the insertion allele frequency of five polymorphic RC-L1s, out of the set of seven, for their presence or absence, did not identify an increased number individually or when combined in individuals with the disease. However, we did identify reduced levels of methylation of RC-L1s in the motor cortex of those individuals with both familial and sporadic ALS compared to control brains. The changes to the regulation of the loci encompassing these RC-L1s demonstrated tissue specificity and could be related to the disease process.

The hominid-specific retrotransposon SINE-VNTR-Alu (SVA) is a composite element that has contributed to the genetic variation between individuals and influenced genomic structure and function. SVAs are involved in modulating gene expression and splicing patterns, altering mRNA levels and sequences, and have been associated with the development of disease. We evaluated the genome-wide effects of SVAs present in the reference genome on transcript sequence and expression in the CNS of individuals with and without the neurodegenerative disorder Amyotrophic Lateral Sclerosis (ALS). This study identified SVAs in the exons of 179 known transcripts, several of which were expressed in a tissue-specific manner, as well as 92 novel exonisation events occurring in the motor cortex. An analysis of 65 reference genome SVAs polymorphic for their presence/absence in the ALS consortium cohort did not identify any elements that were significantly associated with disease status, age at onset, and survival. However, there were transcripts, such as transferrin and HLA-A, that were differentially expressed between those with or without disease, and expression levels were associated with the genotype of proximal SVAs. This study demonstrates the functional consequences of several SVA elements altering mRNA splicing patterns and expression levels in tissues of the CNS.

Osteophytes are a prominent feature of osteoarthritis (OA) joints and one of the clinical hallmarks of the disease progression. Research on osteophytes is fragmentary and modes of its contribution to OA pathology are obscure.

GTP cyclohydrolase 1, encoded by the GCH1 gene, is an essential enzyme for dopamine production in nigrostriatal cells. Loss-of-function mutations in GCH1 result in severe reduction of dopamine synthesis in nigrostriatal cells and are the most common cause of DOPA-responsive dystonia, a rare disease that classically presents in childhood with generalized dystonia and a dramatic long-lasting response to levodopa. We describe clinical, genetic and nigrostriatal dopaminergic imaging ([(123)I]N-ω-fluoropropyl-2β-carbomethoxy-3β-(4-iodophenyl) tropane single photon computed tomography) findings of four unrelated pedigrees with DOPA-responsive dystonia in which pathogenic GCH1 variants were identified in family members with adult-onset parkinsonism. Dopamine transporter imaging was abnormal in all parkinsonian patients, indicating Parkinson's disease-like nigrostriatal dopaminergic denervation. We subsequently explored the possibility that pathogenic GCH1 variants could contribute to the risk of developing Parkinson's disease, even in the absence of a family history for DOPA-responsive dystonia. The frequency of GCH1 variants was evaluated in whole-exome sequencing data of 1318 cases with Parkinson's disease and 5935 control subjects. Combining cases and controls, we identified a total of 11 different heterozygous GCH1 variants, all at low frequency. This list includes four pathogenic variants previously associated with DOPA-responsive dystonia (Q110X, V204I, K224R and M230I) and seven of undetermined clinical relevance (Q110E, T112A, A120S, D134G, I154V, R198Q and G217V). The frequency of GCH1 variants was significantly higher (Fisher's exact test P-value 0.0001) in cases (10/1318 = 0.75%) than in controls (6/5935 = 0.1%; odds ratio 7.5; 95% confidence interval 2.4-25.3). Our results show that rare GCH1 variants are associated with an increased risk for Parkinson's disease. These findings expand the clinical and biological relevance of GTP cycloydrolase 1 deficiency, suggesting that it not only leads to biochemical striatal dopamine depletion and DOPA-responsive dystonia, but also predisposes to nigrostriatal cell loss. Further insight into GCH1-associated pathogenetic mechanisms will shed light on the role of dopamine metabolism in nigral degeneration and Parkinson's disease.

To test the hypothesis that rs573116164 will have disease-modifying effects in patients with superoxide dismutase 1 (SOD1) familial amyotrophic lateral sclerosis (fALS), we characterized rs573116164 within a cohort of 190 patients with fALS and 560 healthy age-matched controls to assess the variant for association with various measures of disease.