Emily Miyoshi

Biological systems are immensely complex, organized into a multi-scale hierarchy of functional units based on tightly regulated interactions between distinct molecules, cells, organs, and organisms. While experimental methods enable transcriptome-wide measurements across millions of cells, popular bioinformatic tools do not support systems-level analysis. Here we present hdWGCNA, a comprehensive framework for analyzing co-expression networks in high-dimensional transcriptomics data such as single-cell and spatial RNA sequencing (RNA-seq). hdWGCNA provides functions for network inference, gene module identification, gene enrichment analysis, statistical tests, and data visualization. Beyond conventional single-cell RNA-seq, hdWGCNA is capable of performing isoform-level network analysis using long-read single-cell data. We showcase hdWGCNA using data from autism spectrum disorder and Alzheimer's disease brain samples, identifying disease-relevant co-expression network modules. hdWGCNA is directly compatible with Seurat, a widely used R package for single-cell and spatial transcriptomics analysis, and we demonstrate the scalability of hdWGCNA by analyzing a dataset containing nearly 1 million cells.

Microglia are strongly implicated in the development and progression of Alzheimer's disease (AD), yet their impact on pathology and lifespan remains unclear. Here we utilize a CSF1R hypomorphic mouse to generate a model of AD that genetically lacks microglia. The resulting microglial-deficient mice exhibit a profound shift from parenchymal amyloid plaques to cerebral amyloid angiopathy (CAA), which is accompanied by numerous transcriptional changes, greatly increased brain calcification and hemorrhages, and premature lethality. Remarkably, a single injection of wild-type microglia into adult mice repopulates the microglial niche and prevents each of these pathological changes. Taken together, these results indicate the protective functions of microglia in reducing CAA, blood-brain barrier dysfunction, and brain calcification. To further understand the clinical implications of these findings, human AD tissue and iPSC-microglia were examined, providing evidence that microglia phagocytose calcium crystals, and this process is impaired by loss of the AD risk gene, TREM2.

Alzheimer's disease (AD) is a devastating neurological disorder characterized by changes in cell-type proportions and consequently marked alterations of the transcriptome. Here we use a data-driven systems biology meta-analytical approach across three human AD cohorts, encompassing six cortical brain regions, and integrate with multi-scale datasets comprising of DNA methylation, histone acetylation, transcriptome- and genome-wide association studies and quantitative trait loci to further characterize the genetic architecture of AD. We perform co-expression network analysis across more than 1200 human brain samples, identifying robust AD-associated dysregulation of the transcriptome, unaltered in normal human aging. We assess the cell-type specificity of AD gene co-expression changes and estimate cell-type proportion changes in human AD by integrating co-expression modules with single-cell transcriptome data generated from 27 321 nuclei from human postmortem prefrontal cortical tissue. We also show that genetic variants of AD are enriched in a microglial AD-associated module and identify key transcription factors regulating co-expressed modules. Additionally, we validate our results in multiple published human AD gene expression datasets, which can be easily accessed using our online resource (https://swaruplab.bio.uci.edu/consensusAD).

The complexity of affected brain regions and cell types is a challenge for Huntington's disease (HD) treatment. Here we use single nucleus RNA sequencing to investigate molecular pathology in the cortex and striatum from R6/2 mice and human HD post-mortem tissue. We identify cell type-specific and -agnostic signatures suggesting oligodendrocytes (OLs) and oligodendrocyte precursors (OPCs) are arrested in intermediate maturation states. OL-lineage regulators OLIG1 and OLIG2 are negatively correlated with CAG length in human OPCs, and ATACseq analysis of HD mouse NeuN-negative cells shows decreased accessibility regulated by OL maturation genes. The data implicates glucose and lipid metabolism in abnormal cell maturation and identify PRKCE and Thiamine Pyrophosphokinase 1 (TPK1) as central genes. Thiamine/biotin treatment of R6/1 HD mice to compensate for TPK1 dysregulation restores OL maturation and rescues neuronal pathology. Our insights into HD OL pathology spans multiple brain regions and link OL maturation deficits to abnormal thiamine metabolism.