Mapping medically relevant RNA isoform diversity in the aged human frontal cortex with deep long-read RNA-seq.

Determining whether the RNA isoforms from medically relevant genes have distinct functions could facilitate direct targeting of RNA isoforms for disease treatment. Here, as a step toward this goal for neurological diseases, we sequenced 12 postmortem, aged human frontal cortices (6 Alzheimer disease cases and 6 controls; 50% female) using one Oxford Nanopore PromethION flow cell per sample. We identified 1,917 medically relevant genes expressing multiple isoforms in the frontal cortex where 1,018 had multiple isoforms with different protein-coding sequences. Of these 1,018 genes, 57 are implicated in brain-related diseases including major depression, schizophrenia, Parkinson's disease and Alzheimer disease. Our study also uncovered 53 new RNA isoforms in medically relevant genes, including several where the new isoform was one of the most highly expressed for that gene. We also reported on five mitochondrially encoded, spliced RNA isoforms. We found 99 differentially expressed RNA isoforms between cases with Alzheimer disease and controls.

Surveying the landscape of RNA isoform diversity and expression across 9 GTEx tissues using long-read sequencing data.

Even though alternative RNA splicing was discovered nearly 50 years ago (1977), we still understand very little about most isoforms arising from a single gene, including in which tissues they are expressed and if their functions differ. Human gene annotations suggest remarkable transcriptional complexity, with approximately 252,798 distinct RNA isoform annotations from 62,710 gene bodies (Ensembl v109; 2023), emphasizing the need to understand their biological effects. For example, 256 gene bodies have ≥50 annotated isoforms and 30 have ≥100, where one protein-coding gene (MAPK10) even has 192 distinct RNA isoform annotations. Whether such isoform diversity results from biological redundancy or spurious alternative splicing (i.e., noise), or whether individual isoforms have specialized functions (even if subtle) remains a mystery for most genes. Recent studies by Aguzzoli-Heberle et al., Leung et al., and Glinos et al. demonstrated long-read RNAseq enables improved RNA isoform quantification for essentially any tissue, cell type, or biological condition (e.g., disease, development, aging, etc.), making it possible to better assess individual isoform expression and function. While each study provided important discoveries related to RNA isoform diversity, deeper exploration is needed. We sought to quantify and characterize real isoform usage across tissues (compared to annotations). We used long-read RNAseq data from 58 GTEx samples across nine tissues (three brain, two heart, muscle, lung, liver, and cultured fibroblasts) generated by Glinos et al. and found considerable isoform diversity within and across tissues. Cerebellar hemisphere was the most transcriptionally complex tissue (22,522 distinct isoforms; 3,726 unique); liver was least diverse (12,435 distinct isoforms; 1,039 unique). We highlight gene clusters exhibiting high tissue-specific isoform diversity per tissue (e.g., TPM1 expresses 19 in heart's atrial appendage). We also validated 447 of the 700 new isoforms discovered by Aguzzoli-Heberle et al. and found that 88 were expressed in all nine tissues, while 58 were specific to a single tissue. This study represents a broad survey of the RNA isoform landscape, demonstrating isoform diversity across nine tissues and emphasizes the need to better understand how individual isoforms from a single gene body contribute to human health and disease.

Using deep long-read RNAseq in Alzheimer's disease brain to assess medical relevance of RNA isoform diversity.

Due to alternative splicing, human protein-coding genes average over eight RNA isoforms, resulting in nearly four distinct protein coding sequences per gene. Long-read RNAseq (IsoSeq) enables more accurate quantification of isoforms, shedding light on their specific roles. To assess the medical relevance of measuring RNA isoform expression, we sequenced 12 aged human frontal cortices (6 Alzheimer's disease cases and 6 controls; 50% female) using one Oxford Nanopore PromethION flow cell per sample. Our study uncovered 53 new high-confidence RNA isoforms in medically relevant genes, including several where the new isoform was one of the most highly expressed for that gene. Specific examples include WDR4 (61%; microcephaly), MYL3 (44%; hypertrophic cardiomyopathy), and MTHFS (25%; major depression, schizophrenia, bipolar disorder). Other notable genes with new high-confidence isoforms include CPLX2 (10%; schizophrenia, epilepsy) and MAOB (9%; targeted for Parkinson's disease treatment). We identified 1,917 medically relevant genes expressing multiple isoforms in human frontal cortex, where 1,018 had multiple isoforms with different protein coding sequences, demonstrating the need to better understand how individual isoforms from a single gene body are involved in human health and disease, if at all. Exactly 98 of the 1,917 genes are implicated in brain-related diseases, including Alzheimer's disease genes such as APP (Aß precursor protein; five), MAPT (tau protein; four), and BIN1 (eight). As proof of concept, we also found 99 differentially expressed RNA isoforms between Alzheimer's cases and controls, despite the genes themselves not exhibiting differential expression. Our findings highlight the significant knowledge gaps in RNA isoform diversity and their medical relevance. Deep long-read RNA sequencing will be necessary going forward to fully comprehend the medical relevance of individual isoforms for a "single" gene.

The Polygenic Risk Score Knowledge Base offers a centralized online repository for calculating and contextualizing polygenic risk scores.

The process of identifying suitable genome-wide association (GWA) studies and formatting the data to calculate multiple polygenic risk scores on a single genome can be laborious. Here, we present a centralized polygenic risk score calculator currently containing over 250,000 genetic variant associations from the NHGRI-EBI GWAS Catalog for users to easily calculate sample-specific polygenic risk scores with comparable results to other available tools. Polygenic risk scores are calculated either online through the Polygenic Risk Score Knowledge Base (PRSKB; https://prs.byu.edu ) or via a command-line interface. We report study-specific polygenic risk scores across the UK Biobank, 1000 Genomes, and the Alzheimer's Disease Neuroimaging Initiative (ADNI), contextualize computed scores, and identify potentially confounding genetic risk factors in ADNI. We introduce a streamlined analysis tool and web interface to calculate and contextualize polygenic risk scores across various studies, which we anticipate will facilitate a wider adaptation of polygenic risk scores in future disease research.

Transcriptional response to VZV infection is modulated by RNA polymerase III in lung epithelial cell lines.

Ancestral RNA polymerase III (Pol III) is a multi-subunit polymerase responsible for transcription of short non-coding RNA, such as double-stranded short interspersed nuclear elements (SINEs). Although SINE ncRNAs are generally transcriptionally repressed, they can be induced in response to viral infections and can stimulate immune signaling pathways. Indeed, mutations in RNA Pol III have been associated with poor antiviral interferon response following infection with varicella zoster virus (VZV). In this study, we probed the role of Pol III transcripts in the detection and initial immune response to VZV by characterizing the transcriptional response following VZV infection of wild type A549 lung epithelial cells as well as A549 cells lacking specific RNA sensors MAVS and TLR3, or interferon-stimulated genes RNase L and PKR in presence or absence of functional RNA Pol III. Multiple components of the antiviral sensing and interferon signaling pathways were involved in restricting VZV replication in lung epithelial cells thus suggesting an innate defense system with built-in redundancy. In addition, RNA Pol III silencing altered the antiviral transcriptional program indicating that it plays an essential role in the sensing of VZV infection.

Web-Based Protein Interactions Calculator Identifies Likely Proteome Coevolution with Alzheimer's Disease-Associated Proteins.

Protein-protein functional interactions arise from either transitory or permanent biomolecular associations and often lead to the coevolution of the interacting residues. Although mutual information has traditionally been used to identify coevolving residues within the same protein, its application between coevolving proteins remains largely uncharacterized. Therefore, we developed the Protein Interactions Calculator (PIC) to efficiently identify coevolving residues between two protein sequences using mutual information. We verified the algorithm using 2102 known human protein interactions and 233 known bacterial protein interactions, with a respective 1975 and 252 non-interacting protein controls. The average PIC score for known human protein interactions was 4.5 times higher than non-interacting proteins (p = 1.03 × 10-108) and 1.94 times higher in bacteria (p = 1.22 × 10-35). We then used the PIC scores to determine the probability that two proteins interact. Using those probabilities, we paired 37 Alzheimer's disease-associated proteins with 8608 other proteins and determined the likelihood that each pair interacts, which we report through a web interface. The PIC had significantly higher sensitivity and residue-specific resolution not available in other algorithms. Therefore, we propose that the PIC can be used to prioritize potential protein interactions, which can lead to a better understanding of biological processes and additional therapeutic targets belonging to protein interaction groups.

The Ramp Atlas: facilitating tissue and cell-specific ramp sequence analyses through an intuitive web interface.

Ramp sequences occur when the average translational efficiency of codons near the 5' end of highly expressed genes is significantly lower than the rest of the gene sequence, which counterintuitively increases translational efficiency by decreasing downstream ribosomal collisions. Here, we show that the relative codon adaptiveness within different tissues changes the existence of a ramp sequence without altering the underlying genetic code. We present the first comprehensive analysis of tissue and cell type-specific ramp sequences and report 3108 genes with ramp sequences that change between tissues and cell types, which corresponds with increased gene expression within those tissues and cells. The Ramp Atlas (https://ramps.byu.edu/) allows researchers to query precomputed ramp sequences in 18 388 genes across 62 tissues and 66 cell types and calculate tissue-specific ramp sequences from user-uploaded FASTA files through an intuitive web interface. We used The Ramp Atlas to identify seven SARS-CoV-2 genes and seven human SARS-CoV-2 entry factor genes with tissue-specific ramp sequences that may help explain viral proliferation within those tissues. We anticipate that The Ramp Atlas will facilitate personalized and creative tissue-specific ramp sequence analyses for both human and viral genes that will increase our ability to utilize this often-overlooked regulatory region.

Curated variation benchmarks for challenging medically relevant autosomal genes.

The repetitive nature and complexity of some medically relevant genes poses a challenge for their accurate analysis in a clinical setting. The Genome in a Bottle Consortium has provided variant benchmark sets, but these exclude nearly 400 medically relevant genes due to their repetitiveness or polymorphic complexity. Here, we characterize 273 of these 395 challenging autosomal genes using a haplotype-resolved whole-genome assembly. This curated benchmark reports over 17,000 single-nucleotide variations, 3,600 insertions and deletions and 200 structural variations each for human genome reference GRCh37 and GRCh38 across HG002. We show that false duplications in either GRCh37 or GRCh38 result in reference-specific, missed variants for short- and long-read technologies in medically relevant genes, including CBS, CRYAA and KCNE1. When masking these false duplications, variant recall can improve from 8% to 100%. Forming benchmarks from a haplotype-resolved whole-genome assembly may become a prototype for future benchmarks covering the whole genome.

Pairwise Correlation Analysis of the Alzheimer's Disease Neuroimaging Initiative (ADNI) Dataset Reveals Significant Feature Correlation.

The Alzheimer's Disease Neuroimaging Initiative (ADNI) contains extensive patient measurements (e.g., magnetic resonance imaging [MRI], biometrics, RNA expression, etc.) from Alzheimer's disease (AD) cases and controls that have recently been used by machine learning algorithms to evaluate AD onset and progression. While using a variety of biomarkers is essential to AD research, highly correlated input features can significantly decrease machine learning model generalizability and performance. Additionally, redundant features unnecessarily increase computational time and resources necessary to train predictive models. Therefore, we used 49,288 biomarkers and 793,600 extracted MRI features to assess feature correlation within the ADNI dataset to determine the extent to which this issue might impact large scale analyses using these data. We found that 93.457% of biomarkers, 92.549% of the gene expression values, and 100% of MRI features were strongly correlated with at least one other feature in ADNI based on our Bonferroni corrected α (p-value ≤ 1.40754 × 10-13). We provide a comprehensive mapping of all ADNI biomarkers to highly correlated features within the dataset. Additionally, we show that significant correlation within the ADNI dataset should be resolved before performing bulk data analyses, and we provide recommendations to address these issues. We anticipate that these recommendations and resources will help guide researchers utilizing the ADNI dataset to increase model performance and reduce the cost and complexity of their analyses.

Long-read targeted sequencing uncovers clinicopathological associations for C9orf72-linked diseases.

To examine the length of a hexanucleotide expansion in C9orf72, which represents the most frequent genetic cause of frontotemporal lobar degeneration and motor neuron disease, we employed a targeted amplification-free long-read sequencing technology: No-Amp sequencing. In our cross-sectional study, we assessed cerebellar tissue from 28 well-characterized C9orf72 expansion carriers. We obtained 3507 on-target circular consensus sequencing reads, of which 814 bridged the C9orf72 repeat expansion (23%). Importantly, we observed a significant correlation between expansion sizes obtained using No-Amp sequencing and Southern blotting (P = 5.0 × 10-4). Interestingly, we also detected a significant survival advantage for individuals with smaller expansions (P = 0.004). Additionally, we uncovered that smaller expansions were significantly associated with higher levels of C9orf72 transcripts containing intron 1b (P = 0.003), poly(GP) proteins (P = 1.3 × 10- 5), and poly(GA) proteins (P = 0.005). Thorough examination of the composition of the expansion revealed that its GC content was extremely high (median: 100%) and that it was mainly composed of GGGGCC repeats (median: 96%), suggesting that expanded C9orf72 repeats are quite pure. Taken together, our findings demonstrate that No-Amp sequencing is a powerful tool that enables the discovery of relevant clinicopathological associations, highlighting the important role played by the cerebellar size of the expanded repeat in C9orf72-linked diseases.

Systematic analysis of dark and camouflaged genes reveals disease-relevant genes hiding in plain sight.

BACKGROUND: The human genome contains "dark" gene regions that cannot be adequately assembled or aligned using standard short-read sequencing technologies, preventing researchers from identifying mutations within these gene regions that may be relevant to human disease. Here, we identify regions with few mappable reads that we call dark by depth, and others that have ambiguous alignment, called camouflaged. We assess how well long-read or linked-read technologies resolve these regions. RESULTS: Based on standard whole-genome Illumina sequencing data, we identify 36,794 dark regions in 6054 gene bodies from pathways important to human health, development, and reproduction. Of these gene bodies, 8.7% are completely dark and 35.2% are ≥ 5% dark. We identify dark regions that are present in protein-coding exons across 748 genes. Linked-read or long-read sequencing technologies from 10x Genomics, PacBio, and Oxford Nanopore Technologies reduce dark protein-coding regions to approximately 50.5%, 35.6%, and 9.6%, respectively. We present an algorithm to resolve most camouflaged regions and apply it to the Alzheimer's Disease Sequencing Project. We rescue a rare ten-nucleotide frameshift deletion in CR1, a top Alzheimer's disease gene, found in disease cases but not in controls. CONCLUSIONS: While we could not formally assess the association of the CR1 frameshift mutation with Alzheimer's disease due to insufficient sample-size, we believe it merits investigating in a larger cohort. There remain thousands of potentially important genomic regions overlooked by short-read sequencing that are largely resolved by long-read technologies.

Common DNA Variants Accurately Rank an Individual of Extreme Height.

Polygenic scores (or genetic risk scores) quantify the aggregate of small effects from many common genetic loci that have been associated with a trait through genome-wide association. Polygenic scores were first used successfully in schizophrenia and have since been applied to multiple phenotypes including multiple sclerosis, rheumatoid arthritis, and height. Because human height is an easily-measured and complex polygenic trait, polygenic height scores provide exciting insights into the predictability of aggregate common variant effect on the phenotype. Shawn Bradley is an extremely tall former professional basketball player from Brigham Young University and the National Basketball Association (NBA), measuring 2.29 meters (7'6″, 99.99999th percentile for height) tall, with no known medical conditions. Here, we present a case where a rare combination of common SNPs in one individual results in an extremely high polygenic height score that is correlated with an extreme phenotype. While polygenic scores are not clinically significant in the average case, our findings suggest that for extreme phenotypes, polygenic scores may be more successful for the prediction of individuals.

Long-read sequencing across the C9orf72 'GGGGCC' repeat expansion: implications for clinical use and genetic discovery efforts in human disease.

BACKGROUND: Many neurodegenerative diseases are caused by nucleotide repeat expansions, but most expansions, like the C9orf72 'GGGGCC' (G4C2) repeat that causes approximately 5-7% of all amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) cases, are too long to sequence using short-read sequencing technologies. It is unclear whether long-read sequencing technologies can traverse these long, challenging repeat expansions. Here, we demonstrate that two long-read sequencing technologies, Pacific Biosciences' (PacBio) and Oxford Nanopore Technologies' (ONT), can sequence through disease-causing repeats cloned into plasmids, including the FTD/ALS-causing G4C2 repeat expansion. We also report the first long-read sequencing data characterizing the C9orf72 G4C2 repeat expansion at the nucleotide level in two symptomatic expansion carriers using PacBio whole-genome sequencing and a no-amplification (No-Amp) targeted approach based on CRISPR/Cas9. RESULTS: Both the PacBio and ONT platforms successfully sequenced through the repeat expansions in plasmids. Throughput on the MinION was a challenge for whole-genome sequencing; we were unable to attain reads covering the human C9orf72 repeat expansion using 15 flow cells. We obtained 8× coverage across the C9orf72 locus using the PacBio Sequel, accurately reporting the unexpanded allele at eight repeats, and reading through the entire expansion with 1324 repeats (7941 nucleotides). Using the No-Amp targeted approach, we attained > 800× coverage and were able to identify the unexpanded allele, closely estimate expansion size, and assess nucleotide content in a single experiment. We estimate the individual's repeat region was > 99% G4C2 content, though we cannot rule out small interruptions. CONCLUSIONS: Our findings indicate that long-read sequencing is well suited to characterizing known repeat expansions, and for discovering new disease-causing, disease-modifying, or risk-modifying repeat expansions that have gone undetected with conventional short-read sequencing. The PacBio No-Amp targeted approach may have future potential in clinical and genetic counseling environments. Larger and deeper long-read sequencing studies in C9orf72 expansion carriers will be important to determine heterogeneity and whether the repeats are interrupted by non-G4C2 content, potentially mitigating or modifying disease course or age of onset, as interruptions are known to do in other repeat-expansion disorders. These results have broad implications across all diseases where the genetic etiology remains unclear.

Microglial translational profiling reveals a convergent APOE pathway from aging, amyloid, and tau.

Alzheimer's disease (AD) is an age-associated neurodegenerative disease characterized by amyloidosis, tauopathy, and activation of microglia, the brain resident innate immune cells. We show that a RiboTag translational profiling approach can bypass biases due to cellular enrichment/cell sorting. Using this approach in models of amyloidosis, tauopathy, and aging, we revealed a common set of alterations and identified a central APOE-driven network that converged on CCL3 and CCL4 across all conditions. Notably, aged females demonstrated a significant exacerbation of many of these shared transcripts in this APOE network, revealing a potential mechanism for increased AD susceptibility in females. This study has broad implications for microglial transcriptomic approaches and provides new insights into microglial pathways associated with different pathological aspects of aging and AD.

Poly(GR) impairs protein translation and stress granule dynamics in C9orf72-associated frontotemporal dementia and amyotrophic lateral sclerosis.

The major genetic cause of frontotemporal dementia (FTD) and amyotrophic lateral sclerosis (ALS) is a C9orf72 G4C2 repeat expansion1,2. Proposed mechanisms by which the expansion causes c9FTD/ALS include toxicity from repeat-containing RNA and from dipeptide repeat proteins translated from these transcripts. To investigate the contribution of poly(GR) dipeptide repeat proteins to c9FTD/ALS pathogenesis in a mammalian in vivo model, we generated mice that expressed GFP-(GR)100 in the brain. GFP-(GR)100 mice developed age-dependent neurodegeneration, brain atrophy, and motor and memory deficits through the accumulation of diffuse, cytoplasmic poly(GR). Poly(GR) co-localized with ribosomal subunits and the translation initiation factor eIF3η in GFP-(GR)100 mice and, of importance, in c9FTD/ALS patients. Combined with the differential expression of ribosome-associated genes in GFP-(GR)100 mice, these findings demonstrate poly(GR)-mediated ribosomal distress. Indeed, poly(GR) inhibited canonical and non-canonical protein translation in HEK293T cells, and also induced the formation of stress granules and delayed their disassembly. These data suggest that poly(GR) contributes to c9FTD/ALS by impairing protein translation and stress granule dynamics, consequently causing chronic cellular stress and preventing cells from mounting an effective stress response. Decreasing poly(GR) and/or interrupting interactions between poly(GR) and ribosomal and stress granule-associated proteins may thus represent potential therapeutic strategies to restore homeostasis.

Corticobasal degeneration with TDP-43 pathology presenting with progressive supranuclear palsy syndrome: a distinct clinicopathologic subtype.

Corticobasal degeneration (CBD) is a clinically heterogeneous tauopathy, which has overlapping clinicopathologic and genetic characteristics with progressive supranuclear palsy (PSP). This study aimed to elucidate whether transactive response DNA-binding protein of 43 kDa (TDP-43) pathology contributes to clinicopathologic heterogeneity of CBD. Paraffin-embedded sections of the midbrain, pons, subthalamic nucleus, and basal forebrain from 187 autopsy-confirmed CBD cases were screened with immunohistochemistry for phospho-TDP-43. In cases with TDP-43 pathology, additional brain regions (i.e., precentral, cingulate, and superior frontal gyri, hippocampus, medulla, and cerebellum) were immunostained. Hierarchical clustering analysis was performed based on the topographical distribution and severity of TDP-43 pathology, and clinicopathologic and genetic features were compared between the clusters. TDP-43 pathology was observed in 45% of CBD cases, most frequently in midbrain tegmentum (80% of TDP-43-positive cases), followed by subthalamic nucleus (69%). TDP-43-positive CBD was divided into TDP-limited (52%) and TDP-severe (48%) by hierarchical clustering analysis. TDP-severe patients were more likely to have been diagnosed clinically as PSP compared to TDP-limited and TDP-negative patients (80 vs 32 vs 30%, P < 0.001). The presence of downward gaze palsy was the strongest factor for the antemortem diagnosis of PSP, and severe TDP-43 pathology in the midbrain tectum was strongly associated with downward gaze palsy. In addition, tau burden in the olivopontocerebellar system was significantly greater in TDP-positive than TDP-negative CBD. Genetic analyses revealed that MAPT H1/H1 genotype frequency was significantly lower in TDP-severe than in TDP-negative and TDP-limited CBD (65 vs 89 vs 91%, P < 0.001). The homozygous minor allele frequencies in GRN rs5848 and TMEM106B rs3173615 were not significantly different between the three groups. In conclusion, the present study indicates that CBD with severe TDP-43 pathology is a distinct clinicopathologic subtype of CBD, characterized by PSP-like clinical presentations, severe tau pathology in the olivopontocerebellar system, and low frequency of MAPT H1 haplotype.

An Epigenetic Spin to ALS and FTD.

Amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) are two devastating and lethal neurodegenerative diseases seen comorbidly in up to 15% of patients. Despite several decades of research, no effective treatment or disease-modifying strategies have been developed. We now understand more than before about the genetics and biology behind ALS and FTD, but the genetic etiology for the majority of patients is still unknown and the phenotypic variability observed across patients, even those carrying the same mutation, is enigmatic. Additionally, susceptibility factors leading to neuronal vulnerability in specific central nervous system regions involved in disease are yet to be identified. As the inherited but dynamic epigenome acts as a cell-specific interface between the inherited fixed genome and both cell-intrinsic mechanisms and environmental input, adaptive epigenetic changes might contribute to the ALS/FTD aspects we still struggle to comprehend. This chapter summarizes our current understanding of basic epigenetic mechanisms, how they relate to ALS and FTD, and their potential as therapeutic targets. A clear understanding of the biological mechanisms driving these two currently incurable diseases is urgent-well-needed therapeutic strategies need to be developed soon. Disease-specific epigenetic changes have already been observed in patients and these might be central to this endeavor.

Correction to: Linkage, whole genome sequence, and biological data implicate variants in RAB10 in Alzheimer's disease resilience.

CORRECTION: The original version of this article [1] unfortunately contained a typographical error. The 'Alzheimer's Disease Neuroimaging Initiative' was erroneously included as 'Alzheimer's Disease Neuroimaging Initative' in the author list of the article.

Linkage, whole genome sequence, and biological data implicate variants in RAB10 in Alzheimer's disease resilience.

BACKGROUND: While age and the APOE ε4 allele are major risk factors for Alzheimer's disease (AD), a small percentage of individuals with these risk factors exhibit AD resilience by living well beyond 75 years of age without any clinical symptoms of cognitive decline. METHODS: We used over 200 "AD resilient" individuals and an innovative, pedigree-based approach to identify genetic variants that segregate with AD resilience. First, we performed linkage analyses in pedigrees with resilient individuals and a statistical excess of AD deaths. Second, we used whole genome sequences to identify candidate SNPs in significant linkage regions. Third, we replicated SNPs from the linkage peaks that reduced risk for AD in an independent dataset and in a gene-based test. Finally, we experimentally characterized replicated SNPs. RESULTS: Rs142787485 in RAB10 confers significant protection against AD (p value = 0.0184, odds ratio = 0.5853). Moreover, we replicated this association in an independent series of unrelated individuals (p value = 0.028, odds ratio = 0.69) and used a gene-based test to confirm a role for RAB10 variants in modifying AD risk (p value = 0.002). Experimentally, we demonstrated that knockdown of RAB10 resulted in a significant decrease in Aß42 (p value = 0.0003) and in the Aß42/Aß40 ratio (p value = 0.0001) in neuroblastoma cells. We also found that RAB10 expression is significantly elevated in human AD brains (p value = 0.04). CONCLUSIONS: Our results suggest that RAB10 could be a promising therapeutic target for AD prevention. In addition, our gene discovery approach can be expanded and adapted to other phenotypes, thus serving as a model for future efforts to identify rare variants for AD and other complex human diseases.

Conserved DNA methylation combined with differential frontal cortex and cerebellar expression distinguishes C9orf72-associated and sporadic ALS, and implicates SERPINA1 in disease.

We previously found C9orf72-associated (c9ALS) and sporadic amyotrophic lateral sclerosis (sALS) brain transcriptomes comprise thousands of defects, among which, some are likely key contributors to ALS pathogenesis. We have now generated complementary methylome data and combine these two data sets to perform a comprehensive "multi-omic" analysis to clarify the molecular mechanisms initiating RNA misregulation in ALS. We found that c9ALS and sALS patients have generally distinct but overlapping methylome profiles, and that the c9ALS- and sALS-affected genes and pathways have similar biological functions, indicating conserved pathobiology in disease. Our results strongly implicate SERPINA1 in both C9orf72 repeat expansion carriers and non-carriers, where expression levels are greatly increased in both patient groups across the frontal cortex and cerebellum. SERPINA1 expression is particularly pronounced in C9orf72 repeat expansion carriers for both brain regions, where SERPINA1 levels are strictly down regulated across most human tissues, including the brain, except liver and blood, and are not measurable in E18 mouse brain. The altered biological networks we identified contain critical molecular players known to contribute to ALS pathology, which also interact with SERPINA1. Our comprehensive combined methylation and transcription study identifies new genes and highlights that direct genetic and epigenetic changes contribute to c9ALS and sALS pathogenesis.