The nutrient sensor CRTC & Sarcalumenin / Thinman represent a new pathway in cardiac hypertrophy.

Obesity and type 2 diabetes are at epidemic levels and a significant proportion of these patients are diagnosed with left ventricular hypertrophy. CREB R egulated T ranscription C o-activator ( CRTC ) is a key regulator of metabolism in mammalian hepatocytes, where it is activated by calcineurin (CaN) to increase expression of gluconeogenic genes. CaN is known its role in pathological cardiac hypertrophy, however, a role for CRTC in the heart has not been identified. In Drosophila , CRTC null mutants have little body fat and exhibit severe cardiac restriction, myofibrillar disorganization, cardiac fibrosis and tachycardia, all hallmarks of heart disease. Cardiac-specific knockdown of CRTC , or its coactivator CREBb , mimicked the reduced body fat and heart defects of CRTC null mutants. Comparative gene expression in CRTC loss- or gain-of-function fly hearts revealed contra-regulation of genes involved in glucose, fatty acid, and amino acid metabolism, suggesting that CRTC also acts as a metabolic switch in the heart. Among the contra-regulated genes with conserved CREB binding sites, we identified the fly ortholog of Sarcalumenin, which is a Ca 2+ -binding protein in the sarcoplasmic reticulum. Cardiac knockdown recapitulated the loss of CRTC cardiac restriction and fibrotic phenotypes, suggesting it is a downstream effector of CRTC we named thinman ( tmn ). Importantly, cardiac overexpression of either CaN or CRTC in flies caused hypertrophy that was reversed in a CRTC mutant background, suggesting CRTC mediates hypertrophy downstream of CaN, perhaps as an alternative to NFAT. CRTC novel role in the heart is likely conserved in vertebrates as knockdown in zebrafish also caused cardiac restriction, as in fl ies. These data suggest that CRTC is involved in myocardial cell maintenance and that CaN-CRTC- Sarcalumenin/ tmn signaling represents a novel and conserved pathway underlying cardiac hypertrophy.

Mitochondrial MICOS complex genes, implicated in hypoplastic left heart syndrome, maintain cardiac contractility and actomyosin integrity.

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart disease (CHD) with a likely oligogenic etiology, but our understanding of the genetic complexities and pathogenic mechanisms leading to HLHS is limited. We performed whole genome sequencing (WGS) on 183 HLHS patient-parent trios to identify candidate genes, which were functionally tested in the Drosophila heart model. Bioinformatic analysis of WGS data from an index family of a HLHS proband born to consanguineous parents prioritized 9 candidate genes with rare, predicted damaging homozygous variants. Of them, cardiac-specific knockdown (KD) of mitochondrial MICOS complex subunit dCHCHD3/6 resulted in drastically compromised heart contractility, diminished levels of sarcomeric actin and myosin, reduced cardiac ATP levels, and mitochondrial fission-fusion defects. These defects were similar to those inflicted by cardiac KD of ATP synthase subunits of the electron transport chain (ETC), consistent with the MICOS complex's role in maintaining cristae morphology and ETC assembly. Five additional HLHS probands harbored rare, predicted damaging variants in CHCHD3 or CHCHD6. Hypothesizing an oligogenic basis for HLHS, we tested 60 additional prioritized candidate genes from these patients for genetic interactions with CHCHD3/6 in sensitized fly hearts. Moderate KD of CHCHD3/6 in combination with Cdk12 (activator of RNA polymerase II), RNF149 (goliath, E3 ubiquitin ligase), or SPTBN1 (ß-Spectrin, scaffolding protein) caused synergistic heart defects, suggesting the likely involvement of diverse pathways in HLHS. Further elucidation of novel candidate genes and genetic interactions of potentially disease-contributing pathways is expected to lead to a better understanding of HLHS and other CHDs.

Multiplatform modeling of atrial fibrillation identifies phospholamban as a central regulator of cardiac rhythm.

Atrial fibrillation (AF) is a common and genetically inheritable form of cardiac arrhythmia; however, it is currently not known how these genetic predispositions contribute to the initiation and/or maintenance of AF-associated phenotypes. One major barrier to progress is the lack of experimental systems to investigate the effects of gene function on rhythm parameters in models with human atrial and whole-organ relevance. Here, we assembled a multi-model platform enabling high-throughput characterization of the effects of gene function on action potential duration and rhythm parameters using human induced pluripotent stem cell-derived atrial-like cardiomyocytes and a Drosophila heart model, and validation of the findings using computational models of human adult atrial myocytes and tissue. As proof of concept, we screened 20 AF-associated genes and identified phospholamban loss of function as a top conserved hit that shortens action potential duration and increases the incidence of arrhythmia phenotypes upon stress. Mechanistically, our study reveals that phospholamban regulates rhythm homeostasis by functionally interacting with L-type Ca2+ channels and NCX. In summary, our study illustrates how a multi-model system approach paves the way for the discovery and molecular delineation of gene regulatory networks controlling atrial rhythm with application to AF.

Trajectory reconstruction identifies dysregulation of perinatal maturation programs in pluripotent stem cell-derived cardiomyocytes.

A limitation in the application of pluripotent stem cell-derived cardiomyocytes (PSC-CMs) is the failure of these cells to achieve full functional maturity. The mechanisms by which directed differentiation differs from endogenous development, leading to consequent PSC-CM maturation arrest, remain unclear. Here, we generate a single-cell RNA sequencing (scRNA-seq) reference of mouse in vivo CM maturation with extensive sampling of previously difficult-to-isolate perinatal time periods. We subsequently generate isogenic embryonic stem cells to create an in vitro scRNA-seq reference of PSC-CM-directed differentiation. Through trajectory reconstruction, we identify an endogenous perinatal maturation program that is poorly recapitulated in vitro. By comparison with published human datasets, we identify a network of nine transcription factors (TFs) whose targets are consistently dysregulated in PSC-CMs across species. Notably, these TFs are only partially activated in common ex vivo approaches to engineer PSC-CM maturation. Our study can be leveraged toward improving the clinical viability of PSC-CMs.

Conserved transcription factors promote cell fate stability and restrict reprogramming potential in differentiated cells.

Defining the mechanisms safeguarding cell fate identity in differentiated cells is crucial to improve 1) - our understanding of how differentiation is maintained in healthy tissues or altered in a disease state, and 2) - our ability to use cell fate reprogramming for regenerative purposes. Here, using a genome-wide transcription factor screen followed by validation steps in a variety of reprogramming assays (cardiac, neural and iPSC in fibroblasts and endothelial cells), we identified a set of four transcription factors (ATF7IP, JUNB, SP7, and ZNF207 [AJSZ]) that robustly opposes cell fate reprogramming in both lineage and cell type independent manners. Mechanistically, our integrated multi-omics approach (ChIP, ATAC and RNA-seq) revealed that AJSZ oppose cell fate reprogramming by 1) - maintaining chromatin enriched for reprogramming TF motifs in a closed state and 2) - downregulating genes required for reprogramming. Finally, KD of AJSZ in combination with MGT overexpression, significantly reduced scar size and improved heart function by 50%, as compared to MGT alone post-myocardial infarction. Collectively, our study suggests that inhibition of barrier to reprogramming mechanisms represents a promising therapeutic avenue to improve adult organ function post-injury.

Genetic architecture of natural variation of cardiac performance from flies to humans.

Deciphering the genetic architecture of human cardiac disorders is of fundamental importance but their underlying complexity is a major hurdle. We investigated the natural variation of cardiac performance in the sequenced inbred lines of the Drosophila Genetic Reference Panel (DGRP). Genome-wide associations studies (GWAS) identified genetic networks associated with natural variation of cardiac traits which were used to gain insights as to the molecular and cellular processes affected. Non-coding variants that we identified were used to map potential regulatory non-coding regions, which in turn were employed to predict transcription factors (TFs) binding sites. Cognate TFs, many of which themselves bear polymorphisms associated with variations of cardiac performance, were also validated by heart-specific knockdown. Additionally, we showed that the natural variations associated with variability in cardiac performance affect a set of genes overlapping those associated with average traits but through different variants in the same genes. Furthermore, we showed that phenotypic variability was also associated with natural variation of gene regulatory networks. More importantly, we documented correlations between genes associated with cardiac phenotypes in both flies and humans, which supports a conserved genetic architecture regulating adult cardiac function from arthropods to mammals. Specifically, roles for PAX9 and EGR2 in the regulation of the cardiac rhythm were established in both models, illustrating that the characteristics of natural variations in cardiac function identified in Drosophila can accelerate discovery in humans.

Nascent polypeptide-Associated Complex and Signal Recognition Particle have cardiac-specific roles in heart development and remodeling.

Establishing a catalog of Congenital Heart Disease (CHD) genes and identifying functional networks would improve our understanding of its oligogenic underpinnings. Our studies identified protein biogenesis cofactors Nascent polypeptide-Associated Complex (NAC) and Signal-Recognition-Particle (SRP) as disease candidates and novel regulators of cardiac differentiation and morphogenesis. Knockdown (KD) of the alpha- (Nacα) or beta-subunit (bicaudal, bic) of NAC in the developing Drosophila heart disrupted cardiac developmental remodeling resulting in a fly with no heart. Heart loss was rescued by combined KD of Nacα with the posterior patterning Hox gene Abd-B. Consistent with a central role for this interaction in cardiogenesis, KD of Nacα in cardiac progenitors derived from human iPSCs impaired cardiac differentiation while co-KD with human HOXC12 and HOXD12 rescued this phenotype. Our data suggest that Nacα KD preprograms cardioblasts in the embryo for abortive remodeling later during metamorphosis, as Nacα KD during translation-intensive larval growth or pupal remodeling only causes moderate heart defects. KD of SRP subunits in the developing fly heart produced phenotypes that targeted specific segments and cell types, again suggesting cardiac-specific and spatially regulated activities. Together, we demonstrated directed function for NAC and SRP in heart development, and that regulation of NAC function depends on Hox genes.

Atrial Fibrillation Genomics: Discovery and Translation.

PURPOSE OF REVIEW: Our understanding of the fundamental cellular and molecular factors leading to atrial fibrillation (AF) remains stagnant despite significant advancement in ablation and device technologies. Diagnosis and prevention strategies fall behind that of treatment, but expanding knowledge in AF genetics holds the potential to drive progress. We aim to review how an understanding of the genetic contributions to AF can guide an approach to individualized risk stratification and novel avenues in drug discovery. RECENT FINDINGS: Rare familial forms of AF identified monogenic contributions to the development of AF. Genome-wide association studies (GWAS) further identified single-nucleotide polymorphisms (SNPs) suggesting polygenic and multiplex nature of this common disease. Polygenic risk scores accounting for the multitude of associated SNPs that each confer mildly elevated risk have been developed to translate genetic information into clinical practice, though shortcomings remain. Additionally, novel laboratory methods have been empowered by recent genetic findings to enhance drug discovery efforts. AF is increasingly recognized as a disease with a significant genetic component. With expanding sequencing technologies and accessibility, polygenic risk scores can help identify high risk individuals. Advancement in digital health tools, artificial intelligence and machine learning based on standard electrocardiograms, and genomic driven drug discovery may be integrated to deliver a sophisticated level of precision medicine in this modern era of emphasis on prevention. Randomized, prospective studies to demonstrate clinical benefits of these available tools are needed to validate this approach.

PGC1/PPAR drive cardiomyocyte maturation at single cell level via YAP1 and SF3B2.

Cardiomyocytes undergo significant structural and functional changes after birth, and these fundamental processes are essential for the heart to pump blood to the growing body. However, due to the challenges of isolating single postnatal/adult myocytes, how individual newborn cardiomyocytes acquire multiple aspects of the mature phenotype remains poorly understood. Here we implement large-particle sorting and analyze single myocytes from neonatal to adult hearts. Early myocytes exhibit wide-ranging transcriptomic and size heterogeneity that is maintained until adulthood with a continuous transcriptomic shift. Gene regulatory network analysis followed by mosaic gene deletion reveals that peroxisome proliferator-activated receptor coactivator-1 signaling, which is active in vivo but inactive in pluripotent stem cell-derived cardiomyocytes, mediates the shift. This signaling simultaneously regulates key aspects of cardiomyocyte maturation through previously unrecognized proteins, including YAP1 and SF3B2. Our study provides a single-cell roadmap of heterogeneous transitions coupled to cellular features and identifies a multifaceted regulator controlling cardiomyocyte maturation.

Small-molecule probe reveals a kinase cascade that links stress signaling to TCF/LEF and Wnt responsiveness.

Wnt signaling plays a central role in tissue maintenance and cancer. Wnt activates downstream genes through ß-catenin, which interacts with TCF/LEF transcription factors. A major question is how this signaling is coordinated relative to tissue organization and renewal. We used a recently described class of small molecules that binds tubulin to reveal a molecular cascade linking stress signaling through ATM, HIPK2, and p53 to the regulation of TCF/LEF transcriptional activity. These data suggest a mechanism by which mitotic and genotoxic stress can indirectly modulate Wnt responsiveness to exert coherent control over cell shape and renewal. These findings have implications for understanding tissue morphogenesis and small-molecule anticancer therapeutics.

Patient-specific genomics and cross-species functional analysis implicate LRP2 in hypoplastic left heart syndrome.

Congenital heart diseases (CHDs), including hypoplastic left heart syndrome (HLHS), are genetically complex and poorly understood. Here, a multidisciplinary platform was established to functionally evaluate novel CHD gene candidates, based on whole-genome and iPSC RNA sequencing of a HLHS family-trio. Filtering for rare variants and altered expression in proband iPSCs prioritized 10 candidates. siRNA/RNAi-mediated knockdown in healthy human iPSC-derived cardiomyocytes (hiPSC-CM) and in developing Drosophila and zebrafish hearts revealed that LDL receptor-related protein LRP2 is required for cardiomyocyte proliferation and differentiation. Consistent with hypoplastic heart defects, compared to patents the proband's iPSC-CMs exhibited reduced proliferation. Interestingly, rare, predicted-damaging LRP2 variants were enriched in a HLHS cohort; however, understanding their contribution to HLHS requires further investigation. Collectively, we have established a multi-species high-throughput platform to rapidly evaluate candidate genes and their interactions during heart development, which are crucial first steps toward deciphering oligogenic underpinnings of CHDs, including hypoplastic left hearts.

Dietary Emulsifier Sodium Stearoyl Lactylate Alters Gut Microbiota in vitro and Inhibits Bacterial Butyrate Producers.

Dietary emulsifiers are widely used in industrially processed foods, although the effects of these food additives on human gut microbiota are not well studied. Here, we investigated the effects of five different emulsifiers [glycerol monoacetate, glycerol monostearate, glycerol monooleate, propylene glycol monostearate, and sodium stearoyl lactylate (SSL)] on fecal microbiota in vitro. We found that 0.025% (w/v) of SSL reduced the relative abundance of the bacterial class Clostridia and others. The relative abundance of the families Clostridiaceae, Lachnospiraceae, and Ruminococcaceae was substantially reduced whereas that of Bacteroidaceae and Enterobacteriaceae was increased. Given the marked impact of SSL on Clostridia, we used genome reconstruction to predict community-wide production of short-chain fatty acids, which were experimentally assessed by GC-MS analysis. SSL significantly reduced concentrations of butyrate, and increased concentrations of propionate compared to control cultures. The presence of SSL increased lipopolysaccharide, LPS and flagellin in cultured communities, thereby enhancing the proinflammatory potential of SSL-selected bacterial communities.

Silencing of CCR4-NOT complex subunits affects heart structure and function.

The identification of genetic variants that predispose individuals to cardiovascular disease and a better understanding of their targets would be highly advantageous. Genome-wide association studies have identified variants that associate with QT-interval length (a measure of myocardial repolarization). Three of the strongest associating variants (single-nucleotide polymorphisms) are located in the putative promotor region of CNOT1, a gene encoding the central CNOT1 subunit of CCR4-NOT: a multifunctional, conserved complex regulating gene expression and mRNA stability and turnover. We isolated the minimum fragment of the CNOT1 promoter containing all three variants from individuals homozygous for the QT risk alleles and demonstrated that the haplotype associating with longer QT interval caused reduced reporter expression in a cardiac cell line, suggesting that reduced CNOT1 expression might contribute to abnormal QT intervals. Systematic siRNA-mediated knockdown of CCR4-NOT components in human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) revealed that silencing CNOT1 and other CCR4-NOT genes reduced their proliferative capacity. Silencing CNOT7 also shortened action potential duration. Furthermore, the cardiac-specific knockdown of Drosophila orthologs of CCR4-NOT genes in vivo (CNOT1/Not1 and CNOT7/8/Pop2) was either lethal or resulted in dilated cardiomyopathy, reduced contractility or a propensity for arrhythmia. Silencing CNOT2/Not2, CNOT4/Not4 and CNOT6/6L/twin also affected cardiac chamber size and contractility. Developmental studies suggested that CNOT1/Not1 and CNOT7/8/Pop2 are required during cardiac remodeling from larval to adult stages. To summarize, we have demonstrated how disease-associated genes identified by GWAS can be investigated by combining human cardiomyocyte cell-based and whole-organism in vivo heart models. Our results also suggest a potential link of CNOT1 and CNOT7/8 to QT alterations and further establish a crucial role of the CCR4-NOT complex in heart development and function.This article has an associated First Person interview with the first author of the paper.

AlleleProfileR: A versatile tool to identify and profile sequence variants in edited genomes.

Gene editing strategies, such as zinc-finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and clustered regularly interspaced short palindromic repeat/Cas9 (CRISPR/Cas9), are revolutionizing biology. However, quantitative and sensitive detection of targeted mutations are required to evaluate and quantify the genome editing outcomes. Here we present AlleleProfileR, a new analysis tool, written in a combination of R and C++, with the ability to batch process the sequence analysis of large and complex genome editing experiments, including the recently developed base editing technologies.

Model system identification of novel congenital heart disease gene candidates: focus on RPL13.

Genetics is a significant factor contributing to congenital heart disease (CHD), but our understanding of the genetic players and networks involved in CHD pathogenesis is limited. Here, we searched for de novo copy number variations (CNVs) in a cohort of 167 CHD patients to identify DNA segments containing potential pathogenic genes. Our search focused on new candidate disease genes within 19 deleted de novo CNVs, which did not cover known CHD genes. For this study, we developed an integrated high-throughput phenotypical platform to probe for defects in cardiogenesis and cardiac output in human induced pluripotent stem cell (iPSC)-derived multipotent cardiac progenitor (MCPs) cells and, in parallel, in the Drosophila in vivo heart model. Notably, knockdown (KD) in MCPs of RPL13, a ribosomal gene and SON, an RNA splicing cofactor, reduced proliferation and differentiation of cardiomyocytes, while increasing fibroblasts. In the fly, heart-specific RpL13 KD, predominantly at embryonic stages, resulted in a striking 'no heart' phenotype. KD of Son and Pdss2, among others, caused structural and functional defects, including reduced or abolished contractility, respectively. In summary, using a combination of human genetics and cardiac model systems, we identified new genes as candidates for causing human CHD, with particular emphasis on ribosomal genes, such as RPL13. This powerful, novel approach of combining cardiac phenotyping in human MCPs and in the in vivo Drosophila heart at high throughput will allow for testing large numbers of CHD candidates, based on patient genomic data, and for building upon existing genetic networks involved in heart development and disease.

Generation of First Heart Field-like Cardiac Progenitors and Ventricular-like Cardiomyocytes from Human Pluripotent Stem Cells.

The generation of large amounts of functional human pluripotent stem cells-derived cardiac progenitors and cardiomyocytes of defined heart field origin is a pre-requisite for cell-based cardiac therapies and disease modeling. We have recently shown that Id genes are both necessary and sufficient to specify first heart field progenitors during vertebrate development. This differentiation protocol leverages these findings and uses Id1 overexpression in combination with Activin A as potent specifying cues to produce first heart field-like (FHF-L) progenitors. Importantly, resulting progenitors efficiently differentiate (~70-90%) into ventricular-like cardiomyocytes. Here we describe a detailed method to 1) generate Id1-overexpressing hPSCs and 2) differentiate scalable quantities of cryopreservable FHF-L progenitors and ventricular-like cardiomyocytes.

Acellular therapeutic approach for heart failure: in vitro production of extracellular vesicles from human cardiovascular progenitors.

Aims: We have shown that extracellular vesicles (EVs) secreted by embryonic stem cell-derived cardiovascular progenitor cells (Pg) recapitulate the therapeutic effects of their parent cells in a mouse model of chronic heart failure (CHF). Our objectives are to investigate whether EV released by more readily available cell sources are therapeutic, whether their effectiveness is influenced by the differentiation state of the secreting cell, and through which mechanisms they act. Methods and results: The total EV secreted by human induced pluripotent stem cell-derived cardiovascular progenitors (iPSC-Pg) and human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CM) were isolated by ultracentrifugation and characterized by Nanoparticle Tracking Analysis, western blot, and cryo-electron microscopy. In vitro bioactivity assays were used to evaluate their cellular effects. Cell and EV microRNA (miRNA) content were assessed by miRNA array. Myocardial infarction was induced in 199 nude mice. Three weeks later, mice with left ventricular ejection fraction (LVEF) ≤ 45% received transcutaneous echo-guided injections of iPSC-CM (1.4 × 106, n = 19), iPSC-Pg (1.4 × 106, n = 17), total EV secreted by 1.4 × 106 iPSC-Pg (n = 19), or phosphate-buffered saline (control, n = 17) into the peri-infarct myocardium. Seven weeks later, hearts were evaluated by echocardiography, histology, and gene expression profiling, blinded to treatment group. In vitro, EV were internalized by target cells, increased cell survival, cell proliferation, and endothelial cell migration in a dose-dependent manner and stimulated tube formation. Extracellular vesicles were rich in miRNAs and most of the 16 highly abundant, evolutionarily conserved miRNAs are associated with tissue-repair pathways. In vivo, EV outperformed cell injections, significantly improving cardiac function through decreased left ventricular volumes (left ventricular end systolic volume: -11%, P < 0.001; left ventricular end diastolic volume: -4%, P = 0.002), and increased LVEF (+14%, P < 0.0001) relative to baseline values. Gene profiling revealed that EV-treated hearts were enriched for tissue reparative pathways. Conclusion: Extracellular vesicles secreted by iPSC-Pg are effective in the treatment of CHF, possibly, in part, through their specific miRNA signature and the associated stimulation of distinct cardioprotective pathways. The processing and regulatory advantages of EV could make them effective substitutes for cell transplantation.

An Automated Platform for Assessment of Congenital and Drug-Induced Arrhythmia with hiPSC-Derived Cardiomyocytes.

The ability to produce unlimited numbers of human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs) harboring disease and patient-specific gene variants creates a new paradigm for modeling congenital heart diseases (CHDs) and predicting proarrhythmic liabilities of drug candidates. However, a major roadblock to implementing hiPSC-CM technology in drug discovery is that conventional methods for monitoring action potential (AP) kinetics and arrhythmia phenotypes in vitro have been too costly or technically challenging to execute in high throughput. Herein, we describe the first large-scale, fully automated and statistically robust analysis of AP kinetics and drug-induced proarrhythmia in hiPSC-CMs. The platform combines the optical recording of a small molecule fluorescent voltage sensing probe (VoltageFluor2.1.Cl), an automated high throughput microscope and automated image analysis to rapidly generate physiological measurements of cardiomyocytes (CMs). The technique can be readily adapted on any high content imager to study hiPSC-CM physiology and predict the proarrhythmic effects of drug candidates.

Id genes are essential for early heart formation.

Deciphering the fundamental mechanisms controlling cardiac specification is critical for our understanding of how heart formation is initiated during embryonic development and for applying stem cell biology to regenerative medicine and disease modeling. Using systematic and unbiased functional screening approaches, we discovered that the Id family of helix-loop-helix proteins is both necessary and sufficient to direct cardiac mesoderm formation in frog embryos and human embryonic stem cells. Mechanistically, Id proteins specify cardiac cell fate by repressing two inhibitors of cardiogenic mesoderm formation-Tcf3 and Foxa2-and activating inducers Evx1, Grrp1, and Mesp1. Most importantly, CRISPR/Cas9-mediated ablation of the entire Id (Id1-4) family in mouse embryos leads to failure of anterior cardiac progenitor specification and the development of heartless embryos. Thus, Id proteins play a central and evolutionarily conserved role during heart formation and provide a novel means to efficiently produce cardiovascular progenitors for regenerative medicine and drug discovery applications.

Early molecular events during retinoic acid induced differentiation of neuromesodermal progenitors.

Bipotent neuromesodermal progenitors (NMPs) residing in the caudal epiblast drive coordinated body axis extension by generating both posterior neuroectoderm and presomitic mesoderm. Retinoic acid (RA) is required for body axis extension, however the early molecular response to RA signaling is poorly defined, as is its relationship to NMP biology. As endogenous RA is first seen near the time when NMPs appear, we used WNT/FGF agonists to differentiate embryonic stem cells to NMPs which were then treated with a short 2-h pulse of 25 nM RA or 1 µM RA followed by RNA-seq transcriptome analysis. Differential expression analysis of this dataset indicated that treatment with 25 nM RA, but not 1 µM RA, provided physiologically relevant findings. The 25 nM RA dataset yielded a cohort of previously known caudal RA target genes including Fgf8 (repressed) and Sox2 (activated), plus novel early RA signaling targets with nearby conserved RA response elements. Importantly, validation of top-ranked genes in vivo using RA-deficient Raldh2-/- embryos identified novel examples of RA activation (Nkx1-2, Zfp503, Zfp703, Gbx2, Fgf15, Nt5e) or RA repression (Id1) of genes expressed in the NMP niche or progeny. These findings provide evidence for early instructive and permissive roles of RA in controlling differentiation of NMPs to neural and mesodermal lineages.