We have developed a high-throughput amplification method for generating robust gene expression profiles using single cell or low RNA inputs.

Progressive myoclonus epilepsy of Unverricht-Lundborg type (EPM1, OMIM254800) is an autosomal recessive neurodegenerative disorder characterized by stimulus-sensitive and action-activated myoclonus, tonic-clonic epileptic seizures, and ataxia. Loss-of-function mutations in the gene encoding the cysteine protease inhibitor cystatin B (CSTB) underlie EPM1. The deficiency of CSTB in mice (Cstb-/- mice) generates a phenotype resembling the symptoms of EPM1 patients and is accompanied by microglial activation at two weeks of age and an upregulation of immune system-associated genes in the cerebellum at one month of age. To shed light on molecular pathways and processes linked to CSTB deficiency in microglia we characterized the transcriptome of cultured Cstb-/- mouse microglia using microarray hybridization and RNA sequencing (RNA-seq). The gene expression profiles obtained with these two techniques were in good accordance and not polarized to either pro- or anti-inflammatory status. In Cstb-/- microglia, altogether 184 genes were differentially expressed. Of these, 33 genes were identified by both methods. Several interferon-regulated genes were weaker expressed in Cstb-/- microglia compared to control. This was confirmed by quantitative real-time PCR of the transcripts Irf7 and Stat1. Subsequently, we explored the biological context of CSTB deficiency in microglia more deeply by functional enrichment and canonical pathway analysis. This uncovered a potential role for CSTB in chemotaxis, antigen-presentation, and in immune- and defense response-associated processes by altering JAK-STAT pathway signaling. These data support and expand the previously suggested involvement of inflammatory processes to the disease pathogenesis of EPM1 and connect CSTB deficiency in microglia to altered expression of interferon-regulated genes.

The cell cycle progression in mouse embryonic stem cells (mESCs) is controlled by ion fluxes that alter cell volume [1]. This suggests that ion fluxes might control dynamic changes in morphology over the cell cycle, such as rounding up of the cell at mitosis. However, specific channels regulating such dynamic changes and the possible interactions with actomyosin complex have not been clearly identified. Following RNAseq transcriptome analysis of cell cycle sorted mESCs, we found that expression of the K(+) ion channel Erg1 peaked in G1 cell cycle phase, which was confirmed by immunostaining. Inhibition of Erg channel activity caused loss of G1 phase cells via non-apoptotic cell death. Cells first lost the ability of membrane blebbing, a typical feature of cultured embryonic stem cells. Continued Erg inhibition further increased cell volume and the cell eventually ruptured. In addition, atomic force measurements on live cells revealed a decreased cortical stiffness after treatment, suggesting alterations in actomyosin organization. When the intracellular osmotic pressure was experimentally decreased by hypertonic solution or block of K(+) ion import via the Na, K-ATPase, cell viability was restored and cells acquired normal volume and blebbing activity. Our results suggest that Erg channels have a critical function in K(+) ion homeostasis of mESCs over the cell cycle, and that cell death following Erg inhibition is a consequence of the inability to regulate cell volume.

The cell cycle coordinates core functions such as replication and cell division. However, cell-cycle-regulated transcription in the control of non-core functions, such as cell identity maintenance through specific transcription factors (TFs) and signalling pathways remains unclear. Here, we provide a resource consisting of mapped transcriptomes in unsynchronized HeLa and U2OS cancer cells sorted for cell cycle phase by Fucci reporter expression. We developed a novel algorithm for data analysis that enables efficient visualization and data comparisons and identified cell cycle synchronization of Notch signalling and TFs associated with development. Furthermore, the cell cycle synchronizes with the circadian clock, providing a possible link between developmental transcriptional networks and the cell cycle. In conclusion we find that cell cycle synchronized transcriptional patterns are temporally compartmentalized and more complex than previously anticipated, involving genes, which control cell identity and development.

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.

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.

Infertility is a worldwide concern that can be treated with in vitro fertilization (IVF). Improvements in IVF and infertility treatment depend largely on better understanding of the molecular mechanisms for human preimplantation development. Several large-scale studies have been conducted to identify gene expression patterns for the first five days of human development, and many functional studies utilize mouse as a model system. We have identified genes of possible importance for this time period by analyzing human microarray data and available data from online databases. We selected 70 candidate genes for human preimplantation development and investigated their expression in the early mouse development from oocyte to the 8-cell stage. Maternally loaded genes expectedly decreased in expression during development both in human and mouse. We discovered that 25 significantly upregulated genes after fertilization in human included 13 genes whose orthologs in mouse behaved differently and mimicked the expression profile of maternally expressed genes. Our findings highlight many significant differences in gene expression patterns during mouse and human preimplantation development. We also describe four cancer-testis antigen families that are also highly expressed in human embryos: PRAME, SSX, GAGE and MAGEA.