Molecular convergence of risk variants for congenital heart defects leveraging a regulatory map of the human fetal heart
Authors: Ma X. Rosa , Conley Stephanie D. , Kosicki Michael , Bredikhin Danila , Cui Ran , Tran Steven , Sheth Maya U. , Qiu Wei-Lin , Chen Sijie , Kundu Soumya , Kang Helen Y. , Amgalan Dulguun , Munger Chad J. , Duan Lauren , Dang Katherine , Rubio Oriane Matthys , Kany Shinwan , Zamirpour Siavash , DePaolo John , Padmanabhan Arun ,   Birth Defects Research Laboratory, Olgin Jeffrey , Damrauer Scott , Andersson Robin , Gu Mingxia , Priest James R. , Quertermous Thomas , Qiu Xiaojie , Rabinovitch Marlene , Visel Axel , Pennacchio Len , Kundaje Anshul , Glass Ian A. , Gifford Casey A. , Pirruccello James P. , Goodyer William R. , Engreitz Jesse M. 
Published: November 2024

Congenital heart defects (CHD) arise in part due to inherited genetic variants that alter genes and noncoding regulatory elements in the human genome. These variants are thought to act during fetal development to influence the formation of different heart structures. However, identifying the genes, pathways, and cell types that mediate these effects has been challenging due to the immense diversity of cell types involved in heart development as well as the superimposed complexities of interpreting noncoding sequences. As such, understanding the molecular functions of both noncoding and coding variants remains paramount to our fundamental understanding of cardiac development and CHD. Here, we created a gene regulation map of the healthy human fetal heart across developmental time, and applied it to interpret the functions of variants associated with CHD and quantitative cardiac traits. We collected single-cell multiomic data from 734,000 single cells sampled from 41 fetal hearts spanning post-conception weeks 6 to 22, enabling the construction of gene regulation maps in 90 cardiac cell types and states, including rare populations of cardiac conduction cells. Through an unbiased analysis of all 90 cell types, we find that both rare coding variants associated with CHD and common noncoding variants associated with valve traits converge to affect valvular interstitial cells (VICs). VICs are enriched for high expression of known CHD genes previously identified through mapping of rare coding variants. Eight CHD genes, as well as other genes in similar molecular pathways, are linked to common noncoding variants associated with other valve diseases or traits via enhancers in VICs. In addition, certain common noncoding variants impact enhancers with activities highly specific to particular subanatomic structures in the heart, illuminating how such variants can impact specific aspects of heart structure and function. Together, these results implicate new enhancers, genes, and cell types in the genetic etiology of CHD, identify molecular convergence of common noncoding and rare coding variants on VICs, and suggest a more expansive view of the cell types instrumental in genetic risk for CHD, beyond the working cardiomyocyte. This regulatory map of the human fetal heart will provide a foundational resource for understanding cardiac development, interpreting genetic variants associated with heart disease, and discovering targets for cell-type specific therapies.