VEGF-C/VEGFR-3 signaling plays a central role in lymphatic development, regulating the budding of lymphatic progenitor cells from embryonic veins and maintaining the expression of PROX1 during later developmental stages. However, how VEGFR-3 activation translates into target gene expression is still not completely understood. We used cap analysis of gene expression (CAGE) RNA sequencing to characterize the transcriptional changes invoked by VEGF-C in LECs and to identify the transcription factors (TFs) involved. We found that MAFB, a TF involved in differentiation of various cell types, is rapidly induced and activated by VEGF-C. MAFB induced expression of PROX1 as well as other TFs and markers of differentiated LECs, indicating a role in the maintenance of the mature LEC phenotype. Correspondingly, E14.5 Mafb(-/-) embryos showed impaired lymphatic patterning in the skin. This suggests that MAFB is an important TF involved in lymphangiogenesis.

Big leaps in science happen when scientists from different backgrounds interact. In the past 15 years, the FANTOM Consortium has brought together scientists from different fields to analyze and interpret genomic data produced with novel technologies, including mouse full-length cDNAs and, more recently, expression profiling at single-nucleotide resolution by cap-analysis gene expression. The FANTOM Consortium has provided the most comprehensive mouse cDNA collection for functional studies and extensive maps of the human and mouse transcriptome comprising promoters, enhancers, as well as the network of their regulatory interactions. More importantly, serendipitous observations of the FANTOM dataset led us to realize that the mammalian genome is pervasively transcribed, even from retrotransposon elements, which were previously considered junk DNA. The majority of products from the mammalian genome are long non-coding RNAs (lncRNAs), including sense-antisense, intergenic, and enhancer RNAs. While the biological function has been elucidated for some lncRNAs, more than 98 % of them remain without a known function. We argue that large-scale studies are urgently needed to address the functional role of lncRNAs.

Thrombospondin-2 (TSP2) is an anti-angiogenic matricellular protein that inhibits tumor growth and angiogenesis. Tumor-associated blood vascular endothelial cells (BECs) were isolated from human invasive bladder cancers and from matched normal bladder tissue by immuno-laser capture microdissection. Exon expression profiling analyses revealed a particularly high expression of a short TSP2 transcript containing only the last 9 (3') exons of the full-length TSP2 transcript. Using 5' and 3' RACE (rapid amplification of cDNA ends) and Sanger sequencing, we confirmed the existence of the shorter transcript of TSP2 (sTSP2) and determined its sequence which completely lacked the anti-angiogenic thrombospondin type 1 repeats domain. The largest open reading frame predicted within the transcript comprises 209 amino acids and matches almost completely the C-terminal lectin domain of full-length TSP2. We produced recombinant sTSP2 and found that unlike the full-length TSP2, sTSP2 did not inhibit vascular endothelial growth factor-A-induced proliferation of cultured human BECs, but in contrast when combined with TSP2 blocked the inhibitory effects of TSP2 on BEC proliferation. In vivo studies with stably transfected A431 squamous cell carcinoma cells revealed that full-length TSP2, but not sTSP2, inhibited tumor growth and angiogenesis. This study reveals that the transcriptional program of tumor stromal cells can change to transcribe a new version of an endogenous angiogenesis inhibitor that has lost its anti-angiogenic activity.

The human body consists of 37 trillion single cells represented by over 50 organs that are stitched together to make us who we are, yet we still have very little understanding about the basic units of our body: what cell types and states make up our organs both compositionally and spatially. Previous efforts to profile a wide range of human cell types have been attempted by the FANTOM and GTEx consortia. Now, with the advancement in genomic technologies, profiling the human body at single-cell resolution is possible and will generate an unprecedented wealth of data that will accelerate basic and clinical research with tangible applications to future medicine. To date, several major organs have been profiled, but the challenges lie in ways to integrate single-cell genomics data in a meaningful way. In recent years, several consortia have begun to introduce harmonization and equity in data collection and analysis. Herein, we introduce existing and nascent single-cell genomics consortia, and present benefits to necessitate single-cell genomic consortia in a regional environment to achieve the universal human cell reference dataset.

Recent studies have revealed the importance of long noncoding RNAs (lncRNAs) as tissue-specific regulators of gene expression. There is ample evidence that distinct types of vasculature undergo tight transcriptional control to preserve their structure, identity, and functions. We determine a comprehensive map of lineage-specific lncRNAs in human dermal lymphatic and blood vascular endothelial cells (LECs and BECs), combining RNA-Seq and CAGE-Seq. Subsequent antisense oligonucleotide-knockdown transcriptomic profiling of two LEC- and two BEC-specific lncRNAs identifies LETR1 as a critical gatekeeper of the global LEC transcriptome. Deep RNA-DNA, RNA-protein interaction studies, and phenotype rescue analyses reveal that LETR1 is a nuclear trans-acting lncRNA modulating, via key epigenetic factors, the expression of essential target genes, including KLF4 and SEMA3C, governing the growth and migratory ability of LECs. Together, our study provides several lines of evidence supporting the intriguing concept that every cell type expresses precise lncRNA signatures to control lineage-specific regulatory programs.

The human brain is a complex, three-dimensional structure. To better recapitulate brain complexity, recent efforts have focused on the development of human-specific midbrain organoids. Human iPSC-derived midbrain organoids consist of differentiated and functional neurons, which contain active synapses, as well as astrocytes and oligodendrocytes. However, the absence of microglia, with their ability to remodel neuronal networks and phagocytose apoptotic cells and debris, represents a major disadvantage for the current midbrain organoid systems. Additionally, neuroinflammation-related disease modeling is not possible in the absence of microglia. So far, no studies about the effects of human iPSC-derived microglia on midbrain organoid neural cells have been published. Here we describe an approach to derive microglia from human iPSCs and integrate them into iPSC-derived midbrain organoids. Using single nuclear RNA Sequencing, we provide a detailed characterization of microglia in midbrain organoids as well as the influence of their presence on the other cells of the organoids. Furthermore, we describe the effects that microglia have on cell death and oxidative stress-related gene expression. Finally, we show that microglia in midbrain organoids affect synaptic remodeling and increase neuronal excitability. Altogether, we show a more suitable system to further investigate brain development, as well as neurodegenerative diseases and neuroinflammation.

Efficient induction of the otic placode, the developmental origin of the inner ear from human pluripotent stem cells (hPSCs), provides a robust platform for otic development and sensorineural hearing loss modelling. Nevertheless, there remains a limited capacity of otic lineage specification from hPSCs by stepwise differentiation methods, since the critical factors for successful otic cell differentiation have not been thoroughly investigated. In this study, we developed a novel differentiation system involving the use of a three-dimensional (3D) floating culture with signalling factors for generating otic cell lineages via stepwise differentiation of hPSCs.

Cytotoxic CD4+ T cells (CD4-CTLs) show the presence of cytolytic granules, which include the enzymes granzyme and perforin. The cells have a pathogenic and protective role in various diseases, including cancer, viral infection, and autoimmune disease. In mice, cytotoxic CD4+ T cells express CD8αα+ and reside in the intestine (mouse CD4+CTLs; mCD4-CTLs). The population of cytotoxic CD4+ T cells in the human intestine is currently unknown. Moreover, it is unclear how cytotoxic CD4 T cells change in patients with inflammatory bowel disease (IBD). Here, we aimed to identify cytotoxic CD4+ T cells in the human intestine and analyze the characteristics of the population in patients with IBD using single-cell RNA-seq (scRNA-seq). In CD4+ T cells, granzyme and perforin expression was high in humanMAIT (hMAIT) cells and hCD4+CD8A+ T cell cluster. Both CD4 and CD8A were expressed in hTreg, hMAIT, and hCD4+CD8A+ T cell clusters. Next we performed fast gene set enrichment analysis to identify cell populations that showed homology to mCD4CTLs. The analysis identified the hCD4+CD8A+ T cell cluster (hCTL-like population; hCD4-CTL) similar to mouse CTLs. The percentage of CD4+CD8A+ T cells among the total CD4+ T cells in the inflamed intestine of the patients with Crohn's disease was significantly reduced compared with that in the noninflamed intestine of the patients. In summary, we identified cytotoxic CD4+CD8+ T cells in the small intestine of humans. The integration of the mouse and human sc-RNA-seq data analysis highlight an approach to identify human cell populations related to mouse cell populations, which may help determine the functional properties of several human cell populations in mice.

Tumor-associated macrophages (TAMs) have been implicated in the promotion of breast cancer growth and metastasis, and a strong infiltration by TAMs has been associated with estrogen receptor (ER)-negative tumors and poor prognosis. However, the molecular mechanisms behind these observations are unclear. We investigated macrophage activation in response to co-culture with several breast cancer cell lines (T47D, MCF-7, BT-474, SKBR-3, Cal-51 and MDA-MB-231) and found that high granulocyte colony-stimulating factor (G-CSF) secretion by the triple-negative breast cancer (TNBC) cell line MDA-MB-231 gave rise to immunosuppressive HLA-DRlo macrophages that promoted migration of breast cancer cells via secretion of TGF-α. In human breast cancer samples (n = 548), G-CSF was highly expressed in TNBC (p < 0.001) and associated with CD163+ macrophages (p < 0.0001), poorer overall survival (OS) (p = 0.021) and significantly increased numbers of TGF-α+ cells. While G-CSF blockade in the 4T1 mammary tumor model promoted maturation of MHCIIhi blood monocytes and TAMs and significantly reduced lung metastasis, anti-CSF-1R treatment promoted MHCIIloF4/80hiMRhi anti-inflammatory TAMs and enhanced lung metastasis in the presence of high G-CSF levels. Combined anti-G-CSF and anti-CSF-1R therapy significantly increased lymph node metastases, possibly via depletion of the so-called "gate-keeper" subcapsular sinus macrophages. These results indicate that G-CSF promotes the anti-inflammatory phenotype of tumor-induced macrophages when CSF-1R is inhibited and therefore caution against the use of M-CSF/CSF-1R targeting agents in tumors with high G-CSF expression.

Adipogenesis is the process of cell differentiation by which mesenchymal stem cells (MSC) become adipocytes. Investigating the transcriptional regulatory process during adipogenesis may provide strategies to prevent obesity and other metabolic disorders. In recent years, numerous zinc finger proteins (ZFPs) have been implicated in regulating differentiation and cell fate determination. To investigate the regulatory role of ZFPs involved in adipogenesis, we performed genome-wide microarray expression profiling of an adipogenesis time series. Particularly focusing on the transiently responsive ZFPs, we identified and characterized the functional role of ZNF395 in adipogenesis. A systematic ablation of the ZNF395 transcript during adipogenesis revealed 40% reduction of adipocytes when compared to control. Furthermore, the number of adipocytes as well as the expression of key adipocyte markers were greatly induced when MSC were co-transduced with ZNF395 and PPARG2. To further elucidate the functional role of ZNF395 during adipogenesis, we attempted to trans-differentiate human dermal fibroblasts with PPARG2. The test remarkably revealed that ZNF395 in conjunction with PPARG2 greatly induced adipogenesis from dermal fibroblasts when compared to PPARG2 alone. These loss and gain of function experiments firmly establish that ZNF395 coordinate the transcriptional regulatory pathway with PPARG2, which may be necessary for the genesis of adipocytes.

The lymphatic vascular system maintains tissue fluid homeostasis, helps mediate afferent immune responses, and promotes cancer metastasis. To address the role microRNAs (miRNAs) play in the development and function of the lymphatic vascular system, we defined the in vitro miRNA expression profiles of primary human lymphatic endothelial cells (LECs) and blood vascular endothelial cells (BVECs) and identified four BVEC signature and two LEC signature miRNAs. Their vascular lineage-specific expression patterns were confirmed in vivo by quantitative real-time PCR and in situ hybridization. Functional characterization of the BVEC signature miRNA miR-31 identified a novel BVEC-specific posttranscriptional regulatory mechanism that inhibits the expression of lymphatic lineage-specific transcripts in vitro. We demonstrate that suppression of lymphatic differentiation is partially mediated via direct repression of PROX1, a transcription factor that functions as a master regulator of lymphatic lineage-specific differentiation. Finally, in vivo studies of Xenopus and zebrafish demonstrated that gain of miR-31 function impaired venous sprouting and lymphatic vascular development, thus highlighting the importance of miR-31 as a negative regulator of lymphatic development. Collectively, our findings identify miR-31 is a potent regulator of vascular lineage-specific differentiation and development in vertebrates.

Transcriptional Regulatory Networks (TRNs) coordinate multiple transcription factors (TFs) in concert to maintain tissue homeostasis and cellular function. The re-establishment of target cell TRNs has been previously implicated in direct trans-differentiation studies where the newly introduced TFs switch on a set of key regulatory factors to induce de novo expression and function. However, the extent to which TRNs in starting cell types, such as dermal fibroblasts, protect cells from undergoing cellular reprogramming remains largely unexplored. In order to identify TFs specific to maintaining the fibroblast state, we performed systematic knockdown of 18 fibroblast-enriched TFs and analyzed differential mRNA expression against the same 18 genes, building a Matrix-RNAi. The resulting expression matrix revealed seven highly interconnected TFs. Interestingly, suppressing four out of seven TFs generated lipid droplets and induced PPARG and CEBPA expression in the presence of adipocyte-inducing medium only, while negative control knockdown cells maintained fibroblastic character in the same induction regime. Global gene expression analyses further revealed that the knockdown-induced adipocytes expressed genes associated with lipid metabolism and significantly suppressed fibroblast genes. Overall, this study reveals the critical role of the TRN in protecting cells against aberrant reprogramming, and demonstrates the vulnerability of donor cell's TRNs, offering a novel strategy to induce transgene-free trans-differentiations.

Hydroxymethylcytosines (hmC), one of several reported cytosine modifications, was recently found to be enriched in embryonic stem cells and neuronal cells, and thought to play an important role in regulating gene expression and cell specification. However, unlike methylcytosines (mC), the fate of hmC beyond DNA replication is not well understood. Here, to monitor the status of hmC during DNA replication, we prepared a stable episomal vector-based monitoring system called MoCEV in 293T cells. The MoCEV system containing fully hydroxymethylated-cytosine fragments revealed a significant modification towards mC after several rounds of DNA replication. Strikingly this modification was specifically observed at the CpG sites (71.9% of cytosines), whereas only 1.1% of modified cytosines were detected at the non-CpG sites. Since the unmodified MoCEV did not undergo any DNA methylation during cell division, the results strongly suggest that somatic cells undergo hmC to mC specifically at the CpG sites during cell division.

Single-cell transcriptomic profiling is a powerful tool to explore cellular heterogeneity. However, most of these methods focus on the 3'-end of polyadenylated transcripts and provide only a partial view of the transcriptome. We introduce C1 CAGE, a method for the detection of transcript 5'-ends with an original sample multiplexing strategy in the C1TM microfluidic system. We first quantifiy the performance of C1 CAGE and find it as accurate and sensitive as other methods in the C1 system. We then use it to profile promoter and enhancer activities in the cellular response to TGF-β of lung cancer cells and discover subpopulations of cells differing in their response. We also describe enhancer RNA dynamics revealing transcriptional bursts in subsets of cells with transcripts arising from either strand in a mutually exclusive manner, validated using single molecule fluorescence in situ hybridization.

MAFB is a transcription factor involved in the terminal differentiation of several cell types, including macrophages and keratinocytes. MAFB is also expressed in lymphatic endothelial cells (LECs) and is upregulated by VEGF-C/VEGFR-3 signaling. Recent studies have revealed that MAFB regulates several genes involved in lymphatic differentiation and that global Mafb knockout mice show defects in patterning of lymphatic vessels during embryogenesis. However, it has remained unknown whether this effect is LEC-intrinsic and whether MAFB might also be involved in postnatal lymphangiogenesis. We established conditional, lymphatic-specific Mafb knockout mice and found comparable lymphatic patterning defects during embryogenesis as in the global MAFB knockout. Lymphatic MAFB deficiency resulted in increased lymphatic branching in the diaphragm at P7, but had no major effect on lymphatic patterning or function in healthy adult mice. By contrast, tumor-induced lymphangiogenesis was enhanced in mice lacking lymphatic MAFB. Together, these data reveal that LEC-expressed MAFB is involved in lymphatic vascular morphogenesis during embryonic and postnatal development as well as in pathological conditions. Therefore, MAFB could represent a target for therapeutic modulation of lymphangiogenesis.

Long noncoding RNAs (lncRNAs) constitute the majority of transcripts in the mammalian genomes, and yet, their functions remain largely unknown. As part of the FANTOM6 project, we systematically knocked down the expression of 285 lncRNAs in human dermal fibroblasts and quantified cellular growth, morphological changes, and transcriptomic responses using Capped Analysis of Gene Expression (CAGE). Antisense oligonucleotides targeting the same lncRNAs exhibited global concordance, and the molecular phenotype, measured by CAGE, recapitulated the observed cellular phenotypes while providing additional insights on the affected genes and pathways. Here, we disseminate the largest-to-date lncRNA knockdown data set with molecular phenotyping (over 1000 CAGE deep-sequencing libraries) for further exploration and highlight functional roles for ZNF213-AS1 and lnc-KHDC3L-2.

The brainstem is a posterior region of the brain, composed of three parts, midbrain, pons, and medulla oblongata. It is critical in controlling heartbeat, blood pressure, and respiration, all of which are life-sustaining functions, and therefore, damages to or disorders of the brainstem can be lethal. Brain organoids derived from human pluripotent stem cells (hPSCs) recapitulate the course of human brain development and are expected to be useful for medical research on central nervous system disorders. However, existing organoid models are limited in the extent hPSCs recapitulate human brain development and hence are not able to fully elucidate the diseases affecting various components of the brain such as brainstem. Here, we developed a method to generate human brainstem organoids (hBSOs), containing midbrain/hindbrain progenitors, noradrenergic and cholinergic neurons, dopaminergic neurons, and neural crest lineage cells. Single-cell RNA sequence (scRNA-seq) analysis, together with evidence from proteomics and electrophysiology, revealed that the cellular population in these organoids was similar to that of the human brainstem, which raises the possibility of making use of hBSOs in investigating central nervous system disorders affecting brainstem and in efficient drug screenings.

Classical transection studies suggest that, in addition to the hypothalamus, the brainstem is essential for non-rapid eye movement (NREM) sleep. The circuits underlying this function, however, have remained largely unknown. We identified a circuit distributed in the midbrain, pons, and medulla that promotes NREM sleep in mice. We focused on the sublaterodorsal tegmentum, an area implicated in dual regulation of REM and NREM sleep. Transcriptomic and genetic analyses revealed that neurons positive for the neuropeptide neurotensin promote NREM sleep. Further analyses identified downstream NREM sleep-promoting neurons in the dorsal deep mesencephalic nucleus, the lateral part of the periaqueductal gray, and the medial vestibular nucleus that were also neurotensinergic. Infusion of neurotensin into the fourth ventricle induced NREM sleep-like cortical activity, whereas mice deficient for neurotensin exhibited increased REM sleep, implicating the involvement of the neuropeptide itself. These findings identify a widely distributed NREM sleep-regulating circuit in the brainstem with a common molecular property.

Human multipotent mesenchymal stromal/stem cells (MSCs) have been utilized in cell therapy for various diseases and their clinical applications are expected to increase in the future. However, the variation in MSC-based product quality due to the MSC heterogeneity has resulted in significant constraints in the clinical utility of MSCs. Therefore, we hypothesized that it might be important to identify and ensure/enrich suitable cell subpopulations for therapies using MSC-based products. In this study, we aimed to identify functional cell subpopulations to predict the efficacy of angiogenic therapy using bone marrow-derived MSCs (BM-MSCs). To assess its angiogenic potency, we observed various levels of vascular endothelial growth factor (VEGF) secretion among 11 donor-derived BM-MSC lines under in vitro ischemic culture conditions. Next, by clarifying the heterogeneity of BM-MSCs using single-cell RNA-sequencing analysis, we identified a functional cell subpopulation that contributed to the overall VEGF production in BM-MSC lines under ischemic conditions. We also found that leucine-rich repeat-containing 75A (LRRC75A) was more highly expressed in this cell subpopulation than in the others. Importantly, knockdown of LRRC75A using small interfering RNA resulted in significant inhibition of VEGF secretion in ischemic BM-MSCs, indicating that LRRC75A regulates VEGF secretion under ischemic conditions. Therefore, LRRC75A may be a useful biomarker to identify cell subpopulations that contribute to the angiogenic effects of BM-MSCs. Our work provides evidence that a strategy based on single-cell transcriptome profiles is effective for identifying functional cell subpopulations in heterogeneous MSC-based products.

Cellular reprogramming is driven by a defined set of transcription factors; however, the regulatory logic that underlies cell-type specification and diversification remains elusive. Single-cell RNA-seq provides unprecedented coverage to measure dynamic molecular changes at the single-cell resolution. Here, we multiplex and ectopically express 20 pro-neuronal transcription factors in human dermal fibroblasts and demonstrate a widespread diversification of neurons based on cell morphology and canonical neuronal marker expressions. Single-cell RNA-seq analysis reveals diverse and distinct neuronal subtypes, including reprogramming processes that strongly correlate with the developing brain. Gene mapping of 20 exogenous pro-neuronal transcription factors further unveiled key determinants responsible for neuronal lineage specification and a regulatory logic dictating neuronal diversification, including glutamatergic and cholinergic neurons. The multiplex scRNA-seq approach is a robust and scalable approach to elucidate lineage and cellular specification across various biological systems.