Multi-chamber cardioids unravel human heart development and cardiac defects.

The number one cause of human fetal death are defects in heart development. Because the human embryonic heart is inaccessible and the impacts of mutations, drugs, and environmental factors on the specialized functions of different heart compartments are not captured by in vitro models, determining the underlying causes is difficult. Here, we established a human cardioid platform that recapitulates the development of all major embryonic heart compartments, including right and left ventricles, atria, outflow tract, and atrioventricular canal. By leveraging 2D and 3D differentiation, we efficiently generated progenitor subsets with distinct first, anterior, and posterior second heart field identities. This advance enabled the reproducible generation of cardioids with compartment-specific in vivo-like gene expression profiles, morphologies, and functions. We used this platform to unravel the ontogeny of signal and contraction propagation between interacting heart chambers and dissect how mutations, teratogens, and drugs cause compartment-specific defects in the developing human heart.

Spatially organized cellular communities form the developing human heart.

The heart, which is the first organ to develop, is highly dependent on its form to function1,2. However, how diverse cardiac cell types spatially coordinate to create the complex morphological structures that are crucial for heart function remains unclear. Here we integrated single-cell RNA-sequencing with high-resolution multiplexed error-robust fluorescence in situ hybridization to resolve the identity of the cardiac cell types that develop the human heart. This approach also provided a spatial mapping of individual cells that enables illumination of their organization into cellular communities that form distinct cardiac structures. We discovered that many of these cardiac cell types further specified into subpopulations exclusive to specific communities, which support their specialization according to the cellular ecosystem and anatomical region. In particular, ventricular cardiomyocyte subpopulations displayed an unexpected complex laminar organization across the ventricular wall and formed, with other cell subpopulations, several cellular communities. Interrogating cell-cell interactions within these communities using in vivo conditional genetic mouse models and in vitro human pluripotent stem cell systems revealed multicellular signalling pathways that orchestrate the spatial organization of cardiac cell subpopulations during ventricular wall morphogenesis. These detailed findings into the cellular social interactions and specialization of cardiac cell types constructing and remodelling the human heart offer new insights into structural heart diseases and the engineering of complex multicellular tissues for human heart repair.

A wearable cardiac ultrasound imager.

Continuous imaging of cardiac functions is highly desirable for the assessment of long-term cardiovascular health, detection of acute cardiac dysfunction and clinical management of critically ill or surgical patients1-4. However, conventional non-invasive approaches to image the cardiac function cannot provide continuous measurements owing to device bulkiness5-11, and existing wearable cardiac devices can only capture signals on the skin12-16. Here we report a wearable ultrasonic device for continuous, real-time and direct cardiac function assessment. We introduce innovations in device design and material fabrication that improve the mechanical coupling between the device and human skin, allowing the left ventricle to be examined from different views during motion. We also develop a deep learning model that automatically extracts the left ventricular volume from the continuous image recording, yielding waveforms of key cardiac performance indices such as stroke volume, cardiac output and ejection fraction. This technology enables dynamic wearable monitoring of cardiac performance with substantially improved accuracy in various environments.

Phospholipase Cε hydrolyzes perinuclear phosphatidylinositol 4-phosphate to regulate cardiac hypertrophy.

Phospholipase Cε (PLCε) is a multifunctional enzyme implicated in cardiovascular, pancreatic, and inflammatory functions. Here we show that conditional deletion of PLCε in mouse cardiac myocytes protects from stress-induced pathological hypertrophy. PLCε small interfering RNA (siRNA) in ventricular myocytes decreases endothelin-1 (ET-1)-dependent elevation of nuclear calcium and activation of nuclear protein kinase D (PKD). PLCε scaffolded to muscle-specific A kinase-anchoring protein (mAKAP), along with PKCε and PKD, localizes these components at or near the nuclear envelope, and this complex is required for nuclear PKD activation. Phosphatidylinositol 4-phosphate (PI4P) is identified as a perinuclear substrate in the Golgi apparatus for mAKAP-scaffolded PLCε. We conclude that perinuclear PLCε, scaffolded to mAKAP in cardiac myocytes, responds to hypertrophic stimuli to generate diacylglycerol (DAG) from PI4P in the Golgi apparatus, in close proximity to the nuclear envelope, to regulate activation of nuclear PKD and hypertrophic signaling pathways.

Single-nucleus profiling of human dilated and hypertrophic cardiomyopathy.

Heart failure encompasses a heterogeneous set of clinical features that converge on impaired cardiac contractile function1,2 and presents a growing public health concern. Previous work has highlighted changes in both transcription and protein expression in failing hearts3,4, but may overlook molecular changes in less prevalent cell types. Here we identify extensive molecular alterations in failing hearts at single-cell resolution by performing single-nucleus RNA sequencing of nearly 600,000 nuclei in left ventricle samples from 11 hearts with dilated cardiomyopathy and 15 hearts with hypertrophic cardiomyopathy as well as 16 non-failing hearts. The transcriptional profiles of dilated or hypertrophic cardiomyopathy hearts broadly converged at the tissue and cell-type level. Further, a subset of hearts from patients with cardiomyopathy harbour a unique population of activated fibroblasts that is almost entirely absent from non-failing samples. We performed a CRISPR-knockout screen in primary human cardiac fibroblasts to evaluate this fibrotic cell state transition; knockout of genes associated with fibroblast transition resulted in a reduction of myofibroblast cell-state transition upon TGFß1 stimulation for a subset of genes. Our results provide insights into the transcriptional diversity of the human heart in health and disease as well as new potential therapeutic targets and biomarkers for heart failure.

Cardiomyocyte orientation recovery at micrometer scale reveals long-axis fiber continuum in heart walls.

Coordinated cardiomyocyte contraction drives the mammalian heart to beat and circulate blood. No consensus model of cardiomyocyte geometrical arrangement exists, due to the limited spatial resolution of whole heart imaging methods and the piecemeal nature of studies based on histological sections. By combining microscopy and computer vision, we produced the first-ever three-dimensional cardiomyocyte orientation reconstruction across mouse ventricular walls at the micrometer scale, representing a gain of three orders of magnitude in spatial resolution. We recovered a cardiomyocyte arrangement aligned to the long-axis direction of the outer ventricular walls. This cellular network lies in a thin shell and forms a continuum with longitudinally arranged cardiomyocytes in the inner walls, with a complex geometry at the apex. Our reconstruction methods can be applied at fine spatial scales to further understanding of heart wall electrical function and mechanics, and set the stage for the study of micron-scale fiber remodeling in heart disease.

Llgl1 mediates timely epicardial emergence and establishment of an apical laminin sheath around the trabeculating cardiac ventricle.

During heart development, the embryonic ventricle becomes enveloped by the epicardium, which adheres to the outer apical surface of the heart. This is concomitant with onset of ventricular trabeculation, where a subset of cardiomyocytes lose apicobasal polarity and delaminate basally from the ventricular wall. Llgl1 regulates the formation of apical cell junctions and apicobasal polarity, and we investigated its role in ventricular wall maturation. We found that llgl1 mutant zebrafish embryos exhibit aberrant apical extrusion of ventricular cardiomyocytes. While investigating apical cardiomyocyte extrusion, we identified a basal-to-apical shift in laminin deposition from the internal to the external ventricular wall. We find that epicardial cells express several laminin subunits as they adhere to the ventricle, and that the epicardium is required for laminin deposition on the ventricular surface. In llgl1 mutants, timely establishment of the epicardial layer is disrupted due to delayed emergence of epicardial cells, resulting in delayed apical deposition of laminin on the ventricular surface. Together, our analyses reveal an unexpected role for Llgl1 in correct timing of epicardial development, supporting integrity of the ventricular myocardial wall.

Mechanical forces remodel the cardiac extracellular matrix during zebrafish development.

The cardiac extracellular matrix (cECM) is fundamental for organ morphogenesis and maturation, during which time it undergoes remodeling, yet little is known about whether mechanical forces generated by the heartbeat regulate this remodeling process. Using zebrafish as a model and focusing on stages when cardiac valves and trabeculae form, we found that altering cardiac contraction impairs cECM remodeling. Longitudinal volumetric quantifications in wild-type animals revealed region-specific dynamics: cECM volume decreases in the atrium but not in the ventricle or atrioventricular canal. Reducing cardiac contraction resulted in opposite effects on the ventricular and atrial ECM, whereas increasing the heart rate affected the ventricular ECM but had no effect on the atrial ECM, together indicating that mechanical forces regulate the cECM in a chamber-specific manner. Among the ECM remodelers highly expressed during cardiac morphogenesis, we found one that was upregulated in non-contractile hearts, namely tissue inhibitor of matrix metalloproteinase 2 (timp2). Loss- and gain-of-function analyses of timp2 revealed its crucial role in cECM remodeling. Altogether, our results indicate that mechanical forces control cECM remodeling in part through timp2 downregulation.

Fueling the heartbeat: Dynamic regulation of intracellular ATP during excitation-contraction coupling in ventricular myocytes.

The heart beats approximately 100,000 times per day in humans, imposing substantial energetic demands on cardiac muscle. Adenosine triphosphate (ATP) is an essential energy source for normal function of cardiac muscle during each beat, as it powers ion transport, intracellular Ca2+ handling, and actin-myosin cross-bridge cycling. Despite this, the impact of excitation-contraction coupling on the intracellular ATP concentration ([ATP]i) in myocytes is poorly understood. Here, we conducted real-time measurements of [ATP]i in ventricular myocytes using a genetically encoded ATP fluorescent reporter. Our data reveal rapid beat-to-beat variations in [ATP]i. Notably, diastolic [ATP]i was <1 mM, which is eightfold to 10-fold lower than previously estimated. Accordingly, ATP-sensitive K+ (KATP) channels were active at physiological [ATP]i. Cells exhibited two distinct types of ATP fluctuations during an action potential: net increases (Mode 1) or decreases (Mode 2) in [ATP]i. Mode 1 [ATP]i increases necessitated Ca2+ entry and release from the sarcoplasmic reticulum (SR) and were associated with increases in mitochondrial Ca2+. By contrast, decreases in mitochondrial Ca2+ accompanied Mode 2 [ATP]i decreases. Down-regulation of the protein mitofusin 2 reduced the magnitude of [ATP]i fluctuations, indicating that SR-mitochondrial coupling plays a crucial role in the dynamic control of ATP levels. Activation of ß-adrenergic receptors decreased [ATP]i, underscoring the energetic impact of this signaling pathway. Finally, our work suggests that cross-bridge cycling is the largest consumer of ATP in a ventricular myocyte during an action potential. These findings provide insights into the energetic demands of EC coupling and highlight the dynamic nature of ATP concentrations in cardiac muscle.

Emerging Technologies in Cardiac Pacing.

Cardiac pacing to treat bradyarrhythmias has evolved in recent decades. Recognition that a substantial proportion of pacemaker-dependent patients can develop heart failure due to electrical and mechanical dyssynchrony from traditional right ventricular apical pacing has led to development of more physiologic pacing methods that better mimic normal cardiac conduction and provide synchronized ventricular contraction. Conventional biventricular pacing has been shown to benefit patients with heart failure and conduction system disease but can be limited by scarring and fibrosis. His bundle pacing and left bundle branch area pacing are novel techniques that can provide more physiologic ventricular activation as an alternative to conventional or biventricular pacing. Leadless pacing has emerged as another alternative pacing technique to overcome limitations in conventional transvenous pacemaker systems. Our objective is to review the evolution of cardiac pacing and explore these new advances in pacing strategies.

Parallels between oncogene-driven cardiac hyperplasia and heart regeneration in zebrafish.

The human heart is poorly regenerative and cardiac tumors are extremely rare. Whether the adult zebrafish myocardium is responsive to oncogene overexpression and how this condition affects its intrinsic regenerative capacity remains unknown. Here, we have established a strategy of inducible and reversible expression of HRASG12V in zebrafish cardiomyocytes. This approach stimulated a hyperplastic cardiac enlargement within 16 days. The phenotype was suppressed by rapamycin-mediated inhibition of TOR signaling. As TOR signaling is also required for heart restoration after cryoinjury, we compared transcriptomes of hyperplastic and regenerating ventricles. Both conditions were associated with upregulation of cardiomyocyte dedifferentiation and proliferation factors, as well as with similar microenvironmental responses, such as deposition of nonfibrillar Collagen XII and recruitment of immune cells. Among the differentially expressed genes, many proteasome and cell-cycle regulators were upregulated only in oncogene-expressing hearts. Preconditioning of the heart with short-term oncogene expression accelerated cardiac regeneration after cryoinjury, revealing a beneficial synergism between both programs. Identification of the molecular bases underlying the interplay between detrimental hyperplasia and advantageous regeneration provides new insights into cardiac plasticity in adult zebrafish.

3D Imaging Reveals Complex Microvascular Remodeling in the Right Ventricle in Pulmonary Hypertension.

BACKGROUND: Pathogenic concepts of right ventricular (RV) failure in pulmonary arterial hypertension focus on a critical loss of microvasculature. However, the methods underpinning prior studies did not take into account the 3-dimensional (3D) aspects of cardiac tissue, making accurate quantification difficult. We applied deep-tissue imaging to the pressure-overloaded RV to uncover the 3D properties of the microvascular network and determine whether deficient microvascular adaptation contributes to RV failure. METHODS: Heart sections measuring 250-µm-thick were obtained from mice after pulmonary artery banding (PAB) or debanding PAB surgery and properties of the RV microvascular network were assessed using 3D imaging and quantification. Human heart tissues harvested at the time of transplantation from pulmonary arterial hypertension cases were compared with tissues from control cases with normal RV function. RESULTS: Longitudinal 3D assessment of PAB mouse hearts uncovered complex microvascular remodeling characterized by tortuous, shorter, thicker, highly branched vessels, and overall preserved microvascular density. This remodeling process was reversible in debanding PAB mice in which the RV function recovers over time. The remodeled microvasculature tightly wrapped around the hypertrophied cardiomyocytes to maintain a stable contact surface to cardiomyocytes as an adaptation to RV pressure overload, even in end-stage RV failure. However, microvasculature-cardiomyocyte contact was impaired in areas with interstitial fibrosis where cardiomyocytes displayed signs of hypoxia. Similar to PAB animals, microvascular density in the RV was preserved in patients with end-stage pulmonary arterial hypertension, and microvascular architectural changes appeared to vary by etiology, with patients with pulmonary veno-occlusive disease displaying a lack of microvascular complexity with uniformly short segments. CONCLUSIONS: 3D deep tissue imaging of the failing RV in PAB mice, pulmonary hypertension rats, and patients with pulmonary arterial hypertension reveals complex microvascular changes to preserve the microvascular density and maintain a stable microvascular-cardiomyocyte contact. Our studies provide a novel framework to understand microvascular adaptation in the pressure-overloaded RV that focuses on cell-cell interaction and goes beyond the concept of capillary rarefaction.

Cells of the adult human heart.

Cardiovascular disease is the leading cause of death worldwide. Advanced insights into disease mechanisms and therapeutic strategies require a deeper understanding of the molecular processes involved in the healthy heart. Knowledge of the full repertoire of cardiac cells and their gene expression profiles is a fundamental first step in this endeavour. Here, using state-of-the-art analyses of large-scale single-cell and single-nucleus transcriptomes, we characterize six anatomical adult heart regions. Our results highlight the cellular heterogeneity of cardiomyocytes, pericytes and fibroblasts, and reveal distinct atrial and ventricular subsets of cells with diverse developmental origins and specialized properties. We define the complexity of the cardiac vasculature and its changes along the arterio-venous axis. In the immune compartment, we identify cardiac-resident macrophages with inflammatory and protective transcriptional signatures. Furthermore, analyses of cell-to-cell interactions highlight different networks of macrophages, fibroblasts and cardiomyocytes between atria and ventricles that are distinct from those of skeletal muscle. Our human cardiac cell atlas improves our understanding of the human heart and provides a valuable reference for future studies.

Hey2 enhancer activity defines unipotent progenitors for left ventricular cardiomyocytes in juxta-cardiac field of early mouse embryo.

The cardiac crescent is the first structure of the heart and contains progenitor cells of the first heart field, which primarily differentiate into left ventricular cardiomyocytes. The interface between the forming cardiac crescent and extraembryonic tissue is known as the juxta-cardiac field (JCF), and progenitor cells in this heart field contribute to the myocardium of the left ventricle and atrioventricular canal as well as the epicardium. However, it is unclear whether there are progenitor cells that differentiate specifically into left ventricular cardiomyocytes. We have previously demonstrated that an enhancer of the gene encoding the Hey2 bHLH transcriptional repressor is activated in the ventricular myocardium during mouse embryonic development. In this study, we aimed to investigate the characteristics of cardiomyocyte progenitor cells and their cell lineages by analyzing Hey2 enhancer activity at the earliest stages of heart formation. We found that the Hey2 enhancer initiated its activity prior to cardiomyocyte differentiation within the JCF. Hey2 enhancer-active cells were present rostrally to the Tbx5-expressing region at the early phase of cardiac crescent formation and differentiated exclusively into left ventricular cardiomyocytes in a lineage distinct from the Tbx5-positive lineage. By the late phase of cardiac crescent formation, Hey2 enhancer activity became significantly overlapped with Tbx5 expression in cells that contribute to the left ventricular myocardium. Our study reveals that a population of unipotent progenitor cells for left ventricular cardiomyocytes emerge in the JCF, providing further insight into the mode of cell type diversification during early cardiac development.

Proteomics and phosphoproteomics of failing human left ventricle identifies dilated cardiomyopathy-associated phosphorylation of CTNNA3.

The prognosis and treatment outcomes of heart failure (HF) patients rely heavily on disease etiology, yet the majority of underlying signaling mechanisms are complex and not fully elucidated. Phosphorylation is a major point of protein regulation with rapid and profound effects on the function and activity of protein networks. Currently, there is a lack of comprehensive proteomic and phosphoproteomic studies examining cardiac tissue from HF patients with either dilated dilated cardiomyopathy (DCM) or ischemic cardiomyopathy (ICM). Here, we used a combined proteomic and phosphoproteomic approach to identify and quantify more than 5,000 total proteins with greater than 13,000 corresponding phosphorylation sites across explanted left ventricle (LV) tissue samples, including HF patients with DCM vs. nonfailing controls (NFC), and left ventricular infarct vs. noninfarct, and periinfarct vs. noninfarct regions of HF patients with ICM. Each pair-wise comparison revealed unique global proteomic and phosphoproteomic profiles with both shared and etiology-specific perturbations. With this approach, we identified a DCM-associated hyperphosphorylation cluster in the cardiomyocyte intercalated disc (ICD) protein, αT-catenin (CTNNA3). We demonstrate using both ex vivo isolated cardiomyocytes and in vivo using an AAV9-mediated overexpression mouse model, that CTNNA3 phosphorylation at these residues plays a key role in maintaining protein localization at the cardiomyocyte ICD to regulate conductance and cell-cell adhesion. Collectively, this integrative proteomic/phosphoproteomic approach identifies region- and etiology-associated signaling pathways in human HF and describes a role for CTNNA3 phosphorylation in the pathophysiology of DCM.

Dgcr8 functions in the secondary heart field for outflow tract and right ventricle development in mammals.

The DGCR8 gene, encoding a critical miRNA processing protein, maps within the hemizygous region in patients with 22q11.2 deletion syndrome. Most patients have malformations of the cardiac outflow tract that is derived in part from the anterior second heart field (aSHF) mesoderm. To understand the function of Dgcr8 in the aSHF, we inactivated it in mice using Mef2c-AHF-Cre. Inactivation resulted in a fully penetrant persistent truncus arteriosus and a hypoplastic right ventricle leading to lethality by E14.5. To understand the molecular mechanism for this phenotype, we performed gene expression profiling of the aSHF and the cardiac outflow tract with right ventricle in conditional null versus normal mouse littermates at stage E9.5 prior to morphology changes. We identified dysregulation of mRNA gene expression, of which some are relevant to cardiogenesis. Many pri-miRNA genes were strongly increased in expression in mutant embryos along with reduced expression of mature miRNA genes. We further examined the individual, mature miRNAs that were decreased in expression along with pri-miRNAs that were accumulated that could be direct effects due to loss of Dgcr8. Among these genes, were miR-1a, miR-133a, miR-134, miR143 and miR145a, which have known functions in heart development. These early mRNA and miRNA changes may in part, explain the first steps that lead to the resulting phenotype in Dgcr8 aSHF conditional mutant embryos.

The biochemically defined super relaxed state of myosin-A paradox.

The biochemical SRX (super-relaxed) state of myosin has been defined as a low ATPase activity state. This state can conserve energy when the myosin is not recruited for muscle contraction. The SRX state has been correlated with a structurally defined ordered (versus disordered) state of muscle thick filaments. The two states may be linked via a common interacting head motif (IHM) where the two heads of heavy meromyosin (HMM), or myosin, fold back onto each other and form additional contacts with S2 and the thick filament. Experimental observations of the SRX, IHM, and the ordered form of thick filaments, however, do not always agree, and result in a series of unresolved paradoxes. To address these paradoxes, we have reexamined the biochemical measurements of the SRX state for porcine cardiac HMM. In our hands, the commonly employed mantATP displacement assay was unable to quantify the population of the SRX state with all data fitting very well by a single exponential. We further show that mavacamten inhibits the basal ATPases of both porcine ventricle HMM and S1 (Ki, 0.32 and 1.76 µM respectively) while dATP activates HMM cooperatively without any evidence of an SRX state. A combination of our experimental observations and theories suggests that the displacement of mantATP in purified proteins is not a reliable assay to quantify the SRX population. This means that while the structurally defined IHM and ordered thick filaments clearly exist, great care must be employed when using the mantATP displacement assay.

Mineralocorticoid Receptor Antagonism Prevents Aortic Plaque Progression and Reduces Left Ventricular Mass and Fibrosis in Patients With Type 2 Diabetes and Chronic Kidney Disease: The MAGMA Trial.

BACKGROUND: Persistent mineralocorticoid receptor activation is a pathologic response in type 2 diabetes and chronic kidney disease. Whereas mineralocorticoid receptor antagonists are beneficial in reducing cardiovascular complications, direct mechanistic pathways for these effects in humans are lacking. METHODS: The MAGMA trial (Mineralocorticoid Receptor Antagonism Clinical Evaluation in Atherosclerosis) was a randomized, double-blind, placebo-controlled trial in patients with high-risk type 2 diabetes with chronic kidney disease (not receiving dialysis) on maximum tolerated renin-angiotensin system blockade. The primary end point was change in thoracic aortic wall volume, expressed as absolute or percent value (ΔTWV or ΔPWV), using 3T magnetic resonance imaging at 12 months. Secondary end points were changes in left ventricle (LV) mass; LV fibrosis, measured as a change in myocardial native T1; and 24-hour ambulatory and central aortic blood pressures. Tertiary end points included plasma proteomic changes in 7596 plasma proteins using an aptamer-based assay. RESULTS: A total of 79 patients were randomized to placebo (n=42) or 25 mg of spironolactone daily (n=37). After a modified intent-to-treat, including available baseline data of study end points, patients who completed the trial protocol were included in the final analyses. At the 12-month follow-up, the average change in PWV was 7.1±10.7% in the placebo group and 0.87±10.0% in the spironolactone group (P=0.028), and ΔTWV was 1.2±1.7 cm3 in the placebo group and 0.037±1.9 cm3 in the spironolactone group (P=0.022). Change in LV mass was 3.1±8.4 g in the placebo group and -5.8±8.4 g in the spironolactone group (P=0.001). Changes in LV T1 values were significantly different between the placebo and spironolactone groups (26.0±41.9 ms in the placebo group versus a decrease of -10.1±36.3 ms in the spironolactone group; P=6.33×10-4). Mediation analysis revealed that the spironolactone effect on thoracic aortic wall volume and myocardial mass remained significant after adjustment for ambulatory and central blood pressures. Proteomic analysis revealed a dominant effect of spironolactone on pathways involving oxidative stress, inflammation, and leukocyte activation. CONCLUSIONS: Among patients with diabetes with moderate to severe chronic kidney disease at elevated cardiovascular risk, treatment with spironolactone prevented progression of aortic wall volume and resulted in regression of LV mass and favorable alterations in native T1, suggesting amelioration of left-ventricular fibrosis. REGISTRATION: URL: https://www.clinicaltrials.gov; Unique identifier: NCT02169089.

Redefining tissue specificity of genetic regulation of gene expression in the presence of allelic heterogeneity.

Uncovering the functional impact of genetic variation on gene expression is important in understanding tissue biology and the pathogenesis of complex traits. Despite large efforts to map expression quantitative trait loci (eQTLs) across many human tissues, our ability to translate those findings to understanding human disease has been incomplete, and the majority of disease loci are not explained by association with expression of a target gene. Cell-type specificity and the presence of multiple independent causal variants for many eQTLs are potential confounders contributing to the apparent discrepancy with disease loci. In this study, we investigate the tissue specificity of genetic effects on gene expression and the overlap with disease loci while considering the presence of multiple causal variants within and across tissues. We find evidence of pervasive tissue specificity of eQTLs, often masked by linkage disequilibrium that misleads traditional meta-analytic approaches. We propose CAFEH (colocalization and fine-mapping in the presence of allelic heterogeneity), a Bayesian method that integrates genetic association data across multiple traits, incorporating linkage disequilibrium to identify causal variants. CAFEH outperforms previous approaches in colocalization and fine-mapping. Using CAFEH, we show that genes with highly tissue-specific genetic effects are under greater selection, enriched in differentiation and developmental processes, and more likely to be involved in human disease. Last, we demonstrate that CAFEH can efficiently leverage the widespread allelic heterogeneity in genetic regulation of gene expression to prioritize the target tissue in genome-wide association complex trait loci, thereby improving our ability to interpret complex trait genetics.

Dissecting mechanisms of chamber-specific cardiac differentiation and its perturbation following retinoic acid exposure.

The specification of distinct cardiac lineages occurs before chamber formation and acquisition of bona fide atrial or ventricular identity. However, the mechanisms underlying these early specification events remain poorly understood. Here, we performed single cell analysis at the murine cardiac crescent, primitive heart tube and heart tube stages to uncover the transcriptional mechanisms underlying formation of atrial and ventricular cells. We find that progression towards differentiated cardiomyocytes occurs primarily based on heart field progenitor identity, and that progenitors contribute to ventricular or atrial identity through distinct differentiation mechanisms. We identify new candidate markers that define such differentiation processes and examine their expression dynamics using computational lineage trajectory methods. We further show that exposure to exogenous retinoic acid causes defects in ventricular chamber size, dysregulation in FGF signaling and a shunt in differentiation towards orthogonal lineages. Retinoic acid also causes defects in cell-cycle exit resulting in formation of hypomorphic ventricles. Collectively, our data identify, at a single cell level, distinct lineage trajectories during cardiac specification and differentiation, and the precise effects of manipulating cardiac progenitor patterning via retinoic acid signaling.