Deciphering the role of a SINE-VNTR-Alu retrotransposon polymorphism as a biomarker of Parkinson's disease progression.

SINE-VNTR-Alu (SVA) retrotransposons are transposable elements which represent a source of genetic variation. We previously demonstrated that the presence/absence of a human-specific SVA, termed SVA_67, correlated with the progression of Parkinson's disease (PD). In the present study, we demonstrate that SVA_67 acts as expression quantitative trait loci, thereby exhibiting a strong regulatory effect across the genome using whole genome and transcriptomic data from the Parkinson's progression markers initiative cohort. We further show that SVA_67 is polymorphic for its variable number tandem repeat domain which correlates with both regulatory properties in a luciferase reporter gene assay in vitro and differential expression of multiple genes in vivo. Additionally, this variation's utility as a biomarker is reflected in a correlation with a number of PD progression markers. These experiments highlight the plethora of transcriptomic and phenotypic changes associated with SVA_67 polymorphism which should be considered when investigating the missing heritability of neurodegenerative diseases.

Exploring SVA Insertion Polymorphisms in Shaping Differential Gene Expressions in the Central Nervous System.

Transposable elements (TEs) are repetitive elements which make up around 45% of the human genome. A class of TEs, known as SINE-VNTR-Alu (SVA), demonstrate the capacity to mobilise throughout the genome, resulting in SVA polymorphisms for their presence or absence within the population. Although studies have previously highlighted the involvement of TEs within neurodegenerative diseases, such as Parkinson's disease and amyotrophic lateral sclerosis (ALS), the exact mechanism has yet to be identified. In this study, we used whole-genome sequencing and RNA sequencing data of ALS patients and healthy controls from the New York Genome Centre ALS Consortium to elucidate the influence of reference SVA elements on gene expressions genome-wide within central nervous system (CNS) tissues. To investigate this, we applied a matrix expression quantitative trait loci analysis and demonstrate that reference SVA insertion polymorphisms can significantly modulate the expression of numerous genes, preferentially in the trans position and in a tissue-specific manner. We also highlight that SVAs significantly regulate mitochondrial genes as well as genes within the HLA and MAPT loci, previously associated within neurodegenerative diseases. In conclusion, this study continues to bring to light the effects of polymorphic SVAs on gene regulation and further highlights the importance of TEs within disease pathology.

Rare and common genetic determinants of mitochondrial function determine severity but not risk of amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease involving selective vulnerability of energy-intensive motor neurons (MNs). It has been unclear whether mitochondrial function is an upstream driver or a downstream modifier of neurotoxicity. We separated upstream genetic determinants of mitochondrial function, including genetic variation within the mitochondrial genome or autosomes; from downstream changeable factors including mitochondrial DNA copy number (mtCN). Across three cohorts including 6,437 ALS patients, we discovered that a set of mitochondrial haplotypes, chosen because they are linked to measurements of mitochondrial function, are a determinant of ALS survival following disease onset, but do not modify ALS risk. One particular haplotype appeared to be neuroprotective and was significantly over-represented in two cohorts of long-surviving ALS patients. Causal inference for mitochondrial function was achievable using mitochondrial haplotypes, but not autosomal SNPs in traditional Mendelian randomization (MR). Furthermore, rare loss-of-function genetic variants within, and reduced MN expression of, ACADM and DNA2 lead to ∼50 % shorter ALS survival; both proteins are implicated in mitochondrial function. Both mtCN and cellular vulnerability are linked to DNA2 function in ALS patient-derived neurons. Finally, MtCN responds dynamically to the onset of ALS independently of mitochondrial haplotype, and is correlated with disease severity. We conclude that, based on the genetic measures we have employed, mitochondrial function is a therapeutic target for amelioration of disease severity but not prevention of ALS.

A SINE-VNTR-Alu at the LRIG2 locus is associated with proximal and distal gene expression in CRISPR and population models.

SINE-VNTR-Alu (SVA) retrotransposons represent mobile regulatory elements that have the potential to influence the surrounding genome when they insert into a locus. Evolutionarily recent mobilisation has resulted in loci in the human genome where a given retrotransposon might be observed to be present or absent, termed a retrotransposon insertion polymorphism (RIP). We previously observed that an SVA RIP ~ 2 kb upstream of LRIG2 on chromosome 1, the 'LRIG2 SVA', was associated with differences in local gene expression and methylation, and that the two were correlated. Here, we have used CRISPR-mediated deletion of the LRIG2 SVA in a cell line model to validate that presence of the retrotransposon is directly affecting local expression and provide evidence that is suggestive of a modest role for the SVA in modulating nearby methylation. Additionally, in leveraging an available Hi-C dataset we observed that the LRIG2 SVA was also involved in long-range chromatin interactions with a cluster of genes ~ 300 kb away, and that expression of these genes was to varying degrees associated with dosage of the SVA in both CRISPR cell line and population models. Altogether, these data support a regulatory role for SVAs in the modulation of gene expression, with the latter potentially involving chromatin looping, consistent with the model that RIPs may contribute to interpersonal differences in transcriptional networks.

Unsupervised machine learning identifies distinct ALS molecular subtypes in post-mortem motor cortex and blood expression data.

Amyotrophic lateral sclerosis (ALS) displays considerable clinical and genetic heterogeneity. Machine learning approaches have previously been utilised for patient stratification in ALS as they can disentangle complex disease landscapes. However, lack of independent validation in different populations and tissue samples have greatly limited their use in clinical and research settings. We overcame these issues by performing hierarchical clustering on the 5000 most variably expressed autosomal genes from motor cortex expression data of people with sporadic ALS from the KCL BrainBank (N = 112). Three molecular phenotypes linked to ALS pathogenesis were identified: synaptic and neuropeptide signalling, oxidative stress and apoptosis, and neuroinflammation. Cluster validation was achieved by applying linear discriminant analysis models to cases from TargetALS US motor cortex (N = 93), as well as Italian (N = 15) and Dutch (N = 397) blood expression datasets, for which there was a high assignment probability (80-90%) for each molecular subtype. The ALS and motor cortex specificity of the expression signatures were tested by mapping KCL BrainBank controls (N = 59), and occipital cortex (N = 45) and cerebellum (N = 123) samples from TargetALS to each cluster, before constructing case-control and motor cortex-region logistic regression classifiers. We found that the signatures were not only able to distinguish people with ALS from controls (AUC 0.88 ± 0.10), but also reflect the motor cortex-based disease process, as there was perfect discrimination between motor cortex and the other brain regions. Cell types known to be involved in the biological processes of each molecular phenotype were found in higher proportions, reinforcing their biological interpretation. Phenotype analysis revealed distinct cluster-related outcomes in both motor cortex datasets, relating to disease onset and progression-related measures. Our results support the hypothesis that different mechanisms underpin ALS pathogenesis in subgroups of patients and demonstrate potential for the development of personalised treatment approaches. Our method is available for the scientific and clinical community at https://alsgeclustering.er.kcl.ac.uk .

Regulation of expression quantitative trait loci by SVA retrotransposons within the major histocompatibility complex.

Genomic and transcriptomic studies of expression quantitative trait loci (eQTL) revealed that SINE-VNTR-Alu (SVA) retrotransposon insertion polymorphisms (RIPs) within human genomes markedly affect the co-expression of many coding and noncoding genes by coordinated regulatory processes. This study examined the polymorphic SVA modulation of gene co-expression within the major histocompatibility complex (MHC) genomic region where more than 160 coding genes are involved in innate and adaptive immunity. We characterized the modulation of SVA RIPs utilizing the genomic and transcriptomic sequencing data obtained from whole blood of 1266 individuals in the Parkinson's Progression Markers Initiative (PPMI) cohort that included an analysis of human leukocyte antigen (HLA)-A regulation in a subpopulation of the cohort. The regulatory properties of eight SVAs located within the class I and class II MHC regions were associated with differential co-expression of 71 different genes within and 75 genes outside the MHC region. Some of the same genes were affected by two or more different SVA. Five SVA are annotated in the human genomic reference sequence GRCh38.p14/hg38, whereas the other three were novel insertions within individuals. We also examined and found distinct structural effects (long and short variants and the CT internal variants) for one of the SVA (R_SVA_24) insertions on the differential expression of the HLA-A gene within a subpopulation (550 individuals) of the PPMI cohort. This is the first time that many HLA and non-HLA genes (multilocus expression units) and splicing mechanisms have been shown to be regulated by eight structurally polymorphic SVA within the MHC genomic region by applying precise statistical analysis of RNA data derived from the blood samples of a human cohort population. This study shows that SVA within the MHC region are important regulators or rheostats of gene co-expression that might have potential roles in diversity, health, and disease.

Transcriptomic profiling of cerebrospinal fluid identifies ALS pathway enrichment and RNA biomarkers in MND individuals.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder and the most common form of motor neurone disease (MND) which is characterized by the damage and death of motor neurons in the brain and spinal cord of affected individuals. Due to the heterogeneity of the disease, a better understanding of the interaction between genetics and biochemistry with the identification of biomarkers is crucial for therapy development. In this study, we used cerebrospinal fluid (CSF) RNA-sequencing data from the New York Genome Center (NYGC) ALS Consortium and analyzed differential gene expression between 47 MND individuals and 29 healthy controls. Pathway analysis showed that the affected genes are enriched in many pathways associated with ALS, including nucleocytoplasmic transport, autophagy, and apoptosis. Moreover, we assessed differential expression on both gene- and transcript-based levels and demonstrate that the expression of previously identified potential biomarkers, including CAPG, CCL3, and MAP2, was significantly higher in MND individuals. Ultimately, this study highlights the transcriptomic composition of CSF which enables insights into changes in the brain in ALS and therefore increases the confidence in the use of CSF for biomarker development.

NeuroBooster Array: A Genome-Wide Genotyping Platform to Study Neurological Disorders Across Diverse Populations.

Genome-wide genotyping platforms have the capacity to capture genetic variation across different populations, but there have been disparities in the representation of population-dependent genetic diversity. The motivation for pursuing this endeavor was to create a comprehensive genome-wide array capable of encompassing a wide range of neuro-specific content for the Global Parkinson's Genetics Program (GP2) and the Center for Alzheimer's and Related Dementias (CARD). CARD aims to increase diversity in genetic studies, using this array as a tool to foster inclusivity. GP2 is the first supported resource project of the Aligning Science Across Parkinson's (ASAP) initiative that aims to support a collaborative global effort aimed at significantly accelerating the discovery of genetic factors contributing to Parkinson's disease and atypical parkinsonism by generating genome-wide data for over 200,000 individuals in a multi-ancestry context. Here, we present the Illumina NeuroBooster array (NBA), a novel, high-throughput and cost-effective custom-designed content platform to screen for genetic variation in neurological disorders across diverse populations. The NBA contains a backbone of 1,914,934 variants (Infinium Global Diversity Array) complemented with custom content of 95,273 variants implicated in over 70 neurological conditions or traits with potential neurological complications. Furthermore, the platform includes over 10,000 tagging variants to facilitate imputation and analyses of neurodegenerative disease-related GWAS loci across diverse populations. The NBA can identify low frequency variants and accurately impute over 15 million common variants from the latest release of the TOPMed Imputation Server as of August 2023 (reference of over 300 million variants and 90,000 participants). We envisage this valuable tool will standardize genetic studies in neurological disorders across different ancestral groups, allowing researchers to perform genetic research inclusively and at a global scale.

CRISPR deletion of a SINE-VNTR-Alu (SVA_67) retrotransposon demonstrates its ability to differentially modulate gene expression at the MAPT locus.

Background: SINE-VNTR-Alu (SVA) retrotransposons are hominid-specific elements which have been shown to play important roles in processes such as chromatin structure remodelling and regulation of gene expression demonstrating that these repetitive elements exert regulatory functions. We have previously shown that the presence or absence of a specific SVA element, termed SVA_67, was associated with differential expression of several genes at the MAPT locus, a locus associated with Parkinson's Disease (PD) and frontotemporal dementia. However, we were not able to demonstrate that causation of differential gene expression was directed by the SVA due to lack of functional validation. Methods: We performed CRISPR to delete SVA_67 in the HEK293 cell line. Quantification of target gene expression was performed using qPCR to assess the effects on expression in response to the deletion of SVA_67. Differences between CRISPR edit and control cell lines were analysed using two-tailed t-test with a minimum 95% confidence interval to determine statistical significance. Results: In this study, we provide data highlighting the SVA-specific effect on differential gene expression. We demonstrate that the hemizygous deletion of the endogenous SVA_67 in CRISPR edited cell lines was associated with differential expression of several genes at the MAPT locus associated with neurodegenerative diseases including KANSL1, MAPT and LRRC37A. Discussion: This data is consistent with our previous bioinformatic work of differential gene expression analysis using transcriptomic data from the Parkinson's Progression Markers Initiative (PPMI) cohort. As SVAs have regulatory influences on gene expression, and insertion polymorphisms contribute to interpersonal differences in expression patterns, these results highlight the potential contribution of these elements to complex diseases with potentially many genetic components, such as PD.

Reference LINE-1 insertion polymorphisms correlate with Parkinson's disease progression and differential transcript expression in the PPMI cohort.

Long interspersed nuclear element-1 (LINE-1/L1) retrotransposons make up 17% of the human genome. They represent one class of transposable elements with the capacity to both mobilize autonomously and in trans via the mobilization of other elements, primarily Alu and SVA elements. Reference LINE-1 elements are, by definition, found in the reference genome, however, due to the polymorphic nature of these elements, variation for presence or absence is present within the population. We used a combination of clinical and transcriptomic data from the Parkinson's Progression Markers Initiative (PPMI) and applied matrix expression quantitative trait loci analysis and linear mixed-effects models involving 114 clinical, biochemical and imaging data from the PPMI cohort to elucidate the role of reference LINE-1 insertion polymorphism on both gene expression genome-wide and progression of Parkinson's disease (PD). We demonstrate that most LINE-1 insertion polymorphisms are capable of regulating gene expression, preferentially in trans, including previously identified PD risk loci. In addition, we show that 70 LINE-1 elements were associated with longitudinal changes of at least one PD progression marker, including ipsilateral count density ratio and UPDRS scores which are indicators of degeneration and severity. In conclusion, this study highlights the effect of the polymorphic nature of LINE-1 retrotransposons on gene regulation and progression of PD which underlines the importance of analyzing transposable elements within complex diseases.

A Genome-Wide Screen for the Exonisation of Reference SINE-VNTR-Alus and Their Expression in CNS Tissues of Individuals with Amyotrophic Lateral Sclerosis.

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.

Aspirin modulates production of pro-inflammatory and pro-resolving mediators in endothelial cells.

Endothelial cells synthesize biochemical signals to coordinate a response to insults, resolve inflammation and restore barrier integrity. Vascular cells release a variety of vasoactive bioactive lipid metabolites during the inflammatory response and produce pro-resolving mediators (e.g., Lipoxin A4, LXA4) in cooperation with leukocytes and platelets to bring a halt to inflammation. Aspirin, used in a variety of cardiovascular and pro-thrombotic disorders (e.g., atherosclerosis, angina, preeclampsia), potently inhibits proinflammatory eicosanoid formation. Moreover, aspirin stimulates the synthesis of pro-resolving lipid mediators (SPM), so-called Aspirin-Triggered Lipoxins (ATL). We demonstrate that cytokines stimulated a time- and dose-dependent increase in PGI2 (6-ketoPGF1α) and PGE2 formation that is blocked by aspirin. Eicosanoid production was caused by cytokine-induced expression of cyclooxygenase-2 (COX-2). We also detected increased production of pro-resolving LXA4 in cytokine-stimulated endothelial cells. The R-enantiomer of LXA4, 15-epi-LXA4, was enhanced by aspirin, but only in the presence of cytokine challenge, indicating dependence on COX-2 expression. In contrast to previous reports, we detected arachidonate 5-lipoxygenase (ALOX5) mRNA expression and its cognate protein (5-lipoxygenase, 5-LOX), suggesting that endothelial cells possess the enzymatic machinery necessary to synthesize both pro-inflammatory and pro-resolving lipid mediators independent of added leukocytes or platelets. Finally, we observed that, endothelial cells produced LTB4 in the absence of leukocytes. These results indicate that endothelial cells produce both pro-inflammatory and pro-resolving lipid mediators in the absence of other cell types and aspirin exerts pleiotropic actions influencing both COX and LOX pathways.

Genome-Wide Analysis of Structural Variants in Parkinson Disease.

OBJECTIVE: Identification of genetic risk factors for Parkinson disease (PD) has to date been primarily limited to the study of single nucleotide variants, which only represent a small fraction of the genetic variation in the human genome. Consequently, causal variants for most PD risk are not known. Here we focused on structural variants (SVs), which represent a major source of genetic variation in the human genome. We aimed to discover SVs associated with PD risk by performing the first large-scale characterization of SVs in PD. METHODS: We leveraged a recently developed computational pipeline to detect and genotype SVs from 7,772 Illumina short-read whole genome sequencing samples. Using this set of SV variants, we performed a genome-wide association study using 2,585 cases and 2,779 controls and identified SVs associated with PD risk. Furthermore, to validate the presence of these variants, we generated a subset of matched whole-genome long-read sequencing data. RESULTS: We genotyped and tested 3,154 common SVs, representing over 412 million nucleotides of previously uncatalogued genetic variation. Using long-read sequencing data, we validated the presence of three novel deletion SVs that are associated with risk of PD from our initial association analysis, including a 2 kb intronic deletion within the gene LRRN4. INTERPRETATION: We identified three SVs associated with genetic risk of PD. This study represents the most comprehensive assessment of the contribution of SVs to the genetic risk of PD to date. ANN NEUROL 2023;93:1012-1022.

RetroSnake: A modular pipeline to detect human endogenous retroviruses in genome sequencing data.

Human endogenous retroviruses (HERVs) integrated into the human genome as a result of ancient exogenous infections and currently comprise ∼8% of our genome. The members of the most recently acquired HERV family, HERV-Ks, still retain the potential to produce viral molecules and have been linked to a wide range of diseases including cancer and neurodegeneration. Although a range of tools for HERV detection in NGS data exist, most of them lack wet lab validation and they do not cover all steps of the analysis. Here, we describe RetroSnake, an end-to-end, modular, computationally efficient, and customizable pipeline for the discovery of HERVs in short-read NGS data. RetroSnake is based on an extensively wet-lab validated protocol, it covers all steps of the analysis from raw data to the generation of annotated results presented as an interactive html file, and it is easy to use by life scientists without substantial computational training. Availability and implementation: The Pipeline and an extensive documentation are available on GitHub.

Src Family Kinases Facilitate the Crosstalk between CGRP and Cytokines in Sensitizing Trigeminal Ganglion via Transmitting CGRP Receptor/PKA Pathway.

The communication between calcitonin gene-related peptide (CGRP) and cytokines plays a prominent role in maintaining trigeminal ganglion (TG) and trigeminovascular sensitization. However, the underlying regulatory mechanism is elusive. In this study, we explored the hypothesis that Src family kinases (SFKs) activity facilitates the crosstalk between CGRP and cytokines in sensitizing TG. Mouse TG tissue culture was performed to study CGRP release by enzyme-linked immunosorbent assay, cytokine release by multiplex assay, cytokine gene expression by quantitative polymerase chain reaction, and phosphorylated SFKs level by western blot. The results demonstrated that a SFKs activator, pYEEI (YGRKKRRQRRREPQY(PO3H2)EEIPIYL) alone, did not alter CGRP release or the inflammatory cytokine interleukin-1ß (IL-1ß) gene expression in the mouse TG. In contrast, a SFKs inhibitor, saracatinib, restored CGRP release, the inflammatory cytokines IL-1ß, C-X-C motif ligand 1, C-C motif ligand 2 (CCL2) release, and IL-1ß, CCL2 gene expression when the mouse TG was pre-sensitized with hydrogen peroxide and CGRP respectively. Consistently with this, the phosphorylated SFKs level was increased by both hydrogen peroxide and CGRP in the mouse TG, which was reduced by a CGRP receptor inhibitor BIBN4096 and a protein kinase A (PKA) inhibitor PKI (14-22) Amide. The present study demonstrates that SFKs activity plays a pivotal role in facilitating the crosstalk between CGRP and cytokines by transmitting CGRP receptor/PKA signaling to potentiate TG sensitization and ultimately trigeminovascular sensitization.

A polymorphic transcriptional regulatory domain in the amyotrophic lateral sclerosis risk gene CFAP410 correlates with differential isoform expression.

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.

Non-reference genome transposable elements (TEs) have a significant impact on the progression of the Parkinson's disease.

The pathophysiology of Parkinson's disease (PD) is a complex process of the interaction between genetic and environmental factors. Studies on the genetic component of PD have predominantly focused on single nucleotide polymorphisms (SNPs) using a cross-sectional case-control design in large genome-wide association studies. This approach while giving insight into a significant portion of the genetics of PD does not fully account for all the genetic components resulting in missing heritability. In this study, we approached this problem by focusing on the non-reference genome transposable elements (TEs) and their impact on the progression of PD using a longitudinal study design within the Parkinson's progression markers initiative (PPMI) cohort. We analyzed 2886 Alu repeats, 360 LINE1 and 128 SINE-VNTR-Alus (SVAs) that were called from the whole-genome sequence data which are not within the reference genome. The presence or absence of these non-reference TE variants is known as a retrotransposon insertion polymorphism, and measuring this polymorphism describes the impact of TEs on the traits. The variations for the presence or absence of the non-reference TE elements were modeled to align with the changes in the 114 outcome measures during the five-year follow-up period of the PPMI cohort. Linear mixed-effects models were used, and many TEs were found to have a highly significant effect on the longitudinal changes in the clinically important PD outcomes such as UPDRS subscale II, UPDRS total scores, and modified Schwab and England ADL scale. In addition, the progression of several imaging and functional measures, including the Caudate/Putamen ratio and levodopa equivalent daily dose (LEDD) were also significantly affected by the TEs. In conclusion, this study identified the overwhelming effect of the non-reference TEs on the progression of PD and is a good example of the impact the variations in the "junk DNA" have on complex diseases.

Characterisation of retrotransposon insertion polymorphisms in whole genome sequencing data from individuals with amyotrophic lateral sclerosis.

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.

Characterisation of the Function of a SINE-VNTR-Alu Retrotransposon to Modulate Isoform Expression at the MAPT Locus.

SINE-VNTR-Alu retrotransposons represent one class of transposable elements which contribute to the regulation and evolution of the primate genome and have the potential to be involved in genetic instability and disease progression. However, these polymorphic elements have not been extensively analysed when addressing the missing heritability of neurodegenerative diseases, including Parkinson's disease (PD) and amyotrophic lateral sclerosis (ALS). SVA_67, a retrotransposon insertion polymorphism, is located in a 1.8 Mb region of high linkage disequilibrium, called the MAPT locus, which is known to contribute to increased risk of developing PD, frontotemporal dementia and other tauopathies. To investigate the role of SVA_67 in directing differential gene expression at this locus, we characterised the impact of SVA_67 allele dosage on isoform expression of several genes in the MAPT locus using the datasets from both the Parkinson's Progression Markers Initiative and New York Genome Center Consortium Target ALS cohort. The Parkinson's data was from gene expression in the blood and the ALS data from a variety of CNS regions and allowed us to demonstrate that SVA_67 presence or absence correlated with both isoform- and tissue-specific expression of multiple genes at this locus. This study highlights the importance of addressing SVA polymorphism in disease genetics to gain insight into a better understanding of the role of these regulatory domains to a variety of neurodegenerative diseases.

Longitudinal intronic RNA-Seq analysis of Parkinson's disease patients reveals disease-specific nascent transcription.

Transcriptomic studies usually focus on either gene or exon-based annotations, and only limited experiments have reported changes in reads mapping to introns. The analysis of intronic reads allows the detection of nascent transcription that is not influenced by steady-state RNA levels and provides information on actively transcribed genes. Here, we describe substantial intronic transcriptional changes in Parkinson's disease (PD) patients compared to healthy controls (CO) at two different timepoints; at the time of diagnosis (BL) and three years later (V08). We used blood RNA-Seq data from the Parkinson's Progression Markers Initiative (PPMI) cohort and identified significantly changed transcription of intronic reads only in PD patients during this follow-up period. In CO subjects, only nine transcripts demonstrated differentially expressed introns between visits. However, in PD patients, 4873 transcripts had differentially expressed introns at visit V08 compared to BL, many of them in genes previously associated with neurodegenerative diseases, such as LRRK2, C9orf72, LGALS3, KANSL1AS1, and ALS2. In addition, at the time of diagnosis (BL visit), we identified 836 transcripts (e.g. SNCA, DNAJC19, PRRG4) and at visit V08, 2184 transcripts (e.g. PINK1, GBA, ALS2, PLEKHM1) with differential intronic expression specific to PD patients. In contrast, reads mapping to exonic regions demonstrated little variation indicating highly specific changes only in intronic transcription. Our study demonstrated that PD is characterized by substantial changes in the nascent transcription, and description of these changes could help to understand the molecular pathology underpinning this disease.