We propose to study gene regulatory noise in differentiated cardiomyocytes using single cell technology and use these functional data to identidy novel CVD risk loci. Genome-wide association studies (GWAS) have identified many variants associated with cardiovascular- related diseases, some of which are novel. However, similar to studies of other common diseases, these identified risk-associated variants fail to explain a significant portion of the genetic heritability of cardiovascular disease (CVD). Moreover, many associated variants are non-coding without obvious function, though putatively, these are involved in gene regulation. By combining results of GWAS with functional genomic data (for example, eQTL mapping), one can identify variants that influence molecular functions and are also associated with disease risk. Individual cells are expected to tolerate uncertainties in the form of both external and internal perturbations arising from variable environments or mutations. This is especially critical in the context of cell fate transitions during differentiation. It has long been recognized that robustness is an inherent property of all biological systems and is strongly favored by evolution. Depending on their different roles, the regulation of subsets of genes is required to be particularly robust in order to maintain the phenotype or the identity of a cell. Many dynamic physiological processes must also be robust, and as a result, loss or grain of robustness is associated with certain clinically relevant phenotypes and complex genetic disease. In particular, inter-individual variation in penetrance of a mutation associated with a disease, or the variability in response to drug, can also be explained at times by different degrees of robustness. Despite the importance of robustness and the regulation of noise as mechanisms that maintain high fitness, we still have a relatively poor understanding of how robustness is achieved and how is noise being regulated at the molecular level. To take first steps towards understanding how robustness is regulated in humans, and to identify loci where mutations can affect robustness and underlie CVD risk, we will use single cell technology to collect gene expression data from differentiated cardiomyocytes. We will use a detailed time course design and high- resolution single cell gene expression data. Our approach will allow us to map variance QTLs (robustness QTLs) in addition to the more standard expression QTLs. Robustness QTLs may be of particular importance as contributors to CVD risk, a diseases in which threshold effects are predominant. Specifically, we will propose to collect single cell RNA-seq throughout cardiomyocyte differentiation of 70 Hutterite individuals (Aim 1), map eQTLs, as well as genetic loci associated with inter-individual variation in gene expression robustness (Aim 2), and Integrate eQTL and robustness QTL mapping with GWAS results to identify variants associated with CVD risk and with CVD-related phenotypes (Aim 3).

Public Health Relevance

The goals of this study are to characterize regulatory robustness in single-cell gene expression data from differentiated cardiomyocytes, identify genomic loci in which genetic variation is associated with inter-individual differences in regulatory noise (robustness QTLs), and use them to identify novel loci that are associated with CVD risk.

National Institute of Health (NIH)
National Heart, Lung, and Blood Institute (NHLBI)
Research Project (R01)
Project #
Application #
Study Section
Genetics of Health and Disease Study Section (GHD)
Program Officer
Papanicolaou, George
Project Start
Project End
Budget Start
Budget End
Support Year
Fiscal Year
Total Cost
Indirect Cost
University of Chicago
Schools of Medicine
United States
Zip Code
Banovich, Nicholas E; Li, Yang I; Raj, Anil et al. (2018) Impact of regulatory variation across human iPSCs and differentiated cells. Genome Res 28:122-131
Knowles, David A; Burrows, Courtney K; Blischak, John D et al. (2018) Determining the genetic basis of anthracycline-cardiotoxicity by molecular response QTL mapping in induced cardiomyocytes. Elife 7:
Engelmann, Brett W; Hsiao, Chiaowen Joyce; Blischak, John D et al. (2018) A Methodological Assessment and Characterization of Genetically-Driven Variation in Three Human Phosphoproteomes. Sci Rep 8:12106
Tung, Po-Yuan; Blischak, John D; Hsiao, Chiaowen Joyce et al. (2017) Batch effects and the effective design of single-cell gene expression studies. Sci Rep 7:39921
Cusanovich, Darren A; Caliskan, Minal; Billstrand, Christine et al. (2016) Integrated analyses of gene expression and genetic association studies in a founder population. Hum Mol Genet 25:2104-2112
Burrows, Courtney K; Banovich, Nicholas E; Pavlovic, Bryan J et al. (2016) Genetic Variation, Not Cell Type of Origin, Underlies the Majority of Identifiable Regulatory Differences in iPSCs. PLoS Genet 12:e1005793
Thomas, Samantha M; Kagan, Courtney; Pavlovic, Bryan J et al. (2015) Reprogramming LCLs to iPSCs Results in Recovery of Donor-Specific Gene Expression Signature. PLoS Genet 11:e1005216
Zhou, Xiang; Stephens, Matthew (2014) Efficient multivariate linear mixed model algorithms for genome-wide association studies. Nat Methods 11:407-9
Zhou, Xiang; Carbonetto, Peter; Stephens, Matthew (2013) Polygenic modeling with bayesian sparse linear mixed models. PLoS Genet 9:e1003264
Mizrahi-Man, Orna; Davenport, Emily R; Gilad, Yoav (2013) Taxonomic classification of bacterial 16S rRNA genes using short sequencing reads: evaluation of effective study designs. PLoS One 8:e53608

Showing the most recent 10 out of 16 publications