Population-scale tissue transcriptomics maps long non-coding RNAs to complex disease.

Long non-coding RNA (lncRNA) genes have well-established and important impacts on molecular and cellular functions. However, among the thousands of lncRNA genes, it is still a major challenge to identify the subset with disease or trait relevance. To systematically characterize these lncRNA genes, we used Genotype Tissue Expression (GTEx) project v8 genetic and multi-tissue transcriptomic data to profile the expression, genetic regulation, cellular contexts, and trait associations of 14,100 lncRNA genes across 49 tissues for 101 distinct complex genetic traits. Using these approaches, we identified 1,432 lncRNA gene-trait associations, 800 of which were not explained by stronger effects of neighboring protein-coding genes. This included associations between lncRNA quantitative trait loci and inflammatory bowel disease, type 1 and type 2 diabetes, and coronary artery disease, as well as rare variant associations to body mass index.

The kinesin motor KIF1C is a putative transporter of the exon junction complex in neuronal cells.

Neurons critically depend on regulated RNA localization and tight control of spatio-temporal gene expression to maintain their morphological and functional integrity. Mutations in the kinesin motor protein gene KIF1C cause Hereditary Spastic Paraplegia, an autosomal recessive disease leading to predominant degeneration of the long axons of central motoneurons. In this study we aimed to gain insight into the molecular function of KIF1C and understand how KIF1C dysfunction contributes to motoneuron degeneration. We used affinity proteomics in neuronally differentiated neuroblastoma cells (SH-SY5Y) to identify the protein complex associated with KIF1C in neuronal cells; candidate interactions were then validated by immunoprecipitation and mislocalization of putative KIF1C cargoes was studied by immunostainings. We found KIF1C to interact with all core components of the exon junction complex (EJC); expression of mutant KIF1C in neuronal cells leads to loss of the typical localization distally in neurites. Instead, EJC core components accumulate in the pericentrosomal region, here co-localizing with mutant KIF1C. These findings suggest KIF1C as a neuronal transporter of the EJC. Interestingly, the binding of KIF1C to the EJC is RNA-mediated, as treatment with RNAse prior to immunoprecipitation almost completely abolishes the interaction. Silica-based solid-phase extraction of UV-crosslinked RNA-protein complexes furthermore supports direct interaction of KIF1C with RNA, as recently also demonstrated for kinesin heavy chain. Taken together, our findings are consistent with a model where KIF1C transports mRNA in an EJC-bound and therefore transcriptionally silenced state along neurites, thus providing the missing link between the EJC and mRNA localization in neurons.

Candidate-based screening via gene modulation in human neurons and astrocytes implicates FERMT2 in Aß and TAU proteostasis.

Large-scale 'omic' studies investigating the pathophysiological processes that lead to Alzheimer's disease (AD) dementia have identified an increasing number of susceptibility genes, many of which are poorly characterized and have not previously been implicated in AD. Here, we evaluated the utility of human induced pluripotent stem cell-derived neurons and astrocytes as tools to systematically test AD-relevant cellular phenotypes following perturbation of candidate genes identified by genome-wide studies. Lentiviral-mediated delivery of shRNAs was used to modulate expression of 66 genes in astrocytes and 52 genes in induced neurons. Five genes (CNN2, GBA, GSTP1, MINT2 and FERMT2) in neurons and nine genes (CNN2, ITGB1, MINT2, SORL1, VLDLR, NPC1, NPC2, PSAP and SCARB2) in astrocytes significantly altered extracellular amyloid-ß (Aß) levels. Knockdown of AP3M2, CNN2, GSTP1, NPC1, NPC2, PSAP and SORL1 reduced interleukin-6 levels in astrocytes. Only knockdown of FERMT2 led to a reduction in the proportion of TAU that is phosphorylated. Further, CRISPR-Cas9 targeting of FERMT2 in both familial AD (fAD) and fAD-corrected human neurons validated the findings of reduced extracellular Aß. Interestingly, FERMT2 reduction had no effect on the Aß42:40 ratio in corrected neurons and a reduction of phospho-tau, but resulted in an elevation in Aß42:40 ratio and no reduction in phospho-tau in fAD neurons. Taken together, this study has prioritized 15 genes as being involved in contributing to Aß accumulation, phosphorylation of tau and/or cytokine secretion, and, as illustrated with FERMT2, it sets the stage for further cell-type-specific dissection of the role of these genes in AD.

Sensitized mutagenesis screen in Factor V Leiden mice identifies thrombosis suppressor loci.

Factor V Leiden (F5L ) is a common genetic risk factor for venous thromboembolism in humans. We conducted a sensitized N-ethyl-N-nitrosourea (ENU) mutagenesis screen for dominant thrombosuppressor genes based on perinatal lethal thrombosis in mice homozygous for F5L (F5L/L ) and haploinsufficient for tissue factor pathway inhibitor (Tfpi+/- ). F8 deficiency enhanced the survival of F5L/LTfpi+/- mice, demonstrating that F5L/LTfpi+/- lethality is genetically suppressible. ENU-mutagenized F5L/L males and F5L/+Tfpi+/- females were crossed to generate 6,729 progeny, with 98 F5L/LTfpi+/- offspring surviving until weaning. Sixteen lines, referred to as "modifier of Factor 5 Leiden (MF5L1-16)," exhibited transmission of a putative thrombosuppressor to subsequent generations. Linkage analysis in MF5L6 identified a chromosome 3 locus containing the tissue factor gene (F3). Although no ENU-induced F3 mutation was identified, haploinsufficiency for F3 (F3+/- ) suppressed F5L/LTfpi+/- lethality. Whole-exome sequencing in MF5L12 identified an Actr2 gene point mutation (p.R258G) as the sole candidate. Inheritance of this variant is associated with suppression of F5L/LTfpi+/- lethality (P = 1.7 × 10-6), suggesting that Actr2p.R258G is thrombosuppressive. CRISPR/Cas9 experiments to generate an independent Actr2 knockin/knockout demonstrated that Actr2 haploinsufficiency is lethal, supporting a hypomorphic or gain-of-function mechanism of action for Actr2p.R258G Our findings identify F8 and the Tfpi/F3 axis as key regulators in determining thrombosis balance in the setting of F5L and also suggest a role for Actr2 in this process.

NMNAT2:HSP90 Complex Mediates Proteostasis in Proteinopathies.

Nicotinamide mononucleotide adenylyl transferase 2 (NMNAT2) is neuroprotective in numerous preclinical models of neurodegeneration. Here, we show that brain nmnat2 mRNA levels correlate positively with global cognitive function and negatively with AD pathology. In AD brains, NMNAT2 mRNA and protein levels are reduced. NMNAT2 shifts its solubility and colocalizes with aggregated Tau in AD brains, similar to chaperones, which aid in the clearance or refolding of misfolded proteins. Investigating the mechanism of this observation, we discover a novel chaperone function of NMNAT2, independent from its enzymatic activity. NMNAT2 complexes with heat shock protein 90 (HSP90) to refold aggregated protein substrates. NMNAT2's refoldase activity requires a unique C-terminal ATP site, activated in the presence of HSP90. Furthermore, deleting NMNAT2 function increases the vulnerability of cortical neurons to proteotoxic stress and excitotoxicity. Interestingly, NMNAT2 acts as a chaperone to reduce proteotoxic stress, while its enzymatic activity protects neurons from excitotoxicity. Taken together, our data indicate that NMNAT2 exerts its chaperone or enzymatic function in a context-dependent manner to maintain neuronal health.

24-hour rhythms of DNA methylation and their relation with rhythms of RNA expression in the human dorsolateral prefrontal cortex.

Circadian rhythms modulate the biology of many human tissues, including brain tissues, and are driven by a near 24-hour transcriptional feedback loop. These rhythms are paralleled by 24-hour rhythms of large portions of the transcriptome. The role of dynamic DNA methylation in influencing these rhythms is uncertain. While recent work in Neurospora suggests that dynamic site-specific circadian rhythms of DNA methylation may play a role in modulating the fungal molecular clock, such rhythms and their relationship to RNA expression have not, to our knowledge, been elucidated in mammalian tissues, including human brain tissues. We hypothesized that 24-hour rhythms of DNA methylation exist in the human brain, and play a role in driving 24-hour rhythms of RNA expression. We analyzed DNA methylation levels in post-mortem human dorsolateral prefrontal cortex samples from 738 subjects. We assessed for 24-hour rhythmicity of 420,132 DNA methylation sites throughout the genome by considering methylation levels as a function of clock time of death and parameterizing these data using cosine functions. We determined global statistical significance by permutation. We then related rhythms of DNA methylation with rhythms of RNA expression determined by RNA sequencing. We found evidence of significant 24-hour rhythmicity of DNA methylation. Regions near transcription start sites were enriched for high-amplitude rhythmic DNA methylation sites, which were in turn time locked to 24-hour rhythms of RNA expression of nearby genes, with the nadir of methylation preceding peak transcript expression by 1-3 hours. Weak ante-mortem rest-activity rhythms were associated with lower amplitude DNA methylation rhythms as were older age and the presence of Alzheimer's disease. These findings support the hypothesis that 24-hour rhythms of DNA methylation, particularly near transcription start sites, may play a role in driving 24-hour rhythms of gene expression in the human dorsolateral prefrontal cortex, and may be affected by age and Alzheimer's disease.

Genetic epistasis regulates amyloid deposition in resilient aging.

INTRODUCTION: The brain-derived neurotrophic factor (BDNF) interacts with important genetic Alzheimer's disease (AD) risk factors. Specifically, variants within the SORL1 gene determine BDNF's ability to reduce amyloid ß (Aß) in vitro. We sought to test whether functional BDNF variation interacts with SORL1 genotypes to influence expression and downstream AD-related processes in humans. METHODS: We analyzed postmortem brain RNA sequencing and neuropathological data for 441 subjects from the Religious Orders Study/Memory and Aging Project and molecular and structural neuroimaging data for 1285 subjects from the Alzheimer's Disease Neuroimaging Initiative. RESULTS: We found one SORL1 RNA transcript strongly regulated by SORL1-BDNF interactions in elderly without pathological AD and showing stronger associations with diffuse than neuritic Aß plaques. The same SORL1-BDNF interactions also significantly influenced Aß load as measured with [18F]Florbetapir positron emission tomography. DISCUSSION: Our results bridge the gap between risk and resilience factors for AD, demonstrating interdependent roles of established SORL1 and BDNF functional genotypes.

Long runs of homozygosity are enriched for deleterious variation.

Exome sequencing offers the potential to study the population-genomic variables that underlie patterns of deleterious variation. Runs of homozygosity (ROH) are long stretches of consecutive homozygous genotypes probably reflecting segments shared identically by descent as the result of processes such as consanguinity, population size reduction, and natural selection. The relationship between ROH and patterns of predicted deleterious variation can provide insight into the way in which these processes contribute to the maintenance of deleterious variants. Here, we use exome sequencing to examine ROH in relation to the distribution of deleterious variation in 27 individuals of varying levels of apparent inbreeding from 6 human populations. A significantly greater fraction of all genome-wide predicted damaging homozygotes fall in ROH than would be expected from the corresponding fraction of nondamaging homozygotes in ROH (p < 0.001). This pattern is strongest for long ROH (p < 0.05). ROH, and especially long ROH, harbor disproportionately more deleterious homozygotes than would be expected on the basis of the total ROH coverage of the genome and the genomic distribution of nondamaging homozygotes. The results accord with a hypothesis that recent inbreeding, which generates long ROH, enables rare deleterious variants to exist in homozygous form. Thus, just as inbreeding can elevate the occurrence of rare recessive diseases that represent homozygotes for strongly deleterious mutations, inbreeding magnifies the occurrence of mildly deleterious variants as well.

Methylation profiles in peripheral blood CD4+ lymphocytes versus brain: The relation to Alzheimer's disease pathology.

INTRODUCTION: We investigated the change in DNA methylation in peripheral blood CD4+ lymphocytes over time, examined the relation between CD4+ lymphocytes and brain methylation, and compared their associations with AD pathology. METHODS: Genome-wide methylation was measured three times in 41 older persons using Illumina Infinium HumanMethylation450 array. The two CD4+ lymphocytes measures were at study baseline and proximate to death. Brain tissue came from frozen dorsolateral prefrontal cortex. RESULTS: Global methylation features were conserved across tissue. At individual CpG sites, methylation level was concordant between the two CD4+ lymphocytes but more diffuse between CD4+ lymphocytes and brain. Previous associations of brain methylation with neuritic plaques at target methylation sites were not replicated in CD4+ lymphocytes. DISCUSSION: There is no strong evidence of change in CD4+ lymphocytes methylation among older persons over an average of 7.5 years. Methylation associations with AD pathology found in neocortex are not directly reflected in CD4+ lymphocytes.

Dominant mutation of CCDC78 in a unique congenital myopathy with prominent internal nuclei and atypical cores.

Congenital myopathies are clinically and genetically heterogeneous diseases that typically present in childhood with hypotonia and weakness and are most commonly defined by changes observed in muscle biopsy. Approximately 40% of congenital myopathies are currently genetically unresolved. We identified a family with dominantly inherited congenital myopathy characterized by distal weakness and biopsy changes that included core-like areas and increased internalized nuclei. To identify the causative genetic abnormality in this family, we performed linkage analysis followed by whole-exome capture and next-generation sequencing. A splice-acceptor variant in previously uncharacterized CCDC78 was detected in affected individuals and absent in unaffected family members and > 10,000 controls. This variant alters RNA-transcript processing and results in a 222 bp in-frame insertion. CCDC78 is expressed in skeletal muscle, enriched in the perinuclear region and the triad, and found in intracellular aggregates in patient muscle. Modeling of the CCDC78 mutation in zebrafish resulted in changes mirroring the human disease that included altered motor function and abnormal muscle ultrastructure. Using a combination of linkage analysis, next-generation sequencing, and modeling in the zebrafish, we have identified a CCDC78 mutation associated with a unique myopathy with prominent internal nuclei and atypical cores.

Unraveling Axonal Transcriptional Landscapes: Insights from iPSC-Derived Cortical Neurons and Implications for Motor Neuron Degeneration.

Neuronal function and pathology are deeply influenced by the distinct molecular profiles of the axon and soma. Traditional studies have often overlooked these differences due to the technical challenges of compartment specific analysis. In this study, we employ a robust RNA-sequencing (RNA-seq) approach, using microfluidic devices, to generate high-quality axonal transcriptomes from iPSC-derived cortical neurons (CNs). We achieve high specificity of axonal fractions, ensuring sample purity without contamination. Comparative analysis revealed a unique and specific transcriptional landscape in axonal compartments, characterized by diverse transcript types, including protein-coding mRNAs, ribosomal proteins (RPs), mitochondrial-encoded RNAs, and long non-coding RNAs (lncRNAs). Previous works have reported the existence of transcription factors (TFs) in the axon. Here, we detect a subset of previously unreported TFs specific to the axon and indicative of their active participation in transcriptional regulation. To investigate transcripts and pathways essential for central motor neuron (MN) degeneration and maintenance we analyzed KIF1C-knockout (KO) CNs, modeling hereditary spastic paraplegia (HSP), a disorder associated with prominent length-dependent degeneration of central MN axons. We found that several key factors crucial for survival and health were absent in KIF1C-KO axons, highlighting a possible role of these also in other neurodegenerative diseases. Taken together, this study underscores the utility of microfluidic devices in studying compartment-specific transcriptomics in human neuronal models and reveals complex molecular dynamics of axonal biology. The impact of KIF1C on the axonal transcriptome not only deepens our understanding of MN diseases but also presents a promising avenue for exploration of compartment specific disease mechanisms.

SR-TWAS: Leveraging Multiple Reference Panels to Improve TWAS Power by Ensemble Machine Learning.

Multiple reference panels of a given tissue or multiple tissues often exist, and multiple regression methods could be used for training gene expression imputation models for TWAS. To leverage expression imputation models (i.e., base models) trained with multiple reference panels, regression methods, and tissues, we develop a Stacked Regression based TWAS (SR-TWAS) tool which can obtain optimal linear combinations of base models for a given validation transcriptomic dataset. Both simulation and real studies showed that SR-TWAS improved power, due to increased effective training sample sizes and borrowed strength across multiple regression methods and tissues. Leveraging base models across multiple reference panels, tissues, and regression methods, our real application studies identified 6 independent significant risk genes for Alzheimer's disease (AD) dementia for supplementary motor area tissue and 9 independent significant risk genes for Parkinson's disease (PD) for substantia nigra tissue. Relevant biological interpretations were found for these significant risk genes.

Unravelling undiagnosed rare disease cases by HiFi long-read genome sequencing.

Solve-RD is a pan-European rare disease (RD) research program that aims to identify disease-causing genetic variants in previously undiagnosed RD families. We utilised 10-fold coverage HiFi long-read sequencing (LRS) for detecting causative structural variants (SVs), single nucleotide variants (SNVs), insertion-deletions (InDels), and short tandem repeat (STR) expansions in extensively studied RD families without clear molecular diagnoses. Our cohort includes 293 individuals from 114 genetically undiagnosed RD families selected by European Rare Disease Network (ERN) experts. Of these, 21 families were affected by so-called 'unsolvable' syndromes for which genetic causes remain unknown, and 93 families with at least one individual affected by a rare neurological, neuromuscular, or epilepsy disorder without genetic diagnosis despite extensive prior testing. Clinical interpretation and orthogonal validation of variants in known disease genes yielded thirteen novel genetic diagnoses due to de novo and rare inherited SNVs, InDels, SVs, and STR expansions. In an additional four families, we identified a candidate disease-causing SV affecting several genes including an MCF2 / FGF13 fusion and PSMA3 deletion. However, no common genetic cause was identified in any of the 'unsolvable' syndromes. Taken together, we found (likely) disease-causing genetic variants in 13.0% of previously unsolved families and additional candidate disease-causing SVs in another 4.3% of these families. In conclusion, our results demonstrate the added value of HiFi long-read genome sequencing in undiagnosed rare diseases.

Differential gene expression analysis based on linear mixed model corrects false positive inflation for studying quantitative traits.

Differential gene expression (DGE) analysis has been widely employed to identify genes expressed differentially with respect to a trait of interest using RNA sequencing (RNA-Seq) data. Recent RNA-Seq data with large samples pose challenges to existing DGE methods, which were mainly developed for dichotomous traits and small sample sizes. Especially, existing DGE methods are likely to result in inflated false positive rates. To address this gap, we employed a linear mixed model (LMM) that has been widely used in genetic association studies for DGE analysis of quantitative traits. We first applied the LMM method to the discovery RNA-Seq data of dorsolateral prefrontal cortex (DLPFC) tissue (n = 632) with four continuous measures of Alzheimer's Disease (AD) cognitive and neuropathologic traits. The quantile-quantile plots of p-values showed that false positive rates were well calibrated by LMM, whereas other methods not accounting for sample-specific mixed effects led to serious inflation. LMM identified 37 potentially significant genes with differential expression in DLPFC for at least one of the AD traits, 17 of which were replicated in the additional RNA-Seq data of DLPFC, supplemental motor area, spinal cord, and muscle tissues. This application study showed not only well calibrated DGE results by LMM, but also possibly shared gene regulatory mechanisms of AD traits across different relevant tissues.

Dissecting the human leptomeninges at single-cell resolution.

Emerging evidence shows that the meninges conduct essential immune surveillance and immune defense at the brain border, and the dysfunction of meningeal immunity contributes to aging and neurodegeneration. However, no study exists on the molecular properties of cell types within human leptomeninges. Here, we provide single nuclei profiling of dissected postmortem leptomeninges from aged individuals. We detect diverse cell types, including unique meningeal endothelial, mural, and fibroblast subtypes. For immune cells, we show that most T cells express CD8 and bear characteristics of tissue-resident memory T cells. We also identify distinct subtypes of border-associated macrophages (BAMs) that display differential gene expressions from microglia and express risk genes for Alzheimer's Disease (AD), as nominated by genome-wide association studies (GWAS). We discover cell-type-specific differentially expressed genes in individuals with Alzheimer's dementia, particularly in fibroblasts and BAMs. Indeed, when cultured, leptomeningeal cells display the signature of ex vivo AD fibroblasts upon amyloid-ß treatment. We further explore ligand-receptor interactions within the leptomeningeal niche and computationally infer intercellular communications in AD. Thus, our study establishes a molecular map of human leptomeningeal cell types, providing significant insight into the border immune and fibrotic responses in AD.

Inferring protein expression changes from mRNA in Alzheimer's dementia using deep neural networks.

Identifying the molecular systems and proteins that modify the progression of Alzheimer's disease and related dementias (ADRD) is central to drug target selection. However, discordance between mRNA and protein abundance, and the scarcity of proteomic data, has limited our ability to advance candidate targets that are mainly based on gene expression. Therefore, by using a deep neural network that predicts protein abundance from mRNA expression, here we attempt to track the early protein drivers of ADRD. Specifically, by applying the clei2block deep learning model to 1192 brain RNA-seq samples, we identify protein modules and disease-associated expression changes that were not directly observed at the mRNA level. Moreover, pseudo-temporal trajectory inference based on the predicted proteome became more closely correlated with cognitive decline and hippocampal atrophy compared to RNA-based trajectories. This suggests that the predicted changes in protein expression could provide a better molecular representation of ADRD progression. Furthermore, overlaying clinical traits on protein pseudotime trajectory identifies protein modules altered before cognitive impairment. These results demonstrate how our method can be used to identify potential early protein drivers and possible drug targets for treating and/or preventing ADRD.

A cortical immune network map identifies distinct microglial transcriptional programs associated with ß-amyloid and Tau pathologies.

Microglial dysfunction has been proposed as one of the many cellular mechanisms that can contribute to the development of Alzheimer's disease (AD). Here, using a transcriptional network map of the human frontal cortex, we identify five modules of co-expressed genes related to microglia and assess their role in the neuropathologic features of AD in 540 subjects from two cohort studies of brain aging. Two of these transcriptional programs-modules 113 and 114-relate to the accumulation of ß-amyloid, while module 5 relates to tau pathology. We replicate these associations in brain epigenomic data and in two independent datasets. In terms of tau, we propose that module 5, a marker of activated microglia, may lead to tau accumulation and subsequent cognitive decline. We validate our model further by showing that three representative module 5 genes (ACADVL, TRABD, and VASP) encode proteins that are upregulated in activated microglia in AD.

Atlas of RNA editing events affecting protein expression in aged and Alzheimer's disease human brain tissue.

RNA editing is a feature of RNA maturation resulting in the formation of transcripts whose sequence differs from the genome template. Brain RNA editing may be altered in Alzheimer's disease (AD). Here, we analyzed data from 1,865 brain samples covering 9 brain regions from 1,074 unrelated subjects on a transcriptome-wide scale to identify inter-regional differences in RNA editing. We expand the list of known brain editing events by identifying 58,761 previously unreported events. We note that only a small proportion of these editing events are found at the protein level in our proteome-wide validation effort. We also identified the occurrence of editing events associated with AD dementia, neuropathological measures and longitudinal cognitive decline in: SYT11, MCUR1, SOD2, ORAI2, HSDL2, PFKP, and GPRC5B. Thus, we present an extended reference set of brain RNA editing events, identify a subset that are found to be expressed at the protein level, and extend the narrative of transcriptomic perturbation in AD to RNA editing.

Enrichment of SARM1 alleles encoding variants with constitutively hyperactive NADase in patients with ALS and other motor nerve disorders.

SARM1, a protein with critical NADase activity, is a central executioner in a conserved programme of axon degeneration. We report seven rare missense or in-frame microdeletion human SARM1 variant alleles in patients with amyotrophic lateral sclerosis (ALS) or other motor nerve disorders that alter the SARM1 auto-inhibitory ARM domain and constitutively hyperactivate SARM1 NADase activity. The constitutive NADase activity of these seven variants is similar to that of SARM1 lacking the entire ARM domain and greatly exceeds the activity of wild-type SARM1, even in the presence of nicotinamide mononucleotide (NMN), its physiological activator. This rise in constitutive activity alone is enough to promote neuronal degeneration in response to otherwise non-harmful, mild stress. Importantly, these strong gain-of-function alleles are completely patient-specific in the cohorts studied and show a highly significant association with disease at the single gene level. These findings of disease-associated coding variants that alter SARM1 function build on previously reported genome-wide significant association with ALS for a neighbouring, more common SARM1 intragenic single nucleotide polymorphism (SNP) to support a contributory role of SARM1 in these disorders. A broad phenotypic heterogeneity and variable age-of-onset of disease among patients with these alleles also raises intriguing questions about the pathogenic mechanism of hyperactive SARM1 variants.