Experimental and Computational Methods for Allelic Imbalance Analysis from Single-Nucleus RNA-seq Data.

Single-cell RNA-seq (scRNA-seq) is emerging as a powerful tool for understanding gene function across diverse cells. Recently, this has included the use of allele-specific expression (ASE) analysis to better understand how variation in the human genome affects RNA expression at the single-cell level. We reasoned that because intronic reads are more prevalent in single-nucleus RNA-Seq (snRNA-Seq), and introns are under lower purifying selection and thus enriched for genetic variants, that snRNA-seq should facilitate single-cell analysis of ASE. Here we demonstrate how experimental and computational choices can improve the results of allelic imbalance analysis. We explore how experimental choices, such as RNA source, read length, sequencing depth, genotyping, etc., impact the power of ASE-based methods. We developed a new suite of computational tools to process and analyze scRNA-seq and snRNA-seq for ASE. As hypothesized, we extracted more ASE information from reads in intronic regions than those in exonic regions and show how read length can be set to increase power. Additionally, hybrid selection improved our power to detect allelic imbalance in genes of interest. We also explored methods to recover allele-specific isoform expression levels from both long- and short-read snRNA-seq. To further investigate ASE in the context of human disease, we applied our methods to a Parkinson's disease cohort of 94 individuals and show that ASE analysis had more power than eQTL analysis to identify significant SNP/gene pairs in our direct comparison of the two methods. Overall, we provide an end-to-end experimental and computational approach for future studies.

Parkin coregulates glutathione metabolism in adult mammalian brain.

We recently discovered that the expression of PRKN, a young-onset Parkinson disease-linked gene, confers redox homeostasis. To further examine the protective effects of parkin in an oxidative stress model, we first combined the loss of prkn with Sod2 haploinsufficiency in mice. Although adult prkn-/-//Sod2± animals did not develop dopamine cell loss in the S. nigra, they had more reactive oxidative species and a higher concentration of carbonylated proteins in the brain; bi-genic mice also showed a trend for more nitrotyrosinated proteins. Because these redox changes were seen in the cytosol rather than mitochondria, we next explored the thiol network in the context of PRKN expression. We detected a parkin deficiency-associated increase in the ratio of reduced glutathione (GSH) to oxidized glutathione (GSSG) in murine brain, PRKN-linked human cortex and several cell models. This shift resulted from enhanced recycling of GSSG back to GSH via upregulated glutathione reductase activity; it also correlated with altered activities of redox-sensitive enzymes in mitochondria isolated from mouse brain (e.g., aconitase-2; creatine kinase). Intriguingly, human parkin itself showed glutathione-recycling activity in vitro and in cells: For each GSSG dipeptide encountered, parkin regenerated one GSH molecule and was S-glutathionylated by the other (GSSG + P-SH [Formula: see text] GSH + P-S-SG), including at cysteines 59, 95 and 377. Moreover, parkin's S-glutathionylation was reversible by glutaredoxin activity. In summary, we found that PRKN gene expression contributes to the network of available thiols in the cell, including by parkin's participation in glutathione recycling, which involves a reversible, posttranslational modification at select cysteines. Further, parkin's impact on redox homeostasis in the cytosol can affect enzyme activities elsewhere, such as in mitochondria. We posit that antioxidant functions of parkin may explain many of its previously described, protective effects in vertebrates and invertebrates that are unrelated to E3 ligase activity.

Age-associated insolubility of parkin in human midbrain is linked to redox balance and sequestration of reactive dopamine metabolites.

The mechanisms by which parkin protects the adult human brain from Parkinson disease remain incompletely understood. We hypothesized that parkin cysteines participate in redox reactions and that these are reflected in its posttranslational modifications. We found that in post mortem human brain, including in the Substantia nigra, parkin is largely insoluble after age 40 years; this transition is linked to its oxidation, such as at residues Cys95 and Cys253. In mice, oxidative stress induces posttranslational modifications of parkin cysteines that lower its solubility in vivo. Similarly, oxidation of recombinant parkin by hydrogen peroxide (H2O2) promotes its insolubility and aggregate formation, and in exchange leads to the reduction of H2O2. This thiol-based redox activity is diminished by parkin point mutants, e.g., p.C431F and p.G328E. In prkn-null mice, H2O2 levels are increased under oxidative stress conditions, such as acutely by 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine toxin exposure or chronically due to a second, genetic hit; H2O2 levels are also significantly increased in parkin-deficient human brain. In dopamine toxicity studies, wild-type parkin, but not disease-linked mutants, protects human dopaminergic cells, in part through lowering H2O2. Parkin also neutralizes reactive, electrophilic dopamine metabolites via adduct formation, which occurs foremost at the primate-specific residue Cys95. Further, wild-type but not p.C95A-mutant parkin augments melanin formation in vitro. By probing sections of adult, human midbrain from control individuals with epitope-mapped, monoclonal antibodies, we found specific and robust parkin reactivity that co-localizes with neuromelanin pigment, frequently within LAMP-3/CD63+ lysosomes. We conclude that oxidative modifications of parkin cysteines are associated with protective outcomes, which include the reduction of H2O2, conjugation of reactive dopamine metabolites, sequestration of radicals within insoluble aggregates, and increased melanin formation. The loss of these complementary redox effects may augment oxidative stress during ageing in dopamine-producing cells of mutant PRKN allele carriers, thereby enhancing the risk of Parkinson's-linked neurodegeneration.