Genetic mapping of etiologic brain cell types for obesity

Scalable Hypothalamic Arcuate Neuron Differentiation from Human Pluripotent Stem Cells Suitable for Modeling Metabolic and Reproductive Disorders

The hypothalamus, composed of several nuclei, is essential for maintaining our body's homeostasis. The arcuate nucleus (ARC), located in the mediobasal hypothalamus, contains neuronal populations with eminent roles in energy and glucose homeostasis as well as reproduction. These neuronal populations are of great interest for translational research. To fulfill this promise, we used a robotic cell culture platform to provide a scalable and chemically defined approach for differentiating human pluripotent stem cells (hPSCs) into pro-opiomelanocortin (POMC), somatostatin (SST), tyrosine hydroxylase (TH) and gonadotropin-releasing hormone (GnRH) neuronal subpopulations with an ARC-like signature. This robust approach is reproducible across several distinct hPSC lines and exhibits a stepwise induction of key ventral diencephalon and ARC markers in transcriptomic profiling experiments. This is further corroborated by direct comparison to human fetal hypothalamus, and the enriched expression of genes implicated in obesity and type 2 diabetes (T2D). Genome-wide chromatin accessibility profiling by ATAC-seq identified accessible regulatory regions that can be utilized to predict candidate enhancers related to metabolic disorders and hypothalamic development. In depth molecular, cellular, and functional experiments unveiled the responsiveness of the hPSC-derived hypothalamic neurons to hormonal stimuli, such as insulin, neuropeptides including kisspeptin, and incretin mimetic drugs such as Exendin-4, highlighting their potential utility as physiologically relevant cellular models for disease studies. In addition, differential glucose and insulin treatments uncovered adaptability within the generated ARC neurons in the dynamic regulation of POMC and insulin receptors. In summary, the establishment of this model represents a novel, chemically defined, and scalable platform for manufacturing large numbers of hypothalamic arcuate neurons and serves as a valuable resource for modeling metabolic and reproductive disorders.

Mendelian randomization analysis of the brain, cerebrospinal fluid, and plasma proteome identifies potential drug targets for attention deficit hyperactivity disorder

BACKGROUND

The need for new therapeutics for attention deficit hyperactivity disorder (ADHD) is evident. Brain, cerebrospinal fluid (CSF), and plasma protein biomarkers with causal genetic evidence could represent potential drug targets. However, a comprehensive screen of the proteome has not yet been conducted.

METHODS

We employed a three-pronged approach using Mendelian Randomization (MR) and Bayesian colocalization analysis. Firstly, we studied 608 brains, 214 CSF, and 612 plasma proteins as potential causal mediators of ADHD using MR analysis. Secondly, we analysed the consistency of the discovered biomarkers across three distinct subtypes of ADHD: childhood, persistent, and late-diagnosed ADHD. Finally, we extended our analysis to examine the correlation between identified biomarkers and Tourette syndrome and pervasive autism spectrum disorder (ASD), conditions often linked with ADHD. To validate the MR findings, we conducted sensitivity analysis. Additionally, we performed cell type analysis on the human brain to identify risk genes that are notably enriched in various brain cell types.

FINDINGS

After applying Bonferroni correction, we found that the risk of ADHD was increased by brain proteins GMPPB, NAA80, HYI, CISD2, and HYI, TIE1 in CSF and plasma. Proteins GMPPB, NAA80, ICA1L, CISD2, TIE1, and RMDN1 showed overlapped loci with ADHD risk through Bayesian colocalization. Overexpression of GMPPB protein was linked to an increase in the risk for all three ADHD subtypes. While ICA1L provided protection against both ASD and ADHD, CISD2 increased the probability of both disorders. Cell-specific studies revealed that GMPPB, NAA80, ICA1L, and CISD2 were predominantly present on the surface of excitatory-inhibitory neurons.

INTERPRETATION

Our comprehensive MR investigation of the brain, CSF, and plasma proteomes revealed seven proteins with causal connections to ADHD. Particularly, GMPPB and TIE1 emerged as intriguing targets for potential ADHD therapy.

FUNDING

This work was partly funded by the Key R & D Program of Zhejiang (T.L. 2022C03096); the National Natural Science Foundation of China Project (C.Z. 82001413); Postdoctoral Foundation of West China Hospital (C.Z. 2020HXBH163).

Obesity and the cerebral cortex: Underlying neurobiology in mice and humans

Obesity is a major modifiable risk factor for Alzheimer's disease (AD), characterized by progressive atrophy of the cerebral cortex. The neurobiology of obesity contributions to AD is poorly understood. Here we show with in vivo MRI that diet-induced obesity decreases cortical volume in mice, and that higher body adiposity associates with lower cortical volume in humans. Single-nuclei transcriptomics of the mouse cortex reveals that dietary obesity promotes an array of neuron-adverse transcriptional dysregulations, which are mediated by an interplay of excitatory neurons and glial cells, and which involve microglial activation and lowered neuronal capacity for neuritogenesis and maintenance of membrane potential. The transcriptional dysregulations of microglia, more than of other cell types, are like those in AD, as assessed with single-nuclei cortical transcriptomics in a mouse model of AD and two sets of human donors with the disease. Serial two-photon tomography of microglia demonstrates microgliosis throughout the mouse cortex. The spatial pattern of adiposity-cortical volume associations in human cohorts interrogated together with in silico bulk and single-nucleus transcriptomic data from the human cortex implicated microglia (along with other glial cells and subtypes of excitatory neurons), and it correlated positively with the spatial profile of cortical atrophy in patients with mild cognitive impairment and AD. Thus, multi-cell neuron-adverse dysregulations likely contribute to the loss of cortical tissue in obesity. The dysregulations of microglia may be pivotal to the obesity-related risk of AD.

Identifying novel proteins for suicide attempt by integrating proteomes from brain and blood with genome-wide association data

Genome-wide association studies (GWASs) have identified risk loci for suicide attempt (SA), but deciphering how they confer risk for SA remains largely unknown. This study aims to identify the key proteins and gain insights into SA pathogenesis. We integrated data from the brain proteome (N = 376) and blood proteome (N = 35,559) and combined it with the largest SA GWAS summary statistics to date (N = 518,612). A comprehensive set of methods was employed, including Mendelian randomization (MR), Steiger filtering, Bayesian colocalization, proteome‑wide association studies (PWAS), transcript-levels, cell-type specificity, correlation, and protein-protein interaction (PPI) network analysis. Validation was performed using other protein datasets and the SA dataset from FinnGen study. We identified ten proteins (GLRX5, GMPPB, B3GALTL, FUCA2, TTLL12, ADCK1, MMAA, HIBADH, ACP1, DOC2A) associated with SA in brain proteomics. GLRX5, GMPPB, and FUCA2 showed strong colocalization evidence and were supported by PWAS and transcript-level analysis, and were predominantly expressed in glutamatergic neuronal cells. In blood proteomics, one significant protein (PEAR1) and three near-significant proteins (NDE1, EVA1C, B4GALT2) were identified, but lacked colocalization evidence. Moreover, despite the limited correlation between the same protein in brain and blood, the PPI network analysis provided new insights into the interaction between brain and blood in SA. Furthermore, GLRX5 was associated with the GSTP1, the target of Clozapine. The comprehensive analysis provides strong evidence supporting a causal association between three genetically determined brain proteins (GLRX5, GMPPB, and FUCA2) with SA. These findings offer valuable insights into SA's underlying mechanisms and potential therapeutic approaches.

Unraveling NEK4 as a Potential Drug Target in Schizophrenia and Bipolar I Disorder: A Proteomic and Genomic Approach

BACKGROUND AND HYPOTHESIS

Investigating the shared brain protein and genetic components of schizophrenia (SCZ) and bipolar I disorder (BD-I) presents a unique opportunity to understand the underlying pathophysiological processes and pinpoint potential drug targets.

STUDY DESIGN

To identify overlapping susceptibility brain proteins in SCZ and BD-I, we carried out proteome-wide association studies (PWAS) and Mendelian Randomization (MR) by integrating human brain protein quantitative trait loci with large-scale genome-wide association studies for both disorders. We utilized transcriptome-wide association studies (TWAS) to determine the consistency of mRNA-protein dysregulation in both disorders. We applied pleiotropy-informed conditional false discovery rate (pleioFDR) analysis to identify common risk genetic loci for SCZ and BD-I. Additionally, we performed a cell-type-specific analysis in the human brain to detect risk genes notably enriched in distinct brain cell types. The impact of risk gene overexpression on dendritic arborization and axon length in neurons was also examined.

STUDY RESULTS

Our PWAS identified 42 proteins associated with SCZ and 14 with BD-I, among which NEK4, HARS2, SUGP1, and DUS2 were common to both conditions. TWAS and MR analysis verified the significant risk gene NEK4 for both SCZ and BD-I. PleioFDR analysis further supported genetic risk loci associated with NEK4 for both conditions. The cell-type specificity analysis revealed that NEK4 is expressed on the surface of glutamatergic neurons, and its overexpression enhances dendritic arborization and axon length in cultured primary neurons.

CONCLUSIONS

These findings underscore a shared genetic origin for SCZ and BD-I, offering novel insights for potential therapeutic target identification.

Revealing brain cell-stratified causality through dissecting causal variants according to their cell-type-specific effects on gene expression

The human brain has been implicated in the pathogenesis of several complex diseases. Taking advantage of single-cell techniques, genome-wide association studies (GWAS) have taken it a step further and revealed brain cell-type-specific functions for disease loci. However, genetic causal associations inferred by Mendelian randomization (MR) studies usually include all instrumental variables from GWAS, which hampers the understanding of cell-specific causality. Here, we developed an analytical framework, Cell-Stratified MR (csMR), to investigate cell-stratified causality through colocalizing GWAS signals with single-cell eQTL from different brain cells. By applying to obesity-related traits, our results demonstrate the cell-type-specific effects of GWAS variants on gene expression, and indicate the benefits of csMR to identify cell-type-specific causal effect that is often hidden from bulk analyses. We also found csMR valuable to reveal distinct causal pathways between different obesity indicators. These findings suggest the value of our approach to prioritize target cells for extending genetic causation studies.

Prioritization of Kidney Cell Types Highlights Myofibroblast Cells in Regulating Human Blood Pressure

Introduction

Blood pressure (BP) is a highly heritable trait with over 2000 underlying genomic loci identified to date. Although the kidney plays a key role, little is known about specific cell types involved in the genetic regulation of BP.

Methods

Here, we applied stratified linkage disequilibrium score (LDSC) regression to connect BP genome-wide association studies (GWAS) results to specific cell types of the mature human kidney. We used the largest single-stage BP genome-wide analysis to date, including up to 1,028,980 adults of European ancestry, and single-cell transcriptomic data from 14 mature human kidneys, with mean age of 41 years.

Results

Our analyses prioritized myofibroblasts and endothelial cells, among the total of 33 annotated cell type, as specifically involved in BP regulation (P < 0.05/33, i.e., 0.001515). Enrichment of heritability for systolic BP (SBP) was observed in myofibroblast cells in mature human kidney cortex, and enrichment of heritability for diastolic BP (DBP) was observed in descending vasa recta and peritubular capillary endothelial cells as well as stromal myofibroblast cells. The new finding of myofibroblast, the significant cell type for both BP traits, was consistent in 8 replication efforts using 7 sets of independent data, including in human fetal kidney, in East-Asian (EAS) ancestry, using mouse single-cell RNA sequencing (scRNA-seq) data, and when using another prioritization method.

Conclusion

Our findings provide a solid basis for follow-up studies to further identify genes and mechanisms in myofibroblast cells that underlie the regulation of BP.

Genome-wide meta-analyses of restless legs syndrome yield insights into genetic architecture, disease biology and risk prediction

Restless legs syndrome (RLS) affects up to 10% of older adults. Their healthcare is impeded by delayed diagnosis and insufficient treatment. To advance disease prediction and find new entry points for therapy, we performed meta-analyses of genome-wide association studies in 116,647 individuals with RLS (cases) and 1,546,466 controls of European ancestry. The pooled analysis increased the number of risk loci eightfold to 164, including three on chromosome X. Sex-specific meta-analyses revealed largely overlapping genetic predispositions of the sexes (rg = 0.96). Locus annotation prioritized druggable genes such as glutamate receptors 1 and 4, and Mendelian randomization indicated RLS as a causal risk factor for diabetes. Machine learning approaches combining genetic and nongenetic information performed best in risk prediction (area under the curve (AUC) = 0.82-0.91). In summary, we identified targets for drug development and repurposing, prioritized potential causal relationships between RLS and relevant comorbidities and risk factors for follow-up and provided evidence that nonlinear interactions are likely relevant to RLS risk prediction.

Obesity affects brain cortex gene expression in an APOE genotype and sex dependent manner

OBJECTIVE

Obesity is the top modifiable risk factor for Alzheimer's disease. We hypothesized that high fat diet (HFD)-induced obesity alters brain transcriptomics in APOE-genotype and sex dependent manners. Here, we investigated interactions between HFD, APOE, and sex, using a knock-in mouse model of the human APOE3 and APOE4 alleles.

METHODS

Six-month-old APOE3-TR and APOE4-TR mice were treated with either HFD or control chow. After 4 months, total RNA was extracted from the cerebral cortices and analyzed by poly-A enriched RNA sequencing on the Illumina platform.

RESULTS

Female mice demonstrated profound HFD-induced transcriptomic changes while there was little to no effect in males. In females, APOE3 brains demonstrated about five times more HFD-induced transcriptomic changes (399 up-regulated and 107 down-regulated genes) compared to APOE4 brains (30 up-regulated and 60 down-regulated). Unsupervised clustering analysis revealed two gene sets that responded to HFD in APOE3 mice but not in APOE4 mice. Pathway analysis demonstrated that HFD in APOE3 mice affected cortical pathways related to feeding behavior, blood circulation, circadian rhythms, extracellular matrix, and cell adhesion.

CONCLUSIONS

Female mice and APOE3 mice have the strongest cortical transcriptomic responses to HFD related to feeding behavior and extracellular matrix remodeling. The relative lack of response of the APOE4 brain to stress associated with obesity may leave it more susceptible to additional stresses that occur with aging and in AD.

Shared genetics linking sociability with the brain's default mode network

The brain's default mode network (DMN) plays a role in social cognition, with altered DMN function being associated with social impairments across various neuropsychiatric disorders. In the present study, we examined the genetic relationship between sociability and DMN-related resting-state functional magnetic resonance imaging (rs-fMRI) traits. To this end, we used genome-wide association summary statistics for sociability and 31 activity and 64 connectivity DMN-related rs-fMRI traits (N=34,691-342,461). First, we examined global and local genetic correlations between sociability and the rs-fMRI traits. Second, to assess putatively causal relationships between the traits, we conducted bi-directional Mendelian randomisation (MR) analyses. Finally, we prioritised genes influencing both sociability and rs-fMRI traits by combining three methods: gene-expression eQTL MR analyses, the CELLECT framework using single-nucleus RNA-seq data, and network propagation in the context of a protein-protein interaction network. Significant local genetic correlations were found between sociability and two rs-fMRI traits, one representing spontaneous activity within the temporal cortex, the other representing connectivity between the frontal/cingulate and angular/temporal cortices. Sociability affected 12 rs-fMRI traits when allowing for weakly correlated genetic instruments. Combing all three methods for gene prioritisation, we defined 17 highly prioritised genes, with DRD2 and LINGO1 showing the most robust evidence across all analyses. By integrating genetic and transcriptomics data, our gene prioritisation strategy may serve as a blueprint for future studies. The prioritised genes could be explored as potential biomarkers for social dysfunction in the context of neuropsychiatric disorders and as drug target genes.

Cell type-specific network analysis in Diversity Outbred mice identifies genes potentially responsible for human bone mineral density GWAS associations

Genome-wide association studies (GWASs) have identified many sources of genetic variation associated with bone mineral density (BMD), a clinical predictor of fracture risk and osteoporosis. Aside from the identification of causal genes, other difficult challenges to informing GWAS include characterizing the roles of predicted causal genes in disease and providing additional functional context, such as the cell type predictions or biological pathways in which causal genes operate. Leveraging single-cell transcriptomics (scRNA-seq) can assist in informing BMD GWAS by linking disease-associated variants to genes and providing a cell type context for which these causal genes drive disease. Here, we use large-scale scRNA-seq data from bone marrow-derived stromal cells cultured under osteogenic conditions (BMSC-OBs) from Diversity Outbred (DO) mice to generate cell type-specific networks and contextualize BMD GWAS-implicated genes. Using trajectories inferred from the scRNA-seq data, we identify networks enriched with genes that exhibit the most dynamic changes in expression across trajectories. We discover 21 network driver genes, which are likely to be causal for human BMD GWAS associations that colocalize with expression/splicing quantitative trait loci (eQTL/sQTL). These driver genes, including Fgfrl1 and Tpx2, along with their associated networks, are predicted to be novel regulators of BMD via their roles in the differentiation of mesenchymal lineage cells. In this work, we showcase the use of single-cell transcriptomics from mouse bone-relevant cells to inform human BMD GWAS and prioritize genetic targets with potential causal roles in the development of osteoporosis.

Integrated Single-Cell Transcriptomic Atlas of Human Kidney Endothelial Cells

Key Points

We created a comprehensive reference atlas of normal human kidney endothelial cells. We confirmed that endothelial cell types in the human kidney were also highly conserved in the mouse kidney.

Background

Kidney endothelial cells are exposed to different microenvironmental conditions that support specific physiologic processes. However, the heterogeneity of human kidney endothelial cells has not yet been systematically described.

Methods

We reprocessed and integrated seven human kidney control single-cell/single-nucleus RNA sequencing datasets of >200,000 kidney cells in the same process.

Results

We identified five major cell types, 29,992 of which were endothelial cells. Endothelial cell reclustering identified seven subgroups that differed in molecular characteristics and physiologic functions. Mapping new data to a normal kidney endothelial cell atlas allows rapid data annotation and analysis. We confirmed that endothelial cell types in the human kidney were also highly conserved in the mouse kidney and identified endothelial marker genes that were conserved in humans and mice, as well as differentially expressed genes between corresponding subpopulations. Furthermore, combined analysis of single-cell transcriptome data with public genome-wide association study data showed a significant enrichment of endothelial cells, especially arterial endothelial cells, in BP heritability. Finally, we identified M1 and M12 from coexpression networks in endothelial cells that may be deeply involved in BP regulation.

Conclusions

We created a comprehensive reference atlas of normal human kidney endothelial cells that provides the molecular foundation for understanding how the identity and function of kidney endothelial cells are altered in disease, aging, and between species. Finally, we provide a publicly accessible online tool to explore the datasets described in this work (https://vascularmap.jiangnan.edu.cn ).

Hypothalamic GABAergic Neurons Expressing Cellular Retinoic Acid Binding Protein 1 (CRABP1) Are Sensitive to Metabolic Status and Liraglutide in Male Mice

INTRODUCTION

Owing to their privileged anatomical location, neurons of the arcuate nucleus of the hypothalamus (ARC) play critical roles in sensing and responding to metabolic signals such as leptin and glucagon-like peptide 1 (GLP-1). In addition to the well-known proopiomelanocortin (POMC)- and agouti-related peptide (AgRP)-expressing neurons, subpopulations of GABAergic neurons are emerging as key regulators of energy balance. However, the precise identity of these metabolic neurons is still elusive. Here, we identified and characterized the molecular signature of a novel population of GABAergic neurons of the ARC expressing Cellular retinoic acid binding protein 1 (Crabp1).

METHODS

Using a combination of immunohistochemistry and in situ hybridization techniques, we investigated the expression of Crabp1 across the mouse brain and characterized the molecular identity of Crabp1ARC neurons. We also determined whether Crabp1ARC neurons are sensitive to fasting, leptin, and GLP1R agonism by assessing cFOS immunoreactivity as a marker of neuronal activity.

RESULTS

Crabp1ARC neurons represent a novel GABAergic neuronal population robustly enriched in the ARC and are distinct from the prototypical melanocortin neurons. Crabp1ARC neurons overlap with three subpopulations of yet uncharacterized ARC neurons expressing Htr3b, Tbx19, and Tmem215. Notably, Crabp1ARC neurons express receptors for metabolic hormones and their activity is modulated by the nutritional state and GLP1R agonism.

CONCLUSION

Crabp1ARC neurons represent a novel heterogeneous population of GABAergic neurons sensitive to metabolic status.

GWAS-informed data integration and non-coding CRISPRi screen illuminate genetic etiology of bone mineral density

Over 1,100 independent signals have been identified with genome-wide association studies (GWAS) for bone mineral density (BMD), a key risk factor for mortality-increasing fragility fractures; however, the effector gene(s) for most remain unknown. Informed by a variant-to-gene mapping strategy implicating 89 non-coding elements predicted to regulate osteoblast gene expression at BMD GWAS loci, we executed a single-cell CRISPRi screen in human fetal osteoblast 1.19 cells (hFOBs). The BMD relevance of hFOBs was supported by heritability enrichment from cross-cell type stratified LD-score regression involving 98 cell types grouped into 15 tissues. 24 genes showed perturbation in the screen, with four (ARID5B, CC2D1B, EIF4G2, and NCOA3) exhibiting consistent effects upon siRNA knockdown on three measures of osteoblast maturation and mineralization. Lastly, additional heritability enrichments, genetic correlations, and multi-trait fine-mapping revealed that many BMD GWAS signals are pleiotropic and likely mediate their effects via non-bone tissues that warrant attention in future screens.

Light controls mesophyll-specific post-transcriptional splicing of photoregulatory genes by AtPRMT5

Light plays a central role in plant growth and development, providing an energy source and governing various aspects of plant morphology. Previous study showed that many polyadenylated full-length RNA molecules within the nucleus contain unspliced introns (post-transcriptionally spliced introns, PTS introns), which may play a role in rapidly responding to changes in environmental signals. However, the mechanism underlying post-transcriptional regulation during initial light exposure of young, etiolated seedlings remains elusive. In this study, we used FLEP-seq2, a Nanopore-based sequencing technique, to analyze nuclear RNAs in Arabidopsis (Arabidopsis thaliana) seedlings under different light conditions and found numerous light-responsive PTS introns. We also used single-nucleus RNA sequencing (snRNA-seq) to profile transcripts in single nucleus and investigate the distribution of light-responsive PTS introns across distinct cell types. We established that light-induced PTS introns are predominant in mesophyll cells during seedling de-etiolation following exposure of etiolated seedlings to light. We further demonstrated the involvement of the splicing-related factor A. thaliana PROTEIN ARGININE METHYLTRANSFERASE 5 (AtPRMT5), working in concert with the E3 ubiquitin ligase CONSTITUTIVE PHOTOMORPHOGENIC 1 (COP1), a critical repressor of light signaling pathways. We showed that these two proteins orchestrate light-induced PTS events in mesophyll cells and facilitate chloroplast development, photosynthesis, and morphogenesis in response to ever-changing light conditions. These findings provide crucial insights into the intricate mechanisms underlying plant acclimation to light at the cell-type level.

3D genomic features across >50 diverse cell types reveal insights into the genomic architecture of childhood obesity

The prevalence of childhood obesity is increasing worldwide, along with the associated common comorbidities of type 2 diabetes and cardiovascular disease in later life. Motivated by evidence for a strong genetic component, our prior genome-wide association study (GWAS) efforts for childhood obesity revealed 19 independent signals for the trait; however, the mechanism of action of these loci remains to be elucidated. To molecularly characterize these childhood obesity loci we sought to determine the underlying causal variants and the corresponding effector genes within diverse cellular contexts. Integrating childhood obesity GWAS summary statistics with our existing 3D genomic datasets for 57 human cell types, consisting of high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq, and RNA-seq, we applied stratified LD score regression and calculated the proportion of genome-wide SNP heritability attributable to cell type-specific features, revealing pancreatic alpha cell enrichment as the most statistically significant. Subsequent chromatin contact-based fine-mapping was carried out for genome-wide significant childhood obesity loci and their linkage disequilibrium proxies to implicate effector genes, yielded the most abundant number of candidate variants and target genes at the BDNF, ADCY3, TMEM18 and FTO loci in skeletal muscle myotubes and the pancreatic beta-cell line, EndoC-BH1. One novel implicated effector gene, ALKAL2 - an inflammation-responsive gene in nerve nociceptors - was observed at the key TMEM18 locus across multiple immune cell types. Interestingly, this observation was also supported through colocalization analysis using expression quantitative trait loci (eQTL) derived from the Genotype-Tissue Expression (GTEx) dataset, supporting an inflammatory and neurologic component to the pathogenesis of childhood obesity. Our comprehensive appraisal of 3D genomic datasets generated in a myriad of different cell types provides genomic insights into pediatric obesity pathogenesis.

Identification of novel proteins for sleep apnea by integrating genome-wide association data and human brain proteomes

BACKGROUND

Sleep apnea is regarded as a significant global public health issue. The relationship between sleep apnea and nervous system diseases is intricate, yet the precise mechanism remains unclear.

METHODS

In this study, we conducted a comprehensive analysis integrating the human brain proteome and transcriptome with sleep apnea genome-wide association study (GWAS), employing genome-wide association study (PWAS), transcriptome-wide association study (TWAS), Mendelian randomization (MR), and colocalization analysis to identify brain proteins associated with sleep apnea.

RESULTS

The discovery PWAS identified six genes (CNNM2, XRCC6, C3orf18, CSDC2, SQRDL, and DGUOK) whose altered protein abundances in the brain were found to be associated with sleep apnea. The independent confirmatory PWAS successfully replicated four out of these six genes (CNNM2, C3orf18, CSDC2, and SQRDL). The transcriptome level TWAS analysis further confirmed two out of the four genes (C3orf18 and CSDC2). The subsequent two-sample Mendelian randomization provided compelling causal evidence supporting the association of C3orf18, CSDC2, CNNM2, and SQRDL with sleep apnea. The co-localization analysis further supported the association between CSDC2 and sleep apnea (posterior probability of hypothesis 4 = 0.75).

CONCLUSIONS

In summary, the integration of brain proteomic and transcriptomic data provided multifaceted evidence supporting causal relationships between four specific brain proteins (CSDC2, C3orf18, CNNM2, and SQRDL) and sleep apnea. Our findings provide new insights into the molecular basis of sleep apnea in the brain, promising to advance understanding of its pathogenesis in future research.

Connecting genomic results for psychiatric disorders to human brain cell types and regions reveals convergence with functional connectivity

Understanding the temporal and spatial brain locations etiological for psychiatric disorders is essential for targeted neurobiological research. Integration of genomic insights from genome-wide association studies with single-cell transcriptomics is a powerful approach although past efforts have necessarily relied on mouse atlases. Leveraging a comprehensive atlas of the adult human brain, we prioritized cell types via the enrichment of SNP-heritabilities for brain diseases, disorders, and traits, progressing from individual cell types to brain regions. Our findings highlight specific neuronal clusters significantly enriched for the SNP-heritabilities for schizophrenia, bipolar disorder, and major depressive disorder along with intelligence, education, and neuroticism. Extrapolation of cell-type results to brain regions reveals important patterns for schizophrenia with distinct subregions in the hippocampus and amygdala exhibiting the highest significance. Cerebral cortical regions display similar enrichments despite the known prefrontal dysfunction in those with schizophrenia highlighting the importance of subcortical connectivity. Using functional MRI connectivity from cases with schizophrenia and neurotypical controls, we identified brain networks that distinguished cases from controls that also confirmed involvement of the central and lateral amygdala, hippocampal body, and prefrontal cortex. Our findings underscore the value of single-cell transcriptomics in decoding the polygenicity of psychiatric disorders and offer a promising convergence of genomic, transcriptomic, and brain imaging modalities toward common biological targets.

An RNA-seq atlas of mouse brain areas during fasting and diet-induced obesity

Mammalian energy homeostasis is primarilly regulated by the hypothalamus and hindbrain, with the hippocampus, midbrain nuclei, and other regions implicated by evidence from human genetics studies. To understand how these non-canonical brain regions respond to imbalances in energy homeostasis, we performed two experiments examining the effects of different diets in male C57BL6 mice. In our first study, groups of six pair-housed mice were given access to chow, high-fat diet or fasted for 16 hours. In our subsequent study, two groups of 10 mice were single-housed and given access to chow or fasted for 24 h. We recorded food intake for each cage, the change in body weight for each animal, and collected hypothalamus, hippocampus, superior colliculus, inferior colliculus, frontal cortex, and zona incerta-centric samples. We performed bulk RNA sequencing on 185 samples and validated them by a series of quality control assessments including alignment quality and gene expression profiling. We believe these studies capture the transcriptomic effects of acute fasting and high-fat diet in the rodent brain and provide a valuable reference.

Meta-analysis of single-cell RNA-sequencing data for depicting the transcriptomic landscape of chronic obstructive pulmonary disease

Chronic obstructive pulmonary disease (COPD) is a respiratory disease characterized by airflow limitation and chronic inflammation of the lungs that is a leading cause of death worldwide. Since the complete pathological mechanisms at the single-cell level are not fully understood yet, an integrative approach to characterizing the single-cell-resolution landscape of COPD is required. To identify the cell types and mechanisms associated with the development of COPD, we conducted a meta-analysis using three single-cell RNA-sequencing datasets of COPD. Among the 154,011 cells from 16 COPD patients and 18 healthy subjects, 17 distinct cell types were observed. Of the 17 cell types, monocytes, mast cells, and alveolar type 2 cells (AT2 cells) were found to be etiologically implicated in COPD based on genetic and transcriptomic features. The most transcriptomically diversified states of the three etiological cell types showed significant enrichment in immune/inflammatory responses (monocytes and mast cells) and/or mitochondrial dysfunction (monocytes and AT2 cells). We then identified three chemical candidates that may potentially induce COPD by modulating gene expression patterns in the three etiological cell types. Overall, our study suggests the single-cell level mechanisms underlying the pathogenesis of COPD and may provide information on toxic compounds that could be potential risk factors for COPD.

Prediabetic HbA1c and Cortical Atrophy: Underlying Neurobiology

OBJECTIVE

To investigate the relationship between blood glycated hemoglobin (HbA1c) and cerebral cortical thickness (CT) and identify potential cellular mechanisms involved.

RESEARCH DESIGN AND METHODS

A cohort of 30,579 adults age 45 to 81 (mean ± SD: 64 ± 7.5) years with available data on brain MRI and blood HbA1c levels was analyzed. The relationship between HbA1c and CT was probed using independent spatial profiles of cell-specific gene expression. Lastly, a genome-wide association study was conducted on the shared variance between HbA1c and CT.

RESULTS

The HbA1c-CT association was noncontinuous, emerging negatively within the prediabetic range (39.6 mmol/mol). This association was strongest in brain regions with higher expression of genes specific to excitatory neurons and lower expression of genes specific to astrocytes and microglia. A significant locus implicated mitochondrial maintenance and ATP generation.

CONCLUSIONS

Effective glycemia control at prediabetic levels is warranted to preserve brain health and prevent prediabetes-related neurobiologic perturbations.

Central Adiposity Increases Risk of Kidney Stone Disease through Effects on Serum Calcium Concentrations

SIGNIFICANCE STATEMENT

Kidney stone disease is a common disorder with poorly understood pathophysiology. Observational and genetic studies indicate that adiposity is associated with an increased risk of kidney stone disease. However, the relative contribution of general and central adipose depots and the mechanisms by which effects of adiposity on kidney stone disease are mediated have not been defined. Using conventional and genetic epidemiological techniques, we demonstrate that general and central adiposity are independently associated with kidney stone disease. In addition, one mechanism by which central adiposity increases risk of kidney stone disease is by increasing serum calcium concentration. Therapies targeting adipose depots may affect calcium homeostasis and help to prevent kidney stone disease.

BACKGROUND

Kidney stone disease affects approximately 10% of individuals in their lifetime and is frequently recurrent. The disease is linked to obesity, but the mechanisms mediating this association are uncertain.

METHODS

Associations of adiposity and incident kidney stone disease were assessed in the UK Biobank over a mean of 11.6 years/person. Genome-wide association studies and Mendelian randomization (MR) analyses were undertaken in the UK Biobank, FinnGen, and in meta-analyzed cohorts to identify factors that affect kidney stone disease risk.

RESULTS

Observational analyses on UK Biobank data demonstrated that increasing central and general adiposity is independently associated with incident kidney stone formation. Multivariable MR, using meta-analyzed UK Biobank and FinnGen data, established that risk of kidney stone disease increases by approximately 21% per one standard deviation increase in body mass index (BMI, a marker of general adiposity) independent of waist-to-hip ratio (WHR, a marker of central adiposity) and approximately 24% per one standard deviation increase of WHR independent of BMI. Genetic analyses indicate that higher WHR, but not higher BMI, increases risk of kidney stone disease by elevating adjusted serum calcium concentrations (β=0.12 mmol/L); WHR mediates 12%-15% of its effect on kidney stone risk in this way.

CONCLUSIONS

Our study indicates that visceral adipose depots elevate serum calcium concentrations, resulting in increased risk of kidney stone disease. These findings highlight the importance of weight loss in individuals with recurrent kidney stones and suggest that therapies targeting adipose depots may affect calcium homeostasis and contribute to prevention of kidney stone disease.

Deciphering the genetic landscape of obesity: a data-driven approach to identifying plausible causal genes and therapeutic targets

OBJECTIVES

Genome-wide association studies (GWAS) have successfully revealed numerous susceptibility loci for obesity. However, identifying the causal genes, pathways, and tissues/cell types responsible for these associations remains a challenge, and standardized analysis workflows are lacking. Additionally, due to limited treatment options for obesity, there is a need for the development of new pharmacological therapies. This study aimed to address these issues by performing step-wise utilization of knowledgebase for gene prioritization and assessing the potential relevance of key obesity genes as therapeutic targets.

METHODS AND RESULTS

First, we generated a list of 28,787 obesity-associated SNPs from the publicly available GWAS dataset (approximately 800,000 individuals in the GIANT meta-analysis). Then, we prioritized 1372 genes with significant in silico evidence against genomic and transcriptomic data, including transcriptionally regulated genes in the brain from transcriptome-wide association studies. In further narrowing down the gene list, we selected key genes, which we found to be useful for the discovery of potential drug seeds as demonstrated in lipid GWAS separately. We thus identified 74 key genes for obesity, which are highly interconnected and enriched in several biological processes that contribute to obesity, including energy expenditure and homeostasis. Of 74 key genes, 37 had not been reported for the pathophysiology of obesity. Finally, by drug-gene interaction analysis, we detected 23 (of 74) key genes that are potential targets for 78 approved and marketed drugs.

CONCLUSIONS

Our results provide valuable insights into new treatment options for obesity through a data-driven approach that integrates multiple up-to-date knowledgebases.

Integrative single-cell meta-analysis reveals disease-relevant vascular cell states and markers in human atherosclerosis

Coronary artery disease (CAD) is characterized by atherosclerotic plaque formation in the arterial wall. CAD progression involves complex interactions and phenotypic plasticity among vascular and immune cell lineages. Single-cell RNA-seq (scRNA-seq) studies have highlighted lineage-specific transcriptomic signatures, but human cell phenotypes remain controversial. Here, we perform an integrated meta-analysis of 22 scRNA-seq libraries to generate a comprehensive map of human atherosclerosis with 118,578 cells. Besides characterizing granular cell-type diversity and communication, we leverage this atlas to provide insights into smooth muscle cell (SMC) modulation. We integrate genome-wide association study data and uncover a critical role for modulated SMC phenotypes in CAD, myocardial infarction, and coronary calcification. Finally, we identify fibromyocyte/fibrochondrogenic SMC markers (LTBP1 and CRTAC1) as proxies of atherosclerosis progression and validate these through omics and spatial imaging analyses. Altogether, we create a unified atlas of human atherosclerosis informing cell state-specific mechanistic and translational studies of cardiovascular diseases.

Functional interrogation of twenty type 2 diabetes-associated genes using isogenic human embryonic stem cell-derived β-like cells

Genetic studies have identified numerous loci associated with type 2 diabetes (T2D), but the functional roles of many loci remain unexplored. Here, we engineered isogenic knockout human embryonic stem cell lines for 20 genes associated with T2D risk. We examined the impacts of each knockout on β cell differentiation, functions, and survival. We generated gene expression and chromatin accessibility profiles on β cells derived from each knockout line. Analyses of T2D-association signals overlapping HNF4A-dependent ATAC peaks identified a likely causal variant at the FAIM2 T2D-association signal. Additionally, the integrative association analyses identified four genes (CP, RNASE1, PCSK1N, and GSTA2) associated with insulin production, and two genes (TAGLN3 and DHRS2) associated with β cell sensitivity to lipotoxicity. Finally, we leveraged deep ATAC-seq read coverage to assess allele-specific imbalance at variants heterozygous in the parental line and identified a single likely functional variant at each of 23 T2D-association signals.

Time-varying and tissue-dependent effects of adiposity on leptin levels: A Mendelian randomization study

Background

Findings from Mendelian randomization (MR) studies are conventionally interpreted as lifelong effects, which typically do not provide insight into the molecular mechanisms underlying the effect of an exposure on an outcome. In this study, we apply two recently developed MR approaches (known as 'lifecourse' and 'tissue-partitioned' MR) to investigate lifestage-specific effects and tissues of action in the relationship between adiposity and circulating leptin levels.

Methods

Genetic instruments for childhood and adult adiposity were incorporated into a multivariable MR (MVMR) framework to estimate lifestage-specific effects on leptin levels measured during early life (mean age: 10 y) in the Avon Longitudinal Study of Parents and Children and in adulthood (mean age: 55 y) using summary-level data from the deCODE Health study. This was followed by partitioning body mass index (BMI) instruments into those whose effects are putatively mediated by gene expression in either subcutaneous adipose or brain tissues, followed by using MVMR to simultaneously estimate their separate effects on childhood and adult leptin levels.

Results

There was strong evidence that childhood adiposity has a direct effect on leptin levels at age 10 y in the lifecourse (β = 1.10 SD change in leptin levels, 95% CI = 0.90-1.30, p=6 × 10-28), whereas evidence of an indirect effect was found on adulthood leptin along the causal pathway involving adulthood body size (β = 0.74, 95% CI = 0.62-0.86, p=1 × 10-33). Tissue-partitioned MR analyses provided evidence to suggest that BMI exerts its effect on leptin levels during both childhood and adulthood via brain tissue-mediated pathways (β = 0.79, 95% CI = 0.22-1.36, p=6 × 10-3 and β = 0.51, 95% CI = 0.32-0.69, p=1 × 10-7, respectively).

Conclusions

Our findings demonstrate the use of lifecourse MR to disentangle direct and indirect effects of early-life exposures on time-varying complex outcomes. Furthermore, by integrating tissue-specific data, we highlight the etiological importance of appetite regulation in the effect of adiposity on leptin levels.

Funding

This work was supported by the Integrative Epidemiology Unit, which receives funding from the UK Medical Research Council and the University of Bristol (MC_UU_00011/1).

Expanding causal genes for Parkinson's disease via multi-omics analysis

Genome‑wide association studies (GWASs) have revealed numerous loci associated with Parkinson's disease (PD). However, some potential causal/risk genes were still not revealed and no etiological therapies are available. To find potential causal genes and explore genetically supported drug targets for PD is urgent. By integrating the expression quantitative trait loci (eQTL) and protein quantitative trait loci (pQTL) datasets from multiple tissues (blood, cerebrospinal fluid (CSF) and brain) and PD GWAS summary statistics, a pipeline combing Mendelian randomization (MR), Steiger filtering analysis, Bayesian colocalization, fine mapping, Protein-protein network and enrichment analysis were applied to identify potential causal genes for PD. As a result, GPNMB displayed a robust causal role for PD at the protein level in the blood, CSF and brain, and transcriptional level in the brain, while the protective role of CD38 (in brain pQTL and eQTL) was also identified. We also found inconsistent roles of DGKQ on PD between protein and mRNA levels. Another 9 proteins (CTSB, ARSA, SEC23IP, CD84, ENTPD1, FCGR2B, BAG3, SNCA, FCGR2A) were associated with the risk for PD based on only a single pQTL after multiple corrections. We also identified some proteins' interactions with known PD causative genes and therapeutic targets. In conclusion, this study suggested GPNMB, CD38, and DGKQ may act in the pathogenesis of PD, but whether the other proteins involved in PD needs more evidence. These findings would help to uncover the genes underlying PD and prioritize targets for future therapeutic interventions.

Single-nucleus transcriptomes reveal spatiotemporal symbiotic perception and early response in Medicago

Establishing legume-rhizobial symbiosis requires precise coordination of complex responses in a time- and cell type-specific manner. Encountering Rhizobium, rapid changes of gene expression levels in host plants occur in the first few hours, which prepare the plants to turn off defence and form a symbiotic relationship with the microbes. Here, we applied single-nucleus RNA sequencing to characterize the roots of Medicago truncatula at 30 min, 6 h and 24 h after nod factor treatment. We found drastic global gene expression reprogramming at 30 min in the epidermis and cortex and most of these changes were restored at 6 h. Moreover, plant defence response genes are activated at 30 min and subsequently suppressed at 6 h in non-meristem cells. Only in the cortical cells but not in other cell types, we found the flavonoid synthase genes required to recruit rhizobia are highly expressed 30 min after inoculation with nod factors. A gene module enriched for symbiotic nitrogen fixation genes showed that MtFER (MtFERONIA) and LYK3 (LysM domain receptor-like kinase 3) share similar responses to symbiotic signals. We further found that MtFER can be phosphorylated by LYK3 and it participates in rhizobial symbiosis. Our results expand our understanding of dynamic spatiotemporal symbiotic responses at the single-cell level.

GWAS of random glucose in 476,326 individuals provide insights into diabetes pathophysiology, complications and treatment stratification

Conventional measurements of fasting and postprandial blood glucose levels investigated in genome-wide association studies (GWAS) cannot capture the effects of DNA variability on 'around the clock' glucoregulatory processes. Here we show that GWAS meta-analysis of glucose measurements under nonstandardized conditions (random glucose (RG)) in 476,326 individuals of diverse ancestries and without diabetes enables locus discovery and innovative pathophysiological observations. We discovered 120 RG loci represented by 150 distinct signals, including 13 with sex-dimorphic effects, two cross-ancestry and seven rare frequency signals. Of these, 44 loci are new for glycemic traits. Regulatory, glycosylation and metagenomic annotations highlight ileum and colon tissues, indicating an underappreciated role of the gastrointestinal tract in controlling blood glucose. Functional follow-up and molecular dynamics simulations of lower frequency coding variants in glucagon-like peptide-1 receptor (GLP1R), a type 2 diabetes treatment target, reveal that optimal selection of GLP-1R agonist therapy will benefit from tailored genetic stratification. We also provide evidence from Mendelian randomization that lung function is modulated by blood glucose and that pulmonary dysfunction is a diabetes complication. Our investigation yields new insights into the biology of glucose regulation, diabetes complications and pathways for treatment stratification.

Single-Cell Transcriptomics of Bone Marrow Stromal Cells in Diversity Outbred Mice: A Model for Population-Level scRNA-Seq Studies

Genome-wide association studies (GWASs) have advanced our understanding of the genetics of osteoporosis; however, the challenge has been converting associations to causal genes. Studies have utilized transcriptomics data to link disease-associated variants to genes, but few population transcriptomics data sets have been generated on bone at the single-cell level. To address this challenge, we profiled the transcriptomes of bone marrow-derived stromal cells (BMSCs) cultured under osteogenic conditions from five diversity outbred (DO) mice using single-cell RNA-seq (scRNA-seq). The goal of the study was to determine if BMSCs could serve as a model to generate cell type-specific transcriptomic profiles of mesenchymal lineage cells from large populations of mice to inform genetic studies. By enriching for mesenchymal lineage cells in vitro, coupled with pooling of multiple samples and downstream genotype deconvolution, we demonstrate the scalability of this model for population-level studies. We demonstrate that dissociation of BMSCs from a heavily mineralized matrix had little effect on viability or their transcriptomic signatures. Furthermore, we show that BMSCs cultured under osteogenic conditions are diverse and consist of cells with characteristics of mesenchymal progenitors, marrow adipogenic lineage precursors (MALPs), osteoblasts, osteocyte-like cells, and immune cells. Importantly, all cells were similar from a transcriptomic perspective to cells isolated in vivo. We employed scRNA-seq analytical tools to confirm the biological identity of profiled cell types. SCENIC was used to reconstruct gene regulatory networks (GRNs), and we observed that cell types show GRNs expected of osteogenic and pre-adipogenic lineage cells. Further, CELLECT analysis showed that osteoblasts, osteocyte-like cells, and MALPs captured a significant component of bone mineral density (BMD) heritability. Together, these data suggest that BMSCs cultured under osteogenic conditions coupled with scRNA-seq can be used as a scalable and biologically informative model to generate cell type-specific transcriptomic profiles of mesenchymal lineage cells in large populations. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Multivariate genome-wide analysis of aging-related traits identifies novel loci and new drug targets for healthy aging

The concept of aging is complex, including many related phenotypes such as healthspan, lifespan, extreme longevity, frailty and epigenetic aging, suggesting shared biological underpinnings; however, aging-related endpoints have been primarily assessed individually. Using data from these traits and multivariate genome-wide association study methods, we modeled their underlying genetic factor ('mvAge'). mvAge (effective n = ~1.9 million participants of European ancestry) identified 52 independent variants in 38 genomic loci. Twenty variants were novel (not reported in input genome-wide association studies). Transcriptomic imputation identified age-relevant genes, including VEGFA and PHB1. Drug-target Mendelian randomization with metformin target genes showed a beneficial impact on mvAge (P value = 8.41 × 10-5). Similarly, genetically proxied thiazolidinediones (P value = 3.50 × 10-10), proprotein convertase subtilisin/kexin 9 inhibition (P value = 1.62 × 10-6), angiopoietin-like protein 4, beta blockers and calcium channel blockers also had beneficial Mendelian randomization estimates. Extending the drug-target Mendelian randomization framework to 3,947 protein-coding genes prioritized 122 targets. Together, these findings will inform future studies aimed at improving healthy aging.

Single-cell transcriptomic profiling of human pancreatic islets reveals genes responsive to glucose exposure over 24 hours

Disruption of pancreatic islet function and glucose homeostasis can lead to the development of sustained hyperglycemia, beta cell glucotoxicity, and ultimately type 2 diabetes (T2D). In this study, we sought to explore the effects of hyperglycemia on human pancreatic islet (HPI) gene expression by exposing HPIs from two donors to low (2.8mM) and high (15.0mM) glucose concentrations over 24 hours, assaying the transcriptome at seven time points using single-cell RNA sequencing (scRNA-seq). We modeled time as both a discrete and continuous variable to determine momentary and longitudinal changes in transcription associated with islet time in culture or glucose exposure. Across all cell types, we identified 1,528 genes associated with time, 1,185 genes associated with glucose exposure, and 845 genes associated with interaction effects between time and glucose. We clustered differentially expressed genes across cell types and found 347 modules of genes with similar expression patterns across time and glucose conditions, including two beta cell modules enriched in genes associated with T2D. Finally, by integrating genomic features from this study and genetic summary statistics for T2D and related traits, we nominate 363 candidate effector genes that may underlie genetic associations for T2D and related traits.

Identifying novel proteins underlying loneliness by integrating GWAS summary data with human brain proteomes

Enduring loneliness is associated with mental disorders and physical diseases. Although genome-wide association studies (GWAS) have identified risk loci associated with loneliness, how these loci confer the risk remains largely unknown. In the current study, we aimed to investigate key proteins underlying loneliness in the brain by integrating human brain proteomes and transcriptomes with loneliness GWAS to perform a discovery proteome-wide association study (PWAS), followed by a confirmatory PWAS, transcriptome-wide association analysis (TWAS), Mendelian randomization (MR), Steigering filtering analysis and Bayesian colocalization analysis. Moreover, given the fact that loneliness is associated with mental disorders, we explored the shared genetic architecture between loneliness and mental disorders. Totally, we identified 18 genes to be associated with loneliness via their cis-regulated brain protein abundance. Eleven of the 18 genes (61.1%) were replicated in the confirmatory PWAS, and mRNA levels of 4 genes were further validated to be associated with loneliness.MR and genetic colocalization analysis further confirmed that the increased protein abundance of ALDH2 and ICA1L was protective against loneliness, while the increased protein abundance of GPX1 was a risk for developing loneliness. Furthermore, we found genetic correlations, bidirectional causal associations and overlapping phenotype-associated protein profiles between loneliness and mental disorders including major depression and schizophrenia. In summary, our findings provided clues about the brain-related molecular basis underlying loneliness, which warrants further investigation.

Adipogenic and SWAT cells separate from a common progenitor in human brown and white adipose depots

Adipocyte function is a major determinant of metabolic disease, warranting investigations of regulating mechanisms. We show at single-cell resolution that progenitor cells from four human brown and white adipose depots separate into two main cell fates, an adipogenic and a structural branch, developing from a common progenitor. The adipogenic gene signature contains mitochondrial activity genes, and associates with genome-wide association study traits for fat distribution. Based on an extracellular matrix and developmental gene signature, we name the structural branch of cells structural Wnt-regulated adipose tissue-resident (SWAT) cells. When stripped from adipogenic cells, SWAT cells display a multipotent phenotype by reverting towards progenitor state or differentiating into new adipogenic cells, dependent on media. Label transfer algorithms recapitulate the cell types in human adipose tissue datasets. In conclusion, we provide a differentiation map of human adipocytes and define the multipotent SWAT cell, providing a new perspective on adipose tissue regulation.

Cross-disorder genetic analysis of immune diseases reveals distinct gene associations that converge on common pathways

Genome-wide association studies (GWAS) have mapped thousands of susceptibility loci associated with immune-mediated diseases. To assess the extent of the genetic sharing across nine immune-mediated diseases we apply genomic structural equation modelling to GWAS data from European populations. We identify three disease groups: gastrointestinal tract diseases, rheumatic and systemic diseases, and allergic diseases. Although loci associated with the disease groups are highly specific, they converge on perturbing the same pathways. Finally, we test for colocalization between loci and single-cell eQTLs derived from peripheral blood mononuclear cells. We identify the causal route by which 46 loci predispose to three disease groups and find evidence for eight genes being candidates for drug repurposing. Taken together, here we show that different constellations of diseases have distinct patterns of genetic associations, but that associated loci converge on perturbing different nodes in T cell activation and signalling pathways.

Functional interrogation of twenty type 2 diabetes-associated genes using isogenic hESC-derived β-like cells

Genetic studies have identified numerous loci associated with type 2 diabetes (T2D), but the functional role of many loci has remained unexplored. In this study, we engineered isogenic knockout human embryonic stem cell (hESC) lines for 20 genes associated with T2D risk. We systematically examined β-cell differentiation, insulin production and secretion, and survival. We performed RNA-seq and ATAC-seq on hESC-β cells from each knockout line. Analyses of T2D GWAS signals overlapping with HNF4A-dependent ATAC peaks identified a specific SNP as a likely causal variant. In addition, we performed integrative association analyses and identified four genes ( CP, RNASE1, PCSK1N and GSTA2 ) associated with insulin production, and two genes ( TAGLN3 and DHRS2 ) associated with sensitivity to lipotoxicity. Finally, we leveraged deep ATAC-seq read coverage to assess allele-specific imbalance at variants heterozygous in the parental hESC line, to identify a single likely functional variant at each of 23 T2D GWAS signals.

Brain transcriptome-wide association study implicates novel risk genes underlying schizophrenia risk

BACKGROUND

To identify risk genes whose expression are regulated by the reported risk variants and to explore the potential regulatory mechanism in schizophrenia (SCZ).

METHODS

We systematically integrated three independent brain expression quantitative traits (eQTLs) (CommonMind, GTEx, and BrainSeq Phase 2, a total of 1039 individuals) and GWAS data (56 418 cases and 78 818 controls), with the use of transcriptome-wide association study (TWAS). Diffusion magnetic resonance imaging was utilized to quantify the integrity of white matter bundles and determine whether polygenic risk of novel genes linked to brain structure was present in patients with first-episode antipsychotic SCZ.

RESULTS

TWAS showed that eight risk genes (CORO7, DDAH2, DDHD2, ELAC2, GLT8D1, PCDHA8, THOC7, and TYW5) reached transcriptome-wide significance (TWS) level. These findings were confirmed by an independent integrative approach (i.e. Sherlock). We further conducted conditional analyses and identified the potential risk genes that driven the TWAS association signal in each locus. Gene expression analysis showed that several TWS genes (including CORO7, DDAH2, DDHD2, ELAC2, GLT8D1, THOC7 and TYW5) were dysregulated in the dorsolateral prefrontal cortex of SCZ cases compared with controls. TWS genes were mainly expressed on the surface of glutamatergic neurons, GABAergic neurons, and microglia. Finally, SCZ cases had a substantially greater TWS genes-based polygenic risk (PRS) compared to controls, and we showed that fractional anisotropy of the cingulum-hippocampus mediates the influence of TWS genes PRS on SCZ.

CONCLUSIONS

Our findings identified novel SCZ risk genes and highlighted the importance of the TWS genes in frontal-limbic dysfunctions in SCZ, indicating possible therapeutic targets.

Multi-omic underpinnings of epigenetic aging and human longevity

Biological aging is accompanied by increasing morbidity, mortality, and healthcare costs; however, its molecular mechanisms are poorly understood. Here, we use multi-omic methods to integrate genomic, transcriptomic, and metabolomic data and identify biological associations with four measures of epigenetic age acceleration and a human longevity phenotype comprising healthspan, lifespan, and exceptional longevity (multivariate longevity). Using transcriptomic imputation, fine-mapping, and conditional analysis, we identify 22 high confidence associations with epigenetic age acceleration and seven with multivariate longevity. FLOT1, KPNA4, and TMX2 are novel, high confidence genes associated with epigenetic age acceleration. In parallel, cis-instrument Mendelian randomization of the druggable genome associates TPMT and NHLRC1 with epigenetic aging, supporting transcriptomic imputation findings. Metabolomics Mendelian randomization identifies a negative effect of non-high-density lipoprotein cholesterol and associated lipoproteins on multivariate longevity, but not epigenetic age acceleration. Finally, cell-type enrichment analysis implicates immune cells and precursors in epigenetic age acceleration and, more modestly, multivariate longevity. Follow-up Mendelian randomization of immune cell traits suggests lymphocyte subpopulations and lymphocytic surface molecules affect multivariate longevity and epigenetic age acceleration. Our results highlight druggable targets and biological pathways involved in aging and facilitate multi-omic comparisons of epigenetic clocks and human longevity.

Integrated single-nucleus and spatial transcriptomics captures transitional states in soybean nodule maturation

Legumes form symbiosis with rhizobium leading to the development of nitrogen-fixing nodules. By integrating single-nucleus and spatial transcriptomics, we established a cell atlas of soybean nodules and roots. In central infected zones of nodules, we found that uninfected cells specialize into functionally distinct subgroups during nodule development, and revealed a transitional subtype of infected cells with enriched nodulation-related genes. Overall, our results provide a single-cell perspective for understanding rhizobium-legume symbiosis.

Single-cell omics: A new direction for functional genetic research in human diseases and animal models

Over the past decade, with the development of high-throughput single-cell sequencing technology, single-cell omics has been emerged as a powerful tool to understand the molecular basis of cellular mechanisms and refine our knowledge of diverse cell states. They can reveal the heterogeneity at different genetic layers and elucidate their associations by multiple omics analysis, providing a more comprehensive genetic map of biological regulatory networks. In the post-GWAS era, the molecular biological mechanisms influencing human diseases will be further elucidated by single-cell omics. This review mainly summarizes the development and trend of single-cell omics. This involves single-cell omics technologies, single-cell multi-omics technologies, multiple omics data integration methods, applications in various human organs and diseases, classic laboratory cell lines, and animal disease models. The review will reveal some perspectives for elucidating human diseases and constructing animal models.

Identification of MRAP protein family as broad-spectrum GPCR modulators

BACKGROUND

The melanocortin receptor accessory proteins (MRAP1 and MRAP2) are well-known endocrine regulators for the trafficking and signalling of all five melanocortin receptors (MC1R-MC5R). The observation of MRAP2 on regulating several non-melanocortin G protein-coupled receptors (GPCRs) has been sporadically reported, whereas other endogenous GPCR partners of the MRAP protein family are largely unknown.

METHODS

Here, we performed single-cell transcriptome analysis and drew a fine GPCR blueprint and MRAPs-associated network of two major endocrine organs, the hypothalamus and adrenal gland at single-cell resolution. We also integrated multiple bulk RNA-seq profiles and single-cell datasets of human and mouse tissues, and narrowed down a list of 48 GPCRs with strong endogenous co-expression correlation with MRAPs.

RESULTS

36 and 46 metabolic-related GPCRs were consequently identified as novel interacting partners of MRAP1 or MRAP2, respectively. MRAPs exhibited protein-protein interactions and varying pharmacological properties on the surface translocation, constitutive activities and ligand-stimulated downstream signalling of these GPCRs. Knockdown of MRAP2 expression by hypothalamic administration of adeno-associated virus (AAV) packed shRNA stimulated body weight gain in mouse model. Co-injection of corticotropinreleasing factor (CRF), the agonist of corticotropin releasing hormone receptor 1 (CRHR1), suppressed feeding behaviour in a MRAP2-dependent manner.

CONCLUSIONS

Collectively, our study has comprehensively elucidated the complex GPCR networks in two major endocrine organs and redefined the MRAP protein family as broad-spectrum GPCR modulators. MRAP proteins not only serve as a vital endocrine pivot on the regulation of global GPCR activities in vivo that could explain the composite physiological phenotypes of the MRAP2 null murine model but also provide us with new insights of the phenotyping investigation of GPCR-MRAP functional complexes.

scGWAS: landscape of trait-cell type associations by integrating single-cell transcriptomics-wide and genome-wide association studies

BACKGROUND

The rapid accumulation of single-cell RNA sequencing (scRNA-seq) data presents unique opportunities to decode the genetically mediated cell-type specificity in complex diseases. Here, we develop a new method, scGWAS, which effectively leverages scRNA-seq data to achieve two goals: (1) to infer the cell types in which the disease-associated genes manifest and (2) to construct cellular modules which imply disease-specific activation of different processes.

RESULTS

scGWAS only utilizes the average gene expression for each cell type followed by virtual search processes to construct the null distributions of module scores, making it scalable to large scRNA-seq datasets. We demonstrated scGWAS in 40 genome-wide association studies (GWAS) datasets (average sample size N ≈ 154,000) using 18 scRNA-seq datasets from nine major human/mouse tissues (totaling 1.08 million cells) and identified 2533 trait and cell-type associations, each with significant modules for further investigation. The module genes were validated using disease or clinically annotated references from ClinVar, OMIM, and pLI variants.

CONCLUSIONS

We showed that the trait-cell type associations identified by scGWAS, while generally constrained to trait-tissue associations, could recapitulate many well-studied relationships and also reveal novel relationships, providing insights into the unsolved trait-tissue associations. Moreover, in each specific cell type, the associations with different traits were often mediated by different sets of risk genes, implying disease-specific activation of driving processes. In summary, scGWAS is a powerful tool for exploring the genetic basis of complex diseases at the cell type level using single-cell expression data.

Upper cortical layer-driven network impairment in schizophrenia

Schizophrenia is one of the most widespread and complex mental disorders. To characterize the impact of schizophrenia, we performed single-nucleus RNA sequencing (snRNA-seq) of >220,000 neurons from the dorsolateral prefrontal cortex of patients with schizophrenia and matched controls. In addition, >115,000 neurons were analyzed topographically by immunohistochemistry. Compositional analysis of snRNA-seq data revealed a reduction in abundance of GABAergic neurons and a concomitant increase in principal neurons, most pronounced for upper cortical layer subtypes, which was substantiated by histological analysis. Many neuronal subtypes showed extensive transcriptomic changes, the most marked in upper-layer GABAergic neurons, including down-regulation in energy metabolism and up-regulation in neurotransmission. Transcription factor network analysis demonstrated a developmental origin of transcriptomic changes. Last, Visium spatial transcriptomics further corroborated upper-layer neuron vulnerability in schizophrenia. Overall, our results point toward general network impairment within upper cortical layers as a core substrate associated with schizophrenia symptomatology.

Single-Cell Chromatin Accessibility Data Combined with GWAS Improves Detection of Relevant Cell Types in 59 Complex Phenotypes

Several disease risk variants reside on non-coding regions of DNA, particularly in open chromatin regions of specific cell types. Identifying the cell types relevant to complex traits through the integration of chromatin accessibility data and genome-wide association studies (GWAS) data can help to elucidate the mechanisms of these traits. In this study, we created a collection of associations between the combinations of chromatin accessibility data (bulk and single-cell) with an array of 201 complex phenotypes. We integrated the GWAS data of these 201 phenotypes with bulk chromatin accessibility data from 137 cell types measured by DNase-I hypersensitive sequencing and found significant results (FDR adjusted p-value ≤ 0.05) for at least one cell type in 21 complex phenotypes, such as atopic dermatitis, Graves’ disease, and body mass index. With the integration of single-cell chromatin accessibility data measured by an assay for transposase-accessible chromatin with high-throughput sequencing (scATAC-seq), taken from 111 adult and 111 fetal cell types, the resolution of association was magnified, enabling the identification of further cell types. This resulted in the identification of significant correlations (FDR adjusted p-value ≤ 0.05) between 15 categories of single-cell subtypes and 59 phenotypes ranging from autoimmune diseases like Graves’ disease to cardiovascular traits like diastolic/systolic blood pressure.

3D chromatin maps of the human pancreas reveal lineage-specific regulatory architecture of T2D risk

Three-dimensional (3D) chromatin organization maps help dissect cell-type-specific gene regulatory programs. Furthermore, 3D chromatin maps contribute to elucidating the pathogenesis of complex genetic diseases by connecting distal regulatory regions and genetic risk variants to their respective target genes. To understand the cell-type-specific regulatory architecture of diabetes risk, we generated transcriptomic and 3D epigenomic profiles of human pancreatic acinar, alpha, and beta cells using single-cell RNA-seq, single-cell ATAC-seq, and high-resolution Hi-C of sorted cells. Comparisons of these profiles revealed differential A/B (open/closed) chromatin compartmentalization, chromatin looping, and transcriptional factor-mediated control of cell-type-specific gene regulatory programs. We identified a total of 4,750 putative causal-variant-to-target-gene pairs at 194 type 2 diabetes GWAS signals using pancreatic 3D chromatin maps. We found that the connections between candidate causal variants and their putative target effector genes are cell-type stratified and emphasize previously underappreciated roles for alpha and acinar cells in diabetes pathogenesis.

Association Between Brain Structure and Alcohol Use Behaviors in Adults: A Mendelian Randomization and Multiomics Study

IMPORTANCE

Past studies have identified associations between brain macrostructure and alcohol use behaviors. However, identifying directional associations between these phenotypes is difficult due to the limitations of observational studies.

OBJECTIVE

To use mendelian randomization (MR) to identify directional associations between brain structure and alcohol use and elucidate the transcriptomic and cellular underpinnings of identified associations.

DESIGN, SETTING, AND PARTICIPANTS

The main source data comprised summary statistics from population-based and case-control genome-wide association studies (GWAS) of neuroimaging, behavioral, and clinical phenotypes (N = 763 874). Using these data, bidirectional and multivariable MR was performed analyzing associations between brain macrostructure and alcohol use. Downstream transcriptome-wide association studies (TWAS) and cell-type enrichment analyses investigated the biology underlying identified associations. The study approach was data driven and did not test any a priori hypotheses. Data were analyzed August 2021 to May 2022.

MAIN OUTCOMES AND MEASURES

Brain structure phenotypes (global cortical thickness [GCT] and global cortical surface area [GCSA] in 33 709 individuals and left-right subcortical volumes in 19 629 individuals) and alcohol use behaviors (alcoholic drinks per week [DPW] in 537 349 individuals, binge drinking frequency in 143 685 individuals, and alcohol use disorder in 8845 individuals vs 20 657 control individuals [total of 29 502]).

RESULTS

The main bidirectional MR analyses were performed in samples totaling 763 874 individuals, among whom more than 94% were of European ancestry, 52% to 54% were female, and the mean cohort ages were 40 to 63 years. Negative associations were identified between genetically predicted GCT and binge drinking (β, -2.52; 95% CI, -4.13 to -0.91) and DPW (β, -0.88; 95% CI, -1.37 to -0.40) at a false discovery rate (FDR) of 0.05. These associations remained significant in multivariable MR models that accounted for neuropsychiatric phenotypes, substance use, trauma, and neurodegeneration. TWAS of GCT and alcohol use behaviors identified 5 genes at the 17q21.31 locus oppositely associated with GCT and binge drinking or DPW (FDR = 0.05). Cell-type enrichment analyses implicated glutamatergic cortical neurons in alcohol use behaviors.

CONCLUSIONS AND RELEVANCE

The findings in this study show that the associations between GCT and alcohol use may reflect a predispositional influence of GCT and that 17q21.31 genes and glutamatergic cortical neurons may play a role in this association. While replication studies are needed, these findings should enhance the understanding of associations between brain structure and alcohol use.

Genome-wide association analyses of physical activity and sedentary behavior provide insights into underlying mechanisms and roles in disease prevention

Although physical activity and sedentary behavior are moderately heritable, little is known about the mechanisms that influence these traits. Combining data for up to 703,901 individuals from 51 studies in a multi-ancestry meta-analysis of genome-wide association studies yields 99 loci that associate with self-reported moderate-to-vigorous intensity physical activity during leisure time (MVPA), leisure screen time (LST) and/or sedentary behavior at work. Loci associated with LST are enriched for genes whose expression in skeletal muscle is altered by resistance training. A missense variant in ACTN3 makes the alpha-actinin-3 filaments more flexible, resulting in lower maximal force in isolated type IIA muscle fibers, and possibly protection from exercise-induced muscle damage. Finally, Mendelian randomization analyses show that beneficial effects of lower LST and higher MVPA on several risk factors and diseases are mediated or confounded by body mass index (BMI). Our results provide insights into physical activity mechanisms and its role in disease prevention.

Identification of novel proteins for lacunar stroke by integrating genome-wide association data and human brain proteomes

BACKGROUND

Previous genome-wide association studies (GWAS) have identified numerous risk genes for lacunar stroke, but it is challenging to decipher how they confer risk for the disease. We employed an integrative analytical pipeline to efficiently transform genetic associations to identify novel proteins for lacunar stroke.

METHODS

We systematically integrated lacunar stroke genome-wide association study (GWAS) (N=7338) with human brain proteomes (N=376) to perform proteome-wide association studies (PWAS), Mendelian randomization (MR), and Bayesian colocalization. We also used an independent human brain proteomic dataset (N=152) to annotate the new genes.

RESULTS

We found that the protein abundance of seven genes (ICA1L, CAND2, ALDH2, MADD, MRVI1, CSPG4, and PTPN11) in the brain was associated with lacunar stroke. These seven genes were mainly expressed on the surface of glutamatergic neurons, GABAergic neurons, and astrocytes. Three genes (ICA1L, CAND2, ALDH2) were causal in lacunar stroke (P < 0.05/proteins identified for PWAS; posterior probability of hypothesis 4 ≥ 75 % for Bayesian colocalization), and they were linked with lacunar stroke in confirmatory PWAS and independent MR. We also found that ICA1L is related to lacunar stroke at the brain transcriptome level.

CONCLUSIONS

Our present proteomic findings have identified ICA1L, CAND2, and ALDH2 as compelling genes that may give key hints for future functional research and possible therapeutic targets for lacunar stroke.

EPIC: Inferring relevant cell types for complex traits by integrating genome-wide association studies and single-cell RNA sequencing

More than a decade of genome-wide association studies (GWASs) have identified genetic risk variants that are significantly associated with complex traits. Emerging evidence suggests that the function of trait-associated variants likely acts in a tissue- or cell-type-specific fashion. Yet, it remains challenging to prioritize trait-relevant tissues or cell types to elucidate disease etiology. Here, we present EPIC (cEll tyPe enrIChment), a statistical framework that relates large-scale GWAS summary statistics to cell-type-specific gene expression measurements from single-cell RNA sequencing (scRNA-seq). We derive powerful gene-level test statistics for common and rare variants, separately and jointly, and adopt generalized least squares to prioritize trait-relevant cell types while accounting for the correlation structures both within and between genes. Using enrichment of loci associated with four lipid traits in the liver and enrichment of loci associated with three neurological disorders in the brain as ground truths, we show that EPIC outperforms existing methods. We apply our framework to multiple scRNA-seq datasets from different platforms and identify cell types underlying type 2 diabetes and schizophrenia. The enrichment is replicated using independent GWAS and scRNA-seq datasets and further validated using PubMed search and existing bulk case-control testing results.