Peripheral nerve injury leads to various injury-induced responses in sensory neurons including physiological pain, neuronal cell death, and nerve regeneration. In this study, we performed single-cell RNA-sequencing (scRNA-seq) analysis of mouse nonpeptidergic nociceptors (NP), peptidergic nociceptors (PEP), and large myelinated sensory neurons (LM) under both control and injury conditions at 3 days after sciatic nerve transection (SNT). After performing principle component and weighted gene co-expression network analysis, we categorized dorsal root ganglion (DRG) neurons into different subtypes and discovered co-regulated injury-response genes including novel regeneration associated genes (RAGs) in association with neuronal development, protein translation and cytoplasm transportation. In addition, we found significant up-regulation of the genes associated with cell death such as Pdcd2 in a subset of NP neurons after axotomy, implicating their actions in neuronal cell death upon nerve injury. Our study revealed the distinctive and sustained heterogeneity of transcriptomic responses to injury at single neuron level, implicating the involvement of different gene regulatory networks in nerve regeneration, neuronal cell death and neuropathy in different population of DRG neurons.

The development of the mammalian retina is a complicated process involving the generation of distinct types of neurons from retinal progenitor cells (RPCs) in a spatiotemporal-specific manner. The progression of RPCs during retinogenesis includes RPC proliferation, cell-fate commitment, and specific neuronal differentiation. In this study, by performing single-cell RNA sequencing of cells isolated from human embryonic stem cell (hESC)-derived 3D retinal organoids, we successfully deconstructed the temporal progression of RPCs during early human retinogenesis. We identified two distinctive subtypes of RPCs with unique molecular profiles, namely multipotent RPCs and neurogenic RPCs. We found that genes related to the Notch and Wnt signaling pathways, as well as chromatin remodeling, were dynamically regulated during RPC commitment. Interestingly, our analysis identified that CCND1, a G₁-phase cell-cycle regulator, was coexpressed with ASCL1 in a cell-cycle-independent manner. Temporally controlled overexpression of CCND1 in retinal organoids demonstrated a role for CCND1 in promoting early retinal neurogenesis. Together, our results revealed critical pathways and novel genes in early retinogenesis of humans.

Background

Single-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.

Results

Here, 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.

Conclusions

Our method can be used to detect transcriptome, methylome, and single nucleotide polymorphism information within single cells to dissect the mechanisms of epigenetic gene regulation.

Circulating cell-free DNAs (cfDNAs) are fragmented DNA molecules released into the blood by cells. Previous studies have suggested that mitochondria-originated cfDNA fragments (mt-cfDNAs) in cancer patients are more fragmented than those from healthy controls. However, it is still unknown where these short mt-cfDNAs originate, and whether the length of mt-cfDNAs can be correlated with tumor burden and cancer progression. In this study, we first performed whole-genome sequencing analysis (WGS) of cfDNAs from a human tumor cell line-xenotransplantation mouse model and found that mt-cfDNAs released from transplanted tumor cells were shorter than the mouse counterpart. We next analyzed blood cfDNA samples from hepatocellular carcinoma and prostate cancer patients and found that mt-cfDNA lengths were inversely related to tumor size as well as the concentration of circulating tumor DNA. Our study suggested that monitoring the size of mt-cfDNAs in cancer patients would be a useful way to estimate tumor burden and cancer progression.