Neighbor-specific gene expression revealed from physically interacting cells during mouse embryonic development.

Development of multicellular organisms is orchestrated by persistent cell-cell communication between neighboring partners. Direct interaction between different cell types can induce molecular signals that dictate lineage specification and cell fate decisions. Current single-cell RNA-seq technology cannot adequately analyze cell-cell contact-dependent gene expression, mainly due to the loss of spatial information. To overcome this obstacle and resolve cell-cell contact-specific gene expression during embryogenesis, we performed RNA sequencing of physically interacting cells (PIC-seq) and assessed them alongside similar single-cell transcriptomes derived from developing mouse embryos between embryonic day (E) 7.5 and E9.5. Analysis of the PIC-seq data identified gene expression signatures that were dependent on the presence of specific neighboring cell types. Our computational predictions, validated experimentally, demonstrated that neural progenitor (NP) cells upregulate Lhx5 and Nkx2-1 genes, when exclusively interacting with definitive endoderm (DE) cells. Moreover, there was a reciprocal impact on the transcriptome of DE cells, as they tend to upregulate Rax and Gsc when in contact with NP cells. Using individual cell transcriptome data, we formulated a means of computationally predicting the impact of one cell type on the transcriptome of its neighboring cell types. We have further developed a distinctive spatial-t-distributed stochastic neighboring embedding to display the pseudospatial distribution of cells in a 2-dimensional space. In summary, we describe an innovative approach to study contact-specific gene regulation during embryogenesis.

VeTra: a tool for trajectory inference based on RNA velocity.

MOTIVATION: Trajectory inference (TI) for single cell RNA sequencing (scRNAseq) data is a powerful approach to interpret dynamic cellular processes such as cell cycle and development. Still, however, accurate inference of trajectory is challenging. Recent development of RNA velocity provides an approach to visualize cell state transition without relying on prior knowledge. RESULTS: To perform TI and group cells based on RNA velocity we developed VeTra. By applying cosine similarity and merging weakly connected components, VeTra identifies cell groups from the direction of cell transition. Besides, VeTra suggests key regulators from the inferred trajectory. VeTra is a useful tool for TI and subsequent analysis. AVAILABILITY AND IMPLEMENTATION: The Vetra is available at https://github.com/wgzgithub/VeTra. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Cardiac fibroblast sub-types in vitro reflect pathological cardiac remodeling in vivo.

Many heart diseases are associated with fibrosis, but it is unclear whether different types of heart disease correlate with different subtypes of activated fibroblasts and to which extent such diversity is modeled during in vitro activation of primary cardiac fibroblasts. Analyzing the expression of 82 fibrosis related genes in 65 heart failure (HF) patients, we identified a panel of 12 genes clearly distinguishing HF patients better from healthy controls than measurement of the collagen-related hydroxyproline content. A subcluster enriched in ischemic HF was recognized, but not for diabetes, high BMI, or gender. Single-cell transcriptomic analysis of in vitro activated mouse cardiac fibroblasts distinguished 6 subpopulations, including a contractile Acta2high precursor population, which was predicted by time trajectory analysis to develop into Acta2low subpopulations with high production of extracellular matrix molecules. The 12 gene profile identified in HF patients showed highest similarity to the fibroblast subset with the strongest expression of extracellular matrix molecules. Population markers identified were furthermore able to clearly cluster different disease stages in a murine model for myocardial infarct. These data suggest that major features of cardiac fibroblast activation in heart failure patients, in murine heart disease models, and in cell culture of primary murine cardiac fibroblast are shared.