Identifying disease-critical cell types and cellular processes by integrating single-cell RNA-sequencing and human genetics

Uncovering cell-type-specific immunomodulatory variants and molecular phenotypes in COVID-19 using structurally resolved protein networks

Immunomodulatory variants that lead to the loss or gain of specific protein interactions often manifest only as organismal phenotypes in infectious disease. Here, we propose a network-based approach to integrate genetic variation with a structurally resolved human protein interactome network to prioritize immunomodulatory variants in COVID-19. We find that, in addition to variants that pass genome-wide significance thresholds, variants at the interface of specific protein-protein interactions, even though they do not meet genome-wide thresholds, are equally immunomodulatory. The integration of these variants with single-cell epigenomic and transcriptomic data prioritizes myeloid and T cell subsets as the most affected by these variants across both the peripheral blood and the lung compartments. Of particular interest is a common coding variant that disrupts the OAS1-PRMT6 interaction and affects downstream interferon signaling. Critically, our framework is generalizable across infectious disease contexts and can be used to implicate immunomodulatory variants that do not meet genome-wide significance thresholds.

Differential expression and co-expression reveal cell types relevant to genetic disorder phenotypes

MOTIVATION

Knowledge of the specific cell types affected by genetic alterations in rare diseases is crucial for advancing diagnostics and treatments. Despite significant progress, the cell types involved in the majority of rare disease manifestations remain largely unknown. In this study, we integrated scRNA-seq data from non-diseased samples with known genetic disorder genes and phenotypic information to predict the specific cell types disrupted by pathogenic mutations for 482 disease phenotypes.

RESULTS

We found significant phenotype-cell type associations focusing on differential expression and co-expression mechanisms. Our analysis revealed that 13% of the associations documented in the literature were captured through differential expression, while 42% were elucidated through co-expression analysis, also uncovering potential new associations. These findings underscore the critical role of cellular context in disease manifestation and highlight the potential of single-cell data for the development of cell-aware diagnostics and targeted therapies for rare diseases.

AVAILABILITY AND IMPLEMENTATION

All code generated in this work is available at https://github.com/SergioAlias/sc-coex.

scPerb: Predict single-cell perturbation via style transfer-based variational autoencoder

INTRODUCTION

Traditional methods for obtaining cellular responses after perturbation are usually labor-intensive and costly, especially when working with multiple different experimental conditions. Therefore, accurate prediction of cellular responses to perturbations is of great importance in computational biology. Existing methodologies, such as graph-based approaches, vector arithmetic, and neural networks, either mix perturbation-related variances with cell-type-specific patterns or implicitly distinguish them within black-box models.

OBJECTIVES

This study aims to introduce and demonstrate a novel framework, scPerb, which explicitly extracts perturbation-related variances and transfers them from unperturbed to perturbed cells to accurately predict the effect of perturbation in single-cell level.

METHODS

scPerb utilizes a style transfer strategy by incorporating a style encoder into the architecture of a variational autoencoder. The style encoder captures the differences in latent representations between unperturbed and perturbed cells, enabling accurate prediction of post-perturbation gene expression data.

RESULTS

Comprehensive comparisons with existing methods demonstrate that scPerb delivers improved performance and higher accuracy in predicting cellular responses to perturbations. Notably, scPerb outperforms other methods across multiple datasets, achieving superior R2 values of 0.98, 0.98, and 0.96 on three benchmarking datasets.

CONCLUSION

scPerb offers a significant advancement in predicting cellular responses by effectively separating and transferring perturbation-related variances. This framework not only enhances prediction accuracy but also provides a robust tool for computational biology, addressing the limitations of current methodologies.

Interpretable GWAS by linking clinical phenotypes to quantifiable immune repertoire components

Bridging the gap between genotype and phenotype in GWAS studies is challenging. A multitude of genetic variants have been associated with immune-related diseases, including cancer, yet the interpretability of most variants remains low. Here, we investigate the quantitative components in the T cell receptor (TCR) repertoire, the frequency of clusters of TCR sequences predicted to have common antigen specificity, to interpret the genetic associations of diverse human diseases. We first developed a statistical model to predict the TCR components using variants in the TRB and HLA loci. Applying this model to over 300,000 individuals in the UK Biobank data, we identified 2309 associations between TCR abundances and various immune diseases. TCR clusters predicted to be pathogenic for autoimmune diseases were significantly enriched for predicted autoantigen-specificity. Moreover, four TCR clusters were associated with better outcomes in distinct cancers, where conventional GWAS cannot identify any significant locus. Collectively, our results highlight the integral role of adaptive immune responses in explaining the associations between genotype and phenotype.

Genes with differential expression across ancestries are enriched in ancestry-specific disease effects likely due to gene-by-environment interactions

Multi-ancestry genome-wide association studies (GWASs) have highlighted the existence of variants with ancestry-specific effect sizes. Understanding where and why these ancestry-specific effects occur is fundamental to understanding the genetic basis of human diseases and complex traits. Here, we characterized genes differentially expressed across ancestries (ancDE genes) at the cell-type level by leveraging single-cell RNA-sequencing data in peripheral blood mononuclear cells for 21 individuals with East Asian (EAS) ancestry and 23 individuals with European (EUR) ancestry (172,385 cells); then, we tested whether variants surrounding those genes were enriched in disease variants with ancestry-specific effect sizes by leveraging ancestry-matched GWASs of 31 diseases and complex traits (average n ∼ 90,000 and ∼ 267,000 in EAS and EUR, respectively). We observed that ancDE genes tended to be cell-type specific and enriched in genes interacting with the environment and in variants with ancestry-specific disease effect sizes, which suggests cell-type-specific, gene-by-environment interactions shared between regulatory and disease architectures. Finally, we illustrated how different environments might have led to ancestry-specific myeloid cell leukemia 1 (MCL1) expression in B cells and ancestry-specific allele effect sizes in lymphocyte count GWASs for variants surrounding MCL1. Our results imply that large single-cell and GWAS datasets from diverse ancestries are required to improve our understanding of human diseases.

Fine mapping of candidate effector genes for heart rate

An elevated resting heart rate (RHR) is associated with increased cardiovascular mortality. Genome-wide association studies (GWAS) have identified > 350 loci. Uniquely, in this study we applied genetic fine-mapping leveraging tissue specific chromatin segmentation and colocalization analyses to identify causal variants and candidate effector genes for RHR. We used RHR GWAS summary statistics from 388,237 individuals of European ancestry from UK Biobank and performed fine mapping using publicly available genomic annotation datasets. High-confidence causal variants (accounting for > 75% posterior probability) were identified, and we collated candidate effector genes using a multi-omics approach that combined evidence from colocalisation with molecular quantitative trait loci (QTLs), and long-range chromatin interaction analyses. Finally, we performed druggability analyses to investigate drug repurposing opportunities. The fine mapping pipeline indicated 442 distinct RHR signals. For 90 signals, a single variant was identified as a high-confidence causal variant, of which 22 were annotated as missense. In trait-relevant tissues, 39 signals colocalised with cis-expression QTLs (eQTLs), 3 with cis-protein QTLs (pQTLs), and 75 had promoter interactions via Hi-C. In total, 262 candidate genes were highlighted (79% had promoter interactions, 15% had a colocalised eQTL, 8% had a missense variant and 1% had a colocalised pQTL), and, for the first time, enrichment in nervous system pathways. Druggability analyses highlighted ACHE, CALCRL, MYT1 and TDP1 as potential targets. Our genetic fine-mapping pipeline prioritised 262 candidate genes for RHR that warrant further investigation in functional studies, and we provide potential therapeutic targets to reduce RHR and cardiovascular mortality.

Oral delivery of sodium alginate/chitosan bilayer microgels loaded with Lactobacillus rhamnosus GG for targeted therapy of ulcerative colitis

Probiotics regulate intestinal flora balance and enhance the intestinal barrier, which is useful in preventing and treating colitis. However, they have strict storage requirements. In addition, they degrade in a strongly acidic environment, resulting in a significant decrease in their activity when used as microbial agents. Lactobacillus rhamnosus GG (LGG) was loaded into acid-resistant and colon-targeting double-layer microgels. The inner layer consists of guar gum (GG) and low methoxyl pectin (LMP), which can achieve retention and degradation in the colon. To achieve colon localization, the outer layer was composed of chitosan (CS) and sodium alginate (SA). The formulation demonstrated favorable bio-responses across various pH conditions in vitro and sustained release of LGG in the colon lesions. Bare LGG survival decreased by 52.2 % in simulated gastric juice (pH 1.2) for 2 h, whereas that of encapsulated LGG decreased by 18.5 %. In the DSS-induced inflammatory model, LGG-loaded microgel significantly alleviated UC symptoms in mice and reduced inflammatory factor levels in the colon. Encapsulation of LGG improved its stability in acidic conditions, thus increasing its content at the colon lesions and reducing pathogenic bacteria. These findings provide an experimental basis and a technical reference for developing and applying probiotic microgel preparations.

An efficient, not-only-linear correlation coefficient based on clustering

Identifying meaningful patterns in data is crucial for understanding complex biological processes, particularly in transcriptomics, where genes with correlated expression often share functions or contribute to disease mechanisms. Traditional correlation coefficients, which primarily capture linear relationships, may overlook important nonlinear patterns. We introduce the clustermatch correlation coefficient (CCC), a not-only-linear coefficient that utilizes clustering to efficiently detect both linear and nonlinear associations. CCC outperforms standard methods by revealing biologically meaningful patterns that linear-only coefficients miss and is faster than state-of-the-art coefficients such as the maximal information coefficient. When applied to human gene expression data from genotype-tissue expression (GTEx), CCC identified robust linear relationships and nonlinear patterns, such as sex-specific differences, that are undetectable by standard methods. Highly ranked gene pairs were enriched for interactions in integrated networks built from protein-protein interactions, transcription factor regulation, and chemical and genetic perturbations, suggesting that CCC can detect functional relationships missed by linear-only approaches. CCC is a highly efficient, next-generation, not-only-linear correlation coefficient for genome-scale data. A record of this paper's transparent peer review process is included in the supplemental information.

Identifying genetic variants associated with chromatin looping and genome function

Here we present a comprehensive HiChIP dataset on naïve CD4 T cells (nCD4) from 30 donors and identify QTLs that associate with genotype-dependent and/or allele-specific variation of HiChIP contacts defining loops between active regulatory regions (iQTLs). We observe a substantial overlap between iQTLs and previously defined eQTLs and histone QTLs, and an enrichment for fine-mapped QTLs and GWAS variants. Furthermore, we describe a distinct subset of nCD4 iQTLs, for which the significant variation of chromatin contacts in nCD4 are translated into significant eQTL trends in CD4 T cell memory subsets. Finally, we define connectivity-QTLs as iQTLs that are significantly associated with concordant genotype-dependent changes in chromatin contacts over a broad genomic region (e.g., GWAS SNP in the RNASET2 locus). Our results demonstrate the importance of chromatin contacts as a complementary modality for QTL mapping and their power in identifying previously uncharacterized QTLs linked to cell-specific gene expression and connectivity.

Schizophrenia genomics: genetic complexity and functional insights

Determining the causes of schizophrenia has been a notoriously intractable problem, resistant to a multitude of investigative approaches over centuries. In recent decades, genomic studies have delivered hundreds of robust findings that implicate nearly 300 common genetic variants (via genome-wide association studies) and more than 20 rare variants (via whole-exome sequencing and copy number variant studies) as risk factors for schizophrenia. In parallel, functional genomic and neurobiological studies have provided exceptionally detailed information about the cellular composition of the brain and its interconnections in neurotypical individuals and, increasingly, in those with schizophrenia. Taken together, these results suggest unexpected complexity in the mechanisms that drive schizophrenia, pointing to the involvement of ensembles of genes (polygenicity) rather than single-gene causation. In this Review, we describe what we now know about the genetics of schizophrenia and consider the neurobiological implications of this information.

Recent developments and new directions in the use of natural products for the treatment of inflammatory bowel disease

BACKGROUND

Inflammatory bowel disease (IBD) represents a significant global health challenge, and there is an urgent need to explore novel therapeutic interventions. Natural products have demonstrated highly promising effectiveness in the treatment of IBD.

PURPOSE

This study systematically reviews the latest research advancements in leveraging natural products for IBD treatment.

METHODS

This manuscript strictly adheres to the PRISMA guidelines. Relevant literature on the effects of natural products on IBD was retrieved from the PubMed, Web of Science and Cochrane Library databases using the search terms "natural product," "inflammatory bowel disease," "colitis," "metagenomics", "target identification", "drug delivery systems", "polyphenols," "alkaloids," "terpenoids," and so on. The retrieved data were then systematically summarized and reviewed.

RESULTS

This review assessed the different effects of various natural products, such as polyphenols, alkaloids, terpenoids, quinones, and others, in the treatment of IBD. While these natural products offer promising avenues for IBD management, they also face challenges in terms of clinical translation and drug discovery. The advent of metagenomics, single-cell sequencing, target identification techniques, drug delivery systems, and other cutting-edge technologies heralds a new era in overcoming these challenges.

CONCLUSION

This paper provides an overview of current research progress in utilizing natural products for the treatment of IBD, exploring how contemporary technological innovations can aid in discovering and harnessing bioactive natural products for the treatment of IBD.

scParser: sparse representation learning for scalable single-cell RNA sequencing data analysis

The rapid rise in the availability and scale of scRNA-seq data needs scalable methods for integrative analysis. Though many methods for data integration have been developed, few focus on understanding the heterogeneous effects of biological conditions across different cell populations in integrative analysis. Our proposed scalable approach, scParser, models the heterogeneous effects from biological conditions, which unveils the key mechanisms by which gene expression contributes to phenotypes. Notably, the extended scParser pinpoints biological processes in cell subpopulations that contribute to disease pathogenesis. scParser achieves favorable performance in cell clustering compared to state-of-the-art methods and has a broad and diverse applicability.

A general principle governing neuronal evolution reveals a human-accelerated neuron type potentially underlying the high prevalence of autism in humans

The remarkable ability of a single genome sequence to encode a diverse collection of distinct cell types, including the thousands of cell types found in the mammalian brain, is a key characteristic of multicellular life. While it has been observed that some cell types are far more evolutionarily conserved than others, the factors driving these differences in evolutionary rate remain unknown. Here, we hypothesized that highly abundant neuronal cell types may be under greater selective constraint than rarer neuronal types, leading to variation in their rates of evolution. To test this, we leveraged recently published cross-species single-nucleus RNA-sequencing datasets from three distinct regions of the mammalian neocortex. We found a strikingly consistent relationship where more abundant neuronal subtypes show greater gene expression conservation between species, which replicated across three independent datasets covering >106 neurons from six species. Based on this principle, we discovered that the most abundant type of neocortical neurons-layer 2/3 intratelencephalic excitatory neurons-has evolved exceptionally quickly in the human lineage compared to other apes. Surprisingly, this accelerated evolution was accompanied by the dramatic down-regulation of autism-associated genes, which was likely driven by polygenic positive selection specific to the human lineage. In sum, we introduce a general principle governing neuronal evolution and suggest that the exceptionally high prevalence of autism in humans may be a direct result of natural selection for lower expression of a suite of genes that conferred a fitness benefit to our ancestors while also rendering an abundant class of neurons more sensitive to perturbation.

Enhanced In Situ Spatial Proteomics by Effective Combination of MALDI Imaging and LC-MS/MS

Highly specialized cells are fundamental for the proper functioning of complex organs. Variations in cell-type-specific gene expression and protein composition have been linked to a variety of diseases. Investigation of the distinctive molecular makeup of these cells within tissues is therefore critical in biomedical research. Although several technologies have emerged as valuable tools to address this cellular heterogeneity, most workflows lack sufficient in situ resolution and are associated with high costs and extremely long analysis times. Here, we present a combination of experimental and computational approaches that allows a more comprehensive investigation of molecular heterogeneity within tissues than by either shotgun LC-MS/MS or MALDI imaging alone. We applied our pipeline to the mouse brain, which contains a wide variety of cell types that not only perform unique functions but also exhibit varying sensitivities to insults. We explored the distinct neuronal populations within the hippocampus, a brain region crucial for learning and memory that is involved in various neurological disorders. As an example, we identified the groups of proteins distinguishing the neuronal populations of the dentate gyrus (DG) and the cornu ammonis (CA) in the same brain section. Most of the annotated proteins matched the regional enrichment of their transcripts, thereby validating the method. As the method is highly reproducible, the identification of individual masses through the combination of MALDI-IMS and LC-MS/MS methods can be used for the much faster and more precise interpretation of MALDI-IMS measurements only. This greatly speeds up spatial proteomic analyses and allows the detection of local protein variations within the same population of cells. The method's general applicability has the potential to be used to investigate different biological conditions and tissues and a much higher throughput than other techniques making it a promising approach for clinical routine applications.

Scalable genetic screening for regulatory circuits using compressed Perturb-seq

Pooled CRISPR screens with single-cell RNA sequencing readout (Perturb-seq) have emerged as a key technique in functional genomics, but they are limited in scale by cost and combinatorial complexity. In this study, we modified the design of Perturb-seq by incorporating algorithms applied to random, low-dimensional observations. Compressed Perturb-seq measures multiple random perturbations per cell or multiple cells per droplet and computationally decompresses these measurements by leveraging the sparse structure of regulatory circuits. Applied to 598 genes in the immune response to bacterial lipopolysaccharide, compressed Perturb-seq achieves the same accuracy as conventional Perturb-seq with an order of magnitude cost reduction and greater power to learn genetic interactions. We identified known and novel regulators of immune responses and uncovered evolutionarily constrained genes with downstream targets enriched for immune disease heritability, including many missed by existing genome-wide association studies. Our framework enables new scales of interrogation for a foundational method in functional genomics.

Genetic variation across and within individuals

Germline variation and somatic mutation are intricately connected and together shape human traits and disease risks. Germline variants are present from conception, but they vary between individuals and accumulate over generations. By contrast, somatic mutations accumulate throughout life in a mosaic manner within an individual due to intrinsic and extrinsic sources of mutations and selection pressures acting on cells. Recent advancements, such as improved detection methods and increased resources for association studies, have drastically expanded our ability to investigate germline and somatic genetic variation and compare underlying mutational processes. A better understanding of the similarities and differences in the types, rates and patterns of germline and somatic variants, as well as their interplay, will help elucidate the mechanisms underlying their distinct yet interlinked roles in human health and biology.

RNA Sequencing in Disease Diagnosis

RNA sequencing (RNA-seq) enables the accurate measurement of multiple transcriptomic phenotypes for modeling the impacts of disease variants. Advances in technologies, experimental protocols, and analysis strategies are rapidly expanding the application of RNA-seq to identify disease biomarkers, tissue- and cell-type-specific impacts, and the spatial localization of disease-associated mechanisms. Ongoing international efforts to construct biobank-scale transcriptomic repositories with matched genomic data across diverse population groups are further increasing the utility of RNA-seq approaches by providing large-scale normative reference resources. The availability of these resources, combined with improved computational analysis pipelines, has enabled the detection of aberrant transcriptomic phenotypes underlying rare diseases. Further expansion of these resources, across both somatic and developmental tissues, is expected to soon provide unprecedented insights to resolve disease origin, mechanism of action, and causal gene contributions, suggesting the continued high utility of RNA-seq in disease diagnosis.

Non-coding variants impact cis-regulatory coordination in a cell type-specific manner

BACKGROUND

Interactions among cis-regulatory elements (CREs) play a crucial role in gene regulation. Various approaches have been developed to map these interactions genome-wide, including those relying on interindividual epigenomic variation to identify groups of covariable regulatory elements, referred to as chromatin modules (CMs). While CM mapping allows to investigate the relationship between chromatin modularity and gene expression, the computational principles used for CM identification vary in their application and outcomes.

RESULTS

We comprehensively evaluate and streamline existing CM mapping tools and present guidelines for optimal utilization of epigenome data from a diverse population of individuals to assess regulatory coordination across the human genome. We showcase the effectiveness of our recommended practices by analyzing distinct cell types and demonstrate cell type specificity of CRE interactions in CMs and their relevance for gene expression. Integration of genotype information revealed that many non-coding disease-associated variants affect the activity of CMs in a cell type-specific manner by affecting the binding of cell type-specific transcription factors. We provide example cases that illustrate in detail how CMs can be used to deconstruct GWAS loci, assess variable expression of cell surface receptors in immune cells, and reveal how genetic variation can impact the expression of prognostic markers in chronic lymphocytic leukemia.

CONCLUSIONS

Our study presents an optimal strategy for CM mapping and reveals how CMs capture the coordination of CREs and its impact on gene expression. Non-coding genetic variants can disrupt this coordination, and we highlight how this may lead to disease predisposition in a cell type-specific manner.

Sequential immunotherapy: towards cures for autoimmunity

Despite major progress in the treatment of autoimmune diseases in the past two decades, most therapies do not cure disease and can be associated with increased risk of infection through broad suppression of the immune system. However, advances in understanding the causes of autoimmune disease and clinical data from novel therapeutic modalities such as chimeric antigen receptor T cell therapies provide evidence that it may be possible to re-establish immune homeostasis and, potentially, prolong remission or even cure autoimmune diseases. Here, we propose a 'sequential immunotherapy' framework for immune system modulation to help achieve this ambitious goal. This framework encompasses three steps: controlling inflammation; resetting the immune system through elimination of pathogenic immune memory cells; and promoting and maintaining immune homeostasis via immune regulatory agents and tissue repair. We discuss existing drugs and those in development for each of the three steps. We also highlight the importance of causal human biology in identifying and prioritizing novel immunotherapeutic strategies as well as informing their application in specific patient subsets, enabling precision medicine approaches that have the potential to transform clinical care.

A Novel Macrophage Subpopulation Conveys Increased Genetic Risk of Coronary Artery Disease

BACKGROUND

Coronary artery disease (CAD), the leading cause of death worldwide, is influenced by both environmental and genetic factors. Although over 250 genetic risk loci have been identified through genome-wide association studies, the specific causal variants and their regulatory mechanisms are still largely unknown, particularly in disease-relevant cell types such as macrophages.

METHODS

We utilized single-cell RNA-seq and single-cell multiomics approaches in primary human monocyte-derived macrophages to explore the transcriptional regulatory network involved in a critical pathogenic event of coronary atherosclerosis-the formation of lipid-laden foam cells. The relative genetic contribution to CAD was assessed by partitioning disease heritability across different macrophage subpopulations. Meta-analysis of single-cell RNA-seq data sets from 38 human atherosclerotic samples was conducted to provide high-resolution cross-referencing to macrophage subpopulations in vivo.

RESULTS

We identified 18 782 cis-regulatory elements by jointly profiling the gene expression and chromatin accessibility of >5000 macrophages. Integration with CAD genome-wide association study data prioritized 121 CAD-related genetic variants and 56 candidate causal genes. We showed that CAD heritability was not uniformly distributed and was particularly enriched in the gene programs of a novel CD52-hi lipid-handling macrophage subpopulation. These CD52-hi macrophages displayed significantly less lipoprotein accumulation and were also found in human atherosclerotic plaques. We investigated the cis-regulatory effect of a risk variant rs10488763 on FDX1, implicating the recruitment of AP-1 and C/EBP-β in the causal mechanisms at this locus.

CONCLUSIONS

Our results provide genetic evidence of the divergent roles of macrophage subsets in atherogenesis and highlight lipid-handling macrophages as a key subpopulation through which genetic variants operate to influence disease. These findings provide an unbiased framework for functional fine-mapping of genome-wide association study results using single-cell multiomics and offer new insights into the genotype-environment interactions underlying atherosclerotic disease.

Fine-mapping causal tissues and genes at disease-associated loci

Heritable diseases often manifest in a highly tissue-specific manner, with different disease loci mediated by genes in distinct tissues or cell types. We propose Tissue-Gene Fine-Mapping (TGFM), a fine-mapping method that infers the posterior probability (PIP) for each gene-tissue pair to mediate a disease locus by analyzing GWAS summary statistics (and in-sample LD) and leveraging eQTL data from diverse tissues to build cis-predicted expression models; TGFM also assigns PIPs to causal variants that are not mediated by gene expression in assayed genes and tissues. TGFM accounts for both co-regulation across genes and tissues and LD between SNPs (generalizing existing fine-mapping methods), and incorporates genome-wide estimates of each tissue's contribution to disease as tissue-level priors. TGFM was well-calibrated and moderately well-powered in simulations; unlike previous methods, TGFM was able to attain correct calibration by modeling uncertainty in cis-predicted expression models. We applied TGFM to 45 UK Biobank diseases/traits (average N = 316K) using eQTL data from 38 GTEx tissues. TGFM identified an average of 147 PIP > 0.5 causal genetic elements per disease/trait, of which 11% were gene-tissue pairs. Implicated gene-tissue pairs were concentrated in known disease-critical tissues, and causal genes were strongly enriched in disease-relevant gene sets. Causal gene-tissue pairs identified by TGFM recapitulated known biology (e.g., TPO-thyroid for Hypothyroidism), but also included biologically plausible novel findings (e.g., SLC20A2-artery aorta for Diastolic blood pressure). Further application of TGFM to single-cell eQTL data from 9 cell types in peripheral blood mononuclear cells (PBMC), analyzed jointly with GTEx tissues, identified 30 additional causal gene-PBMC cell type pairs at PIP > 0.5-primarily for autoimmune disease and blood cell traits, including the biologically plausible example of CD52 in classical monocyte cells for Monocyte count. In conclusion, TGFM is a robust and powerful method for fine-mapping causal tissues and genes at disease-associated loci.

Machine learning integrative approaches to advance computational immunology

The study of immunology, traditionally reliant on proteomics to evaluate individual immune cells, has been revolutionized by single-cell RNA sequencing. Computational immunologists play a crucial role in analysing these datasets, moving beyond traditional protein marker identification to encompass a more detailed view of cellular phenotypes and their functional roles. Recent technological advancements allow the simultaneous measurements of multiple cellular components-transcriptome, proteome, chromatin, epigenetic modifications and metabolites-within single cells, including in spatial contexts within tissues. This has led to the generation of complex multiscale datasets that can include multimodal measurements from the same cells or a mix of paired and unpaired modalities. Modern machine learning (ML) techniques allow for the integration of multiple "omics" data without the need for extensive independent modelling of each modality. This review focuses on recent advancements in ML integrative approaches applied to immunological studies. We highlight the importance of these methods in creating a unified representation of multiscale data collections, particularly for single-cell and spatial profiling technologies. Finally, we discuss the challenges of these holistic approaches and how they will be instrumental in the development of a common coordinate framework for multiscale studies, thereby accelerating research and enabling discoveries in the computational immunology field.

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.

Microfluidic device for the high-throughput and selective encapsulation of single target cells

Droplet-based microfluidic technologies for encapsulating single cells have rapidly evolved into powerful tools for single-cell analysis. In conventional passive single-cell encapsulation techniques, because cells arrive randomly at the droplet generation section, to encapsulate only a single cell with high precision, the average number of cells per droplet has to be decreased by reducing the average frequency at which cells arrive relative to the droplet generation rate. Therefore, the encapsulation efficiency for a given droplet generation rate is very low. Additionally, cell sorting operations are required prior to the encapsulation of target cells for specific cell type analysis. To address these challenges, we developed a cell encapsulation technology with a cell sorting function using a microfluidic chip. The microfluidic chip is equipped with an optical detection section to detect the optical information of cells and a sorting section to encapsulate cells into droplets by controlling a piezo element, enabling active encapsulation of only the single target cells. For a particle population including both targeted and non-targeted particles arriving at an average frequency of up to 6000 particles per s, with an average number of particles per droplet of 0.45, our device maintained a high purity above 97.9% for the single-target-particle droplets and achieved an outstanding throughput, encapsulating up to 2900 single target particles per s. The proposed encapsulation technology surpasses the encapsulation efficiency of conventional techniques, provides high efficiency and flexibility for single-cell research, and shows excellent potential for various applications in single-cell analysis.

Inferring causal cell types of human diseases and risk variants from candidate regulatory elements

The heritability of human diseases is extremely enriched in candidate regulatory elements (cRE) from disease-relevant cell types. Critical next steps are to infer which and how many cell types are truly causal for a disease (after accounting for co-regulation across cell types), and to understand how individual variants impact disease risk through single or multiple causal cell types. Here, we propose CT-FM and CT-FM-SNP, two methods that leverage cell-type-specific cREs to fine-map causal cell types for a trait and for its candidate causal variants, respectively. We applied CT-FM to 63 GWAS summary statistics (average N = 417K) using nearly one thousand cRE annotations, primarily coming from ENCODE4. CT-FM inferred 81 causal cell types with corresponding SNP-annotations explaining a high fraction of trait SNP-heritability (~2/3 of the SNP-heritability explained by existing cREs), identified 16 traits with multiple causal cell types, highlighted cell-disease relationships consistent with known biology, and uncovered previously unexplored cellular mechanisms in psychiatric and immune-related diseases. Finally, we applied CT-FM-SNP to 39 UK Biobank traits and predicted high confidence causal cell types for 2,798 candidate causal non-coding SNPs. Our results suggest that most SNPs impact a phenotype through a single cell type, and that pleiotropic SNPs target different cell types depending on the phenotype context. Altogether, CT-FM and CT-FM-SNP shed light on how genetic variants act collectively and individually at the cellular level to impact disease risk.

Functional genomics implicates natural killer cells in the pathogenesis of ankylosing spondylitis

Objective

Multiple lines of evidence indicate that ankylosing spondylitis (AS) is a lymphocyte-driven disease. However, which lymphocyte populations are critical in AS pathogenesis is not known. In this study, we aimed to identify the key cell types mediating the genetic risk in AS using an unbiased functional genomics approach.

Methods

We integrated genome-wide association study (GWAS) data with epigenomic and transcriptomic datasets of human immune cells. To quantify enrichment of cell type-specific open chromatin or gene expression in AS risk loci, we used three published methods that have successfully identified relevant cell types in other diseases. We performed co-localization analyses between GWAS risk loci and genetic variants associated with gene expression (eQTL) to find putative target genes.

Results

Natural killer (NK) cell-specific open chromatin regions are significantly enriched in heritability for AS, compared to other immune cell types such as T cells, B cells, and monocytes. This finding was consistent between two AS GWAS. Using RNA-seq data, we validated that genes in AS risk loci are enriched in NK cell-specific gene expression. Using the human Space-Time Gut Cell Atlas, we also found significant upregulation of AS-associated genes predominantly in NK cells. Co-localization analysis revealed four AS risk loci affecting regulation of candidate target genes in NK cells: two known loci, ERAP1 and TNFRSF1A, and two under-studied loci, ENTR1 (aka SDCCAG3) and B3GNT2.

Conclusion

Our findings suggest that NK cells may play a crucial role in AS development and highlight four putative target genes for functional follow-up in NK cells.

Disentangling associations between complex traits and cell types with seismic

Integrating single-cell RNA sequencing (scRNA-seq) with Genome-Wide Association Studies (GWAS) can help reveal GWAS-associated cell types, furthering our understanding of the cell-type-specific biological processes underlying complex traits and disease. However, current methods have technical limitations that hinder them from making systematic, scalable, interpretable disease-cell-type associations. In order to rapidly and accurately pinpoint associations, we develop a novel framework, seismic, which characterizes cell types using a new specificity score. We compare seismic with alternative methods across over 1,000 cell type characterizations at different granularities and 28 traits, demonstrating that seismic both corroborates findings and identifies trait-relevant cell groups which are not apparent through other methodologies. Furthermore, as part of the seismic framework, the specific genes driving cell type-trait associations can easily be accessed and analyzed, enabling further biological insights. The advantages of seismic are particularly salient in neurodegenerative diseases such as Parkinson's and Alzheimer's, where disease pathology has not only cell-specific manifestations, but also brain region-specific differences. Interestingly, a case study of Alzheimer's disease reveals the importance of considering GWAS endpoints, as studies relying on clinical diagnoses consistently identify microglial associations, while GWAS with a tau biomarker endpoint reveals neuronal associations. In general, seismic is a computationally efficient, powerful, and interpretable approach for identifying associations between complex traits and cell type-specific expression.

Landscape of human organoids: Ideal model in clinics and research

In the last decade, organoid research has entered a golden era, signifying a pivotal shift in the biomedical landscape. The year 2023 marked a milestone with the publication of thousands of papers in this arena, reflecting exponential growth. However, amid this burgeoning expansion, a comprehensive and accurate overview of the field has been conspicuously absent. Our review is intended to bridge this gap, providing a panoramic view of the rapidly evolving organoid landscape. We meticulously analyze the organoid field from eight distinctive vantage points, harnessing our rich experience in academic research, industrial application, and clinical practice. We present a deep exploration of the advances in organoid technology, underpinned by our long-standing involvement in this arena. Our narrative traverses the historical genesis of organoids and their transformative impact across various biomedical sectors, including oncology, toxicology, and drug development. We delve into the synergy between organoids and avant-garde technologies such as synthetic biology and single-cell omics and discuss their pivotal role in tailoring personalized medicine, enhancing high-throughput drug screening, and constructing physiologically pertinent disease models. Our comprehensive analysis and reflective discourse provide a deep dive into the existing landscape and emerging trends in organoid technology. We spotlight technological innovations, methodological evolution, and the broadening spectrum of applications, emphasizing the revolutionary influence of organoids in personalized medicine, oncology, drug discovery, and other fields. Looking ahead, we cautiously anticipate future developments in the field of organoid research, especially its potential implications for personalized patient care, new avenues of drug discovery, and clinical research. We trust that our comprehensive review will be an asset for researchers, clinicians, and patients with keen interest in personalized medical strategies. We offer a broad view of the present and prospective capabilities of organoid technology, encompassing a wide range of current and future applications. In summary, in this review we attempt a comprehensive exploration of the organoid field. We offer reflections, summaries, and projections that might be useful for current researchers and clinicians, and we hope to contribute to shaping the evolving trajectory of this dynamic and rapidly advancing field.

Unraveling the genetic architecture of blood unfolded p-53 among non-demented elderlies: novel candidate genes for early Alzheimer's disease

BACKGROUND

Alzheimer's disease (AD) is a heritable neurodegenerative disease whose long asymptomatic phase makes the early diagnosis of it pivotal. Blood U-p53 has recently emerged as a superior predictive biomarker for AD in the early stages. We hypothesized that genetic variants associated with blood U-p53 could reveal novel loci and pathways involved in the early stages of AD.

RESULTS

We performed a blood U-p53 Genome-wide association study (GWAS) on 484 healthy and mild cognitively impaired subjects from the ADNI cohort using 612,843 Single nucleotide polymorphisms (SNPs). We performed a pathway analysis and prioritized candidate genes using an AD single-cell gene program. We fine-mapped the intergenic SNPs by leveraging a cell-type-specific enhancer-to-gene linking strategy using a brain single-cell multimodal dataset. We validated the candidate genes in an independent brain single-cell RNA-seq and the ADNI blood transcriptome datasets. The rs279686 between AASS and FEZF1 genes was the most significant SNP (p-value = 4.82 × 10-7). Suggestive pathways were related to the immune and nervous systems. Twenty-three candidate genes were prioritized at 27 suggestive loci. Fine-mapping of 5 intergenic loci yielded nine cell-specific candidate genes. Finally, 15 genes were validated in the independent single-cell RNA-seq dataset, and five were validated in the ADNI blood transcriptome dataset.

CONCLUSIONS

We underlined the importance of performing a GWAS on an early-stage biomarker of AD and leveraging functional omics datasets for pinpointing causal genes in AD. Our study prioritized nine genes (SORCS1, KIF5C, TMEFF2, TMEM63C, HLA-E, ATAT1, TUBB, ARID1B, and RUNX1) strongly implicated in the early stages of AD.

Inferring gene regulatory networks from single-cell multiome data using atlas-scale external data

Existing methods for gene regulatory network (GRN) inference rely on gene expression data alone or on lower resolution bulk data. Despite the recent integration of chromatin accessibility and RNA sequencing data, learning complex mechanisms from limited independent data points still presents a daunting challenge. Here we present LINGER (Lifelong neural network for gene regulation), a machine-learning method to infer GRNs from single-cell paired gene expression and chromatin accessibility data. LINGER incorporates atlas-scale external bulk data across diverse cellular contexts and prior knowledge of transcription factor motifs as a manifold regularization. LINGER achieves a fourfold to sevenfold relative increase in accuracy over existing methods and reveals a complex regulatory landscape of genome-wide association studies, enabling enhanced interpretation of disease-associated variants and genes. Following the GRN inference from reference single-cell multiome data, LINGER enables the estimation of transcription factor activity solely from bulk or single-cell gene expression data, leveraging the abundance of available gene expression data to identify driver regulators from case-control studies.

Integrated single-cell and bulk RNA sequencing reveals CREM is involved in the pathogenesis of ulcerative colitis

Background

Ulcerative colitis (UC) is an inflammatory bowel disease characterized by persistent colonic inflammation. Here, we performed a systematic analysis to gain better insights into UC pathogenesis.

Methods

We analyzed two UC-related datasets extracted from the gene expression omnibus database using several bioinformatics tools. The primary cell types and key subgroups of primary cells associated with UC and differentially expressed genes (DEGs) between UC and control samples were identified. The molecular regulation of the key genes was also predicted. The gene ontology and Kyoto encyclopedia of genes and genomes enrichment analyses of marker genes of key cell subgroups and model genes were performed. The expression of key enriched genes was validated in 10 clinical samples using real-time quantitative polymerase chain reaction (RT-qPCR).

Results

Monocytes were identified as the major cell type. Ten differentially expressed marker genes were obtained by intersecting the 3121 DEGs, 38 marker genes in major cell types, and 104 marker genes in key cell subgroups. Four essential genes, associated with immune response, were obtained using support vector machine recursive feature elimination and least absolute shrinkage and selection operator analyses. The four essential genes were highly expressed in Cluster 0 during differentiation. Validation of the four key genes in colonic mucosal biopsy specimens from 10 normal and 10 UC patients revealed that CREM was highly expressed in both the lesion-free sites and lesion sites colonic mucosa of UC patients compared with normal adults.

Conclusions

We identified CREM involved in UC pathogenesis, which is expected to provide a new therapeutic target for UC.

Genome-wide association study identifies 30 obsessive-compulsive disorder associated loci

Obsessive-compulsive disorder (OCD) affects ~1% of the population and exhibits a high SNP-heritability, yet previous genome-wide association studies (GWAS) have provided limited information on the genetic etiology and underlying biological mechanisms of the disorder. We conducted a GWAS meta-analysis combining 53,660 OCD cases and 2,044,417 controls from 28 European-ancestry cohorts revealing 30 independent genome-wide significant SNPs and a SNP-based heritability of 6.7%. Separate GWAS for clinical, biobank, comorbid, and self-report sub-groups found no evidence of sample ascertainment impacting our results. Functional and positional QTL gene-based approaches identified 249 significant candidate risk genes for OCD, of which 25 were identified as putatively causal, highlighting WDR6, DALRD3, CTNND1 and genes in the MHC region. Tissue and single-cell enrichment analyses highlighted hippocampal and cortical excitatory neurons, along with D1- and D2-type dopamine receptor-containing medium spiny neurons, as playing a role in OCD risk. OCD displayed significant genetic correlations with 65 out of 112 examined phenotypes. Notably, it showed positive genetic correlations with all included psychiatric phenotypes, in particular anxiety, depression, anorexia nervosa, and Tourette syndrome, and negative correlations with a subset of the included autoimmune disorders, educational attainment, and body mass index.. This study marks a significant step toward unraveling its genetic landscape and advances understanding of OCD genetics, providing a foundation for future interventions to address this debilitating disorder.

The cells of the sensory epithelium, and not the stria vascularis, are the main cochlear cells related to the genetic pathogenesis of age-related hearing loss

Age-related hearing loss (ARHL) is a major health concern among the elderly population. It is hoped that increasing our understanding of its underlying pathophysiological processes will lead to the development of novel therapies. Recent genome-wide association studies (GWASs) discovered a few dozen genetic variants in association with elevated risk for ARHL. Integrated analysis of GWAS results and transcriptomics data is a powerful approach for elucidating specific cell types that are involved in disease pathogenesis. Intriguingly, recent studies that applied such bioinformatics approaches to ARHL resulted in disagreeing findings as for the key cell types that are most strongly linked to the genetic pathogenesis of ARHL. These conflicting studies pointed either to cochlear sensory epithelial or to stria vascularis cells as the cell types most prominently involved in the genetic basis of ARHL. Seeking to resolve this discrepancy, we integrated the analysis of four ARHL GWAS datasets with four independent inner-ear single-cell RNA-sequencing datasets. Our analysis clearly points to the cochlear sensory epithelial cells as the key cells for the genetic predisposition to ARHL. We also explain the limitation of the bioinformatics analysis performed by previous studies that led to missing the enrichment for ARHL GWAS signal in sensory epithelial cells. Collectively, we show that cochlear epithelial cells, not stria vascularis cells, are the main inner-ear cells related to the genetic pathogenesis of ARHL.

Current insights and future perspectives of ultraviolet radiation (UV) exposure: Friends and foes to the skin and beyond the skin

Ultraviolet (UV) radiation is ubiquitous in the environment, which has been classified as an established human carcinogen. As the largest and outermost organ of the body, direct exposure of skin to sunlight or UV radiation can result in sunburn, inflammation, photo-immunosuppression, photoaging and even skin cancers. To date, there are tactics to protect the skin by preventing UV radiation and reducing the amount of UV radiation to the skin. Nevertheless, deciphering the essential regulatory mechanisms may pave the way for therapeutic interventions against UV-induced skin disorders. Additionally, UV light is considered beneficial for specific skin-related conditions in medical UV therapy. Recent evidence indicates that the biological effects of UV exposure extend beyond the skin and include the treatment of inflammatory diseases, solid tumors and certain abnormal behaviors. This review mainly focuses on the effects of UV on the skin. Moreover, novel findings of the biological effects of UV in other organs and systems are also summarized. Nevertheless, the mechanisms through which UV affects the human organism remain to be fully elucidated to achieve a more comprehensive understanding of its biological effects.

High-resolution omics of vascular ageing and inflammatory pathways in neurodegeneration

High-resolution omics, particularly single-cell and spatial transcriptomic profiling, are rapidly enhancing our comprehension of the normal molecular diversity of gliovascular cells, as well as their age-related changes that contribute to neurodegeneration. With more omic profiling studies being conducted, it is becoming increasingly essential to synthesise valuable information from the rapidly accumulating findings. In this review, we present an overview of the molecular features of neurovascular and glial cells that have been recently discovered through omic profiling, with a focus on those that have potentially significant functional implications and/or show cross-species differences between human and mouse, and that are linked to vascular deficits and inflammatory pathways in ageing and neurodegenerative disorders. Additionally, we highlight the translational applications of omic profiling, and discuss omic-based strategies to accelerate biomarker discovery and facilitate disease course-modifying therapeutics development for neurodegenerative conditions.

Integration of human organoids single-cell transcriptomic profiles and human genetics repurposes critical cell type-specific drug targets for severe COVID-19

Human organoids recapitulate the cell type diversity and function of their primary organs holding tremendous potentials for basic and translational research. Advances in single-cell RNA sequencing (scRNA-seq) technology and genome-wide association study (GWAS) have accelerated the biological and therapeutic interpretation of trait-relevant cell types or states. Here, we constructed a computational framework to integrate atlas-level organoid scRNA-seq data, GWAS summary statistics, expression quantitative trait loci, and gene-drug interaction data for distinguishing critical cell populations and drug targets relevant to coronavirus disease 2019 (COVID-19) severity. We found that 39 cell types across eight kinds of organoids were significantly associated with COVID-19 outcomes. Notably, subset of lung mesenchymal stem cells increased proximity with fibroblasts predisposed to repair COVID-19-damaged lung tissue. Brain endothelial cell subset exhibited significant associations with severe COVID-19, and this cell subset showed a notable increase in cell-to-cell interactions with other brain cell types, including microglia. We repurposed 33 druggable genes, including IFNAR2, TYK2, and VIPR2, and their interacting drugs for COVID-19 in a cell-type-specific manner. Overall, our results showcase that host genetic determinants have cellular-specific contribution to COVID-19 severity, and identification of cell type-specific drug targets may facilitate to develop effective therapeutics for treating severe COVID-19 and its complications.

Single-cell transcriptome landscape of circulating CD4+ T cell populations in autoimmune diseases

CD4+ T cells are key mediators of various autoimmune diseases; however, their role in disease progression remains unclear due to cellular heterogeneity. Here, we evaluated CD4+ T cell subpopulations using decomposition-based transcriptome characterization and canonical clustering strategies. This approach identified 12 independent gene programs governing whole CD4+ T cell heterogeneity, which can explain the ambiguity of canonical clustering. In addition, we performed a meta-analysis using public single-cell datasets of over 1.8 million peripheral CD4+ T cells from 953 individuals by projecting cells onto the reference and cataloging cell frequency and qualitative alterations of the populations in 20 diseases. The analyses revealed that the 12 transcriptional programs were useful in characterizing each autoimmune disease and predicting its clinical status. Moreover, genetic variants associated with autoimmune diseases showed disease-specific enrichment within the 12 gene programs. The results collectively provide a landscape of single-cell transcriptomes of CD4+ T cell subpopulations involved in autoimmune disease.

The genetic basis of autoimmunity seen through the lens of T cell functional traits

Autoimmune disease heritability is enriched in T cell-specific regulatory regions of the genome. Modern-day T cell datasets now enable association studies between single nucleotide polymorphisms (SNPs) and a myriad of molecular phenotypes, including chromatin accessibility, gene expression, transcriptional programs, T cell antigen receptor (TCR) amino acid usage, and cell state abundances. Such studies have identified hundreds of quantitative trait loci (QTLs) in T cells that colocalize with genetic risk for autoimmune disease. The key challenge facing immunologists today lies in synthesizing these results toward a unified understanding of the autoimmune T cell: which genes, cell states, and antigens drive tissue destruction?

Integrating leiomyoma genetics, epigenomics, and single-cell transcriptomics reveals causal genetic variants, genes, and cell types

Uterine fibroids (UF), that can disrupt normal uterine function and cause significant physical and psychological health problems, are observed in nearly 70% of women of reproductive age. Although heritable genetics is a significant risk factor, specific genetic variations and gene targets causally associated with UF are poorly understood. Here, we performed a meta-analysis on existing fibroid genome-wide association studies (GWAS) and integrated the identified risk loci and potentially causal single nucleotide polymorphisms (SNPs) with epigenomics, transcriptomics, 3D chromatin organization from diverse cell types as well as primary UF patient's samples. This integrative analysis identifies 24 UF-associated risk loci that potentially target 394 genes, of which 168 are differentially expressed in UF tumors. Critically, integrating this data with single-cell gene expression data from UF patients reveales the causal cell types with aberrant expression of these target genes. Lastly, CRISPR-based epigenetic repression (dCas9-KRAB) or activation (dCas9-p300) in a UF disease-relevant cell type further refines and narrows down the potential gene targets. Our findings and the methodological approach indicate the effectiveness of integrating multi-omics data with locus-specific epigenetic editing approaches for identifying gene- and celt type-targets of disease-relevant risk loci.

Predicting cell-type specific disease genes of diabetes with the biological network

Type 2 diabetes (T2D) is a chronic condition that can lead to significant harm, such as heart disease, kidney disease, nerve damage, and blindness. Although T2D-related genes have been identified through Genome-wide association studies (GWAS) and various computational methods, the biological mechanism of T2D at the cell type level remains unclear. Exploring cell type-specific genes related to T2D is essential to understand the cellular mechanisms underlying the disease. To address this issue, we introduce DiGCellNet (predicting Disease Genes with Cell type specificity based on biological Networks), a model that integrates graph convolutional network (GCN) and multi-task learning (MTL) to predict T2D-associated cell type-specific genes based on the biological network. Our work represents the first attempt to predict cell type-specific disease genes using GCN and MTL. We evaluate our approach by predicting genes specific to four cell types and demonstrate that the proposed DiGCellNet outperforms other models that combine node embeddings with traditional machine learning algorithms. Moreover, DiGCellNet successfully identifies CALM1 as a gene specific to beta cell type in T2D cases, and this association is confirmed using an independent dataset. The code is available at https://github.com/23AIBox/23AIBox-DiGCellNet.

Convergence of coronary artery disease genes onto endothelial cell programs

Linking variants from genome-wide association studies (GWAS) to underlying mechanisms of disease remains a challenge1-3. For some diseases, a successful strategy has been to look for cases in which multiple GWAS loci contain genes that act in the same biological pathway1-6. However, our knowledge of which genes act in which pathways is incomplete, particularly for cell-type-specific pathways or understudied genes. Here we introduce a method to connect GWAS variants to functions. This method links variants to genes using epigenomics data, links genes to pathways de novo using Perturb-seq and integrates these data to identify convergence of GWAS loci onto pathways. We apply this approach to study the role of endothelial cells in genetic risk for coronary artery disease (CAD), and discover 43 CAD GWAS signals that converge on the cerebral cavernous malformation (CCM) signalling pathway. Two regulators of this pathway, CCM2 and TLNRD1, are each linked to a CAD risk variant, regulate other CAD risk genes and affect atheroprotective processes in endothelial cells. These results suggest a model whereby CAD risk is driven in part by the convergence of causal genes onto a particular transcriptional pathway in endothelial cells. They highlight shared genes between common and rare vascular diseases (CAD and CCM), and identify TLNRD1 as a new, previously uncharacterized member of the CCM signalling pathway. This approach will be widely useful for linking variants to functions for other common polygenic diseases.

SURGE: uncovering context-specific genetic-regulation of gene expression from single-cell RNA sequencing using latent-factor models

Genetic regulation of gene expression is a complex process, with genetic effects known to vary across cellular contexts such as cell types and environmental conditions. We developed SURGE, a method for unsupervised discovery of context-specific expression quantitative trait loci (eQTLs) from single-cell transcriptomic data. This allows discovery of the contexts or cell types modulating genetic regulation without prior knowledge. Applied to peripheral blood single-cell eQTL data, SURGE contexts capture continuous representations of distinct cell types and groupings of biologically related cell types. We demonstrate the disease-relevance of SURGE context-specific eQTLs using colocalization analysis and stratified LD-score regression.

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.

Leveraging single-cell ATAC-seq and RNA-seq to identify disease-critical fetal and adult brain cell types

Prioritizing disease-critical cell types by integrating genome-wide association studies (GWAS) with functional data is a fundamental goal. Single-cell chromatin accessibility (scATAC-seq) and gene expression (scRNA-seq) have characterized cell types at high resolution, and studies integrating GWAS with scRNA-seq have shown promise, but studies integrating GWAS with scATAC-seq have been limited. Here, we identify disease-critical fetal and adult brain cell types by integrating GWAS summary statistics from 28 brain-related diseases/traits (average N = 298 K) with 3.2 million scATAC-seq and scRNA-seq profiles from 83 cell types. We identified disease-critical fetal (respectively adult) brain cell types for 22 (respectively 23) of 28 traits using scATAC-seq, and for 8 (respectively 17) of 28 traits using scRNA-seq. Significant scATAC-seq enrichments included fetal photoreceptor cells for major depressive disorder, fetal ganglion cells for BMI, fetal astrocytes for ADHD, and adult VGLUT2 excitatory neurons for schizophrenia. Our findings improve our understanding of brain-related diseases/traits and inform future analyses.

Affected cell types for hundreds of Mendelian diseases revealed by analysis of human and mouse single-cell data

Mendelian diseases tend to manifest clinically in certain tissues, yet their affected cell types typically remain elusive. Single-cell expression studies showed that overexpression of disease-associated genes may point to the affected cell types. Here, we developed a method that infers disease-affected cell types from the preferential expression of disease-associated genes in cell types (PrEDiCT). We applied PrEDiCT to single-cell expression data of six human tissues, to infer the cell types affected in Mendelian diseases. Overall, we inferred the likely affected cell types for 328 diseases. We corroborated our findings by literature text-mining, expert validation, and recapitulation in mouse corresponding tissues. Based on these findings, we explored characteristics of disease-affected cell types, showed that diseases manifesting in multiple tissues tend to affect similar cell types, and highlighted cases where gene functions could be used to refine inference. Together, these findings expand the molecular understanding of disease mechanisms and cellular vulnerability.

Integrating genetic regulation and single-cell expression with GWAS prioritizes causal genes and cell types for glaucoma

Primary open-angle glaucoma (POAG), characterized by retinal ganglion cell death, is a leading cause of irreversible blindness worldwide. However, its molecular and cellular causes are not well understood. Elevated intraocular pressure (IOP) is a major risk factor, but many patients have normal IOP. Colocalization and Mendelian randomization analysis of >240 POAG and IOP genome-wide association study (GWAS) loci and overlapping expression and splicing quantitative trait loci (e/sQTLs) in 49 GTEx tissues and retina prioritizes causal genes for 60% of loci. These genes are enriched in pathways implicated in extracellular matrix organization, cell adhesion, and vascular development. Analysis of single-nucleus RNA-seq of glaucoma-relevant eye tissues reveals that the POAG and IOP colocalizing genes and genome-wide associations are enriched in specific cell types in the aqueous outflow pathways, retina, optic nerve head, peripapillary sclera, and choroid. This study nominates IOP-dependent and independent regulatory mechanisms, genes, and cell types that may contribute to POAG pathogenesis.

Imaging and AI based chromatin biomarkers for diagnosis and therapy evaluation from liquid biopsies

Multiple genomic and proteomic studies have suggested that peripheral blood mononuclear cells (PBMCs) respond to tumor secretomes and thus could provide possible avenues for tumor prognosis and treatment evaluation. We hypothesized that the chromatin organization of PBMCs obtained from liquid biopsies, which integrates secretome signals with gene expression programs, provides efficient biomarkers to characterize tumor signals and the efficacy of proton therapy in tumor patients. Here, we show that chromatin imaging of PBMCs combined with machine learning methods provides such robust and predictive chromatin biomarkers. We show that such chromatin biomarkers enable the classification of 10 healthy and 10 pan-tumor patients. Furthermore, we extended our pipeline to assess the tumor types and states of 30 tumor patients undergoing (proton) radiation therapy. We show that our pipeline can thereby accurately distinguish between three tumor groups with up to 89% accuracy and enables the monitoring of the treatment effects. Collectively, we show the potential of chromatin biomarkers for cancer diagnostics and therapy evaluation.

Large-scale single-cell RNA sequencing atlases of human immune cells across lifespan: Possibilities and challenges

Single-cell RNA sequencing technologies have successfully been leveraged for immunological insights into human prenatal, pediatric, and adult tissues. These single-cell studies have led to breakthroughs in our understanding of stem, myeloid, and lymphoid cell function. Computational analysis of fetal hematopoietic tissues has uncovered trajectories for T- and B-cell differentiation across multiple organ sites, and how these trajectories might be dysregulated in fetal and pediatric health and disease. As we enter the age of large-scale, multiomic, and integrative single-cell meta-analysis, we assess the advances and challenges of large-scale data generation, analysis, and reanalysis, and data dissemination for a broad range of scientific and clinical communities. We discuss Findable, Accessible, Interoperable, and Reusable data sharing and unified cell ontology languages as strategic areas for progress of the field in the near future. We also reflect on the trend toward deployment of multiomic and spatial genomic platforms within single-cell RNA sequencing projects, and the crucial role these data types will assume in the immediate future toward creation of comprehensive and rich single-cell atlases. We demonstrate using our recent studies of human prenatal and adult hematopoietic tissues the importance of interdisciplinary and collaborative working in science to reveal biological insights in parallel with technological and computational innovations.

Dynamic regulatory elements in single-cell multimodal data implicate key immune cell states enriched for autoimmune disease heritability

In autoimmune diseases such as rheumatoid arthritis, the immune system attacks the body's own cells. Developing a precise understanding of the cell states where noncoding autoimmune risk variants impart causal mechanisms is critical to developing curative therapies. Here, to identify noncoding regions with accessible chromatin that associate with cell-state-defining gene expression patterns, we leveraged multimodal single-nucleus RNA and assay for transposase-accessible chromatin (ATAC) sequencing data across 28,674 cells from the inflamed synovial tissue of 12 donors. Specifically, we used a multivariate Poisson model to predict peak accessibility from single-nucleus RNA sequencing principal components. For 14 autoimmune diseases, we discovered that cell-state-dependent ('dynamic') chromatin accessibility peaks in immune cell types were enriched for heritability, compared with cell-state-invariant ('cs-invariant') peaks. These dynamic peaks marked regulatory elements associated with T peripheral helper, regulatory T, dendritic and STAT1+CXCL10+ myeloid cell states. We argue that dynamic regulatory elements can help identify precise cell states enriched for disease-critical genetic variation.