Single-cell RNA-sequencing of differentiating iPS cells reveals dynamic genetic effects on gene expression

xQTLatlas: a comprehensive resource for human cellular-resolution multi-omics genetic regulatory landscape

Understanding how genetic variants influence molecular phenotypes in different cellular contexts is crucial for elucidating the molecular and cellular mechanisms behind complex traits, which in turn has spurred significant advances in research into molecular quantitative trait locus (xQTL) at the cellular level. With the rapid proliferation of data, there is a critical need for a comprehensive and accessible platform to integrate this information. To meet this need, we developed xQTLatlas (http://www.hitxqtl.org.cn/), a database that provides a multi-omics genetic regulatory landscape at cellular resolution. xQTLatlas compiles xQTL summary statistics from 151 cell types and 339 cell states across 55 human tissues. It organizes these data into 20 xQTL types, based on four distinct discovery strategies, and spans 13 molecular phenotypes. Each entry in xQTLatlas is meticulously annotated with comprehensive metadata, including the origin of the tissue, cell type, cell state and the QTL discovery strategies utilized. Additionally, xQTLatlas features multiscale data exploration tools and a suite of interactive visualizations, facilitating in-depth analysis of cell-level xQTL. xQTLatlas provides a valuable resource for deepening our understanding of the impact of functional variants on molecular phenotypes in different cellular environments, thereby facilitating extensive research efforts.

Individual variation in the emergence of anterior-to-posterior neural fates from human pluripotent stem cells

Variability between human pluripotent stem cell (hPSC) lines remains a challenge and opportunity in biomedicine. In this study, hPSC lines from multiple donors were differentiated toward neuroectoderm and mesendoderm lineages. We revealed dynamic transcriptomic patterns that delineate the emergence of these lineages, which were conserved across lines, along with individual line-specific transcriptional signatures that were invariant throughout differentiation. These transcriptomic signatures predicted an antagonism between SOX21-driven forebrain fates and retinoic acid-induced hindbrain fates. Replicate lines and paired adult tissue demonstrated the stability of these line-specific transcriptomic traits. We show that this transcriptomic variation in lineage bias had both genetic and epigenetic origins, aligned with the anterior-to-posterior structure of early mammalian development, and was present across a large collection of hPSC lines. These findings contribute to developing systematic analyses of PSCs to define the origin and consequences of variation in the early events orchestrating individual human development.

Single-cell eQTL mapping in yeast reveals a tradeoff between growth and reproduction

Expression quantitative trait loci (eQTLs) provide a key bridge between noncoding DNA sequence variants and organismal traits. The effects of eQTLs can differ among tissues, cell types, and cellular states, but these differences are obscured by gene expression measurements in bulk populations. We developed a one-pot approach to map eQTLs in Saccharomyces cerevisiae by single-cell RNA sequencing (scRNA-seq) and applied it to over 100,000 single cells from three crosses. We used scRNA-seq data to genotype each cell, measure gene expression, and classify the cells by cell-cycle stage. We mapped thousands of local and distant eQTLs and identified interactions between eQTL effects and cell-cycle stages. We took advantage of single-cell expression information to identify hundreds of genes with allele-specific effects on expression noise. We used cell-cycle stage classification to map 20 loci that influence cell-cycle progression. One of these loci influenced the expression of genes involved in the mating response. We showed that the effects of this locus arise from a common variant (W82R) in the gene GPA1 , which encodes a signaling protein that negatively regulates the mating pathway. The 82R allele increases mating efficiency at the cost of slower cell-cycle progression and is associated with a higher rate of outcrossing in nature. Our results provide a more granular picture of the effects of genetic variants on gene expression and downstream traits.

Computational methods for allele-specific expression in single cells

Allele-specific expression (ASE) is a powerful signal that can be used to investigate multiple molecular mechanisms, such as cis-regulatory effects and imprinting. Single-cell RNA-sequencing (scRNA-seq) enables ASE characterization at the resolution of individual cells. In this review, we highlight the computational methods for processing and analyzing single-cell ASE data. We first describe a bioinformatics pipeline to obtain ASE counts from raw reads synthesized from previous literature. We then discuss statistical methods for detecting allelic imbalance and its variability across conditions using scRNA-seq data. In addition, we describe other methods that use single-cell ASE to address specific biological questions. Finally, we discuss future directions and emphasize the need for an integrated, optimized bioinformatics pipeline, and further development of statistical methods for different technologies.

Population Diversity at the Single-Cell Level

Population-scale single-cell genomics is a transformative approach for unraveling the intricate links between genetic and cellular variation. This approach is facilitated by cutting-edge experimental methodologies, including the development of high-throughput single-cell multiomics and advances in multiplexed environmental and genetic perturbations. Examining the effects of natural or synthetic genetic variants across cellular contexts provides insights into the mutual influence of genetics and the environment in shaping cellular heterogeneity. The development of computational methodologies further enables detailed quantitative analysis of molecular variation, offering an opportunity to examine the respective roles of stochastic, intercellular, and interindividual variation. Future opportunities lie in leveraging long-read sequencing, refining disease-relevant cellular models, and embracing predictive and generative machine learning models. These advancements hold the potential for a deeper understanding of the genetic architecture of human molecular traits, which in turn has important implications for understanding the genetic causes of human disease.

From genetic associations to genes: methods, applications, and challenges

Genome-wide association studies (GWASs) have identified numerous genetic loci associated with human traits and diseases. However, pinpointing the causal genes remains a challenge, which impedes the translation of GWAS findings into biological insights and medical applications. In this review, we provide an in-depth overview of the methods and technologies used for prioritizing genes from GWAS loci, including gene-based association tests, integrative analysis of GWAS and molecular quantitative trait loci (xQTL) data, linking GWAS variants to target genes through enhancer-gene connection maps, and network-based prioritization. We also outline strategies for generating context-dependent xQTL data and their applications in gene prioritization. We further highlight the potential of gene prioritization in drug repurposing. Lastly, we discuss future challenges and opportunities in this field.

Deciphering the heterogeneity of differentiating hPSC-derived corneal limbal stem cells through single-cell RNA sequencing

A comprehensive understanding of the human pluripotent stem cell (hPSC) differentiation process stands as a prerequisite for the development of hPSC-based therapeutics. In this study, single-cell RNA sequencing (scRNA-seq) was performed to decipher the heterogeneity during differentiation of three hPSC lines toward corneal limbal stem cells (LSCs). The scRNA-seq data revealed nine clusters encompassing the entire differentiation process, among which five followed the anticipated differentiation path of LSCs. The remaining four clusters were previously undescribed cell states that were annotated as either mesodermal-like or undifferentiated subpopulations, and their prevalence was hPSC line dependent. Distinct cluster-specific marker genes identified in this study were confirmed by immunofluorescence analysis and employed to purify hPSC-derived LSCs, which effectively minimized the variation in the line-dependent differentiation efficiency. In summary, scRNA-seq offered molecular insights into the heterogeneity of hPSC-LSC differentiation, allowing a data-driven strategy for consistent and robust generation of LSCs, essential for future advancement toward clinical translation.

Brain Chimeroids reveal individual susceptibility to neurotoxic triggers

Interindividual genetic variation affects the susceptibility to and progression of many diseases1,2. However, efforts to study how individual human brains differ in normal development and disease phenotypes are limited by the paucity of faithful cellular human models, and the difficulty of scaling current systems to represent multiple people. Here we present human brain Chimeroids, a highly reproducible, multidonor human brain cortical organoid model generated by the co-development of cells from a panel of individual donors in a single organoid. By reaggregating cells from multiple single-donor organoids at the neural stem cell or neural progenitor cell stage, we generate Chimeroids in which each donor produces all cell lineages of the cerebral cortex, even when using pluripotent stem cell lines with notable growth biases. We used Chimeroids to investigate interindividual variation in the susceptibility to neurotoxic triggers that exhibit high clinical phenotypic variability: ethanol and the antiepileptic drug valproic acid. Individual donors varied in both the penetrance of the effect on target cell types, and the molecular phenotype within each affected cell type. Our results suggest that human genetic background may be an important mediator of neurotoxin susceptibility and introduce Chimeroids as a scalable system for high-throughput investigation of interindividual variation in processes of brain development and disease.

LMNA -Related Dilated Cardiomyopathy: Single-Cell Transcriptomics during Patient-derived iPSC Differentiation Support Cell type and Lineage-specific Dysregulation of Gene Expression and Development for Cardiomyocytes and Epicardium-Derived Cells with Lamin A/C Haploinsufficiency

LMNA -Related Dilated Cardiomyopathy (DCM) is an autosomal-dominant genetic condition with cardiomyocyte and conduction system dysfunction often resulting in heart failure or sudden death. The condition is caused by mutation in the Lamin A/C ( LMNA ) gene encoding Type-A nuclear lamin proteins involved in nuclear integrity, epigenetic regulation of gene expression, and differentiation. Molecular mechanisms of disease are not completely understood, and there are no definitive treatments to reverse progression or prevent mortality. We investigated possible mechanisms of LMNA -Related DCM using induced pluripotent stem cells derived from a family with a heterozygous LMNA c.357-2A>G splice-site mutation. We differentiated one LMNA mutant iPSC line derived from an affected female (Patient) and two non-mutant iPSC lines derived from her unaffected sister (Control) and conducted single-cell RNA sequencing for 12 samples (4 Patient and 8 Control) across seven time points: Day 0, 2, 4, 9, 16, 19, and 30. Our bioinformatics workflow identified 125,554 cells in raw data and 110,521 (88%) high-quality cells in sequentially processed data. Unsupervised clustering, cell annotation, and trajectory inference found complex heterogeneity: ten main cell types; many possible subtypes; and lineage bifurcation for Cardiac Progenitors to Cardiomyocytes (CM) and Epicardium-Derived Cells (EPDC). Data integration and comparative analyses of Patient and Control cells found cell type and lineage differentially expressed genes (DEG) with enrichment to support pathway dysregulation. Top DEG and enriched pathways included: 10 ZNF genes and RNA polymerase II transcription in Pluripotent cells (PP); BMP4 and TGF Beta/BMP signaling, sarcomere gene subsets and cardiogenesis, CDH2 and EMT in CM; LMNA and epigenetic regulation and DDIT4 and mTORC1 signaling in EPDC. Top DEG also included: XIST and other X-linked genes, six imprinted genes: SNRPN , PWAR6 , NDN , PEG10 , MEG3 , MEG8 , and enriched gene sets in metabolism, proliferation, and homeostasis. We confirmed Lamin A/C haploinsufficiency by allelic expression and Western blot. Our complex Patient-derived iPSC model for Lamin A/C haploinsufficiency in PP, CM, and EPDC provided support for dysregulation of genes and pathways, many previously associated with Lamin A/C defects, such as epigenetic gene expression, signaling, and differentiation. Our findings support disruption of epigenomic developmental programs as proposed in other LMNA disease models. We recognized other factors influencing epigenetics and differentiation; thus, our approach needs improvement to further investigate this mechanism in an iPSC-derived model.

A PLURIPOTENT STEM CELL PLATFORM FOR IN VITRO SYSTEMS GENETICS STUDIES OF MOUSE DEVELOPMENT

The directed differentiation of pluripotent stem cells (PSCs) from panels of genetically diverse individuals is emerging as a powerful experimental system for characterizing the impact of natural genetic variation on developing cell types and tissues. Here, we establish new PSC lines and experimental approaches for modeling embryonic development in a genetically diverse, outbred mouse stock (Diversity Outbred mice). We show that a range of inbred and outbred PSC lines can be stably maintained in the primed pluripotent state (epiblast stem cells -- EpiSCs) and establish the contribution of genetic variation to phenotypic differences in gene regulation and directed differentiation. Using pooled in vitro fertilization, we generate and characterize a genetic reference panel of Diversity Outbred PSCs (n = 230). Finally, we demonstrate the feasibility of pooled culture of Diversity Outbred EpiSCs as "cell villages", which can facilitate the differentiation of large numbers of EpiSC lines for forward genetic screens. These data can complement and inform similar efforts within the stem cell biology and human genetics communities to model the impact of natural genetic variation on phenotypic variation and disease-risk.

Recapitulation of Perturbed Striatal Gene Expression Dynamics of Donors' Brains With Ventral Forebrain Organoids Derived From the Same Individuals With Schizophrenia

OBJECTIVE

Schizophrenia is a brain disorder that originates during neurodevelopment and has complex genetic and environmental etiologies. Despite decades of clinical evidence of altered striatal function in affected patients, studies examining its cellular and molecular mechanisms in humans are limited. To explore neurodevelopmental alterations in the striatum associated with schizophrenia, the authors established a method for the differentiation of induced pluripotent stem cells (iPSCs) into ventral forebrain organoids (VFOs).

METHODS

VFOs were generated from postmortem dural fibroblast-derived iPSCs of four individuals with schizophrenia and four neurotypical control individuals for whom postmortem caudate genotypes and transcriptomic data were profiled in the BrainSeq neurogenomics consortium. Individuals were selected such that the two groups had nonoverlapping schizophrenia polygenic risk scores (PRSs).

RESULTS

Single-cell RNA sequencing analyses of VFOs revealed differences in developmental trajectory between schizophrenia and control individuals in which inhibitory neuronal cells from the patients exhibited accelerated maturation. Furthermore, upregulated genes in inhibitory neurons in schizophrenia VFOs showed a significant overlap with upregulated genes in postmortem caudate tissue of individuals with schizophrenia compared with control individuals, including the donors of the iPSC cohort.

CONCLUSIONS

The findings suggest that striatal neurons derived from high-PRS individuals with schizophrenia carry abnormalities that originated during early brain development and that the VFO model can recapitulate disease-relevant cell type-specific neurodevelopmental phenotypes in a dish.

Integrating single cell expression quantitative trait loci summary statistics to understand complex trait risk genes

Transcriptome-wide association study (TWAS) is a popular approach to dissect the functional consequence of disease associated non-coding variants. Most existing TWAS use bulk tissues and may not have the resolution to reveal cell-type specific target genes. Single-cell expression quantitative trait loci (sc-eQTL) datasets are emerging. The largest bulk- and sc-eQTL datasets are most conveniently available as summary statistics, but have not been broadly utilized in TWAS. Here, we present a new method EXPRESSO (EXpression PREdiction with Summary Statistics Only), to analyze sc-eQTL summary statistics, which also integrates 3D genomic data and epigenomic annotation to prioritize causal variants. EXPRESSO substantially improves existing methods. We apply EXPRESSO to analyze multi-ancestry GWAS datasets for 14 autoimmune diseases. EXPRESSO uniquely identifies 958 novel gene x trait associations, which is 26% more than the second-best method. Among them, 492 are unique to cell type level analysis and missed by TWAS using whole blood. We also develop a cell type aware drug repurposing pipeline, which leverages EXPRESSO results to identify drug compounds that can reverse disease gene expressions in relevant cell types. Our results point to multiple drugs with therapeutic potentials, including metformin for type 1 diabetes, and vitamin K for ulcerative colitis.

Efficient and accurate mixed model association tool for single-cell eQTL analysis

Understanding the genetic basis of gene expression can help us understand the molecular underpinnings of human traits and disease. Expression quantitative trait locus (eQTL) mapping can help in studying this relationship but have been shown to be very cell-type specific, motivating the use of single-cell RNA sequencing and single-cell eQTLs to obtain a more granular view of genetic regulation. Current methods for single-cell eQTL mapping either rely on the "pseudobulk" approach and traditional pipelines for bulk transcriptomics or do not scale well to large datasets. Here, we propose SAIGE-QTL, a robust and scalable tool that can directly map eQTLs using single-cell profiles without needing aggregation at the pseudobulk level. Additionally, SAIGE-QTL allows for testing the effects of less frequent/rare genetic variation through set-based tests, which is traditionally excluded from eQTL mapping studies. We evaluate the performance of SAIGE-QTL on both real and simulated data and demonstrate the improved power for eQTL mapping over existing pipelines.

Cell-type and dynamic state govern genetic regulation of gene expression in heterogeneous differentiating cultures

Identifying the molecular effects of human genetic variation across cellular contexts is crucial for understanding the mechanisms underlying disease-associated loci, yet many cell-types and developmental stages remain underexplored. Here we harnessed the potential of heterogeneous differentiating cultures ( HDCs ), an in vitro system in which pluripotent cells asynchronously differentiate into a broad spectrum of cell-types. We generated HDCs for 53 human donors and collected single-cell RNA-sequencing data from over 900,000 cells. We identified expression quantitative trait loci in 29 cell-types and characterized regulatory dynamics across diverse differentiation trajectories. This revealed novel regulatory variants for genes involved in key developmental and disease-related processes while replicating known effects from primary tissues, and dynamic regulatory effects associated with a range of complex traits.

Single-nucleotide variant calling in single-cell sequencing data with Monopogen

Single-cell omics technologies enable molecular characterization of diverse cell types and states, but how the resulting transcriptional and epigenetic profiles depend on the cell's genetic background remains understudied. We describe Monopogen, a computational tool to detect single-nucleotide variants (SNVs) from single-cell sequencing data. Monopogen leverages linkage disequilibrium from external reference panels to identify germline SNVs and detects putative somatic SNVs using allele cosegregating patterns at the cell population level. It can identify 100 K to 3 M germline SNVs achieving a genotyping accuracy of 95%, together with hundreds of putative somatic SNVs. Monopogen-derived genotypes enable global and local ancestry inference and identification of admixed samples. It identifies variants associated with cardiomyocyte metabolic levels and epigenomic programs. It also improves putative somatic SNV detection that enables clonal lineage tracing in primary human clonal hematopoiesis. Monopogen brings together population genetics, cell lineage tracing and single-cell omics to uncover genetic determinants of cellular processes.

Integrating population genetics, stem cell biology and cellular genomics to study complex human diseases

Human pluripotent stem (hPS) cells can, in theory, be differentiated into any cell type, making them a powerful in vitro model for human biology. Recent technological advances have facilitated large-scale hPS cell studies that allow investigation of the genetic regulation of molecular phenotypes and their contribution to high-order phenotypes such as human disease. Integrating hPS cells with single-cell sequencing makes identifying context-dependent genetic effects during cell development or upon experimental manipulation possible. Here we discuss how the intersection of stem cell biology, population genetics and cellular genomics can help resolve the functional consequences of human genetic variation. We examine the critical challenges of integrating these fields and approaches to scaling them cost-effectively and practically. We highlight two areas of human biology that can particularly benefit from population-scale hPS cell studies, elucidating mechanisms underlying complex disease risk loci and evaluating relationships between common genetic variation and pharmacotherapeutic phenotypes.

Investigating the role of common cis-regulatory variants in modifying penetrance of putatively damaging, inherited variants in severe neurodevelopmental disorders

Recent work has revealed an important role for rare, incompletely penetrant inherited coding variants in neurodevelopmental disorders (NDDs). Additionally, we have previously shown that common variants contribute to risk for rare NDDs. Here, we investigate whether common variants exert their effects by modifying gene expression, using multi-cis-expression quantitative trait loci (cis-eQTL) prediction models. We first performed a transcriptome-wide association study for NDDs using 6987 probands from the Deciphering Developmental Disorders (DDD) study and 9720 controls, and found one gene, RAB2A, that passed multiple testing correction (p = 6.7 × 10-7). We then investigated whether cis-eQTLs modify the penetrance of putatively damaging, rare coding variants inherited by NDD probands from their unaffected parents in a set of 1700 trios. We found no evidence that unaffected parents transmitting putatively damaging coding variants had higher genetically-predicted expression of the variant-harboring gene than their child. In probands carrying putatively damaging variants in constrained genes, the genetically-predicted expression of these genes in blood was lower than in controls (p = 2.7 × 10-3). However, results for proband-control comparisons were inconsistent across different sets of genes, variant filters and tissues. We find limited evidence that common cis-eQTLs modify penetrance of rare coding variants in a large cohort of NDD probands.

Dissecting the biology of feeding and eating disorders

Feeding and eating disorders (FEDs) are heterogenous and characterized by varying patterns of dysregulated eating and weight. Genome-wide association studies (GWASs) are clarifying their underlying biology and their genetic relationship to other psychiatric and metabolic/anthropometric traits. Genetic research on anorexia nervosa (AN) has identified eight significant loci and uncovered genetic correlations implicating both psychiatric and metabolic/anthropometric risk factors. Careful explication of these metabolic contributors may be key to developing effective and enduring treatments for devastating, life-altering, and frequently lethal illnesses. We discuss clinical phenomenology, genomics, phenomics, intestinal microbiota, and functional genomics and propose a path that translates variants to genes, genes to pathways, and pathways to metabolic outcomes to advance the science and eventually treatment of FEDs.

Tissue-specific enhancer-gene maps from multimodal single-cell data identify causal disease alleles

Translating genome-wide association study (GWAS) loci into causal variants and genes requires accurate cell-type-specific enhancer-gene maps from disease-relevant tissues. Building enhancer-gene maps is essential but challenging with current experimental methods in primary human tissues. Here we developed a nonparametric statistical method, SCENT (single-cell enhancer target gene mapping), that models association between enhancer chromatin accessibility and gene expression in single-cell or nucleus multimodal RNA sequencing and ATAC sequencing data. We applied SCENT to 9 multimodal datasets including >120,000 single cells or nuclei and created 23 cell-type-specific enhancer-gene maps. These maps were highly enriched for causal variants in expression quantitative loci and GWAS for 1,143 diseases and traits. We identified likely causal genes for both common and rare diseases and linked somatic mutation hotspots to target genes. We demonstrate that application of SCENT to multimodal data from disease-relevant human tissue enables the scalable construction of accurate cell-type-specific enhancer-gene maps, essential for defining noncoding variant function.

A Multimodal Omics Framework to Empower Target Discovery for Cardiovascular Regeneration

Ischaemic heart disease is a global healthcare challenge with high morbidity and mortality. Early revascularisation in acute myocardial infarction has improved survival; however, limited regenerative capacity and microvascular dysfunction often lead to impaired function and the development of heart failure. New mechanistic insights are required to identify robust targets for the development of novel strategies to promote regeneration. Single-cell RNA sequencing (scRNA-seq) has enabled profiling and analysis of the transcriptomes of individual cells at high resolution. Applications of scRNA-seq have generated single-cell atlases for multiple species, revealed distinct cellular compositions for different regions of the heart, and defined multiple mechanisms involved in myocardial injury-induced regeneration. In this review, we summarise findings from studies of healthy and injured hearts in multiple species and spanning different developmental stages. Based on this transformative technology, we propose a multi-species, multi-omics, meta-analysis framework to drive the discovery of new targets to promote cardiovascular regeneration.

Machine learning nominates the inositol pathway and novel genes in Parkinson's disease

There are 78 loci associated with Parkinson's disease in the most recent genome-wide association study (GWAS), yet the specific genes driving these associations are mostly unknown. Herein, we aimed to nominate the top candidate gene from each Parkinson's disease locus and identify variants and pathways potentially involved in Parkinson's disease. We trained a machine learning model to predict Parkinson's disease-associated genes from GWAS loci using genomic, transcriptomic and epigenomic data from brain tissues and dopaminergic neurons. We nominated candidate genes in each locus and identified novel pathways potentially involved in Parkinson's disease, such as the inositol phosphate biosynthetic pathway (INPP5F, IP6K2, ITPKB and PPIP5K2). Specific common coding variants in SPNS1 and MLX may be involved in Parkinson's disease, and burden tests of rare variants further support that CNIP3, LSM7, NUCKS1 and the polyol/inositol phosphate biosynthetic pathway are associated with the disease. Functional studies are needed to further analyse the involvements of these genes and pathways in Parkinson's disease.

The dynamic genetic determinants of increased transcriptional divergence in spermatids

Cis-genetic effects are key determinants of transcriptional divergence in discrete tissues and cell types. However, how cis- and trans-effects act across continuous trajectories of cellular differentiation in vivo is poorly understood. Here, we quantify allele-specific expression during spermatogenic differentiation at single-cell resolution in an F1 hybrid mouse system, allowing for the comprehensive characterisation of cis- and trans-genetic effects, including their dynamics across cellular differentiation. Collectively, almost half of the genes subject to genetic regulation show evidence for dynamic cis-effects that vary during differentiation. Our system also allows us to robustly identify dynamic trans-effects, which are less pervasive than cis-effects. In aggregate, genetic effects were strongest in round spermatids, which parallels their increased transcriptional divergence we identified between species. Our approach provides a comprehensive quantification of the variability of genetic effects in vivo, and demonstrates a widely applicable strategy to dissect the impact of regulatory variants on gene regulation in dynamic systems.

Characterization of caffeine response regulatory variants in vascular endothelial cells

Genetic variants in gene regulatory sequences can modify gene expression and mediate the molecular response to environmental stimuli. In addition, genotype-environment interactions (GxE) contribute to complex traits such as cardiovascular disease. Caffeine is the most widely consumed stimulant and is known to produce a vascular response. To investigate GxE for caffeine, we treated vascular endothelial cells with caffeine and used a massively parallel reporter assay to measure allelic effects on gene regulation for over 43,000 genetic variants. We identified 665 variants with allelic effects on gene regulation and 6 variants that regulate the gene expression response to caffeine (GxE, false discovery rate [FDR] < 5%). When overlapping our GxE results with expression quantitative trait loci colocalized with coronary artery disease and hypertension, we dissected their regulatory mechanisms and showed a modulatory role for caffeine. Our results demonstrate that massively parallel reporter assay is a powerful approach to identify and molecularly characterize GxE in the specific context of caffeine consumption.

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.

scQTLbase: an integrated human single-cell eQTL database

Genome-wide association studies (GWAS) have identified numerous genetic variants associated with diseases and traits. However, the functional interpretation of these variants remains challenging. Expression quantitative trait loci (eQTLs) have been widely used to identify mutations linked to disease, yet they explain only 20-50% of disease-related variants. Single-cell eQTLs (sc-eQTLs) studies provide an immense opportunity to identify new disease risk genes with expanded eQTL scales and transcriptional regulation at a much finer resolution. However, there is no comprehensive database dedicated to single-cell eQTLs that users can use to search, analyse and visualize them. Therefore, we developed the scQTLbase (http://bioinfo.szbl.ac.cn/scQTLbase), the first integrated human sc-eQTLs portal, featuring 304 datasets spanning 57 cell types and 95 cell states. It contains ∼16 million SNPs significantly associated with cell-type/state gene expression and ∼0.69 million disease-associated sc-eQTLs from 3 333 traits/diseases. In addition, scQTLbase offers sc-eQTL search, gene expression visualization in UMAP plots, a genome browser, and colocalization visualization based on the GWAS dataset of interest. scQTLbase provides a one-stop portal for sc-eQTLs that will significantly advance the discovery of disease susceptibility genes.

Opportunities and tradeoffs in single-cell transcriptomic technologies

Recent technological and algorithmic advances enable single-cell transcriptomic analysis with remarkable depth and breadth. Nonetheless, a persistent challenge is the compromise between the ability to profile high numbers of cells and the achievement of full-length transcript coverage. Currently, the field is progressing and developing new and creative solutions that improve cellular throughput, gene detection sensitivity and full-length transcript capture. Furthermore, long-read sequencing approaches for single-cell transcripts are breaking frontiers that have previously blocked full transcriptome characterization. We here present a comprehensive overview of available options for single-cell transcriptome profiling, highlighting the key advantages and disadvantages of each approach.

Common genetic variation impacts stress response in the brain

To explain why individuals exposed to identical stressors experience divergent clinical outcomes, we determine how molecular encoding of stress modifies genetic risk for brain disorders. Analysis of post-mortem brain (n=304) revealed 8557 stress-interactive expression quantitative trait loci (eQTLs) that dysregulate expression of 915 eGenes in response to stress, and lie in stress-related transcription factor binding sites. Response to stress is robust across experimental paradigms: up to 50% of stress-interactive eGenes validate in glucocorticoid treated hiPSC-derived neurons (n=39 donors). Stress-interactive eGenes show brain region- and cell type-specificity, and, in post-mortem brain, implicate glial and endothelial mechanisms. Stress dysregulates long-term expression of disorder risk genes in a genotype-dependent manner; stress-interactive transcriptomic imputation uncovered 139 novel genes conferring brain disorder risk only in the context of traumatic stress. Molecular stress-encoding explains individualized responses to traumatic stress; incorporating trauma into genomic studies of brain disorders is likely to improve diagnosis, prognosis, and drug discovery.

Mapping the landscape of lineage-specific dynamic regulation of gene expression using single-cell transcriptomics and application to genetics of complex disease

Single-cell transcriptome data can provide insights into how genetic variation influences biological processes involved in human biology and disease. However, the identification of gene-level associations in distinct cell types faces several challenges, including the limited reference resource from population scale studies, data sparsity in single-cell RNA sequencing, and the complex cell-state pattern of expression within individual cell types. Here we develop genetic models of cell type specific and cell state adjusted gene expression in mid-brain neurons in the process of specializing from induced pluripotent stem cells. The resulting framework quantifies the dynamics of the genetic regulation of gene expression and estimates its cell type specificity. As an application, we show that the approach detects known and new genes associated with schizophrenia and enables insights into context-dependent disease mechanisms. We provide a genomic resource from a phenome-wide application of our models to more than 1500 phenotypes from the UK Biobank. Using longitudinal genetically determined expression, we implement a predictive causality framework, evaluating the prediction of future values of a target gene expression using prior values of a putative regulatory gene. Collectively, this work demonstrates the insights that can be gained into the molecular underpinnings of diseases by quantifying the genetic control of gene expression at single-cell resolution.

Identifying biomarkers of heterogeneity and transplantation efficacy in retinal pigment epithelial cells

Transplantation of retinal pigment epithelial (RPE) cells holds great promise for patients with retinal degenerative diseases, such as age-related macular degeneration. In-depth characterization of RPE cell product identity and critical quality attributes are needed to enhance efficacy and safety of replacement therapy strategies. Here, we characterized an adult RPE stem cell-derived (RPESC-RPE) cell product using bulk and single-cell RNA sequencing (scRNA-seq), assessing functional cell integration in vitro into a mature RPE monolayer and in vivo efficacy by vision rescue in the Royal College of Surgeons rats. scRNA-seq revealed several distinct subpopulations in the RPESC-RPE product, some with progenitor markers. We identified RPE clusters expressing genes associated with in vivo efficacy and increased cell integration capability. Gene expression analysis revealed lncRNA (TREX) as a predictive marker of in vivo efficacy. TREX knockdown decreased cell integration while overexpression increased integration in vitro and improved vision rescue in the RCS rats.

Induced Pluripotent Stem Cells in Disease Biology and the Evidence for Their In Vitro Utility

Many human phenotypes are impossible to recapitulate in model organisms or immortalized human cell lines. Induced pluripotent stem cells (iPSCs) offer a way to study disease mechanisms in a variety of differentiated cell types while circumventing ethical and practical issues associated with finite tissue sources and postmortem states. Here, we discuss the broad utility of iPSCs in genetic medicine and describe how they are being used to study musculoskeletal, pulmonary, neurologic, and cardiac phenotypes. We summarize the particular challenges presented by each organ system and describe how iPSC models are being used to address them. Finally, we discuss emerging iPSC-derived organoid models and the potential value that they can bring to studies of human disease.

Leveraging Single-Cell Populations to Uncover the Genetic Basis of Complex Traits

The ease and throughput of single-cell genomics have steadily improved, and its current trajectory suggests that surveying single-cell populations will become routine. We discuss the merger of quantitative genetics with single-cell genomics and emphasize how this synergizes with advantages intrinsic to plants. Single-cell population genomics provides increased detection resolution when mapping variants that control molecular traits, including gene expression or chromatin accessibility. Additionally, single-cell population genomics reveals the cell types in which variants act and, when combined with organism-level phenotype measurements, unveils which cellular contexts impact higher-order traits. Emerging technologies, notably multiomics, can facilitate the measurement of both genetic changes and genomic traits in single cells, enabling single-cell genetic experiments. The implementation of single-cell genetics will advance the investigation of the genetic architecture of complex molecular traits and provide new experimental paradigms to study eukaryotic genetics.

Single-cell allele-specific expression analysis reveals dynamic and cell-type-specific regulatory effects

Differential allele-specific expression (ASE) is a powerful tool to study context-specific cis-regulation of gene expression. Such effects can reflect the interaction between genetic or epigenetic factors and a measured context or condition. Single-cell RNA sequencing (scRNA-seq) allows the measurement of ASE at individual-cell resolution, but there is a lack of statistical methods to analyze such data. We present Differential Allelic Expression using Single-Cell data (DAESC), a powerful method for differential ASE analysis using scRNA-seq from multiple individuals, with statistical behavior confirmed through simulation. DAESC accounts for non-independence between cells from the same individual and incorporates implicit haplotype phasing. Application to data from 105 induced pluripotent stem cell (iPSC) lines identifies 657 genes dynamically regulated during endoderm differentiation, with enrichment for changes in chromatin state. Application to a type-2 diabetes dataset identifies several differentially regulated genes between patients and controls in pancreatic endocrine cells. DAESC is a powerful method for single-cell ASE analysis and can uncover novel insights on gene regulation.

MuDCoD: multi-subject community detection in personalized dynamic gene networks from single-cell RNA sequencing

MOTIVATION

With the wide availability of single-cell RNA-seq (scRNA-seq) technology, population-scale scRNA-seq datasets across multiple individuals and time points are emerging. While the initial investigations of these datasets tend to focus on standard analysis of clustering and differential expression, leveraging the power of scRNA-seq data at the personalized dynamic gene co-expression network level has the potential to unlock subject and/or time-specific network-level variation, which is critical for understanding phenotypic differences. Community detection from co-expression networks of multiple time points or conditions has been well-studied; however, none of the existing settings included networks from multiple subjects and multiple time points simultaneously. To address this, we develop Multi-subject Dynamic Community Detection (MuDCoD) for multi-subject community detection in personalized dynamic gene networks from scRNA-seq. MuDCoD builds on the spectral clustering framework and promotes information sharing among the networks of the subjects as well as networks at different time points. It clusters genes in the personalized dynamic gene networks and reveals gene communities that are variable or shared not only across time but also among subjects.

RESULTS

Evaluation and benchmarking of MuDCoD against existing approaches reveal that MuDCoD effectively leverages apparent shared signals among networks of the subjects at individual time points, and performs robustly when there is no or little information sharing among the networks. Applications to population-scale scRNA-seq datasets of human-induced pluripotent stem cells during dopaminergic neuron differentiation and CD4+ T cell activation indicate that MuDCoD enables robust inference for identifying time-varying personalized gene modules. Our results illustrate how personalized dynamic community detection can aid in the exploration of subject-specific biological processes that vary across time.

AVAILABILITY AND IMPLEMENTATION

MuDCoD is publicly available at https://github.com/bo1929/MuDCoD as a Python package. Implementation includes simulation and real-data experiments together with extensive documentation.

Principles and challenges of modeling temporal and spatial omics data

Studies with temporal or spatial resolution are crucial to understand the molecular dynamics and spatial dependencies underlying a biological process or system. With advances in high-throughput omic technologies, time- and space-resolved molecular measurements at scale are increasingly accessible, providing new opportunities to study the role of timing or structure in a wide range of biological questions. At the same time, analyses of the data being generated in the context of spatiotemporal studies entail new challenges that need to be considered, including the need to account for temporal and spatial dependencies and compare them across different scales, biological samples or conditions. In this Review, we provide an overview of common principles and challenges in the analysis of temporal and spatial omics data. We discuss statistical concepts to model temporal and spatial dependencies and highlight opportunities for adapting existing analysis methods to data with temporal and spatial dimensions.

Reliable multiplex generation of pooled induced pluripotent stem cells

Reprogramming somatic cells into pluripotent stem cells (iPSCs) enables the study of systems in vitro. To increase the throughput of reprogramming, we present induction of pluripotency from pooled cells (iPPC)-an efficient, scalable, and reliable reprogramming procedure. Using our deconvolution algorithm that employs pooled sequencing of single-nucleotide polymorphisms (SNPs), we accurately estimated individual donor proportions of the pooled iPSCs. With iPPC, we concurrently reprogrammed over one hundred donor lymphoblastoid cell lines (LCLs) into iPSCs and found strong correlations of individual donors' reprogramming ability across multiple experiments. Individual donors' reprogramming ability remains consistent across both same-day replicates and multiple experimental runs, and the expression of certain immunoglobulin precursor genes may impact reprogramming ability. The pooled iPSCs were also able to differentiate into cerebral organoids. Our procedure enables a multiplex framework of using pooled libraries of donor iPSCs for downstream research and investigation of in vitro phenotypes.

Corneal regeneration strategies: From stem cell therapy to tissue engineered stem cell scaffolds

Corneal epithelial defects and excessive wound healing might lead to severe complications. As stem cells can self-renew infinitely, they are a promising solution for regenerating the corneal epithelium and treating severe corneal epithelial injury. The chemical and biophysical properties of biological scaffolds, such as the amniotic membrane, fibrin, and hydrogels, can provide the necessary signals for stem cell proliferation and differentiation. Multiple researchers have conducted investigations on these scaffolds and evaluated them as potential therapeutic interventions for corneal disorders. These studies have identified various inherent benefits and drawbacks associated with these scaffolds. In this study, we provided a comprehensive overview of the history and use of various stem cells in corneal repair. We mainly discussed biological scaffolds that are used in stem cell transplantation and innovative materials that are under investigation.

Methods and Insights from Single-Cell Expression Quantitative Trait Loci

Recent advancements in single-cell technologies have enabled expression quantitative trait locus (eQTL) analysis across many individuals at single-cell resolution. Compared with bulk RNA sequencing, which averages gene expression across cell types and cell states, single-cell assays capture the transcriptional states of individual cells, including fine-grained, transient, and difficult-to-isolate populations at unprecedented scale and resolution. Single-cell eQTL (sc-eQTL) mapping can identify context-dependent eQTLs that vary with cell states, including some that colocalize with disease variants identified in genome-wide association studies. By uncovering the precise contexts in which these eQTLs act, single-cell approaches can unveil previously hidden regulatory effects and pinpoint important cell states underlying molecular mechanisms of disease. Here, we present an overview of recently deployed experimental designs in sc-eQTL studies. In the process, we consider the influence of study design choices such as cohort, cell states, and ex vivo perturbations. We then discuss current methodologies, modeling approaches, and technical challenges as well as future opportunities and applications.

Omics data integration in computational biology viewed through the prism of machine learning paradigms

Important quantities of biological data can today be acquired to characterize cell types and states, from various sources and using a wide diversity of methods, providing scientists with more and more information to answer challenging biological questions. Unfortunately, working with this amount of data comes at the price of ever-increasing data complexity. This is caused by the multiplication of data types and batch effects, which hinders the joint usage of all available data within common analyses. Data integration describes a set of tasks geared towards embedding several datasets of different origins or modalities into a joint representation that can then be used to carry out downstream analyses. In the last decade, dozens of methods have been proposed to tackle the different facets of the data integration problem, relying on various paradigms. This review introduces the most common data types encountered in computational biology and provides systematic definitions of the data integration problems. We then present how machine learning innovations were leveraged to build effective data integration algorithms, that are widely used today by computational biologists. We discuss the current state of data integration and important pitfalls to consider when working with data integration tools. We eventually detail a set of challenges the field will have to overcome in the coming years.

Methods and applications for single-cell and spatial multi-omics

The joint analysis of the genome, epigenome, transcriptome, proteome and/or metabolome from single cells is transforming our understanding of cell biology in health and disease. In less than a decade, the field has seen tremendous technological revolutions that enable crucial new insights into the interplay between intracellular and intercellular molecular mechanisms that govern development, physiology and pathogenesis. In this Review, we highlight advances in the fast-developing field of single-cell and spatial multi-omics technologies (also known as multimodal omics approaches), and the computational strategies needed to integrate information across these molecular layers. We demonstrate their impact on fundamental cell biology and translational research, discuss current challenges and provide an outlook to the future.

Single-cell genomics meets human genetics

Single-cell genomic technologies are revealing the cellular composition, identities and states in tissues at unprecedented resolution. They have now scaled to the point that it is possible to query samples at the population level, across thousands of individuals. Combining single-cell information with genotype data at this scale provides opportunities to link genetic variation to the cellular processes underpinning key aspects of human biology and disease. This strategy has potential implications for disease diagnosis, risk prediction and development of therapeutic solutions. But, effectively integrating large-scale single-cell genomic data, genetic variation and additional phenotypic data will require advances in data generation and analysis methods. As single-cell genetics begins to emerge as a field in its own right, we review its current state and the challenges and opportunities ahead.

Functional genomics in stem cell models: considerations and applications

Protocols to differentiate human pluripotent stem cells have advanced in terms of cell type specificity and tissue-level complexity over the past 2 decades, which has facilitated human disease modeling in the most relevant cell types. The ability to generate induced PSCs (iPSCs) from patients further enables the study of disease mutations in an appropriate cellular context to reveal the mechanisms that underlie disease etiology and progression. As iPSC-derived disease models have improved in robustness and scale, they have also been adopted more widely for use in drug screens to discover new therapies and therapeutic targets. Advancement in genome editing technologies, in particular the discovery of CRISPR-Cas9, has further allowed for rapid development of iPSCs containing disease-causing mutations. CRISPR-Cas9 technologies have now evolved beyond creating single gene edits, aided by the fusion of inhibitory (CRISPRi) or activation (CRISPRa) domains to a catalytically dead Cas9 protein, enabling inhibition or activation of endogenous gene loci. These tools have been used in CRISPR knockout, CRISPRi, or CRISPRa screens to identify genetic modifiers that synergize or antagonize with disease mutations in a systematic and unbiased manner, resulting in identification of disease mechanisms and discovery of new therapeutic targets to accelerate drug discovery research. However, many technical challenges remain when applying large-scale functional genomics approaches to differentiated PSC populations. Here we review current technologies in the field of iPSC disease modeling and CRISPR-based functional genomics screens and practical considerations for implementation across a range of modalities, applications, and disease areas, as well as explore CRISPR screens that have been performed in iPSC models to-date and the insights and therapies these screens have produced.

Label-free imaging of 3D pluripotent stem cell differentiation dynamics on chip

Massive, parallelized 3D stem cell cultures for engineering in vitro human cell types require imaging methods with high time and spatial resolution to fully exploit technological advances in cell culture technologies. Here, we introduce a large-scale integrated microfluidic chip platform for automated 3D stem cell differentiation. To fully enable dynamic high-content imaging on the chip platform, we developed a label-free deep learning method called Bright2Nuc to predict in silico nuclear staining in 3D from confocal microscopy bright-field images. Bright2Nuc was trained and applied to hundreds of 3D human induced pluripotent stem cell cultures differentiating toward definitive endoderm on a microfluidic platform. Combined with existing image analysis tools, Bright2Nuc segmented individual nuclei from bright-field images, quantified their morphological properties, predicted stem cell differentiation state, and tracked the cells over time. Our methods are available in an open-source pipeline, enabling researchers to upscale image acquisition and phenotyping of 3D cell culture.

A statistical framework to identify cell types whose genetically regulated proportions are associated with complex diseases

Finding disease-relevant tissues and cell types can facilitate the identification and investigation of functional genes and variants. In particular, cell type proportions can serve as potential disease predictive biomarkers. In this manuscript, we introduce a novel statistical framework, cell-type Wide Association Study (cWAS), that integrates genetic data with transcriptomics data to identify cell types whose genetically regulated proportions (GRPs) are disease/trait-associated. On simulated and real GWAS data, cWAS showed good statistical power with newly identified significant GRP associations in disease-associated tissues. More specifically, GRPs of endothelial and myofibroblasts in lung tissue were associated with Idiopathic Pulmonary Fibrosis and Chronic Obstructive Pulmonary Disease, respectively. For breast cancer, the GRP of blood CD8+ T cells was negatively associated with breast cancer (BC) risk as well as survival. Overall, cWAS is a powerful tool to reveal cell types associated with complex diseases mediated by GRPs.

Impaired AIF-CHCHD4 interaction and mitochondrial calcium overload contribute to auditory neuropathy spectrum disorder in patient-iPSC-derived neurons with AIFM1 variant

Auditory neuropathy spectrum disorder (ANSD) is a hearing impairment caused by dysfunction of inner hair cells, ribbon synapses, spiral ganglion neurons and/or the auditory nerve itself. Approximately 1/7000 newborns have abnormal auditory nerve function, accounting for 10%-14% of cases of permanent hearing loss in children. Although we previously identified the AIFM1 c.1265 G > A variant to be associated with ANSD, the mechanism by which ANSD is associated with AIFM1 is poorly understood. We generated induced pluripotent stem cells (iPSCs) from peripheral blood mononuclear cells (PBMCs) via nucleofection with episomal plasmids. The patient-specific iPSCs were edited via CRISPR/Cas9 technology to generate gene-corrected isogenic iPSCs. These iPSCs were further differentiated into neurons via neural stem cells (NSCs). The pathogenic mechanism was explored in these neurons. In patient cells (PBMCs, iPSCs, and neurons), the AIFM1 c.1265 G > A variant caused a novel splicing variant (c.1267-1305del), resulting in AIF p.R422Q and p.423-435del proteins, which impaired AIF dimerization. Such impaired AIF dimerization then weakened the interaction between AIF and coiled-coil-helix-coiled-coil-helix domain-containing protein 4 (CHCHD4). On the one hand, the mitochondrial import of ETC complex subunits was inhibited, subsequently leading to an increased ADP/ATP ratio and elevated ROS levels. On the other hand, MICU1-MICU2 heterodimerization was impaired, leading to mCa2+ overload. Calpain was activated by mCa2+ and subsequently cleaved AIF for its translocation into the nucleus, ultimately resulting in caspase-independent apoptosis. Interestingly, correction of the AIFM1 variant significantly restored the structure and function of AIF, further improving the physiological state of patient-specific iPSC-derived neurons. This study demonstrates that the AIFM1 variant is one of the molecular bases of ANSD. Mitochondrial dysfunction, especially mCa2+ overload, plays a prominent role in ANSD associated with AIFM1. Our findings help elucidate the mechanism of ANSD and may lead to the provision of novel therapies.

Inferring cell diversity in single cell data using consortium-scale epigenetic data as a biological anchor for cell identity

Methods for cell clustering and gene expression from single-cell RNA sequencing (scRNA-seq) data are essential for biological interpretation of cell processes. Here, we present TRIAGE-Cluster which uses genome-wide epigenetic data from diverse bio-samples to identify genes demarcating cell diversity in scRNA-seq data. By integrating patterns of repressive chromatin deposited across diverse cell types with weighted density estimation, TRIAGE-Cluster determines cell type clusters in a 2D UMAP space. We then present TRIAGE-ParseR, a machine learning method which evaluates gene expression rank lists to define gene groups governing the identity and function of cell types. We demonstrate the utility of this two-step approach using atlases of in vivo and in vitro cell diversification and organogenesis. We also provide a web accessible dashboard for analysis and download of data and software. Collectively, genome-wide epigenetic repression provides a versatile strategy to define cell diversity and study gene regulation of scRNA-seq data.

A village in a dish model system for population-scale hiPSC studies

The mechanisms by which DNA alleles contribute to disease risk, drug response, and other human phenotypes are highly context-specific, varying across cell types and different conditions. Human induced pluripotent stem cells are uniquely suited to study these context-dependent effects but cell lines from hundreds or thousands of individuals are required. Village cultures, where multiple induced pluripotent stem lines are cultured and differentiated in a single dish, provide an elegant solution for scaling induced pluripotent stem experiments to the necessary sample sizes required for population-scale studies. Here, we show the utility of village models, demonstrating how cells can be assigned to an induced pluripotent stem line using single-cell sequencing and illustrating that the genetic, epigenetic or induced pluripotent stem line-specific effects explain a large percentage of gene expression variation for many genes. We demonstrate that village methods can effectively detect induced pluripotent stem line-specific effects, including sensitive dynamics of cell states.

Mapping genetic effects on cell type-specific chromatin accessibility and annotating complex immune trait variants using single nucleus ATAC-seq in peripheral blood

Gene regulation is highly cell type-specific and understanding the function of non-coding genetic variants associated with complex traits requires molecular phenotyping at cell type resolution. In this study we performed single nucleus ATAC-seq (snATAC-seq) and genotyping in peripheral blood mononuclear cells from 13 individuals. Clustering chromatin accessibility profiles of 96,002 total nuclei identified 17 immune cell types and sub-types. We mapped chromatin accessibility QTLs (caQTLs) in each immune cell type and sub-type using individuals of European ancestry which identified 6,901 caQTLs at FDR < .10 and 4,220 caQTLs at FDR < .05, including those obscured from assays of bulk tissue such as with divergent effects on different cell types. For 3,941 caQTLs we further annotated putative target genes of variant activity using single cell co-accessibility, and caQTL variants were significantly correlated with the accessibility level of linked gene promoters. We fine-mapped loci associated with 16 complex immune traits and identified immune cell caQTLs at 622 candidate causal variants, including those with cell type-specific effects. At the 6q15 locus associated with type 1 diabetes, in line with previous reports, variant rs72928038 was a naïve CD4+ T cell caQTL linked to BACH2 and we validated the allelic effects of this variant on regulatory activity in Jurkat T cells. These results highlight the utility of snATAC-seq for mapping genetic effects on accessible chromatin in specific cell types.

Mesenchymal properties of iPSC-derived neural progenitors that generate undesired grafts after transplantation

Although neural stem/progenitor cells derived from human induced pluripotent stem cells (hiPSC-NS/PCs) are expected to be a cell source for cell-based therapy, tumorigenesis of hiPSC-NS/PCs is a potential problem for clinical applications. Therefore, to understand the mechanisms of tumorigenicity in NS/PCs, we clarified the cell populations of NS/PCs. We established single cell-derived NS/PC clones (scNS/PCs) from hiPSC-NS/PCs that generated undesired grafts. Additionally, we performed bioassays on scNS/PCs, which classified cell types within parental hiPSC-NS/PCs. Interestingly, we found unique subsets of scNS/PCs, which exhibited the transcriptome signature of mesenchymal lineages. Furthermore, these scNS/PCs expressed both neural (PSA-NCAM) and mesenchymal (CD73 and CD105) markers, and had an osteogenic differentiation capacity. Notably, eliminating CD73⁺ CD105⁺ cells from among parental hiPSC-NS/PCs ensured the quality of hiPSC-NS/PCs. Taken together, the existence of unexpected cell populations among NS/PCs may explain their tumorigenicity leading to potential safety issues of hiPSC-NS/PCs for future regenerative medicine.

Genetic control of the dynamic transcriptional response to immune stimuli and glucocorticoids at single-cell resolution

Synthetic glucocorticoids, such as dexamethasone, have been used as a treatment for many immune conditions, such as asthma and, more recently, severe COVID-19. Single-cell data can capture more fine-grained details on transcriptional variability and dynamics to gain a better understanding of the molecular underpinnings of inter-individual variation in drug response. Here, we used single-cell RNA-seq to study the dynamics of the transcriptional response to glucocorticoids in activated peripheral blood mononuclear cells from 96 African American children. We used novel statistical approaches to calculate a mean-independent measure of gene expression variability and a measure of transcriptional response pseudotime. Using these approaches, we showed that glucocorticoids reverse the effects of immune stimulation on both gene expression mean and variability. Our novel measure of gene expression response dynamics, based on the diagonal linear discriminant analysis, separated individual cells by response status on the basis of their transcriptional profiles and allowed us to identify different dynamic patterns of gene expression along the response pseudotime. We identified genetic variants regulating gene expression mean and variability, including treatment-specific effects, and showed widespread genetic regulation of the transcriptional dynamics of the gene expression response.

Debiased personalized gene coexpression networks for population-scale scRNA-seq data

Population-scale single-cell RNA-seq (scRNA-seq) data sets create unique opportunities for quantifying expression variation across individuals at the gene coexpression network level. Estimation of coexpression networks is well established for bulk RNA-seq; however, single-cell measurements pose novel challenges owing to technical limitations and noise levels of this technology. Gene-gene correlation estimates from scRNA-seq tend to be severely biased toward zero for genes with low and sparse expression. Here, we present Dozer to debias gene-gene correlation estimates from scRNA-seq data sets and accurately quantify network-level variation across individuals. Dozer corrects correlation estimates in the general Poisson measurement model and provides a metric to quantify genes measured with high noise. Computational experiments establish that Dozer estimates are robust to mean expression levels of the genes and the sequencing depths of the data sets. Compared with alternatives, Dozer results in fewer false-positive edges in the coexpression networks, yields more accurate estimates of network centrality measures and modules, and improves the faithfulness of networks estimated from separate batches of the data sets. We showcase unique analyses enabled by Dozer in two population-scale scRNA-seq applications. Coexpression network-based centrality analysis of multiple differentiating human induced pluripotent stem cell (iPSC) lines yields biologically coherent gene groups that are associated with iPSC differentiation efficiency. Application with population-scale scRNA-seq of oligodendrocytes from postmortem human tissues of Alzheimer's disease and controls uniquely reveals coexpression modules of innate immune response with distinct coexpression levels between the diagnoses. Dozer represents an important advance in estimating personalized coexpression networks from scRNA-seq data.