Differential Impact In Vivo of Pf4-ΔCre-Mediated and Gp1ba-ΔCre-Mediated Depletion of Cyclooxygenase-1 in Platelets in Mice.

BACKGROUND: Low-dose aspirin is widely used for the secondary prevention of cardiovascular disease. The beneficial effects of low-dose aspirin are attributable to its inhibition of platelet Cox (cyclooxygenase)-1-derived thromboxane A2. Until recently, the use of the Pf4 (platelet factor 4) Cre has been the only genetic approach to generating megakaryocyte/platelet ablation of Cox-1 in mice. However, Pf4-ΔCre displays ectopic expression outside the megakaryocyte/platelet lineage, especially during inflammation. The use of the Gp1ba (glycoprotein 1bα) Cre promises a more specific, targeted approach. METHODS: To evaluate the role of Cox-1 in platelets, we crossed Pf4-ΔCre or Gp1ba-ΔCre mice with Cox-1flox/flox mice to generate platelet Cox-1-/- mice on normolipidemic and hyperlipidemic (Ldlr-/-; low-density lipoprotein receptor) backgrounds. RESULTS: Ex vivo platelet aggregation induced by arachidonic acid or adenosine diphosphate in platelet-rich plasma was inhibited to a similar extent in Pf4-ΔCre Cox-1-/-/Ldlr-/- and Gp1ba-ΔCre Cox-1-/-/Ldlr-/- mice. In a mouse model of tail injury, Pf4-ΔCre-mediated and Gp1ba-ΔCre-mediated deletions of Cox-1 were similarly efficient in suppressing platelet prostanoid biosynthesis. Experimental thrombogenesis and attendant blood loss were similar in both models. However, the impact on atherogenesis was divergent, being accelerated in the Pf4-ΔCre mice while restrained in the Gp1ba-ΔCres. In the former, accelerated atherogenesis was associated with greater suppression of PGI2 biosynthesis, a reduction in the lipopolysaccharide-evoked capacity to produce PGE2 (prostaglandin E) and PGD2 (prostanglandin D), activation of the inflammasome, elevated plasma levels of IL-1ß (interleukin), reduced plasma levels of HDL-C (high-density lipoprotein receptor-cholesterol), and a reduction in the capacity for reverse cholesterol transport. By contrast, in the latter, plasma HDL-C and α-tocopherol were elevated, and MIP-1α (macrophage inflammatory protein-1α) and MCP-1 (monocyte chemoattractant protein 1) were reduced. CONCLUSIONS: Both approaches to Cox-1 deletion similarly restrain thrombogenesis, but a differential impact on Cox-1-dependent prostanoid formation by the vasculature may contribute to an inflammatory phenotype and accelerated atherogenesis in Pf4-ΔCre mice.

Regularized sequence-context mutational trees capture variation in mutation rates across the human genome.

Germline mutation is the mechanism by which genetic variation in a population is created. Inferences derived from mutation rate models are fundamental to many population genetics methods. Previous models have demonstrated that nucleotides flanking polymorphic sites-the local sequence context-explain variation in the probability that a site is polymorphic. However, limitations to these models exist as the size of the local sequence context window expands. These include a lack of robustness to data sparsity at typical sample sizes, lack of regularization to generate parsimonious models and lack of quantified uncertainty in estimated rates to facilitate comparison between models. To address these limitations, we developed Baymer, a regularized Bayesian hierarchical tree model that captures the heterogeneous effect of sequence contexts on polymorphism probabilities. Baymer implements an adaptive Metropolis-within-Gibbs Markov Chain Monte Carlo sampling scheme to estimate the posterior distributions of sequence-context based probabilities that a site is polymorphic. We show that Baymer accurately infers polymorphism probabilities and well-calibrated posterior distributions, robustly handles data sparsity, appropriately regularizes to return parsimonious models, and scales computationally at least up to 9-mer context windows. We demonstrate application of Baymer in three ways-first, identifying differences in polymorphism probabilities between continental populations in the 1000 Genomes Phase 3 dataset, second, in a sparse data setting to examine the use of polymorphism models as a proxy for de novo mutation probabilities as a function of variant age, sequence context window size, and demographic history, and third, comparing model concordance between different great ape species. We find a shared context-dependent mutation rate architecture underlying our models, enabling a transfer-learning inspired strategy for modeling germline mutations. In summary, Baymer is an accurate polymorphism probability estimation algorithm that automatically adapts to data sparsity at different sequence context levels, thereby making efficient use of the available data.

Tempo: an unsupervised Bayesian algorithm for circadian phase inference in single-cell transcriptomics.

The circadian clock is a 24 h cellular timekeeping mechanism that regulates human physiology. Answering several fundamental questions in circadian biology will require joint measures of single-cell circadian phases and transcriptomes. However, no widespread experimental approaches exist for this purpose. While computational approaches exist to infer cell phase directly from single-cell RNA-sequencing data, existing methods yield poor circadian phase estimates, and do not quantify estimation uncertainty, which is essential for interpretation of results from very sparse single-cell RNA-sequencing data. To address these unmet needs, we introduce Tempo, a Bayesian variational inference approach that incorporates domain knowledge of the clock and quantifies phase estimation uncertainty. Through simulations and analyses of real data, we demonstrate that Tempo yields more accurate estimates of circadian phase than existing methods and provides well-calibrated uncertainty quantifications. Tempo will facilitate large-scale studies of single-cell circadian transcription.

Applications of single-cell genomics and computational strategies to study common disease and population-level variation.

The advent and rapid development of single-cell technologies have made it possible to study cellular heterogeneity at an unprecedented resolution and scale. Cellular heterogeneity underlies phenotypic differences among individuals, and studying cellular heterogeneity is an important step toward our understanding of the disease molecular mechanism. Single-cell technologies offer opportunities to characterize cellular heterogeneity from different angles, but how to link cellular heterogeneity with disease phenotypes requires careful computational analysis. In this article, we will review the current applications of single-cell methods in human disease studies and describe what we have learned so far from existing studies about human genetic variation. As single-cell technologies are becoming widely applicable in human disease studies, population-level studies have become a reality. We will describe how we should go about pursuing and designing these studies, particularly how to select study subjects, how to determine the number of cells to sequence per subject, and the needed sequencing depth per cell. We also discuss computational strategies for the analysis of single-cell data and describe how single-cell data can be integrated with bulk tissue data and data generated from genome-wide association studies. Finally, we point out open problems and future research directions.

Statistical and machine learning methods for spatially resolved transcriptomics with histology.

Recent developments in spatially resolved transcriptomics (SRT) technologies have enabled scientists to get an integrated understanding of cells in their morphological context. Applications of these technologies in diverse tissues and diseases have transformed our views of transcriptional complexity. Most published studies utilized tools developed for single-cell RNA sequencing (scRNA-seq) for data analysis. However, SRT data exhibit different properties from scRNA-seq. To take full advantage of the added dimension on spatial location information in such data, new methods that are tailored for SRT are needed. Additionally, SRT data often have companion high-resolution histology information available. Incorporating histological features in gene expression analysis is an underexplored area. In this review, we will focus on the statistical and machine learning aspects for SRT data analysis and discuss how spatial location and histology information can be integrated with gene expression to advance our understanding of the transcriptional complexity. We also point out open problems and future research directions in this field.

Single-Cell Genomics Reveals a Novel Cell State During Smooth Muscle Cell Phenotypic Switching and Potential Therapeutic Targets for Atherosclerosis in Mouse and Human.

BACKGROUND: Smooth muscle cells (SMCs) play significant roles in atherosclerosis via phenotypic switching, a pathological process in which SMC dedifferentiation, migration, and transdifferentiation into other cell types. Yet how SMCs contribute to the pathophysiology of atherosclerosis remains elusive. METHODS: To reveal the trajectories of SMC transdifferentiation during atherosclerosis and to identify molecular targets for disease therapy, we combined SMC fate mapping and single-cell RNA sequencing of both mouse and human atherosclerotic plaques. We also performed cell biology experiments on isolated SMC-derived cells, conducted integrative human genomics, and used pharmacological studies targeting SMC-derived cells both in vivo and in vitro. RESULTS: We found that SMCs transitioned to an intermediate cell state during atherosclerosis, which was also found in human atherosclerotic plaques of carotid and coronary arteries. SMC-derived intermediate cells, termed "SEM" cells (stem cell, endothelial cell, monocyte), were multipotent and could differentiate into macrophage-like and fibrochondrocyte-like cells, as well as return toward the SMC phenotype. Retinoic acid (RA) signaling was identified as a regulator of SMC to SEM cell transition, and RA signaling was dysregulated in symptomatic human atherosclerosis. Human genomics revealed enrichment of genome-wide association study signals for coronary artery disease in RA signaling target gene loci and correlation between coronary artery disease risk alleles and repressed expression of these genes. Activation of RA signaling by all-trans RA, an anticancer drug for acute promyelocytic leukemia, blocked SMC transition to SEM cells, reduced atherosclerotic burden, and promoted fibrous cap stability. CONCLUSIONS: Integration of cell-specific fate mapping, single-cell genomics, and human genetics adds novel insights into the complexity of SMC biology and reveals regulatory pathways for therapeutic targeting of SMC transitions in atherosclerotic cardiovascular disease.