Xian Adiconis

Synucleinopathies affect ∼15 million people and are classically divided into neuronal (Parkinson’s disease (PD), dementia with Lewy bodies) and glial (multiple system atrophy) disorders. Here we challenge this dichotomy. We functionally fine-map 90 PD GWAS signals across nine cell types in cortex and substantia nigra using disease-context, population-scale single-nucleus eQTL meta-analysis (N = 1,197), bulk brain eQTL analysis (N = 1,182), and Mendelian randomization. A stringent causal framework integrates single-nucleus allelic imbalance (snASE) with orthogonal validation. We identify 125 functional risk genes for 50 loci—nearly doubling supported genes—and assign genes and cell types to over half of GWAS signals. Unexpectedly, 51% of risk genes are regulated in glia, particularly oligodendrocytes and their precursors. Across cell types, risk converges on a shared glial–neuronal vesiculopathy network. These findings uncover a convergent glial-neuronal risk architecture and establish a single-cell atlas for context-aware gene discovery and precision therapeutics for PD.

Combining allele-specific expression (ASE) analysis with single-cell RNA-seq can elucidate how genomic variation affects RNA expression at the single-cell level. We explore how experimental and computational choices impact the power of ASE-based methods and develop a suite of single-cell ASE computational tools. With single-nucleus RNA-Seq, we extract more ASE information from reads in intronic than exonic regions. We show how read length can increase power and that hybrid selection improves power to detect ASE in targeted genes. We apply our methods to a Parkinson’s disease cohort and show that ASE analysis has more power than eQTL analysis.

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.

Inter-individual genetic variation affects the susceptibility to and progression of many diseases. Efforts to study the molecular mechanisms mediating the impact of human genetic variation on normal development and disease phenotypes are limited, however, by the paucity of faithful cellular human models, and the difficulty of scaling current systems to represent multiple individuals. Here, we present human brain “Chimeroids”, a highly reproducible, multi-donor human brain cortical organoid model generated by the co-development of cells from a panel of individual donors in a single organoid, while maintaining fidelity to endogenous tissue. By reaggregating cells from multiple single-donor organoids at the neural stem or committed progenitor cell stage, we generate Chimeroids in which each donor produces all cell lineages of the cerebral cortex, even when using pluripotent stem cell lines with notable growth biases. We leveraged Chimeroids to investigate inter-individual variation in susceptibility to neurotoxic stressors that exhibit high clinical phenotypic variability: ethanol and the anti-epileptic drug valproic acid. Individual donors varied in both the penetrance of the effect on target cell types, and the molecular phenotype within each affected cell type. Our results show that human genetic background may be an important mediator of neurotoxin susceptibility and introduce Chimeroids as a scalable system for high-throughput investigation of the contribution of human genetic variation to brain development and disease.

Aging is a complex process involving transcriptomic changes associated with deterioration across multiple tissues and organs, including the brain. Recent studies using heterochronic parabiosis have shown that various aspects of aging-associated decline are modifiable or even reversible. To better understand how this occurs, we performed single-cell transcriptomic profiling of young and old mouse brains following parabiosis. For each cell type, we catalogued alterations in gene expression, molecular pathways, transcriptional networks, ligand-receptor interactions, and senescence status. Our analyses identified gene signatures demonstrating that heterochronic parabiosis regulates several hallmarks of aging in a cell-type-specific manner. Brain endothelial cells were found to be especially malleable to this intervention, exhibiting dynamic transcriptional changes that affect vascular structure and function. These findings suggest novel strategies for slowing deterioration and driving regeneration in the aging brain through approaches that do not rely on disease-specific mechanisms or actions of individual circulating factors.

Mutations in non-coding cis -regulatory DNA sequences can alter gene expression, organismal phenotype, and fitness. Fitness landscapes, which map DNA sequence to organismal fitness, are a long-standing goal in biology, but have remained elusive because it is challenging to generalize accurately to the vast space of possible sequences using models built on measurements from a limited number of endogenous regulatory sequences. Here, we construct a sequence-to-expression model for such a landscape and use it to decipher principles of cis -regulatory evolution. Using tens of millions of randomly sampled promoter DNA sequences and their measured expression levels in the yeast Sacccharomyces cerevisiae , we construct a deep transformer neural network model that generalizes with exceptional accuracy, and enables sequence design for gene expression engineering. Using our model, we predict and experimentally validate expression divergence under random genetic drift and strong selection weak mutation regimes, show that conflicting expression objectives in different environments constrain expression adaptation, and find that stabilizing selection on gene expression leads to the moderation of regulatory complexity. We present an approach for detecting selective constraint on gene expression using our model and natural sequence variation, and validate it using observed cis -regulatory diversity across 1,011 yeast strains, cross-species RNA-seq from three different clades, and measured expression-to-fitness curves. Finally, we develop a characterization of regulatory evolvability, use it to visualize fitness landscapes in two dimensions, discover evolvability archetypes, quantify the mutational robustness of individual sequences and highlight the mutational robustness of extant natural regulatory sequence populations. Our work provides a general framework that addresses key questions in the evolution of cis -regulatory sequences.

A multitude of single-cell RNA sequencing methods have been developed in recent years, with dramatic advances in scale and power, and enabling major discoveries and large scale cell mapping efforts. However, these methods have not been systematically and comprehensively benchmarked. Here, we directly compare seven methods for single cell and/or single nucleus profiling from three types of samples – cell lines, peripheral blood mononuclear cells and brain tissue – generating 36 libraries in six separate experiments in a single center. To analyze these datasets, we developed and applied scumi, a flexible computational pipeline that can be used for any scRNA-seq method. We evaluated the methods for both basic performance and for their ability to recover known biological information in the samples. Our study will help guide experiments with the methods in this study as well as serve as a benchmark for future studies and for computational algorithm development.

The mammalian brain is complex, with multiple cell types performing a variety of diverse functions, but exactly how the brain is affected with aging remains largely unknown. Here we performed a single-cell transcriptomic analysis of young and old mouse brains. We provide a comprehensive dataset of aging-related genes, pathways and ligand-receptor interactions in nearly all brain cell types. Our analysis identified gene signatures that vary in a coordinated manner across cell types and gene sets that are regulated in a cell type specific manner, even at times in opposite directions. Thus, our data reveals that aging, rather than inducing a universal program drives a distinct transcriptional course in each cell population. These data provide an important resource for the aging community and highlight key molecular processes, including ribosomal biogenesis, underlying aging. We believe that this large-scale dataset, which is publicly accessible online ( aging-mouse-brain ), will facilitate additional discoveries directed towards understanding and modifying the aging process.

Three members of the Puccini genus, P. triticina ( Pt ), P. striiformis f.sp. tritici ( Pst ), and P. graminis f.sp. tritici ( Pgt ), cause the most common and often most significant foliar diseases of wheat. While similar in biology and life cycle, each species is uniquely adapted and specialized. The genomes of Pt and Pst were sequenced and compared to that of Pgt to identify common and distinguishing gene content, to determine gene variation among wheat rust pathogens, other rust fungi and basidiomycetes, and to identify genes of significance for infection. Pt had the largest genome of the three, estimated at 135 Mb with expansion due to mobile elements and repeats encompassing 50.9% of contig bases; by comparison repeats occupy 31.5% for Pst and 36.5% for Pgt . We find all three genomes are highly heterozygous, with Pst (5.97 SNPs/kb) nearly twice the level detected in Pt (2.57 SNPs/kb) and that previously reported for Pgt . Of 1,358 predicted effectors in Pt , 784 were found expressed across diverse life cycle stages including the sexual stage. Comparison to related fungi highlighted the expansion of gene families involved in transcriptional regulation and nucleotide binding, protein modification, and carbohydrate enzyme degradation. Two allelic homeodomain, HD1 and HD2, pairs and three pheromone receptor ( STE3 ) mating-type genes were identified in each dikaryotic Puccinia species. The HD proteins were active in a heterologous Ustilago maydis mating assay and host induced gene silencing of the HD and STE3 alleles reduced wheat host infection.

Background: Recent studies in budding yeast have shown that antisense transcription occurs at many loci. However, the functional role of antisense transcripts has been demonstrated only in a few cases and it has been suggested that most antisense transcripts may result from promiscuous bi-directional transcription in a dense genome. Results: Here, we use strand-specific RNA sequencing to study anti-sense transcription in Saccharomyces cerevisiae. We detect 1,103 putative antisense transcripts expressed in mid-log phase growth, ranging from 39 short transcripts covering only the 3' UTR of sense genes to 145 long transcripts covering the entire sense open reading frame. Many of these antisense transcripts overlap sense genes that are repressed in mid-log phase and are important in stationary phase, stress response, or meiosis. We validate the differential regulation of 67 antisense transcripts and their sense targets in relevant conditions, including nutrient limitation and environmental stresses. Moreover, we show that several antisense transcripts and, in some cases, their differential expression have been conserved across five species of yeast spanning 150 million years of evolution. Divergence in the regulation of antisense transcripts to two respiratory genes coincides with the evolution of respiro-fermentation. Conclusions: Our work provides support for a global and conserved role for antisense transcription in yeast gene regulation.