Tumor-initiating cells are a subpopulation in aggressive cancers that exhibit traits shared with stem cells, including the ability to self-renew and differentiate, commonly referred to as stemness. In addition, such cells are resistant to chemo- and radiation therapy posing a therapeutic challenge. To uncover stemness-associated functions in glioma-initiating cells (GICs), transcriptome profiles were compared to neural stem cells (NSCs) and gene ontology analysis identified an enrichment of Ca2+ signaling genes in NSCs and the more stem-like (NSC-proximal) GICs. Functional analysis in a set of different GIC lines regarding sensitivity to disturbed homeostasis using A23187 and Thapsigargin, revealed that NSC-proximal GICs were more sensitive, corroborating the transcriptome data. Furthermore, Ca2+ drug sensitivity was reduced in GICs after differentiation, with most potent effect in the NSC-proximal GIC, supporting a stemness-associated Ca2+ sensitivity. NSCs and the NSC-proximal GIC line expressed a larger number of ion channels permeable to potassium, sodium and Ca2+. Conversely, a higher number of and higher expression levels of Ca2+ binding genes that may buffer Ca2+, were expressed in NSC-distal GICs. In particular, expression of the AMPA glutamate receptor subunit GRIA1, was found to associate with Ca2+ sensitive NSC-proximal GICs, and decreased as GICs differentiated along with reduced Ca2+ drug sensitivity. The correlation between high expression of Ca2+ channels (such as GRIA1) and sensitivity to Ca2+ drugs was confirmed in an additional nine novel GIC lines. Calcium drug sensitivity also correlated with expression of the NSC markers nestin (NES) and FABP7 (BLBP, brain lipid-binding protein) in this extended analysis. In summary, NSC-associated NES+/FABP7+/GRIA1+ GICs were selectively sensitive to disturbances in Ca2+ homeostasis, providing a potential target mechanism for eradication of an immature population of malignant cells.

Pleomorphic adenoma gene 1 (PLAG1) is a transcription factor involved in cancer and growth. We discovered a de novo DNA motif containing a PLAG1 binding site in the promoters of genes activated during zygotic genome activation (ZGA) in human embryos. This motif was located within an Alu element in a region that was conserved in the murine B1 element. We show that maternally provided Plag1 is needed for timely mouse preimplantation embryo development. Heterozygous mouse embryos lacking maternal Plag1 showed disrupted regulation of 1,089 genes, spent significantly longer time in the 2-cell stage, and started expressing Plag1 ectopically from the paternal allele. The de novo PLAG1 motif was enriched in the promoters of the genes whose activation was delayed in the absence of Plag1. Further, these mouse genes showed a significant overlap with genes upregulated during human ZGA that also contain the motif. By gene ontology, the mouse and human ZGA genes with de novo PLAG1 motifs were involved in ribosome biogenesis and protein synthesis. Collectively, our data suggest that PLAG1 affects embryo development in mice and humans through a conserved DNA motif within Alu/B1 elements located in the promoters of a subset of ZGA genes.

Standard RNAseq methods using bulk RNA and recent single-cell RNAseq methods use DNA barcodes to identify samples and cells, and the barcoded cDNAs are pooled into a library pool before high throughput sequencing. In cases of single-cell and low-input RNAseq methods, the library is further amplified by PCR after the pooling. Preparation of hundreds or more samples for a large study often requires multiple library pools. However, sometimes correlation between expression profiles among the libraries is low and batch effect biases make integration of data between library pools difficult.

Caffeine is a widely consumed psychoactive substance, but little is known about the effects of caffeine stimulation on global gene expression changes in neurons. Here, we conducted gene expression profiling of human neuroepithelial stem cell-derived neurons, stimulated with normal consumption levels of caffeine (3 μM and 10 μM), over a period of 9 h. We found dosage-dependent activation of immediate early genes after 1 h. Neuronal projection development processes were up-regulated and negative regulation of axon extension processes were down-regulated at 3 h. In addition, genes involved in extracellular matrix organization, response for wound healing, and regulation of immune system processes were down-regulated by caffeine at 3 h. This study identified novel genes within the neuronal projection guidance pathways that respond to acute caffeine stimulation and suggests potential mechanisms for the effects of caffeine on neuronal cells.

Reverse transcriptases derived from Moloney Murine Leukemia Virus (MMLV) have an intrinsic terminal transferase activity, which causes the addition of a few non-templated nucleotides at the 3' end of cDNA, with a preference for cytosine. This mechanism can be exploited to make the reverse transcriptase switch template from the RNA molecule to a secondary oligonucleotide during first-strand cDNA synthesis, and thereby to introduce arbitrary barcode or adaptor sequences in the cDNA. Because the mechanism is relatively efficient and occurs in a single reaction, it has recently found use in several protocols for single-cell RNA sequencing. However, the base preference of the terminal transferase activity is not known in detail, which may lead to inefficiencies in template switching when starting from tiny amounts of mRNA. Here, we used fully degenerate oligos to determine the exact base preference at the template switching site up to a distance of ten nucleotides. We found a strong preference for guanosine at the first non-templated nucleotide, with a greatly reduced bias at progressively more distant positions. Based on this result, and a number of careful optimizations, we report conditions for efficient template switching for cDNA amplification from single cells.

Recently, human PAIRED-LIKE homeobox transcription factor (TF) genes were discovered whose expression is limited to the period of embryo genome activation up to the 8-cell stage. One of these TFs is LEUTX, but its importance for human embryogenesis is still subject to debate. We confirmed that human LEUTX acts as a TAATCC-targeting transcriptional activator, like other K50-type PAIRED-LIKE TFs. Phylogenetic comparisons revealed that Leutx proteins are conserved across Placentalia and comprise two conserved domains, the homeodomain, and a Leutx-specific domain containing putative transcriptional activation motifs (9aaTAD). Examination of human genotype resources revealed 116 allelic variants in LEUTX. Twenty-four variants potentially affect function, but they occur only heterozygously at low frequency. One variant affects a DNA-specificity determining residue, mutationally reachable by a one-base transition. In vitro and in silico experiments showed that this LEUTX mutation (alanine to valine at position 54 in the homeodomain) results in a transactivational loss-of-function to a minimal TAATCC-containing promoter and a 36 bp motif enriched in genes involved in embryo genome activation. A compensatory change in residue 47 restores function. The results support the notion that human LEUTX functions as a transcriptional activator important for human embryogenesis.

T cells expressing the cutaneous lymphocyte antigen (CLA) mediate pathogenic inflammation in atopic dermatitis (AD). The molecular alterations contributing to their dysregulation remain unclear. With the aim to elucidate putative altered pathways in AD we profiled DNA methylation levels and miRNA expression in sorted T cell populations (CD4+, CD4+CD45RA+ naïve, CD4+CLA+, and CD8+) from adult AD patients and healthy controls (HC). Skin homing CD4+CLA+ T cells from AD patients showed significant differences in DNA methylation in 40 genes compared to HC (p < 0.05). Reduced DNA methylation levels in the upstream region of the interleukin-13 gene (IL13) in CD4+CLA+ T cells from AD patients correlated with increased IL13 mRNA expression in these cells. Sixteen miRNAs showed differential expression in CD4+CLA+ T cells from AD patients targeting genes in 202 biological processes (p < 0.05). An integrated network analysis of miRNAs and CpG sites identified two communities of strongly interconnected regulatory elements with strong antagonistic behaviours that recapitulated the differences between AD patients and HC. Functional analysis of the genes linked to these communities revealed their association with key cytokine signaling pathways, MAP kinase signaling and protein ubiquitination. Our findings support that epigenetic mechanisms play a role in the pathogenesis of AD by affecting inflammatory signaling molecules in skin homing CD4+CLA+ T cells and uncover putative molecules participating in AD pathways.

Reprogramming of cellular identity using exogenous expression of transcription factors (TFs) is a powerful and exciting tool for tissue engineering, disease modeling, and regenerative medicine. However, generation of desired cell types using this approach is often plagued by inefficiency, slow conversion, and an inability to produce mature functional cells. Here, we show that expression of constitutively active SMAD2/3 significantly improves the efficiency of induced pluripotent stem cell (iPSC) generation by the Yamanaka factors. Mechanistically, SMAD3 interacts with reprogramming factors and co-activators and co-occupies OCT4 target loci during reprogramming. Unexpectedly, active SMAD2/3 also markedly enhances three other TF-mediated direct reprogramming conversions, from B cells to macrophages, myoblasts to adipocytes, and human fibroblasts to neurons, highlighting broad and general roles for SMAD2/3 as cell-reprogramming potentiators. Our results suggest that co-expression of active SMAD2/3 could enhance multiple types of TF-based cell identity conversion and therefore be a powerful tool for cellular engineering.

Several studies have demonstrated that human embryonic stem cells (hESC) can be differentiated into trophoblast-like cells if exposed to bone morphogenic protein 4 (BMP4) and/or inhibitors of fibroblast growth factor 2 (FGF2) and the transforming growth factor beta (TGF-β)/activin/nodal signalling pathways. The goal of this study was to investigate how the inhibitors of these pathways improve the efficiency of hESC differentiation when compared with basic BMP4 treatment. RNA sequencing was used to analyse the effects of all possible inhibitor combinations on the differentiation of hESC into trophoblast-like cells over 12 days. Genes differentially expressed compared with untreated cells were identified at seven time points. Additionally, expression of total human chorionic gonadotrophin (HCG) and its hyperglycosylated form (HCG-H) were determined by immunoassay from cell culture media. We showed that FGF2 inhibition with BMP4 activation up-regulates syncytiotrophoblast-specific genes (CGA, CGB and LGALS16), induces several molecular pathways involved in embryo implantation and triggers HCG-H production. In contrast, inhibition of the TGF-β/activin/nodal pathway decreases the ability of hESC to form trophoblast-like cells. Information about the conditions needed for hESC differentiation toward trophoblast-like cells helps us to find an optimal model for studying the early development of human trophoblasts in normal and in complicated pregnancy.

Inhibitory signaling during natural killer (NK) cell education translates into increased responsiveness to activation; however, the intracellular mechanism for functional tuning by inhibitory receptors remains unclear. Secretory lysosomes are part of the acidic lysosomal compartment that mediates intracellular signalling in several cell types. Here we show that educated NK cells expressing self-MHC specific inhibitory killer cell immunoglobulin-like receptors (KIR) accumulate granzyme B in dense-core secretory lysosomes that converge close to the centrosome. This discrete morphological phenotype is independent of transcriptional programs that regulate effector function, metabolism and lysosomal biogenesis. Meanwhile, interference of signaling from acidic Ca2+ stores in primary NK cells reduces target-specific Ca2+-flux, degranulation and cytokine production. Furthermore, inhibition of PI(3,5)P2 synthesis, or genetic silencing of the PI(3,5)P2-regulated lysosomal Ca2+-channel TRPML1, leads to increased granzyme B and enhanced functional potential, thereby mimicking the educated state. These results indicate an intrinsic role for lysosomal remodeling in NK cell education.

RMRP was the first non-coding nuclear RNA gene implicated in a disease. Its mutations cause cartilage-hair hypoplasia (CHH), an autosomal recessive skeletal dysplasia with growth failure, immunodeficiency, and a high risk for malignancies. This study aimed to gain further insight into the role of RNA Component of Mitochondrial RNA Processing Endoribonuclease (RMRP) in cellular physiology and disease pathogenesis. We combined transcriptome analysis with single-cell analysis using fibroblasts from CHH patients and healthy controls. To directly assess cell cycle progression, we followed CHH fibroblasts by pulse-labeling and time-lapse microscopy. Transcriptome analysis identified 35 significantly upregulated and 130 downregulated genes in CHH fibroblasts. The downregulated genes were significantly connected to the cell cycle. Multiple other pathways, involving regulation of apoptosis, bone and cartilage formation, and lymphocyte function, were also affected, as well as PI3K-Akt signaling. Cell-cycle studies indicated that the CHH cells were delayed specifically in the passage from G2 phase to mitosis. Our findings expand the mechanistic understanding of CHH, indicate possible pathways for therapeutic intervention and add to the limited understanding of the functions of RMRP.

The pathogenesis of type 1 diabetes (T1D) is associated with genetic predisposition and immunological changes during presymptomatic disease. Differences in immune cell subset numbers and phenotypes between T1D patients and healthy controls have been described; however, the role and function of these changes in the pathogenesis is still unclear. Here we aimed to analyze the transcriptomic landscapes of peripheral blood mononuclear cells (PBMCs) during presymptomatic disease.

Oligodendrogenesis in the human central nervous system has been observed mainly at the second trimester of gestation, a much later developmental stage compared to oligodendrogenesis in mice. Here, we characterize the transcriptomic neural diversity in the human forebrain at post-conception weeks (PCW) 8-10. Using single-cell RNA sequencing, we find evidence of the emergence of a first wave of oligodendrocyte lineage cells as early as PCW 8, which we also confirm at the epigenomic level through the use of single-cell ATAC-seq. Using regulatory network inference, we predict key transcriptional events leading to the specification of oligodendrocyte precursor cells (OPCs). Moreover, by profiling the spatial expression of 50 key genes through the use of in situ sequencing (ISS), we identify regions in the human ventral fetal forebrain where oligodendrogenesis first occurs. Our results indicate evolutionary conservation of the first wave of oligodendrogenesis between mice and humans and describe regulatory mechanisms involved in human OPC specification.

Methods to spatially profile the transcriptome are dominated by a trade-off between resolution and throughput. Here we develop a method named Enhanced ELectric Fluorescence in situ Hybridization (EEL FISH) that can rapidly process large tissue samples without compromising spatial resolution. By electrophoretically transferring RNA from a tissue section onto a capture surface, EEL speeds up data acquisition by reducing the amount of imaging needed, while ensuring that RNA molecules move straight down toward the surface, preserving single-cell resolution. We apply EEL on eight entire sagittal sections of the mouse brain and measure the expression patterns of up to 440 genes to reveal complex tissue organization. Moreover, EEL can be used to study challenging human samples by removing autofluorescent lipofuscin, enabling the spatial transcriptome of the human visual cortex to be visualized. We provide full hardware specifications, all protocols and complete software for instrument control, image processing, data analysis and visualization.

Conventional reprogramming methods rely on the ectopic expression of transcription factors to reprogram somatic cells into induced pluripotent stem cells (iPSCs). The forced expression of transcription factors may lead to off-target gene activation and heterogeneous reprogramming, resulting in the emergence of alternative cell types and aberrant iPSCs. Activation of endogenous pluripotency factors by CRISPR activation (CRISPRa) can reduce this heterogeneity. Here, we describe a high-efficiency reprogramming of human somatic cells into iPSCs using optimized CRISPRa. Efficient reprogramming was dependent on the additional targeting of the embryo genome activation-enriched Alu-motif and the miR-302/367 locus. Single-cell transcriptome analysis revealed that the optimized CRISPRa reprogrammed cells more directly and specifically into the pluripotent state when compared to the conventional reprogramming method. These findings support the use of CRISPRa for high-quality pluripotent reprogramming of human cells.

Recent development in inferring RNA velocity from single-cell RNA-seq opens up exciting new vista into developmental lineage and cellular dynamics. However, the estimated velocity only gives a snapshot of how the transcriptome instantaneously changes in individual cells, and it does not provide quantitative predictions and insights about the whole system. In this work, we develop RNA-ODE, a principled computational framework that extends RNA velocity to quantify systems level dynamics and improve single-cell data analysis. We model the gene expression dynamics by an ordinary differential equation (ODE) based formalism. Given a snapshot of gene expression at one time, RNA-ODE is able to predict and extrapolate the expression trajectory of each cell by solving the dynamic equations. Systematic experiments on simulations and on new data from developing brain demonstrate that RNA-ODE substantially improves many aspects of standard single-cell analysis. By leveraging temporal dynamics, RNA-ODE more accurately estimates cell state lineage and pseudo-time compared to previous state-of-the-art methods. It also infers gene regulatory networks and identifies influential genes whose expression changes can decide cell fate. We expect RNA-ODE to be a Swiss army knife that aids many facets of single-cell RNA-seq analysis.

Recent advances in single-cell transcriptomics have illuminated the diverse neuronal and glial cell types within the human brain. However, the regulatory programs governing cell identity and function remain unclear. Using a single-nucleus assay for transposase-accessible chromatin using sequencing (snATAC-seq), we explored open chromatin landscapes across 1.1 million cells in 42 brain regions from three adults. Integrating this data unveiled 107 distinct cell types and their specific utilization of 544,735 candidate cis-regulatory DNA elements (cCREs) in the human genome. Nearly a third of the cCREs demonstrated conservation and chromatin accessibility in the mouse brain cells. We reveal strong links between specific brain cell types and neuropsychiatric disorders including schizophrenia, bipolar disorder, Alzheimer's disease (AD), and major depression, and have developed deep learning models to predict the regulatory roles of noncoding risk variants in these disorders.

The human brain directs complex behaviors, ranging from fine motor skills to abstract intelligence, but the diversity of cell types that support these skills has not been fully described. In this work, we used single-nucleus RNA sequencing to systematically survey cells across the entire adult human brain. We sampled more than three million nuclei from approximately 100 dissections across the forebrain, midbrain, and hindbrain in three postmortem donors. Our analysis identified 461 clusters and 3313 subclusters organized largely according to developmental origins and revealing high diversity in midbrain and hindbrain neurons. Astrocytes and oligodendrocyte-lineage cells also exhibited regional diversity at multiple scales. The transcriptomic census of the entire human brain presented in this work provides a resource for understanding the molecular diversity of the human brain in health and disease.

The effect of vitamin D at the transcriptome level is poorly understood, and furthermore, it is unclear if it differs between obese and normal-weight subjects. The objective of the study was to explore the transcriptome effects of vitamin D supplementation.

No abstract available