Integrative single cell multiomics analysis of human retina indicates a role for hierarchical transcription factors collaboration in genetic effects on gene regulation

Abstract

Abstract Background Systematic characterization of how genetic variation modulates gene regulation in a cell type specific context is essential for understanding complex traits. To address this question, we profiled gene expression and chromatin state of cells from healthy retinae of 20 human donors with a single-cell multiomics approach, and performed genomic sequencing. Results We mapped single-cell eQTL (sc-eQTLs), single-cell caQTL (sc-caQTL), single-cell allelic specific chromatin accessibility (sc-ASCA) and single-cell allelic specific expression (sc-ASE) in major retinal cell types. By integrating these results, we identified and characterized regulatory elements and genetic variants effective on gene regulation in individual cell types. Most of the sc-eQTLs and sc-caQTLs identified show cell type specific effects, while the cis-elements containing the genetic variants with cell type specific effects tend to be accessible in multiple cell types. Furthermore, the transcription factors with binding sites perturbed by genetic variants tend to have higher expression in the cell types, where the variants have effect, than the cell types where the variants do not have effect. Finally, we identified the enriched cell types, candidate causal variants and genes, and cell type specific regulatory mechanism underlying GWAS loci. Conclusions Overall, genetic effects on gene regulation are highly context dependent. Our results suggest that among cell types sharing a similar lineage, cell type dependent genetic effect is primarily driven by trans-factors rather than cell type specific chromatin state of cis-elements. Our findings indicate a role for hierarchical transcription factors collaboration in cell type specific effects of genetic variants on gene regulation.