Mapping the chromatin occupancy of transcription factors (TFs) is a key step in deciphering developmental transcriptional programs. Here we use biotinylated knockin alleles of seven key cardiac TFs (GATA4, NKX2-5, MEF2A, MEF2C, SRF, TBX5, TEAD1) to sensitively and reproducibly map their genome-wide occupancy in the fetal and adult mouse heart. These maps show that TF occupancy is dynamic between developmental stages and that multiple TFs often collaboratively occupy the same chromatin region through indirect cooperativity. Multi-TF regions exhibit features of functional regulatory elements, including evolutionary conservation, chromatin accessibility, and activity in transcriptional enhancer assays. H3K27ac, a feature of many enhancers, incompletely overlaps multi-TF regions, and multi-TF regions lacking H3K27ac retain conservation and enhancer activity. TEAD1 is a core component of the cardiac transcriptional network, co-occupying cardiac regulatory regions and controlling cardiomyocyte-specific gene functions. Our study provides a resource for deciphering the cardiac transcriptional regulatory network and gaining insights into the molecular mechanisms governing heart development.

Piwi-interacting RNAs are small regulatory RNAs with key roles in transposon silencing and regulation of gametogenesis. The production of mature piwi-interacting RNAs requires a critical step of trimming piwi-interacting RNA intermediates to achieve optimally sized piwi-interacting RNAs. The poly(A)-specific ribonuclease family deadenylase PNLDC1 is implicated in piwi-interacting RNA trimming in silkworms. The physiological function of PNLDC1 in mammals remains unknown. Using Pnldc1-deficient mice, here we show that PNLDC1 is required for piwi-interacting RNA biogenesis, transposon silencing, and spermatogenesis. Pnldc1 mutation in mice inhibits piwi-interacting RNA trimming and causes accumulation of untrimmed piwi-interacting RNA intermediates with 3′ end extension, leading to severe reduction of mature piwi-interacting RNAs in the testis. Pnldc1 mutant mice exhibit disrupted LINE1 retrotransposon silencing and defect in spermiogenesis. Together, these results define PNLDC1 as a mammalian piwi-interacting RNA biogenesis factor that protects the germline genome and ensures normal sperm production in mice.

Significance: Many occlusive vascular diseases in humans are largely dependent upon vascular smooth muscle cell (VSMC) phenotypic switching from a contractile to a proliferative phenotype, contributing to the formation of intimal lesions that eventually block the blood flow. Previous studies showed that the long noncoding RNA (lncRNA) NEAT1 is critical for tumorigenesis. In this report, we showed that NEAT1 expression was not only induced in VSMCs during phenotypic switching but functionally was critical for the smooth muscle phenotypic change. Our study demonstrates an unexpected role of the lncRNA NEAT1 in VSMCs and suggests that NEAT1 is a novel therapeutic target for treating occlusive vascular diseases in humans.

Background: Vascular homeostasis is maintained by the differentiated phenotype of vascular smooth muscle cells (VSMCs). The landscape of protein coding genes comprising the transcriptome of differentiated VSMCs has been intensively investigated but many gaps remain including the emerging roles of noncoding genes. Methods: We reanalyzed large-scale, publicly available bulk and single-cell RNA sequencing datasets from multiple tissues and cell types to identify VSMC-enriched long noncoding RNAs. The in vivo expression pattern of a novel smooth muscle cell (SMC)–expressed long noncoding RNA, Carmn (cardiac mesoderm enhancer-associated noncoding RNA), was investigated using a novel Carmn green fluorescent protein knock-in reporter mouse model. Bioinformatics and quantitative real-time polymerase chain reaction analysis were used to assess CARMN expression changes during VSMC phenotypic modulation in human and murine vascular disease models. In vitro, functional assays were performed by knocking down CARMN with antisense oligonucleotides and overexpressing Carmn by adenovirus in human coronary artery SMCs. Carotid artery injury was performed in SMC-specific Carmn knockout mice to assess neointima formation and the therapeutic potential of reversing CARMN loss was tested in a rat carotid artery balloon injury model. The molecular mechanisms underlying CARMN function were investigated using RNA pull-down, RNA immunoprecipitation, and luciferase reporter assays. Results: We identified CARMN , which was initially annotated as the host gene of the MIR143/145 cluster and recently reported to play a role in cardiac differentiation, as a highly abundant and conserved, SMC-specific long noncoding RNA. Analysis of the Carmn GFP knock-in mouse model confirmed that Carmn is transiently expressed in embryonic cardiomyocytes and thereafter becomes restricted to SMCs. We also found that Carmn is transcribed independently of Mir143/145 . CARMN expression is dramatically decreased by vascular disease in humans and murine models and regulates the contractile phenotype of VSMCs in vitro. In vivo, SMC-specific deletion of Carmn significantly exacerbated, whereas overexpression of Carmn markedly attenuated, injury-induced neointima formation in mouse and rat, respectively. Mechanistically, we found that Carmn physically binds to the key transcriptional cofactor myocardin, facilitating its activity and thereby maintaining the contractile phenotype of VSMCs. Conclusions: CARMN is an evolutionarily conserved SMC-specific long noncoding RNA with a previously unappreciated role in maintaining the contractile phenotype of VSMCs and is the first noncoding RNA discovered to interact with myocardin.

The Tibetan cashmere goat (Capra hircus), one of the most ancient breeds in China, has historically been a critical source of meat and cashmere production for local farmers. To adapt to the high-altitude area, extremely harsh climate, and hypoxic environment that the Tibetan cashmere goat lives in, this goat has developed distinct phenotypic traits compared to lowland breeds. However, the genetic components underlying this phenotypic adaptation remain largely unknown. We obtained 118,700 autosomal SNPs through exome sequencing of 330 cashmere goats located at a wide geographic range, including the Tibetan Plateau and low-altitude regions in China. The great majority of SNPs showed low genetic differentiation among populations; however, approximately 2-3 % of the loci showed more genetic differentiation than expected under a selectively neutral model. Together with a combined analysis of high- and low-altitude breeds, we revealed 339 genes potentially under high-altitude selection. Genes associated with cardiovascular system development were significantly enriched in our study. Among these genes, the most evident one was endothelial PAS domain protein 1 (EPAS1), which has been previously reported to be involved in complex oxygen sensing and significantly associated with high-altitude adaptation of human, dog, and grey wolf. The missense mutation Q579L that we identified in EPAS1, which occurs next to the Hypoxia-Inducible Factor-1 (HIF-1) domain, was exclusively enriched in the high-altitude populations. Our study provides insights concerning the population variation in six different cashmere goat populations in China. The variants in cardiovascular system-related genes may explain the observed phenotypic adaptation of the Tibetan cashmere goat.

Background: Fat tail is a unique trait in sheep acquired during domestication. Several genomic analyses have been conducted in sheep breeds from limited geographic origins to identify the genetic factors underlying this trait. Nevertheless, these studies obtained different candidates. The results of these regional studies were easily biased by the breed structures. Results: To minimize the bias and distinguish the true candidates, we used an extended data set of 968 sheep representing 18 fat-tailed breeds and 14 thin-tailed breeds from around the world, and integrated two statistical tests to detect selection signatures, including Genetic Fixation Index (F_ST) and difference of derived allele frequency (ΔDAF). The results showed that platelet derived growth factor D (PDGFD) exhibited the highest genetic differentiation between fat- and thin-tailed sheep breeds. Analysis of sequence variation identified that a 6.8-kb region within the first intron of PDGFD is likely the target of positive selection and contains regulatory mutation(s) in fat-tailed sheep. Histological and gene expression analyses demonstrated that PDGFD expression is associated with maturation and hemostasis of adipocytes. Further retrospective analysis of public transcriptomic datasets revealed that PDGFD expression is down-regulated during adipogenesis in both human and mouse, and is higher in fat tissues of obese individuals than that in lean individuals. Conclusions: These results reveal that PDGFD is the predominant factor for the fat tail phenotype in sheep by contributing to adiopogenesis and maintaining the hemostasis of mature adipocytes. This study provides insights into the selection of fat-tailed sheep and has important application to animal breeding, as well as obesity-related human diseases.