Deconvolution of polygenic risk score in single cells unravels cellular and molecular heterogeneity of complex human diseases.

Polygenic risk scores (PRSs) are commonly used for predicting an individual's genetic risk of complex diseases. Yet, their implication for disease pathogenesis remains largely limited. Here, we introduce scPRS, a geometric deep learning model that constructs single-cell-resolved PRS leveraging reference single-cell chromatin accessibility profiling data to enhance biological discovery as well as disease prediction. Real-world applications across multiple complex diseases, including type 2 diabetes (T2D), hypertrophic cardiomyopathy (HCM), and Alzheimer's disease (AD), showcase the superior prediction power of scPRS compared to traditional PRS methods. Importantly, scPRS not only predicts disease risk but also uncovers disease-relevant cells, such as hormone-high alpha and beta cells for T2D, cardiomyocytes and pericytes for HCM, and astrocytes, microglia and oligodendrocyte progenitor cells for AD. Facilitated by a layered multi-omic analysis, scPRS further identifies cell-type-specific genetic underpinnings, linking disease-associated genetic variants to gene regulation within corresponding cell types. We substantiate the disease relevance of scPRS-prioritized HCM genes and demonstrate that the suppression of these genes in HCM cardiomyocytes is rescued by Mavacamten treatment. Additionally, we establish a novel microglia-specific regulatory relationship between the AD risk variant rs7922621 and its target genes ANXA11 and TSPAN14. We further illustrate the detrimental effects of suppressing these two genes on microglia phagocytosis. Our work provides a multi-tasking, interpretable framework for precise disease prediction and systematic investigation of the genetic, cellular, and molecular basis of complex diseases, laying the methodological foundation for single-cell genetics.

Deep learning modeling of rare noncoding genetic variants in human motor neurons defines CCDC146 as a therapeutic target for ALS.

Amyotrophic lateral sclerosis (ALS) is a fatal and incurable neurodegenerative disease caused by the selective and progressive death of motor neurons (MNs). Understanding the genetic and molecular factors influencing ALS survival is crucial for disease management and therapeutics. In this study, we introduce a deep learning-powered genetic analysis framework to link rare noncoding genetic variants to ALS survival. Using data from human induced pluripotent stem cell (iPSC)-derived MNs, this method prioritizes functional noncoding variants using deep learning, links cis-regulatory elements (CREs) to target genes using epigenomics data, and integrates these data through gene-level burden tests to identify survival-modifying variants, CREs, and genes. We apply this approach to analyze 6,715 ALS genomes, and pinpoint four novel rare noncoding variants associated with survival, including chr7:76,009,472:C>T linked to CCDC146. CRISPR-Cas9 editing of this variant increases CCDC146 expression in iPSC-derived MNs and exacerbates ALS-specific phenotypes, including TDP-43 mislocalization. Suppressing CCDC146 with an antisense oligonucleotide (ASO), showing no toxicity, completely rescues ALS-associated survival defects in neurons derived from sporadic ALS patients and from carriers of the ALS-associated G4C2-repeat expansion within C9ORF72. ASO targeting of CCDC146 may be a broadly effective therapeutic approach for ALS. Our framework provides a generic and powerful approach for studying noncoding genetics of complex human diseases.

KCNJ2 inhibition mitigates mechanical injury in a human brain organoid model of traumatic brain injury.

Traumatic brain injury (TBI) strongly correlates with neurodegenerative disease. However, it remains unclear which neurodegenerative mechanisms are intrinsic to the brain and which strategies most potently mitigate these processes. We developed a high-intensity ultrasound platform to inflict mechanical injury to induced pluripotent stem cell (iPSC)-derived cortical organoids. Mechanically injured organoids elicit classic hallmarks of TBI, including neuronal death, tau phosphorylation, and TDP-43 nuclear egress. We found that deep-layer neurons were particularly vulnerable to injury and that TDP-43 proteinopathy promotes cell death. Injured organoids derived from C9ORF72 amyotrophic lateral sclerosis/frontotemporal dementia (ALS/FTD) patients displayed exacerbated TDP-43 dysfunction. Using genome-wide CRISPR interference screening, we identified a mechanosensory channel, KCNJ2, whose inhibition potently mitigated neurodegenerative processes in vitro and in vivo, including in C9ORF72 ALS/FTD organoids. Thus, targeting KCNJ2 may reduce acute neuronal death after brain injury, and we present a scalable, genetically flexible cerebral organoid model that may enable the identification of additional modifiers of mechanical stress.

CLIP-Seq analysis enables the design of protective ribosomal RNA bait oligonucleotides against C9ORF72 ALS/FTD poly-GR pathophysiology.

Amyotrophic lateral sclerosis and frontotemporal dementia patients with a hexanucleotide repeat expansion in C9ORF72 (C9-HRE) accumulate poly-GR and poly-PR aggregates. The pathogenicity of these arginine-rich dipeptide repeats (R-DPRs) is thought to be driven by their propensity to bind low-complexity domains of multivalent proteins. However, the ability of R-DPRs to bind native RNA and the significance of this interaction remain unclear. Here, we used computational and experimental approaches to characterize the physicochemical properties of R-DPRs and their interaction with RNA. We find that poly-GR predominantly binds ribosomal RNA (rRNA) in cells and exhibits an interaction that is predicted to be energetically stronger than that for associated ribosomal proteins. Critically, modified rRNA "bait" oligonucleotides restore poly-GR-associated ribosomal deficits and ameliorate poly-GR toxicity in patient neurons and Drosophila models. Our work strengthens the hypothesis that ribosomal function is impaired by R-DPRs, highlights a role for direct rRNA binding in mediating ribosomal dysfunction, and presents a strategy for protecting against C9-HRE pathophysiological mechanisms.

Analysis of proteome-wide degradation dynamics in ALS SOD1 iPSC-derived patient neurons reveals disrupted VCP homeostasis.

Mutations in SOD1 cause amyotrophic lateral sclerosis (ALS) through gain-of-function effects, yet the mechanisms by which misfolded mutant SOD1 (mutSOD1) protein impairs human motor neurons (MNs) remain unclear. Here, we use induced-pluripotent-stem-cell-derived MNs coupled to metabolic stable isotope labeling and mass spectrometry to investigate proteome-wide degradation dynamics. We find several proteins, including the ALS-causal valosin-containing protein (VCP), which predominantly acts in proteasome degradation and autophagy, that degrade slower in mutSOD1 relative to isogenic control MNs. The interactome of VCP is altered in mutSOD1 MNs in vitro, while VCP selectively accumulates in the affected motor cortex of ALS-SOD1 patients. Overexpression of VCP rescues mutSOD1 toxicity in MNs in vitro and in a C. elegans model in vivo, in part due to its ability to modulate the degradation of insoluble mutSOD1. Our results demonstrate that VCP contributes to mutSOD1-dependent degeneration, link two distinct ALS-causal genes, and highlight selective protein degradation impairment in ALS pathophysiology.

The C9ORF72 repeat expansion alters neurodevelopment.

Genetic mutations that cause adult-onset neurodegenerative diseases are often expressed during embryonic stages, but it is unclear whether they alter neurodevelopment and how this might influence disease onset. Here, we show that the most common cause of frontotemporal dementia (FTD) and amyotrophic lateral sclerosis (ALS), a repeat expansion in C9ORF72, restricts neural stem cell proliferation and reduces cortical and thalamic size in utero. Surprisingly, a repeat expansion-derived dipeptide repeat protein (DPR) not known to reduce neuronal viability plays a key role in impairing neurodevelopment. Pharmacologically mimicking the effects of the repeat expansion on neurodevelopment increases susceptibility of C9ORF72 mice to motor defects. Thus, the C9ORF72 repeat expansion stunts development of the brain regions prominently affected in C9ORF72 FTD/ALS patients.

Transcriptomic analysis of feeder-free culture system for maintaining naïve-state pluripotency in human pluripotent stem cells.

Background: Human pluripotent stem cells (hPSCs) such as embryonic stem cells (ESCs) and induced pluripotent stem cells (PSCs) have the capacity of self-renewal and multilineage differentiation in vitro. Conventional hPSCs, which are in a primed state, can produce various types of differentiated cells. However, the variability in their degree of pluripotency and differentiation propensities, which is influenced by the inductive methods and culture conditions, limit their availability. Therefore, PSCs in a naïve state are a promising source of PSCs. Methods: We recently developed a culture system for naïve hPSCs using an inhibitor of the NOTCH signaling pathway and a histone H3 methyltransferase disruptor. This culture system requires feeder cells for stably maintaining the naïve hPSCs. We aimed to develop a culture system for hPSCs that could maintain pluripotency under feeder-free conditions. Results: We used two inhibitors to develop an alternative feeder-free culture system to obtain naïve hPSCs. The naïve cells underwent stable cellular proliferation and were positive for naïve stem cell markers; in addition, they could differentiate into the three germ layers. These feeder-free dome-shaped induced pluripotent stem cells (FFDS-iPSCs) have characteristics similar to that of naïve-like PSCs. Conclusions: The naive hPSCs under feeder-free conditions could ensure supply of cells for various applications in regenerative medicine and disease modeling.

Downregulation of Hsp90 and the antimicrobial peptide Mtk suppresses poly(GR)-induced neurotoxicity in C9ORF72-ALS/FTD.

GGGGCC repeat expansion in the C9ORF72 gene is the most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Repeat RNAs can be translated into dipeptide repeat proteins, including poly(GR), whose mechanisms of action remain largely unknown. In an RNA-seq analysis of poly(GR) toxicity in Drosophila, we found that several antimicrobial peptide genes, such as metchnikowin (Mtk), and heat shock protein (Hsp) genes are activated. Mtk knockdown in the fly eye or in all neurons suppresses poly(GR) neurotoxicity. These findings suggest a cell-autonomous role of Mtk in neurodegeneration. Hsp90 knockdown partially rescues both poly(GR) toxicity in flies and neurodegeneration in C9ORF72 motor neurons derived from induced pluripotent stem cells (iPSCs). Topoisomerase II (TopoII) regulates poly(GR)-induced upregulation of Hsp90 and Mtk. TopoII knockdown also suppresses poly(GR) toxicity in Drosophila and improves survival of C9ORF72 iPSC-derived motor neurons. These results suggest potential novel therapeutic targets for C9ORF72-ALS/FTD.

SYF2 suppression mitigates neurodegeneration in models of diverse forms of ALS.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease caused by many diverse genetic etiologies. Although therapeutics that specifically target causal mutations may rescue individual types of ALS, such approaches cannot treat most patients since they have unknown genetic etiology. Thus, there is a critical need for therapeutic strategies that rescue multiple forms of ALS. Here, we combine phenotypic chemical screening on a diverse cohort of ALS patient-derived neurons with bioinformatic analysis of large chemical and genetic perturbational datasets to identify broadly effective genetic targets for ALS. We show that suppressing the gene-encoding, spliceosome-associated factor SYF2 alleviates TDP-43 aggregation and mislocalization, improves TDP-43 activity, and rescues C9ORF72 and causes sporadic ALS neuron survival. Moreover, Syf2 suppression ameliorates neurodegeneration, neuromuscular junction loss, and motor dysfunction in TDP-43 mice. Thus, suppression of spliceosome-associated factors such as SYF2 may be a broadly effective therapeutic approach for ALS.

PIKFYVE inhibition mitigates disease in models of diverse forms of ALS.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease that results from many diverse genetic causes. Although therapeutics specifically targeting known causal mutations may rescue individual types of ALS, these approaches cannot treat most cases since they have unknown genetic etiology. Thus, there is a pressing need for therapeutic strategies that rescue multiple forms of ALS. Here, we show that pharmacological inhibition of PIKFYVE kinase activates an unconventional protein clearance mechanism involving exocytosis of aggregation-prone proteins. Reducing PIKFYVE activity ameliorates ALS pathology and extends survival of animal models and patient-derived motor neurons representing diverse forms of ALS including C9ORF72, TARDBP, FUS, and sporadic. These findings highlight a potential approach for mitigating ALS pathogenesis that does not require stimulating macroautophagy or the ubiquitin-proteosome system.

Radiogenomics of C9orf72 Expansion Carriers Reveals Global Transposable Element Derepression and Enables Prediction of Thalamic Atrophy and Clinical Impairment.

Hexanucleotide repeat expansion (HRE) within C9orf72 is the most common genetic cause of frontotemporal dementia (FTD). Thalamic atrophy occurs in both sporadic and familial FTD but is thought to distinctly affect HRE carriers. Separately, emerging evidence suggests widespread derepression of transposable elements (TEs) in the brain in several neurodegenerative diseases, including C9orf72 HRE-mediated FTD (C9-FTD). Whether TE activation can be measured in peripheral blood and how the reduction in peripheral C9orf72 expression observed in HRE carriers relates to atrophy and clinical impairment remain unknown. We used FreeSurfer software to assess the effects of C9orf72 HRE and clinical diagnosis (n = 78 individuals, male and female) on atrophy of thalamic nuclei. We also generated a novel, human, whole-blood RNA-sequencing dataset to determine the relationships among peripheral C9orf72 expression, TE activation, thalamic atrophy, and clinical severity (n = 114 individuals, male and female). We confirmed global thalamic atrophy and reduced C9orf72 expression in HRE carriers. Moreover, we identified disproportionate atrophy of the right mediodorsal lateral nucleus in HRE carriers and showed that C9orf72 expression associated with clinical severity, independent of thalamic atrophy. Strikingly, we found global peripheral activation of TEs, including the human endogenous LINE-1 element L1HS L1HS levels were associated with atrophy of multiple pulvinar nuclei, a thalamic region implicated in C9-FTD. Integration of peripheral transcriptomic and neuroimaging data from human HRE carriers revealed atrophy of specific thalamic nuclei, demonstrated that C9orf72 levels relate to clinical severity, and identified marked derepression of TEs, including L1HS, which predicted atrophy of FTD-relevant thalamic nuclei.SIGNIFICANCE STATEMENT Pathogenic repeat expansion in C9orf72 is the most frequent genetic cause of FTD and amyotrophic lateral sclerosis (ALS; C9-FTD/ALS). The clinical, neuroimaging, and pathologic features of C9-FTD/ALS are well characterized, whereas the intersections of transcriptomic dysregulation and brain structure remain largely unexplored. Herein, we used a novel radiogenomic approach to examine the relationship between peripheral blood transcriptomics and thalamic atrophy, a neuroimaging feature disproportionately impacted in C9-FTD/ALS. We confirmed reduction of C9orf72 in blood and found broad dysregulation of transposable elements-genetic elements typically repressed in the human genome-in symptomatic C9orf72 expansion carriers, which associated with atrophy of thalamic nuclei relevant to FTD. C9orf72 expression was also associated with clinical severity, suggesting that peripheral C9orf72 levels capture disease-relevant information.

A "multi-omics" analysis of blood-brain barrier and synaptic dysfunction in APOE4 mice.

Apolipoprotein E4 (APOE4), the main susceptibility gene for Alzheimer's disease, leads to blood-brain barrier (BBB) breakdown in humans and mice. Remarkably, BBB dysfunction predicts cognitive decline and precedes synaptic deficits in APOE4 human carriers. How APOE4 affects BBB and synaptic function at a molecular level, however, remains elusive. Using single-nucleus RNA-sequencing and phosphoproteome and proteome analysis, we show that APOE4 compared with APOE3 leads to an early disruption of the BBB transcriptome in 2-3-mo-old APOE4 knock-in mice, followed by dysregulation in protein signaling networks controlling cell junctions, cytoskeleton, clathrin-mediated transport, and translation in brain endothelium, as well as transcription and RNA splicing suggestive of DNA damage in pericytes. Changes in BBB signaling mechanisms paralleled an early, progressive BBB breakdown and loss of pericytes, which preceded postsynaptic interactome disruption and behavioral deficits that developed 2-5 mo later. Thus, dysregulated signaling mechanisms in endothelium and pericytes in APOE4 mice reflect a molecular signature of a progressive BBB failure preceding changes in synaptic function and behavior.

Whole-genome sequencing reveals that variants in the Interleukin 18 Receptor Accessory Protein 3'UTR protect against ALS.

The noncoding genome is substantially larger than the protein-coding genome but has been largely unexplored by genetic association studies. Here, we performed region-based rare variant association analysis of >25,000 variants in untranslated regions of 6,139 amyotrophic lateral sclerosis (ALS) whole genomes and the whole genomes of 70,403 non-ALS controls. We identified interleukin-18 receptor accessory protein (IL18RAP) 3' untranslated region (3'UTR) variants as significantly enriched in non-ALS genomes and associated with a fivefold reduced risk of developing ALS, and this was replicated in an independent cohort. These variants in the IL18RAP 3'UTR reduce mRNA stability and the binding of double-stranded RNA (dsRNA)-binding proteins. Finally, the variants of the IL18RAP 3'UTR confer a survival advantage for motor neurons because they dampen neurotoxicity of human induced pluripotent stem cell (iPSC)-derived microglia bearing an ALS-associated expansion in C9orf72, and this depends on NF-κB signaling. This study reveals genetic variants that protect against ALS by reducing neuroinflammation and emphasizes the importance of noncoding genetic association studies.

Repairing the blood-brain barrier.

Engineered Wnt ligands specifically target blood-brain barrier function.

Engineering skeletal muscle tissues with advanced maturity improves synapse formation with human induced pluripotent stem cell-derived motor neurons.

To develop effective cures for neuromuscular diseases, human-relevant in vitro models of neuromuscular tissues are critically needed to probe disease mechanisms on a cellular and molecular level. However, previous attempts to co-culture motor neurons and skeletal muscle have resulted in relatively immature neuromuscular junctions (NMJs). In this study, NMJs formed by human induced pluripotent stem cell (hiPSC)-derived motor neurons were improved by optimizing the maturity of the co-cultured muscle tissue. First, muscle tissues engineered from the C2C12 mouse myoblast cell line, cryopreserved primary human myoblasts, and freshly isolated primary chick myoblasts on micromolded gelatin hydrogels were compared. After three weeks, only chick muscle tissues remained stably adhered to hydrogels and exhibited progressive increases in myogenic index and stress generation, approaching values generated by native muscle tissue. After three weeks of co-culture with hiPSC-derived motor neurons, engineered chick muscle tissues formed NMJs with increasing co-localization of pre- and postsynaptic markers as well as increased frequency and magnitude of synaptic activity, surpassing structural and functional maturity of previous in vitro models. Engineered chick muscle tissues also demonstrated increased expression of genes related to sarcomere maturation and innervation over time, revealing new insights into the molecular pathways that likely contribute to enhanced NMJ formation. These approaches for engineering advanced neuromuscular tissues with relatively mature NMJs and interrogating their structure and function have many applications in neuromuscular disease modeling and drug development.

ELAVL4, splicing, and glutamatergic dysfunction precede neuron loss in MAPT mutation cerebral organoids.

Frontotemporal dementia (FTD) because of MAPT mutation causes pathological accumulation of tau and glutamatergic cortical neuronal death by unknown mechanisms. We used human induced pluripotent stem cell (iPSC)-derived cerebral organoids expressing tau-V337M and isogenic corrected controls to discover early alterations because of the mutation that precede neurodegeneration. At 2 months, mutant organoids show upregulated expression of MAPT, glutamatergic signaling pathways, and regulators, including the RNA-binding protein ELAVL4, and increased stress granules. Over the following 4 months, mutant organoids accumulate splicing changes, disruption of autophagy function, and build-up of tau and P-tau-S396. By 6 months, tau-V337M organoids show specific loss of glutamatergic neurons as seen in individuals with FTD. Mutant neurons are susceptible to glutamate toxicity, which can be rescued pharmacologically by the PIKFYVE kinase inhibitor apilimod. Our results demonstrate a sequence of events that precede neurodegeneration, revealing molecular pathways associated with glutamate signaling as potential targets for therapeutic intervention in FTD.

p53 is a central regulator driving neurodegeneration caused by C9orf72 poly(PR).

The most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) is a GGGGCC repeat expansion in the C9orf72 gene. We developed a platform to interrogate the chromatin accessibility landscape and transcriptional program within neurons during degeneration. We provide evidence that neurons expressing the dipeptide repeat protein poly(proline-arginine), translated from the C9orf72 repeat expansion, activate a highly specific transcriptional program, exemplified by a single transcription factor, p53. Ablating p53 in mice completely rescued neurons from degeneration and markedly increased survival in a C9orf72 mouse model. p53 reduction also rescued axonal degeneration caused by poly(glycine-arginine), increased survival of C9orf72 ALS/FTD-patient-induced pluripotent stem cell (iPSC)-derived motor neurons, and mitigated neurodegeneration in a C9orf72 fly model. We show that p53 activates a downstream transcriptional program, including Puma, which drives neurodegeneration. These data demonstrate a neurodegenerative mechanism dynamically regulated through transcription-factor-binding events and provide a framework to apply chromatin accessibility and transcription program profiles to neurodegeneration.

The M1311V variant of ATP7A is associated with impaired trafficking and copper homeostasis in models of motor neuron disease.

Disruption in copper homeostasis causes a number of cognitive and motor deficits. Wilson's disease and Menkes disease are neurodevelopmental disorders resulting from mutations in the copper transporters ATP7A and ATP7B, with ATP7A mutations also causing occipital horn syndrome, and distal motor neuropathy. A 65 year old male presenting with brachial amyotrophic diplegia and diagnosed with amyotrophic lateral sclerosis (ALS) was found to harbor a p.Met1311Val (M1311V) substitution variant in ATP7A. ALS is a fatal neurodegenerative disease associated with progressive muscle weakness, synaptic deficits and degeneration of upper and lower motor neurons. To investigate the potential contribution of the ATP7AM1311V variant to neurodegeneration, we obtained and characterized both patient-derived fibroblasts and patient-derived induced pluripotent stem cells differentiated into motor neurons (iPSC-MNs), and compared them to control cell lines. We found reduced localization of ATP7AM1311V to the trans-Golgi network (TGN) at basal copper levels in patient-derived fibroblasts and iPSC-MNs. In addition, redistribution of ATP7AM1311V out of the TGN in response to increased extracellular copper was defective in patient fibroblasts. This manifested in enhanced intracellular copper accumulation and reduced survival of ATP7AM1311V fibroblasts. iPSC-MNs harboring the ATP7AM1311V variant showed decreased dendritic complexity, aberrant spontaneous firing, and decreased survival. Finally, expression of the ATP7AM1311V variant in Drosophila motor neurons resulted in motor deficits. Apilimod, a drug that targets vesicular transport and recently shown to enhance survival of C9orf72-ALS/FTD iPSC-MNs, also increased survival of ATP7AM1311V iPSC-MNs and reduced motor deficits in Drosophila expressing ATP7AM1311V. Taken together, these observations suggest that ATP7AM1311V negatively impacts its role as a copper transporter and impairs several aspects of motor neuron function and morphology.

Protocol for Differentiation of Human iPSCs into Pulmonary Neuroendocrine Cells.

Pulmonary neuroendocrine cells (PNECs) are sensory cells within the lung airway epithelia. Here, we provide a detailed protocol for generating induced PNECs (iPNECs) from human induced pluripotent stem cells (iPSCs). The cellular and molecular profile of iPNECs resembles primary human PNECs. Primary human PNECs are exceedingly rare, comprising only 1% of the adult lung. Therefore, a self-renewing source of patient-specific iPNECs facilitates the creation of reproducible human cellular models to study lung diseases characterized by PNEC dysfunction. For complete details on the use and execution of this protocol, please refer to Hor et al. (2020).