Disease Risk Genes Research Articles

Enhancing disease risk gene discovery by integrating transcription factor-linked trans-variants into transcriptome-wide association analyses.

Transcriptome-wide association studies (TWAS) have been successful in identifying disease susceptibility genes by integrating cis-variants predicted gene expression with genome-wide association studies (GWAS) data. However, trans-variants for predicting gene expression remain largely unexplored. Here, we introduce transTF-TWAS, which incorporates transcription factor (TF)-linked trans-variants to enhance model building for TF downstream target genes. Using data from the Genotype-Tissue Expression project, we predict gene expression and alternative splicing and applied these prediction models to large GWAS datasets for breast, prostate, lung cancers and other diseases. We demonstrate that transTF-TWAS outperforms other existing TWAS approaches in both constructing gene expression prediction models and identifying disease-associated genes, as shown by simulations and real data analysis. Our transTF-TWAS approach significantly contributes to the discovery of disease risk genes. Findings from this study shed new light on several genetically driven key TF regulators and their associated TF-gene regulatory networks underlying disease susceptibility.

TET2 germline variants promote kidney disease by impairing DNA repair and activating cytosolic nucleotide sensors

Genome-wide association studies (GWAS) have identified over 800 loci associated with kidney function, yet the specific genes, variants, and pathways involved remain elusive. By integrating kidney function GWAS with human kidney expression and methylation quantitative trait analyses, we identified Ten-Eleven Translocation (TET) DNA demethylase 2 (TET2) as a novel kidney disease risk gene. Utilizing single-cell chromatin accessibility and CRISPR-based genome editing, we highlight GWAS variants that influence TET2 expression in kidney proximal tubule cells. Experiments using kidney/tubule-specific Tet2 knockout mice indicated its protective role in cisplatin-induced acute kidney injury, as well as in chronic kidney disease and fibrosis induced by unilateral ureteral obstruction or adenine diet. Single-cell gene profiling of kidneys from Tet2 knockout mice and TET2-knockdown tubule cells revealed the altered expression of DNA damage repair and chromosome segregation genes, notably including INO80, another kidney function GWAS target gene itself. Remarkably, both TET2-null and INO80-null cells exhibited an increased accumulation of micronuclei after injury, leading to the activation of cytosolic nucleotide sensor cGAS-STING. Genetic deletion of cGAS or STING in kidney tubules, or pharmacological inhibition of STING, protected TET2-null mice from disease development. In conclusion, our findings highlight TET2 and INO80 as key genes in the pathogenesis of kidney diseases, indicating the importance of DNA damage repair mechanisms.

Multi-scalar data integration decoding risk genes for chronic kidney disease

BackgroundChronic Kidney Disease (CKD) impacts over 10% of the global population, and recent advancements in high-throughput analytical technologies are uncovering the complex physiology underlying this condition. By integrating Genome-Wide Association Studies (GWAS), RNA sequencing (RNA-seq/RNA array), and single-cell RNA sequencing (scRNA-seq) data, our study aimed to explore the genes and cell types relevant to CKD traits.MethodsGWAS summary data for end-stage renal failure (ESRD) and decreased eGFR (CKD) with or without diabetes and (micro)proteinuria were obtained from the GWAS Catalog and the UK Biobank (UKB) database. Two gene Expression Omnibus (GEO) transcriptome datasets were used to establish glomerular and tubular gene expression differences between CKD patients and healthy individuals. Two scRNA-seq datasets were utilized to obtain the expression of key genes at the single-cell level. The expression profile, differentially expressed genes (DEGs), gene-gene interaction, and pathway enrichment were analysed for these CKD risk genes.ResultsA total of 779 distinct SNPs were identified from GWAS across different CKD traits, involving 681 genes. While many of these risk genes are specific to the CKD traits of renal failure, decreased eGFR, and (micro)proteinuria, they share common pathways, including extracellular matrix (ECM). ECM modeling was enriched in upregulated glomerular and tubular DEGs from CKD kidneys compared to healthy controls, with the expression of relevant collagen genes, such as COL1A2, prevalent in fibroblasts/myofibroblasts. Additionally, immune responses, including T cell differentiation, were dysregulated in CKD kidneys. The late podocyte signature gene THSD7A was enriched in podocytes but downregulated in CKD. We also highlighted that the regulated risk genes of CKD are mainly expressed in tubular cells and immune cells in the kidney.ConclusionsOur integrated analysis highlight the genes, pathways, and relevant cell types associational with the pathogenesis of kidney traits, as a basis for further mechanistic studies to understand the pathogenesis of CKD.

The Alzheimer's disease risk gene CD2AP functions in dendritic spines by remodeling F-actin.

CD2AP was identified as a genetic risk factor for late-onset Alzheimer's disease (LOAD). However, it is unclear how CD2AP contributes to LOAD synaptic dysfunction underlying AD memory deficits. We have shown that loss of CD2AP function increases β-amyloid (Aβ) endocytic production, but it is unknown whether it contributes to synapse dysfunction. As CD2AP is an actin-binding protein, it may also function in F-actin-rich dendritic spines, which are the excitatory postsynaptic compartments. Here, we demonstrate that CD2AP colocalizes with F-actin in dendritic spines of primary mouse cortical neurons of both sexes. Cell-autonomous depletion of CD2AP specifically reduces spine density and volume, resulting in a functional decrease in synapse formation and neuronal network activity. Post-synaptic reexpression of CD2AP, but not blocking Aβ-production, is sufficient to rescue spine density. CD2AP overexpression increases spine density, volume, and synapse formation, while a rare LOAD CD2AP mutation induces aberrant F-actin spine-like protrusions without functional synapses. CD2AP controls postsynaptic actin turnover, with the LOAD mutation in CD2AP decreasing F-actin dynamicity. Our data support that CD2AP risk variants could contribute to LOAD synapse dysfunction by disrupting spine formation and growth by deregulating actin dynamics.Significance statement CD2AP is a candidate genetic risk factor of late-onset Alzheimer's disease (LOAD) expressed in neurons with an unknown impact on synapse dysfunction, one of the causal LOAD mechanisms. Our research has revealed CD2AP as a new synaptic protein and established a connection between a LOAD genetic variant in CD2AP and synaptic dysfunction independent of beta-amyloid accumulation. This study suggests an explanation for the CD2AP-mediated predisposition to AD. Furthermore, we have found that controlling CD2AP's impact on spinal F-actin could be a potential target for therapeutic intervention against LOAD.

CRISPR tiling deletion screens reveal functional enhancers of neuropsychiatric risk genes and allelic compensation effects (ACE) on transcription.

Precise transcriptional regulation is critical for cellular function and development, yet the mechanism of this process remains poorly understood for many genes. To gain a deeper understanding of the regulation of neuropsychiatric disease risk genes, we identified a total of 39 functional enhancers for four dosage-sensitive genes, APP, FMR1, MECP2, and SIN3A, using CRISPR tiling deletion screening in human induced pluripotent stem cell (iPSC)-induced excitatory neurons. We found that enhancer annotation provides potential pathological insights into disease-associated copy number variants. More importantly, we discovered that allelic enhancer deletions at SIN3A could be compensated by increased transcriptional activities from the other intact allele. Such allelic compensation effects (ACE) on transcription is stably maintained during differentiation and, once established, cannot be reversed by ectopic SIN3A expression. Further, ACE at SIN3A occurs through dosage sensing by the promoter. Together, our findings unravel a regulatory compensation mechanism that ensures stable and precise transcriptional output for SIN3A, and potentially other dosage-sensitive genes.

Deficiency of SLC26A3 promotes jejunal barrier damage in metabolic disease-susceptible transgenic pigs

Intestinal disorders are common in metabolic syndrome. However, their pathogenesis is still not fully understood. Pig and human intestines are highly similar in terms of associated metabolic processes. Here, we successfully constructed a metabolic disease-susceptible transgenic (TG) Bama pig model by knocking in three humanized disease risk genes with the CRISPR/Cas9 technique to assess its potential as a model for human intestinal diseases and explore the possible pathological mechanisms involved. We found that jejunal barrier integrity was disrupted and that the infiltration of inflammatory cells increased in TG pigs after high-fat and high-sucrose diet (HFHSD) treatment. We revealed significant differences in the transcriptome, associated microbiome profiles and microbial metabolite short-chain fatty acid (SCFA) content of the jejunum of TG pigs. Notably, we found that SLC26A3 was significantly downregulated in TG pigs. Knockdown or overexpression of the SLC26A3 gene in IPEC-J2 cells significantly affected the expression of MUC2, MUC13 and occludin. Furthermore, in vitro experiments further verified that CDX2 directly regulated the expression of SLC26A3. Mechanistically, CDX2 mediated intestinal barrier function by enhancing the expression of SLC26A3 by binding to its promoter region between −1120 and − 1070 bp. TG pigs represent a promising model that provides new insights into preclinical research on human intestinal metabolic diseases associated with metabolic disorders and revealed that SLC26A3 may be a potential therapeutic target for intestinal metabolic diseases.

The Actin and Microtubule Network Regulator WHAMM Is Identified as a Key Kidney Disease Risk Gene

The Actin and Microtubule Network Regulator WHAMM Is Identified as a Key Kidney Disease Risk Gene

A cell and transcriptome atlas of the human arterial vasculature.

Vascular beds show different propensities for different vascular pathologies, yet mechanisms explaining these fundamental differences remain unknown. We sought to build a transcriptomic, cellular, and spatial atlas of human arterial cells across multiple different arterial segments to understand this phenomenon. We found significant cell type-specific segmental heterogeneity. Determinants of arterial identity are predominantly encoded in fibroblasts and smooth muscle cells, and their differentially expressed genes are particularly enriched for vascular disease-associated loci and genes. Adventitial fibroblast-specific heterogeneity in gene expression coincides with numerous vascular disease risk genes, suggesting a previously unrecognized role for this cell type in disease risk. Adult arterial cells from different segments cluster not by anatomical proximity but by embryonic origin, with differentially regulated genes heavily influenced by developmental master regulators. Non-coding transcriptomes across arterial cells contain extensive variation in lnc-RNAs expressed in cell type- and segment-specific patterns, rivaling heterogeneity in protein coding transcriptomes, and show enrichment for non-coding genetic signals for vascular diseases.

Revealing Molecular Patterns of Alzheimer's Disease Risk Gene Expression Signatures in COVID-19 Brains.

Various virus infections are known to predispose to Alzheimer's disease (AD), and a linkage between COVID-19 and AD has been established. COVID-19 infection modulates the gene expression of the genes implicated in progression of AD. Determination of molecular patterns and codon usage and context analysis for the genes that are modulated during COVID-19 infection and are implicated in AD was the target of the study. Our study employed a comprehensive array of research methods, including relative synonymous codon usage, Codon adaptation index analysis, Neutrality and parity analysis, Rare codon analyses, and codon context analysis. This meticulous approach was crucial in determining the molecular patterns present in genes up or downregulated during COVID-19 infection. G/C ending codons were preferred in upregulated genes while not in downregulated genes, and in both gene sets, longer genes have high expressivity. Similarly, T over A nucleotide was preferred, and selection was the major evolutionary force in shaping codon usage in both gene sets. Apart from stops codons, codons CGU - Arg, AUA - Ile, UUA - Leu, UCG - Ser, GUA - Val, and CGA - Arg in upregulated genes, while CUA - Leu, UCG - Ser, and UUA - Leu in downregulated genes were present below the 0.5%. Glutamine-initiated codon pairs have high residual values in upregulated genes. Identical codon pairs GAG-GAG and GUG-GUG were preferred in both gene sets. The shared and unique molecular features in the up- and downregulated gene sets provide insights into the complex interplay between COVID-19 infection and AD. Further studies are required to elucidate the relationship of these molecular patterns with AD pathology.

Cerebral microvascular changes in healthy carriers of the APOE-ɛ4 Alzheimer's disease risk gene.

APOE-ɛ4 is a genetic risk factor for Alzheimer's disease (AD). AD is associated with reduced cerebral blood flow (CBF) and with microvascular changes that limit the transport of oxygen from blood into brain tissue: reduced microvascular cerebral blood volume and high relative transit time heterogeneity (RTH). Healthy APOE-ɛ4 carriers reveal brain regions with elevated CBF compared with carriers of the common ɛ3 allele. Such asymptomatic hyperemia may reflect microvascular dysfunction: a vascular disease entity characterized by suboptimal tissue oxygen uptake, rather than limited blood flow per se. Here, we used perfusion MRI to show that elevated regional CBF is accompanied by reduced capillary blood volume in healthy APOE-ɛ4 carriers (carriers) aged 30-70 years compared with similarly aged APOE-ɛ3 carriers (noncarriers). Younger carriers have elevated hippocampal RTH and more extreme RTH values throughout both white matter (WM) and cortical gray matter (GM) compared with noncarriers. Older carriers have reduced WM CBF and more extreme GM RTH values than noncarriers. Across all groups, lower WM and hippocampal RTH correlate with higher educational attainment, which is associated with lower AD risk. Three days of dietary nitrate supplementation increased carriers' WM CBF but caused older carriers to score worse on two of six aggregate neuropsychological scores. The intervention improved late recall in younger carriers and in noncarriers. The APOE-ɛ4 gene is associated with microvascular changes that may impair tissue oxygen extraction. We speculate that vascular risk factor control is particularly important for APOE-ɛ4 carriers' healthy aging.

Alzheimer's disease risk gene CD2AP is a dose-sensitive determinant of synaptic structure and plasticity.

CD2-Associated protein (CD2AP) is a candidate susceptibility gene for Alzheimer's disease, but its role in the mammalian central nervous system remains largely unknown. We show that CD2AP protein is broadly expressed in the adult mouse brain, including within cortical and hippocampal neurons, where it is detected at pre-synaptic terminals. Deletion of Cd2ap altered dendritic branching and spine density, and impaired ubiquitin-proteasome system activity. Moreover, in mice harboring either one or two copies of a germline Cd2ap null allele, we noted increased paired-pulse facilitation at hippocampal Schaffer-collateral synapses, consistent with a haploinsufficient requirement for pre-synaptic release. Whereas conditional Cd2ap knockout in the brain revealed no gross behavioral deficits in either 3.5- or 12-month-old mice, Cd2ap heterozygous mice demonstrated subtle impairments in discrimination learning using a touchscreen task. Based on unbiased proteomics, partial or complete loss of Cd2ap triggered perturbation of proteins with roles in protein folding, lipid metabolism, proteostasis, and synaptic function. Overall, our results reveal conserved, dose-sensitive requirements for CD2AP in the maintenance of neuronal structure and function, including synaptic homeostasis and plasticity, and inform our understanding of possible cell-type specific mechanisms in Alzheimer's Disease.

Alzheimer's disease genetic risk and changes in brain atrophy and white matter hyperintensities in cognitively unimpaired adults.

Reduced brain volumes and more prominent white matter hyperintensities on MRI scans are commonly observed among older adults without cognitive impairment. However, it remains unclear whether rates of change in these measures among cognitively normal adults differ as a function of genetic risk for late-onset Alzheimer's disease, including APOE-ɛ4, APOE-ɛ2 and Alzheimer's disease polygenic risk scores (AD-PRS), and whether these relationships are influenced by other variables. This longitudinal study examined the trajectories of regional brain volumes and white matter hyperintensities in relationship to APOE genotypes (N = 1541) and AD-PRS (N = 1093) in a harmonized dataset of middle-aged and older individuals with normal cognition at baseline (mean baseline age = 66 years, SD = 9.6) and an average of 5.3 years of MRI follow-up (max = 24 years). Atrophy on volumetric MRI scans was quantified in three ways: (i) a composite score of regions vulnerable to Alzheimer's disease (SPARE-AD); (ii) hippocampal volume; and (iii) a composite score of regions indexing advanced non-Alzheimer's disease-related brain aging (SPARE-BA). Global white matter hyperintensity volumes were derived from fluid attenuated inversion recovery (FLAIR) MRI. Using linear mixed effects models, there was an APOE-ɛ4 gene-dose effect on atrophy in the SPARE-AD composite and hippocampus, with greatest atrophy among ɛ4/ɛ4 carriers, followed by ɛ4 heterozygouts, and lowest among ɛ3 homozygouts and ɛ2/ɛ2 and ɛ2/ɛ3 carriers, who did not differ from one another. The negative associations of APOE-ɛ4 with atrophy were reduced among those with higher education (P < 0.04) and younger baseline ages (P < 0.03). Higher AD-PRS were also associated with greater atrophy in SPARE-AD (P = 0.035) and the hippocampus (P = 0.014), independent of APOE-ɛ4 status. APOE-ɛ2 status (ɛ2/ɛ2 and ɛ2/ɛ3 combined) was not related to baseline levels or atrophy in SPARE-AD, SPARE-BA or the hippocampus, but was related to greater increases in white matter hyperintensities (P = 0.014). Additionally, there was an APOE-ɛ4 × AD-PRS interaction in relation to white matter hyperintensities (P = 0.038), with greater increases in white matter hyperintensities among APOE-ɛ4 carriers with higher AD-PRS. APOE and AD-PRS associations with MRI measures did not differ by sex. These results suggest that APOE-ɛ4 and AD-PRS independently and additively influence longitudinal declines in brain volumes sensitive to Alzheimer's disease and synergistically increase white matter hyperintensity accumulation among cognitively normal individuals. Conversely, APOE-ɛ2 primarily influences white matter hyperintensity accumulation, not brain atrophy. Results are consistent with the view that genetic factors for Alzheimer's disease influence atrophy in a regionally specific manner, likely reflecting preclinical neurodegeneration, and that Alzheimer's disease risk genes contribute to white matter hyperintensity formation.

Effects of Different Combinations of Phytochemical-Rich Fruits and Vegetables on Chronic Disease Risk Markers and Gene Expression Changes: Insights from the MiBLEND Study, a Randomized Trial.

Adequate fruit and vegetable (F and V) intake, as recommended by the World Health Organization (over 400 g/day), is linked to reduced chronic disease risk. However, human intervention trials, especially with whole F and V and in complex combinations, are lacking. The MiBlend Study explored the effects of various phytochemical-rich F and V combinations on chronic disease risk markers, phytochemical absorption, and gene expression in blood. This randomized cross-over study involved participants consuming two of seven different F and V blends for 2 weeks (450 g/day), following a 2-week low F and V intake period (50 g/day). Each blend represented major phytochemical classes (flavonoids, anthocyanins, carotenoids, and glucosinolates) or combinations thereof. Markers of chronic disease risk, including DNA damage, oxidative stress, and retinal microvasculature, were measured. Increasing F and V intake significantly improved plasma antioxidant capacity, DNA damage protection, and retinal arteriolar dilation. Flavonoid-rich, carotenoid-rich, and complex blends notably reduced DNA damage susceptibility. Anthocyanin-rich and carotenoid-rich interventions were most effective in boosting antioxidant capacity, while blends high in flavonoids, especially combined with anthocyanins, significantly improved retinal microvasculature. Gene expression analysis revealed changes in DNA repair, signal transduction, and transcription processes, indicating mechanisms for these health benefits. The study suggests specific F and V blends can provide targeted health improvements, emphasizing the importance of both overall F and V intake and the specific phytochemical composition for personalized preventive strategies.

Dynamic alterations in gene co-expression networks and gene-transcript associations characterize co-morbidities in cocaine use disorder.

Individuals with cocaine use disorder (CUD) who attempt abstinence experience craving and relapse, which poses challenges in treatment. Longitudinal studies linking behavioral manifestations in CUD to the blood transcriptome in living individuals are limited. Therefore, we investigated the connection between drug use behaviors during abstinence with blood transcriptomics. We conducted a comprehensive longitudinal study involving 12 subjects (9 males, 3 females) with CUD and RNA sequencing on blood collected at a drug-free baseline, and 3, 6 & 9 months thereafter. We categorized subjects into 2 responder groups (high-low) based on scores of drug use variables, and 3 responder groups (low-intermediate-high) on days of abstinence. We investigated differential expression and gene-transcript associations across responder groups at each time point. Lastly, we examined genes that are both co-expressed and showed dynamic expression with time. Genes with significant transcript associations between high and. intermediate days of abstinence at 9 months were notably enriched for cannabis use disorder, drinks weekly, and coronary artery disease risk genes. Time-specific gene co-expression analysis prioritized transcripts related to immune processes, cell cycle, RNA-protein synthesis, and second messenger signaling for days of abstinence. We demonstrate that abstinence reflects robust changes in drug use behaviors and the blood transcriptome in CUD. We also highlight the importance of longitudinal studies to capture complex biological processes during abstinence in CUD.

Enhancing Disease Risk Gene Discovery by Integrating Transcription Factor-Linked Trans-located Variants into Transcriptome-Wide Association Analyses.

Transcriptome-wide association studies (TWAS) have been successful in identifying disease susceptibility genes by integrating cis-variants predicted gene expression with genome-wide association studies (GWAS) data. However, trans-located variants for predicting gene expression remain largely unexplored. Here, we introduce transTF-TWAS, which incorporates transcription factor (TF)-linked trans-located variants to enhance model building. Using data from the Genotype-Tissue Expression project, we predict gene expression and alternative splicing and applied these models to large GWAS datasets for breast, prostate, and lung cancers. We demonstrate that transTF-TWAS outperforms other existing TWAS approaches in both constructing gene prediction models and identifying disease-associated genes, as evidenced by simulations and real data analysis. Our transTF-TWAS approach significantly contributes to the discovery of disease risk genes. Findings from this study have shed new light on several genetically driven key regulators and their associated regulatory networks underlying disease susceptibility.

A systems biology-based identification and invivo functional screening of Alzheimer's disease risk genes reveal modulators of memory function.

Genome-wide association studies (GWASs) have uncovered over 75 genomic loci associated with risk for late-onset Alzheimer's disease (LOAD), but identification of the underlying causal genes remains challenging. Studies of induced pluripotent stem cell (iPSC)-derived neurons from LOAD patients have demonstrated the existence of neuronal cell-intrinsic functional defects. Here, we searched for genetic contributions to neuronal dysfunction in LOAD using an integrative systems approach that incorporated multi-evidence-based gene mapping and network-analysis-based prioritization. A systematic perturbation screening of candidate risk genes in Caenorhabditis elegans (C.elegans) revealed that neuronal knockdown of the LOAD risk gene orthologs vha-10 (ATP6V1G2), cmd-1 (CALM3), amph-1 (BIN1), ephx-1 (NGEF), and pho-5 (ACP2) alters short-/intermediate-term memory function, the cognitive domain affected earliest during LOAD progression. These results highlight the impact of LOAD risk genes on evolutionarily conserved memory function, as mediated through neuronal endosomal dysfunction, and identify new targets for further mechanistic interrogation.

Prediction of motor and non-motor Parkinson’s disease symptoms using serum lipidomics and machine learning: a 2-year study

Identifying biological factors which contribute to the clinical progression of heterogeneous motor and non-motor phenotypes in Parkinson’s disease may help to better understand the disease process. Several lipid-related genetic risk factors for Parkinson’s disease have been identified, and the serum lipid signature of Parkinson’s disease patients is significantly distinguishable from controls. However, the extent to which lipid profiles are associated with clinical outcomes remains unclear. Untargeted high-performance liquid chromatography-tandem mass spectrometry identified >900 serum lipids in Parkinson’s disease subjects at baseline (n = 122), and the potential for machine learning models using these lipids to predict motor and non-motor clinical scores after 2 years (n = 67) was assessed. Machine learning models performed best when baseline serum lipids were used to predict the 2-year future Unified Parkinson’s disease rating scale part three (UPDRS III) and Geriatric Depression Scale scores (both normalised root mean square error = 0.7). Feature analysis of machine learning models indicated that species of lysophosphatidylethanolamine, phosphatidylcholine, platelet-activating factor, sphingomyelin, diacylglycerol and triacylglycerol were top predictors of both motor and non-motor scores. Serum lipids were overall more important predictors of clinical outcomes than subject sex, age and mutation status of the Parkinson’s disease risk gene LRRK2. Furthermore, lipids were found to better predict clinical scales than a panel of 27 serum cytokines previously measured in this cohort (The Michael J. Fox Foundation LRRK2 Clinical Cohort Consortium). These results suggest that lipid changes may be associated with clinical phenotypes in Parkinson’s disease.

Altered glia-neuron communication in Alzheimer’s Disease affects WNT, p53, and NFkB Signaling determined by snRNA-seq

BackgroundAlzheimer’s disease is the most common cause of dementia and is characterized by amyloid-β plaques, tau neurofibrillary tangles, and neuronal loss. Although neuronal loss is a primary hallmark of Alzheimer’s disease, it is known that non-neuronal cell populations are ultimately responsible for maintaining brain homeostasis and neuronal health through neuron-glia and glial cell crosstalk. Many signaling pathways have been proposed to be dysregulated in Alzheimer’s disease, including WNT, TGFβ, p53, mTOR, NFkB, and Pi3k/Akt signaling. Here, we predict altered cell-cell communication between glia and neurons.MethodsUsing public snRNA-sequencing data generated from postmortem human prefrontal cortex, we predicted altered cell-cell communication between glia (astrocytes, microglia, oligodendrocytes, and oligodendrocyte progenitor cells) and neurons (excitatory and inhibitory). We confirmed interactions in a second and third independent orthogonal dataset. We determined cell-type-specificity using Jaccard Similarity Index and investigated the downstream effects of altered interactions in inhibitory neurons through gene expression and transcription factor activity analyses of signaling mediators. Finally, we determined changes in pathway activity in inhibitory neurons.ResultsCell-cell communication between glia and neurons is altered in Alzheimer’s disease in a cell-type-specific manner. As expected, ligands are more cell-type-specific than receptors and targets. We identified ligand-receptor pairs in three independent datasets and found involvement of the Alzheimer’s disease risk genes APP and APOE across datasets. Most of the signaling mediators of these interactions were not significantly differentially expressed, however, the mediators that are also transcription factors had differential activity between AD and control. Namely, MYC and TP53, which are associated with WNT and p53 signaling, respectively, had decreased TF activity in Alzheimer’s disease, along with decreased WNT and p53 pathway activity in inhibitory neurons. Additionally, inhibitory neurons had both increased NFkB signaling pathway activity and increased TF activity of NFIL3, an NFkB signaling-associated transcription factor.ConclusionsCell-cell communication between glia and neurons in Alzheimer’s disease is altered in a cell-type-specific manner involving Alzheimer’s disease risk genes. Signaling mediators had altered transcription factor activity suggesting altered glia-neuron interactions may dysregulate signaling pathways including WNT, p53, and NFkB in inhibitory neurons.

Single-cell genomics and regulatory networks for 388 human brains.

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multiomics datasets into a resource comprising >2.8 million nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550,000 cell type-specific regulatory elements and >1.4 million single-cell expression quantitative trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Xenografted human microglia display diverse transcriptomic states in response to Alzheimer’s disease-related amyloid-β pathology

Microglia are central players in Alzheimer’s disease pathology but analyzing microglial states in human brain samples is challenging due to genetic diversity, postmortem delay and admixture of pathologies. To circumvent these issues, here we generated 138,577 single-cell expression profiles of human stem cell-derived microglia xenotransplanted in the brain of the AppNL-G-F model of amyloid pathology and wild-type controls. Xenografted human microglia adopt a disease-associated profile similar to that seen in mouse microglia, but display a more pronounced human leukocyte antigen or HLA state, likely related to antigen presentation in response to amyloid plaques. The human microglial response also involves a pro-inflammatory cytokine/chemokine cytokine response microglia or CRM response to oligomeric Aβ oligomers. Genetic deletion of TREM2 or APOE as well as APOE polymorphisms and TREM2R47H expression in the transplanted microglia modulate these responses differentially. The expression of other Alzheimer’s disease risk genes is differentially regulated across the distinct cell states elicited in response to amyloid pathology. Thus, we have identified multiple transcriptomic cell states adopted by human microglia in a multipronged response to Alzheimer’s disease-related pathology, which should be taken into account in translational studies.