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

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.

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 SCN1A gene, coding for the voltage-gated Na+ channel alpha subunit NaV1.1, is the clinically most relevant epilepsy gene. With the advent of high-throughput next-generation sequencing, clinical laboratories are generating an ever-increasing catalogue of SCN1A variants. Variants are more likely to be classified as pathogenic if they have already been identified previously in a patient with epilepsy. Here, we critically re-evaluate the pathogenicity of this class of variants in a cohort of patients with common epilepsy syndromes and subsequently ask whether a significant fraction of benign variants have been misclassified as pathogenic.

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.

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.

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.

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.

Progressive myoclonus epilepsy of Unverricht-Lundborg type (EPM1) is an autosomal recessively inherited neurodegenerative disease, manifesting with myoclonus, seizures and ataxia, caused by mutations in the cystatin B (CSTB) gene. With the aim of understanding the molecular basis of pathogenetic events in EPM1 we characterized gene expression changes in the cerebella of pre-symptomatic postnatal day 7 (P7) and symptomatic P30 cystatin B -deficient (Cstb(-/-) ) mice, a model for the disease, and in cultured Cstb(-/-) cerebellar granule cells using a pathway-based approach. Differentially expressed genes in P7 cerebella were connected to synaptic function and plasticity, and in cultured cerebellar granule cells, to cell cycle, cytoskeleton, and intracellular transport. In particular, the gene expression data pinpointed alterations in GABAergic pathway. Electrophysiological recordings from Cstb(-/-) cerebellar Purkinje cells revealed a shift of the balance towards decreased inhibition, yet the amount of inhibitory interneurons was not declined in young animals. Instead, we found diminished number of GABAergic terminals and reduced ligand binding to GABAA receptors in Cstb(-/-) cerebellum. These results suggest that alterations in GABAergic signaling could result in reduced inhibition in Cstb(-/-) cerebellum leading to the hyperexcitable phenotype of Cstb(-/-) mice. At P30, the microarray data revealed a marked upregulation of immune and defense response genes, compatible with the previously reported early glial activation that precedes neuronal degeneration. This further implies the role of early-onset neuroinflammation in the pathogenesis of EPM1.

Unverricht-Lundborg type progressive myoclonus epilepsy (EPM1, OMIM 254800) is an autosomal recessive disorder characterized by onset at the age of 6 to 16 years, incapacitating stimulus-sensitive myoclonus and tonic-clonic epileptic seizures. It is caused by mutations in the gene encoding cystatin B. Previously, widespread white matter changes and atrophy has been detected both in adult EPM1 patients and in 6-month-old cystatin B-deficient mice, a mouse model for the EPM1 disease. In order to elucidate the spatiotemporal dynamics of the brain atrophy and white matter changes in EPM1, we conducted longitudinal in vivo magnetic resonance imaging and ex vivo diffusion tensor imaging accompanied with tract-based spatial statistics analysis to compare volumetric changes and fractional anisotropy in the brains of 1 to 6 months of age cystatin B-deficient and control mice. The results reveal progressive but non-uniform volume loss of the cystatin B-deficient mouse brains, indicating that different neuronal populations possess distinct sensitivity to the damage caused by cystatin B deficiency. The diffusion tensor imaging data reveal early and progressive white matter alterations in cystatin B-deficient mice affecting all major tracts. The results also indicate that the white matter damage in the cystatin B-deficient brain is most likely secondary to glial activation and neurodegenerative events rather than a primary result of CSTB deficiency. The data also show that diffusion tensor imaging combined with TBSS analysis provides a feasible approach not only to follow white matter damage in neurodegenerative mouse models but also to detect fractional anisotropy changes related to normal white matter maturation and reorganisation.

Genetic Generalized Epilepsy (GGE) and benign epilepsy with centro-temporal spikes or Rolandic Epilepsy (RE) are common forms of genetic epilepsies. Rare copy number variants have been recognized as important risk factors in brain disorders. We performed a systematic survey of rare deletions affecting protein-coding genes derived from exome data of patients with common forms of genetic epilepsies. We analysed exomes from 390 European patients (196 GGE and 194 RE) and 572 population controls to identify low-frequency genic deletions. We found that 75 (32 GGE and 43 RE) patients out of 390, i.e. ~19%, carried rare genic deletions. In particular, large deletions (>400 kb) represent a higher burden in both GGE and RE syndromes as compared to controls. The detected low-frequency deletions (1) share genes with brain-expressed exons that are under negative selection, (2) overlap with known autism and epilepsy-associated candidate genes, (3) are enriched for CNV intolerant genes recorded by the Exome Aggregation Consortium (ExAC) and (4) coincide with likely disruptive de novo mutations from the NPdenovo database. Employing several knowledge databases, we discuss the most prominent epilepsy candidate genes and their protein-protein networks for GGE and RE.