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.

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.

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.

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.

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.

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.

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.

Locus specific reduction of L1 expression in the cortices of individuals with amyotrophic lateral sclerosis.

The activation and dysregulation of retrotransposons has been identified in the CNS of individuals with the fatal neurodegenerative disorder Amyotrophic lateral sclerosis (ALS). This includes elements from multiple different families and subfamilies of retrotransposons, however there is limited knowledge of the specific loci from which this expression occurs in ALS. The long interspersed element-1 (L1) is the only autonomous retrotransposon in the human genome and members of this family of elements maintain the ability to mobilise. Despite L1s contributing to 17% of the human genome only 80-100 L1s encode the required proteins for mobilisation and are retrotransposition competent. Identifying the specific loci from which L1 expression occurs will inform on the potential functional consequences of their expression, such as the potential for somatic retrotransposition or DNA damage caused by the endonuclease activity of the ORF2 protein of the L1. Here we characterised L1 loci expression using the L1EM tool ( https://github.com/FenyoLab/L1EM ) in RNA sequencing data from 518 samples across four tissues (motor cortex, frontal cortex, cerebellum and cervical spinal cord) in the Target ALS cohort obtained from the New York Genome Center. There was a significant reduction in total intact L1 expression (those that encode functional proteins) in two brain regions of individuals with ALS compared to controls and clustering of the ALS brain regions occurred based on their intact L1 expression profile. Although overall the levels of L1 expression were reduced in ALS/ALS with other neurological disorder (ND) there were individuals in which L1s were expressed at much higher levels than the rest of the ALS/ALSND cohort. Expressed L1 loci were more frequently located in introns compared to those not expressed and the level of L1 expression positively correlated with the expression of the gene in which it was located. Significant differences were observed in the expression profiles of L1s in ALS and specific features of these elements, such as location in the genome and whether or not they are intact, were significantly associated with those that were expressed in the cohort.

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.

Expression Quantitative Trait Loci (eQTLs) Associated with Retrotransposons Demonstrate their Modulatory Effect on the Transcriptome.

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.

Reference SVA insertion polymorphisms are associated with Parkinson's Disease progression and differential gene expression.

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.

Variable number tandem repeats - Their emerging role in sickness and health.

Understanding the mechanisms regulating tissue specific and stimulus inducible regulation is at the heart of understanding human biology and how this translates to wellbeing, the ageing process, and disease progression. Polymorphic DNA variation is superimposed as an extra layer of complexity in such processes which underpin our individuality and are the focus of personalized medicine. This review focuses on the role and action of repetitive DNA, specifically variable number tandem repeats and SINE-VNTR-Alu domains, highlighting their role in modification of gene structure and gene expression in addition to their polymorphic nature being a genetic modifier of disease risk and progression. Although the literature focuses on their role in disease, it illustrates their potential to be major contributors to normal physiological function. To date, these elements have been under-reported in genomic analysis due to the difficulties in their characterization with short read DNA sequencing methods. However, recent advances in long read sequencing methods should resolve these problems allowing for a greater understanding of their contribution to a host of genomic and functional mechanisms underlying physiology and disease.

Transcript Variants of Genes Involved in Neurodegeneration Are Differentially Regulated by the APOE and MAPT Haplotypes.

Genetic variations at the Apolipoprotein E (ApoE) and microtubule-associated protein tau (MAPT) loci have been implicated in multiple neurogenerative diseases, but their exact molecular mechanisms are unclear. In this study, we performed transcript level linear modelling using the blood whole transcriptome data and genotypes of the 570 subjects in the Parkinson's Progression Markers Initiative (PPMI) cohort. ApoE, MAPT haplotypes and two SNPs at the SNCA locus (rs356181, rs3910105) were used to detect expression quantitative trait loci eQTLs associated with the transcriptome and differential usage of transcript isoforms. As a result, we identified 151 genes associated with the genotypic variations, 29 cis and 122 trans eQTL positions. Profound effect with genome-wide significance of ApoE e4 haplotype on the expression of TOMM40 transcripts was identified. This finding potentially explains in part the frequently established genetic association with the APOE e4 haplotypes in neurodegenerative diseases. Moreover, MAPT haplotypes had significant differential impact on 23 transcripts from the 17q21.31 and 17q24.1 loci. MAPT haplotypes had also the largest up-regulating (256) and the largest down-regulating (-178) effect sizes measured as ß values on two different transcripts from the same gene (LRRC37A2). Intronic SNP in the SNCA gene, rs3910105, differentially induced expression of three SNCA isoforms. In conclusion, this study established clear association between well-known haplotypic variance and transcript specific regulation in the blood. APOE e4 and MAPT H1/H2 haplotypic variants are associated with the expression of several genes related to the neurodegeneration.

CRISPR Deletion of a SVA Retrotransposon Demonstrates Function as a cis-Regulatory Element at the TRPV1/TRPV3 Intergenic Region.

SINE-VNTR-Alu (SVA) retrotransposons are a subclass of transposable elements (TEs) that exist only in primate genomes. TE insertions can be co-opted as cis-regulatory elements (CREs); however, the regulatory potential of SVAs has predominantly been demonstrated using bioinformatic approaches and reporter gene assays. The objective of this study was to demonstrate SVA cis-regulatory activity by CRISPR (clustered regularly interspaced short palindromic repeats) deletion and subsequent measurement of direct effects on local gene expression. We identified a region on chromosome 17 that was enriched with human-specific SVAs. Comparative gene expression analysis at this region revealed co-expression of TRPV1 and TRPV3 in multiple human tissues, which was not observed in mouse, highlighting key regulatory differences between the two species. Furthermore, the intergenic region between TRPV1 and TRPV3 coding sequences contained a human specific SVA insertion located upstream of the TRPV3 promoter and downstream of the 3' end of TRPV1, highlighting this SVA as a candidate to study its potential cis-regulatory activity on both genes. Firstly, we generated SVA reporter gene constructs and demonstrated their transcriptional regulatory activity in HEK293 cells. We then devised a dual-targeting CRISPR strategy to facilitate the deletion of this entire SVA sequence and generated edited HEK293 clonal cell lines containing homozygous and heterozygous SVA deletions. In edited homozygous ∆SVA clones, we observed a significant decrease in both TRPV1 and TRPV3 mRNA expression, compared to unedited HEK293. In addition, we also observed an increase in the variability of mRNA expression levels in heterozygous ∆SVA clones. Overall, in edited HEK293 with SVA deletions, we observed a disruption to the co-expression of TRPV1 and TRPV3. Here we provide an example of a human specific SVA with cis-regulatory activity in situ, supporting the role of SVA retrotransposons as contributors to species-specific gene expression.