Keratinocytes (KCs) are the most frequent cells in the epidermis, and they are often isolated and cultured in vitro to study the molecular biology of the skin. Cultured primary cells and various immortalized cells have been frequently used as skin models but their comparability to intact skin has been questioned. Moreover, when analyzing KC transcriptomes, fluctuation of polyA+ RNA content during the KCs' lifecycle has been omitted.

The mammalian cerebral cortex supports cognitive functions such as sensorimotor integration, memory, and social behaviors. Normal brain function relies on a diverse set of differentiated cell types, including neurons, glia, and vasculature. Here, we have used large-scale single-cell RNA sequencing (RNA-seq) to classify cells in the mouse somatosensory cortex and hippocampal CA1 region. We found 47 molecularly distinct subclasses, comprising all known major cell types in the cortex. We identified numerous marker genes, which allowed alignment with known cell types, morphology, and location. We found a layer I interneuron expressing Pax6 and a distinct postmitotic oligodendrocyte subclass marked by Itpr2. Across the diversity of cortical cell types, transcription factors formed a complex, layered regulatory code, suggesting a mechanism for the maintenance of adult cell type identity.

Cell-free DNA has been used for fetal rhesus factor and sex determination, fetal aneuploidy screening, cancer diagnostics and monitoring, and other applications. However current methods of using cell free DNA require amplification, which leads to allelic dropout and bias especially when starting with small amounts of DNA. Here we describe an amplification-free method for sequencing of cell-free DNA, even from low levels of starting material. We evaluated this method in the context of prenatal diagnosis of fetal aneuploidy and compared it with a PCR-based library preparation method as well as a recently described method using unique molecular identifiers (UMI). All methods performed well, however coverage was increased by the amplification-free method and GC-induced bias was reduced by both the amplification-free method and the UMI method. Future diagnostic applications including whole genome sequencing of cell-free DNA will benefit from amplification-free sequencing.

Evidence for distinct human cancer stem cells (CSCs) remains contentious and the degree to which different cancer cells contribute to propagating malignancies in patients remains unexplored. In low- to intermediate-risk myelodysplastic syndromes (MDS), we establish the existence of rare multipotent MDS stem cells (MDS-SCs), and their hierarchical relationship to lineage-restricted MDS progenitors. All identified somatically acquired genetic lesions were backtracked to distinct MDS-SCs, establishing their distinct MDS-propagating function in vivo. In isolated del(5q)-MDS, acquisition of del(5q) preceded diverse recurrent driver mutations. Sequential analysis in del(5q)-MDS revealed genetic evolution in MDS-SCs and MDS-progenitors prior to leukemic transformation. These findings provide definitive evidence for rare human MDS-SCs in vivo, with extensive implications for the targeting of the cells required and sufficient for MDS-propagation.

PAIRED (PRD)-like homeobox genes belong to a class of predicted transcription factor genes. Several of these PRD-like homeobox genes have been predicted in silico from genomic sequence but until recently had no evidence of transcript expression. We found recently that nine PRD-like homeobox genes, ARGFX, CPHX1, CPHX2, DPRX, DUXA, DUXB, NOBOX, TPRX1 and TPRX2, were expressed in human preimplantation embryos. In the current study we characterized these PRD-like homeobox genes in depth and studied their functions as transcription factors. We cloned multiple transcript variants from human embryos and showed that the expression of these genes is specific to embryos and pluripotent stem cells. Overexpression of the genes in human embryonic stem cells confirmed their roles as transcription factors as either activators (CPHX1, CPHX2, ARGFX) or repressors (DPRX, DUXA, TPRX2) with distinct targets that could be explained by the amino acid sequence in homeodomain. Some PRD-like homeodomain transcription factors had high concordance of target genes and showed enrichment for both developmentally important gene sets and a 36 bp DNA recognition motif implicated in Embryo Genome Activation (EGA). Our data implicate a role for these previously uncharacterized PRD-like homeodomain proteins in the regulation of human embryo genome activation and preimplantation embryo development.

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

The number of mammalian transcripts identified by full-length cDNA projects and genome sequencing projects is increasing remarkably. Clustering them into a strictly nonredundant and comprehensive set provides a platform for functional analysis of the transcriptome and proteome, but the quality of the clustering and predictive usefulness have previously required manual curation to identify truncated transcripts and inappropriate clustering of closely related sequences. A Representative Transcript and Protein Sets (RTPS) pipeline was previously designed to identify the nonredundant and comprehensive set of mouse transcripts based on clustering of a large mouse full-length cDNA set (FANTOM2). Here we propose an alternative method that is more robust, requires less manual curation, and is applicable to other organisms in addition to mouse. RTPSs of human, mouse, and rat have been produced by this method and used for validation. Their comprehensiveness and quality are discussed by comparison with other clustering approaches. The RTPSs are available at .

There is a great need to identify factors that could protect pancreatic beta cells against apoptosis or stimulate their replication and thus prevent or reverse the development of diabetes. One potential candidate is mesencephalic astrocyte-derived neurotrophic factor (MANF), an endoplasmic reticulum (ER) stress inducible protein. Manf knockout mice used as a model of diabetes develop the condition because of increased apoptosis and reduced proliferation of beta cells, apparently related to ER stress. Given this novel association between MANF and beta cell death, we studied the potential of MANF to protect human beta cells against experimentally induced ER stress.

Tumor heterogeneity in breast cancer tumors is today widely recognized. Most of the available knowledge in genetic variation however, relates to the primary tumor while metastatic lesions are much less studied. Many studies have revealed marked alterations of standard prognostic and predictive factors during tumor progression. Characterization of paired primary- and metastatic tissues should therefore be fundamental in order to understand mechanisms of tumor progression, clonal relationship to tumor evolution as well as the therapeutic aspects of systemic disease.

Systems biology is increasingly being applied in nanosafety research for observing and predicting the biological perturbations inflicted by exposure to nanoparticles (NPs). In the present study, we used a combined transcriptomics and proteomics approach to assess the responses of human monocytic cells to Au-NPs of two different sizes with three different surface functional groups, i.e., alkyl ammonium bromide, alkyl sodium carboxylate, or poly(ethylene glycol) (PEG)-terminated Au-NPs. Cytotoxicity screening using THP-1 cells revealed a pronounced cytotoxicity for the ammonium-terminated Au-NPs, while no cell death was seen after exposure to the carboxylated or PEG-modified Au-NPs. Moreover, Au-NR3+ NPs, but not the Au-COOH NPs, were found to trigger dose-dependent lethality in vivo in the model organism, Caenorhabditis elegans. RNA sequencing combined with mass spectrometry-based proteomics predicted that the ammonium-modified Au-NPs elicited mitochondrial dysfunction. The latter results were validated by using an array of assays to monitor mitochondrial function. Au-NR3+ NPs were localized in mitochondria of THP-1 cells. Moreover, the cationic Au-NPs triggered autophagy in macrophage-like RFP-GFP-LC3 reporter cells, and cell death was aggravated upon inhibition of autophagy. Taken together, these studies have disclosed mitochondria-dependent effects of cationic Au-NPs resulting in the rapid demise of the cells.

Developmental cell death plays an important role in the construction of functional neural circuits. In vertebrates, the canonical view proposes a selection of the surviving neurons through stochastic competition for target-derived neurotrophic signals, implying an equal potential for neurons to compete. Here we show an alternative cell fitness selection of neurons that is defined by a specific neuronal heterogeneity code. Proprioceptive sensory neurons that will undergo cell death and those that will survive exhibit different molecular signatures that are regulated by retinoic acid and transcription factors, and are independent of the target and neurotrophins. These molecular features are genetically encoded, representing two distinct subgroups of neurons with contrasted functional maturation states and survival outcome. Thus, in this model, a heterogeneous code of intrinsic cell fitness in neighboring neurons provides differential competitive advantage resulting in the selection of cells with higher capacity to survive and functionally integrate into neural networks.

Single-cell RNA-seq has become routine for discovering cell types and revealing cellular diversity, but archived human brain samples still pose a challenge to current high-throughput platforms. We present STRT-seq-2i, an addressable 9600-microwell array platform, combining sampling by limiting dilution or FACS, with imaging and high throughput at competitive cost. We applied the platform to fresh single mouse cortical cells and to frozen post-mortem human cortical nuclei, matching the performance of a previous lower-throughput platform while retaining a high degree of flexibility, potentially also for other high-throughput applications.

Despite recent advances in understanding microbial diversity in skin homeostasis, the relevance of microbial dysbiosis in inflammatory disease is poorly understood. Here we perform a comparative analysis of skin microbial communities coupled to global patterns of cutaneous gene expression in patients with atopic dermatitis or psoriasis. The skin microbiota is analysed by 16S amplicon or whole genome sequencing and the skin transcriptome by microarrays, followed by integration of the data layers. We find that atopic dermatitis and psoriasis can be classified by distinct microbes, which differ from healthy volunteers microbiome composition. Atopic dermatitis is dominated by a single microbe (Staphylococcus aureus), and associated with a disease relevant host transcriptomic signature enriched for skin barrier function, tryptophan metabolism and immune activation. In contrast, psoriasis is characterized by co-occurring communities of microbes with weak associations with disease related gene expression. Our work provides a basis for biomarker discovery and targeted therapies in skin dysbiosis.

CCHCR1 (Coiled-Coil α-Helical Rod protein 1) is a putative psoriasis candidate gene with the risk alleles CCHCR1*WWCC and *Iso3, the latter inhibiting the translation of isoform 1. CCHCR1 was recently shown to be a centrosomal protein, as well as a component of cytoplasmic processing bodies (P-bodies) that regulate mRNA turnover. The function of CCHCR1 has remained unsettled, partly because of the inconsistent findings; it has been shown to play a wide variety of roles in divergent processes, e.g., cell proliferation and steroidogenesis. Here we utilized RNA sequencing (RNAseq) using HEK293 cells overexpressing isoforms 1 or 3 (Iso1, Iso3 cells), in combination with the coding non-risk or risk (*WWCC) haplotype of CCHCR1. Our aim was to study the overall role of CCHCR1 and the effects of its variants.

Long terminal repeat (LTR) retrotransposons are widely distributed across the human genome. They have accumulated through retroviral integration into germline DNA and are latent genetic modules. Active LTR promoters are observed in germline cells; however, little is known about the mechanisms underlying their active transcription in somatic tissues. Here, by integrating our previous transcriptome data set with publicly available data sets, we show that the LTR families MLT2A1 and MLT2A2 are primarily expressed in human four-cell and eight-cell embryos and are also activated in some adult somatic tissues, particularly pineal gland. Three MLT2A elements function as the promoters and first exons of the protein-coding genes ABCE1, COL5A1, and GALNT13 specifically in the pineal gland of humans but not in that of macaques, suggesting that the exaptation of these LTRs as promoters occurred during recent primate evolution. This analysis provides insight into the possible transition from germline insertion to somatic expression of LTR retrotransposons.

Double homeobox 4 (DUX4) is expressed at the early pre-implantation stage in human embryos. Here we show that induced human DUX4 expression substantially alters the chromatin accessibility of non-coding DNA and activates thousands of newly identified transcribed enhancer-like regions, preferentially located within ERVL-MaLR repeat elements. CRISPR activation of transcribed enhancers by C-terminal DUX4 motifs results in the increased expression of target embryonic genome activation (EGA) genes ZSCAN4 and KHDC1P1. We show that DUX4 is markedly enriched in human zygotes, followed by intense nuclear DUX4 localization preceding and coinciding with minor EGA. DUX4 knockdown in human zygotes led to changes in the EGA transcriptome but did not terminate the embryos. We also show that the DUX4 protein interacts with the Mediator complex via the C-terminal KIX binding motif. Our findings contribute to the understanding of DUX4 as a regulator of the non-coding genome.

The adult human brain comprises more than a thousand distinct neuronal and glial cell types, a diversity that emerges during early brain development. To reveal the precise sequence of events during early brain development, we used single-cell RNA sequencing and spatial transcriptomics and uncovered cell states and trajectories in human brains at 5 to 14 postconceptional weeks (pcw). We identified 12 major classes that are organized as ~600 distinct cell states, which map to precise spatial anatomical domains at 5 pcw. We described detailed differentiation trajectories of the human forebrain and midbrain and found a large number of region-specific glioblasts that mature into distinct pre-astrocytes and pre-oligodendrocyte precursor cells. Our findings reveal the establishment of cell types during the first trimester of human brain development.

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.

Traditionally, neuroscientists have defined the identity of neurons by the cells' location, morphology, connectivity and excitability. However, the direct relationship between these parameters and the molecular phenotypes has remained largely unexplored. Here, we present a method for obtaining full transcriptome data from single neocortical pyramidal cells and interneurons after whole-cell patch-clamp recordings in mouse brain slices. In our approach, termed Patch-seq, a patch-clamp stimulus protocol is followed by the aspiration of the entire somatic compartment into the recording pipette, reverse transcription of RNA including addition of unique molecular identifiers, cDNA amplification, Illumina library preparation and sequencing. We show that Patch-seq reveals a close link between electrophysiological characteristics, responses to acute chemical challenges and RNA expression of neurotransmitter receptors and channels. Moreover, it distinguishes neuronal subpopulations that correspond to both well-established and, to our knowledge, hitherto undescribed neuronal subtypes. Our findings demonstrate the ability of Patch-seq to precisely map neuronal subtypes and predict their network contributions in the brain.

The formation of functional connectivity in the nervous system is governed by axon guidance that instructs nerve growth and branching during development, implying a similarity between neuronal subtypes in terms of nerve extension. We demonstrate the molecular mechanism of another layer of complexity in vertebrates by defining a transcriptional program underlying growth differences between positionally different neurons. The rate of axon extension of the early subset of embryonic dorsal root ganglion sensory neurons is encoded in neurons at different axial levels. This code is determined by a segmental pattern of axial levels of Runx family transcription factor Runx3. Runx3 in turn determines transcription levels of genes encoding cytoskeletal proteins involved in axon extension, including Rock1 and Rock2 which have ongoing activities determining axon growth in early sensory neurons and blocking Rock activity reverses axon extension deficits of Runx3(-/-) neurons. Thus, Runx3 acts to regulate positional differences in axon extension properties apparently without affecting nerve guidance and branching, a principle that could be relevant to other parts of the nervous system.