Revealing Dynamic Mechanisms of Cell Fate Decisions From Single-Cell Transcriptomic Data.

Cell fate decisions play a pivotal role in development, but technologies for dissecting them are limited. We developed a multifunction new method, Topographer, to construct a "quantitative" Waddington's landscape of single-cell transcriptomic data. This method is able to identify complex cell-state transition trajectories and to estimate complex cell-type dynamics characterized by fate and transition probabilities. It also infers both marker gene networks and their dynamic changes as well as dynamic characteristics of transcriptional bursting along the cell-state transition trajectories. Applying this method to single-cell RNA-seq data on the differentiation of primary human myoblasts, we not only identified three known cell types, but also estimated both their fate probabilities and transition probabilities among them. We found that the percent of genes expressed in a bursty manner is significantly higher at (or near) the branch point (~97%) than before or after branch (below 80%), and that both gene-gene and cell-cell correlation degrees are apparently lower near the branch point than away from the branching. Topographer allows revealing of cell fate mechanisms in a coherent way at three scales: cell lineage (macroscopic), gene network (mesoscopic), and gene expression (microscopic)

Simultaneous profiling of transcriptome and DNA methylome from a single cell.

BackgroundSingle-cell transcriptome and single-cell methylome technologies have become powerful tools to study RNA and DNA methylation profiles of single cells at a genome-wide scale. A major challenge has been to understand the direct correlation of DNA methylation and gene expression within single-cells. Due to large cell-to-cell variability and the lack of direct measurements of transcriptome and methylome of the same cell, the association is still unclear.ResultsHere, we describe a novel method (scMT-seq) that simultaneously profiles both DNA methylome and transcriptome from the same cell. In sensory neurons, we consistently identify transcriptome and methylome heterogeneity among single cells but the majority of the expression variance is not explained by proximal promoter methylation, with the exception of genes that do not contain CpG islands. By contrast, gene body methylation is positively associated with gene expression for only those genes that contain a CpG island promoter. Furthermore, using single nucleotide polymorphism patterns from our hybrid mouse model, we also find positive correlation of allelic gene body methylation with allelic expression.ConclusionsOur method can be used to detect transcriptome, methylome, and single nucleotide polymorphism information within single cells to dissect the mechanisms of epigenetic gene regulation

Information visualization for DNA microarray data analysis: A critical review

Graphical representation may provide effective means of making sense of the complexity and sheer volume of data produced by DNA microarray experiments that monitor the expression patterns of thousands of genes simultaneously. The ability to use ldquoabstractrdquo graphical representation to draw attention to areas of interest, and more in-depth visualizations to answer focused questions, would enable biologists to move from a large amount of data to particular records they are interested in, and therefore, gain deeper insights in understanding the microarray experiment results. This paper starts by providing some background knowledge of microarray experiments, and then, explains how graphical representation can be applied in general to this problem domain, followed by exploring the role of visualization in gene expression data analysis. Having set the problem scene, the paper then examines various multivariate data visualization techniques that have been applied to microarray data analysis. These techniques are critically reviewed so that the strengths and weaknesses of each technique can be tabulated. Finally, several key problem areas as well as possible solutions to them are discussed as being a source for future work

Subtle changes in chromatin loop contact propensity are associated with differential gene regulation and expression.

While genetic variation at chromatin loops is relevant for human disease, the relationships between contact propensity (the probability that loci at loops physically interact), genetics, and gene regulation are unclear. We quantitatively interrogate these relationships by comparing Hi-C and molecular phenotype data across cell types and haplotypes. While chromatin loops consistently form across different cell types, they have subtle quantitative differences in contact frequency that are associated with larger changes in gene expression and H3K27ac. For the vast majority of loci with quantitative differences in contact frequency across haplotypes, the changes in magnitude are smaller than those across cell types; however, the proportional relationships between contact propensity, gene expression, and H3K27ac are consistent. These findings suggest that subtle changes in contact propensity have a biologically meaningful role in gene regulation and could be a mechanism by which regulatory genetic variants in loop anchors mediate effects on expression

Predicting gene expression in the human malaria parasite Plasmodium falciparum using histone modification, nucleosome positioning, and 3D localization features.

Empirical evidence suggests that the malaria parasite Plasmodium falciparum employs a broad range of mechanisms to regulate gene transcription throughout the organism's complex life cycle. To better understand this regulatory machinery, we assembled a rich collection of genomic and epigenomic data sets, including information about transcription factor (TF) binding motifs, patterns of covalent histone modifications, nucleosome occupancy, GC content, and global 3D genome architecture. We used these data to train machine learning models to discriminate between high-expression and low-expression genes, focusing on three distinct stages of the red blood cell phase of the Plasmodium life cycle. Our results highlight the importance of histone modifications and 3D chromatin architecture in Plasmodium transcriptional regulation and suggest that AP2 transcription factors may play a limited regulatory role, perhaps operating in conjunction with epigenetic factors

Mammalian gene expression variability is explained by underlying cell state.

Gene expression variability in mammalian systems plays an important role in physiological and pathophysiological conditions. This variability can come from differential regulation related to cell state (extrinsic) and allele-specific transcriptional bursting (intrinsic). Yet, the relative contribution of these two distinct sources is unknown. Here, we exploit the qualitative difference in the patterns of covariance between these two sources to quantify their relative contributions to expression variance in mammalian cells. Using multiplexed error robust RNA fluorescent in situ hybridization (MERFISH), we measured the multivariate gene expression distribution of 150 genes related to Ca2+ signaling coupled with the dynamic Ca2+ response of live cells to ATP. We show that after controlling for cellular phenotypic states such as size, cell cycle stage, and Ca2+ response to ATP, the remaining variability is effectively at the Poisson limit for most genes. These findings demonstrate that the majority of expression variability results from cell state differences and that the contribution of transcriptional bursting is relatively minimal

Transcriptomic profiling of tumor-infiltrating CD4 + TIM-3 + T Cells reveals their suppressive, exhausted, and metastatic characteristics in colorectal cancer patients

T cell immunoglobulin mucin-3 (TIM-3) is an immune checkpoint identified as one of the key players in regulating T-cell responses. Studies have shown that TIM-3 is upregulated in the tumor microenvironment (TME). However, the precise role of TIM-3 in colorectal cancer (CRC) TME is yet to be elucidated. We performed phenotypic and molecular characterization of TIM-3+ T cells in the TME and circulation of CRC patients by analyzing tumor tissues (TT, TILs), normal tissues (NT, NILs), and peripheral blood mononuclear cells (PBMC). TIM-3 was upregulated on both CD4+ and CD3+CD4− (CD8+) TILs. CD4+TIM-3+ TILs expressed higher levels of T regulatory cell (Tregs)-signature genes, including FoxP3 and Helios, compared with their TIM-3− counterparts. Transcriptomic and ingenuity pathway analyses showed that TIM-3 potentially activates inflammatory and tumor metastatic pathways. Moreover, NF-κB-mediated transcription factors were upregulated in CD4+TIM-3+ TILs, which could favor proliferation/invasion and induce inflammatory and T-cell exhaustion pathways. In addition, we found that CD4+TIM-3+ TILs potentially support tumor invasion and metastasis, compared with conventional CD4+CD25+ Tregs in the CRC TME. However, functional studies are warranted to support these findings. In conclusion, this study discloses some of the functional pathways of TIM-3+ TILs, which could improve their targeting in more specific therapeutic approaches in CRC patients

Papillary Thyroid Carcinoma Variants are Characterized by Co-dysregulation of Immune and Cancer Associated Genes.

Papillary thyroid carcinoma (PTC) variants exhibit different prognosis, but critical characteristics of PTC variants that contribute to differences in pathogenesis are not well-known. This study aims to characterize dysregulated immune-associated and cancer-associated genes in three PTC subtypes to explore how the interplay between cancer and immune processes causes differential prognosis. RNA-sequencing data from The Cancer Genome Atlas (TCGA) were used to identify dysregulated genes in each variant. The dysregulation profiles of the subtypes were compared using functional pathways clustering and correlations to relevant clinical variables, genomic alterations, and microRNA regulation. We discovered that the dysregulation profiles of classical PTC (CPTC) and the tall cell variant (TCPTC) are similar and are distinct from that of the follicular variant (FVPTC). However, unique cancer or immune-associated genes are associated with clinical variables for each subtype. Cancer-related genes MUC1, FN1, and S100-family members were the most clinically relevant in CPTC, while APLN and IL16, both immune-related, were clinically relevant in FVPTC. RAET-family members, also immune-related, were clinically relevant in TCPTC. Collectively, our data suggest that dysregulation of both cancer and immune associated genes defines the gene expression landscapes of PTC variants, but different cancer or immune related genes may drive the phenotype of each variant