Generation of a Broadly Useful Model for COVID-19 Pathogenesis, Vaccination, and Treatment.

COVID-19, caused by SARS-CoV-2, is a virulent pneumonia, with >4,000,000 confirmed cases worldwide and >290,000 deaths as of May 15, 2020. It is critical that vaccines and therapeutics be developed very rapidly. Mice, the ideal animal for assessing such interventions, are resistant to SARS-CoV-2. Here, we overcome this difficulty by exogenous delivery of human ACE2 with a replication-deficient adenovirus (Ad5-hACE2). Ad5-hACE2-sensitized mice developed pneumonia characterized by weight loss, severe pulmonary pathology, and high-titer virus replication in lungs. Type I interferon, T cells, and, most importantly, signal transducer and activator of transcription 1 (STAT1) are critical for virus clearance and disease resolution in these mice. Ad5-hACE2-transduced mice enabled rapid assessments of a vaccine candidate, of human convalescent plasma, and of two antiviral therapies (poly I:C and remdesivir). In summary, we describe a murine model of broad and immediate utility to investigate COVID-19 pathogenesis and to evaluate new therapies and vaccines.

Transcription factor Hoxb5 reprograms B cells into functional T lymphocytes.

Deletion of master regulators of the B cell lineage reprograms B cells into T cells. Here we found that the transcription factor Hoxb5, which is expressed in uncommitted hematopoietic progenitor cells but is not present in cells committed to the B cell or T cell lineage, was able to reprogram pro-pre-B cells into functional early T cell lineage progenitors. This reprogramming started in the bone marrow and was completed in the thymus and gave rise to T lymphocytes with transcriptomes, hierarchical differentiation, tissue distribution and immunological functions that closely resembled those of their natural counterparts. Hoxb5 repressed B cell 'master genes', activated regulators of T cells and regulated crucial chromatin modifiers in pro-pre-B cells and ultimately drove the B cell fate-to-T cell fate conversion. Our results provide a de novo paradigm for the generation of functional T cells through reprogramming in vivo.

Publisher Correction: Transcription factor Hoxb5 reprograms B cells into functional T lymphocytes.

In the version of this article initially published, some identification of the supplementary information was incorrect. The items originally called Supplementary Tables 1, 2, 3, 4 and 5 should be Source Data Figures 1, 2, 4, 5 and 7, respectively; those originally called Supplementary Tables 6, 7 and 8 should be Supplementary Tables 1, 2 and 3, respectively; and those originally called Source Data Figures 1, 2, 4, 5 and 7 should be Supplementary Tables 4, 5, 6, 7 and 8, respectively. The errors have been corrected in the HTML version of the article.

WNT signalling control by KDM5C during development affects cognition.

Although KDM5C is one of the most frequently mutated genes in X-linked intellectual disability1, the exact mechanisms that lead to cognitive impairment remain unknown. Here we use human patient-derived induced pluripotent stem cells and Kdm5c knockout mice to conduct cellular, transcriptomic, chromatin and behavioural studies. KDM5C is identified as a safeguard to ensure that neurodevelopment occurs at an appropriate timescale, the disruption of which leads to intellectual disability. Specifically, there is a developmental window during which KDM5C directly controls WNT output to regulate the timely transition of primary to intermediate progenitor cells and consequently neurogenesis. Treatment with WNT signalling modulators at specific times reveal that only a transient alteration of the canonical WNT signalling pathway is sufficient to rescue the transcriptomic and chromatin landscapes in patient-derived cells and to induce these changes in wild-type cells. Notably, WNT inhibition during this developmental period also rescues behavioural changes of Kdm5c knockout mice. Conversely, a single injection of WNT3A into the brains of wild-type embryonic mice cause anxiety and memory alterations. Our work identifies KDM5C as a crucial sentinel for neurodevelopment and sheds new light on KDM5C mutation-associated intellectual disability. The results also increase our general understanding of memory and anxiety formation, with the identification of WNT functioning in a transient nature to affect long-lasting cognitive function.

The RNA m6A reader YTHDC1 silences retrotransposons and guards ES cell identity.

The RNA modification N6-methyladenosine (m6A) has critical roles in many biological processes1,2. However, the function of m6A in the early phase of mammalian development remains poorly understood. Here we show that the m6A reader YT521-B homology-domain-containing protein 1 (YTHDC1) is required for the maintenance of mouse embryonic stem (ES) cells in an m6A-dependent manner, and that its deletion initiates cellular reprogramming to a 2C-like state. Mechanistically, YTHDC1 binds to the transcripts of retrotransposons (such as intracisternal A particles, ERVK and LINE1) in mouse ES cells and its depletion results in the reactivation of these silenced retrotransposons, accompanied by a global decrease in SETDB1-mediated trimethylation at lysine 9 of histone H3 (H3K9me3). We further demonstrate that YTHDC1 and its target m6A RNAs act upstream of SETDB1 to repress retrotransposons and Dux, the master inducer of the two-cell stage (2C)-like program. This study reveals an essential role for m6A RNA and YTHDC1 in chromatin modification and retrotransposon repression.

Resolving Cell Fate Decisions during Somatic Cell Reprogramming by Single-Cell RNA-Seq.

Somatic cells can be reprogrammed into induced pluripotent stem cells (iPSCs), which is a highly heterogeneous process. Here we report the cell fate continuum during somatic cell reprogramming at single-cell resolution. We first develop SOT to analyze cell fate continuum from Oct4/Sox2/Klf4- or OSK-mediated reprogramming and show that cells bifurcate into two categories, reprogramming potential (RP) or non-reprogramming (NR). We further show that Klf4 contributes to Cd34+/Fxyd5+/Psca+ keratinocyte-like NR fate and that IFN-γ impedes the final transition to chimera-competent pluripotency along the RP cells. We analyze more than 150,000 single cells from both OSK and chemical reprograming and identify additional NR/RP bifurcation points. Our work reveals a generic bifurcation model for cell fate decisions during somatic cell reprogramming that may be applicable to other systems and inspire further improvements for reprogramming.

Autism patient-derived SHANK2BY29X mutation affects the development of ALDH1A1 negative dopamine neuron.

Autism spectrum disorder (ASD) encompasses a range of neurodevelopmental conditions. Different mutations on a single ASD gene contribute to heterogeneity of disease phenotypes, possibly due to functional diversity of generated isoforms. SHANK2, a causative gene in ASD, demonstrates this phenomenon, but there is a scarcity of tools for studying endogenous SHANK2 proteins in an isoform-specific manner. Here, we report a point mutation on SHANK2, which is found in a patient with autism, located on exon of the SHANK2B transcript variant (NM_133266.5), hereby SHANK2BY29X. This mutation results in an early stop codon and an aberrant splicing event that impacts SHANK2 transcript variants distinctly. Induced pluripotent stem cells (iPSCs) carrying this mutation, from the patient or isogenic editing, fail to differentiate into functional dopamine (DA) neurons, which can be rescued by genetic correction. Available SMART-Seq single-cell data from human midbrain reveals the abundance of SHANK2B transcript in the ALDH1A1 negative DA neurons. We then show that SHANK2BY29X mutation primarily affects SHANK2B expression and ALDH1A1 negative DA neurons in vitro during early neuronal developmental stage. Mice knocked in with the identical mutation exhibit autistic-like behavior, decreased occupancy of ALDH1A1 negative DA neurons and decreased dopamine release in ventral tegmental area (VTA). Our study provides novel insights on a SHANK2 mutation derived from autism patient and highlights SHANK2B significance in ALDH1A1 negative DA neuron.

Systematic calibration of epitranscriptomic maps using a synthetic modification-free RNA library.

Recent years have witnessed rapid progress in the field of epitranscriptomics. Functional interpretation of the epitranscriptome relies on sequencing technologies that determine the location and stoichiometry of various RNA modifications. However, contradictory results have been reported among studies, bringing the biological impacts of certain RNA modifications into doubt. Here, we develop a synthetic RNA library resembling the endogenous transcriptome but without any RNA modification. By incorporating this modification-free RNA library into established mapping techniques as a negative control, we reveal abundant false positives resulting from sequence bias or RNA structure. After calibration, precise and quantitative mapping expands the understanding of two representative modification types, N6-methyladenosine (m6A) and 5-methylcytosine (m5C). We propose that this approach provides a systematic solution for the calibration of various RNA-modification mappings and holds great promise in epitranscriptomic studies.

Relaxed 3D genome conformation facilitates the pluripotent to totipotent-like state transition in embryonic stem cells.

The 3D genome organization is crucial for gene regulation. Although recent studies have revealed a uniquely relaxed genome conformation in totipotent early blastomeres of both fertilized and cloned embryos, how weakened higher-order chromatin structure is functionally linked to totipotency acquisition remains elusive. Using low-input Hi-C, ATAC-seq and ChIP-seq, we systematically examined the dynamics of 3D genome and epigenome during pluripotent to totipotent-like state transition in mouse embryonic stem cells (ESCs). The spontaneously converted 2-cell-embryo-like cells (2CLCs) exhibited more relaxed chromatin architecture compared to ESCs, including global weakening of both enhancer-promoter interactions and TAD insulation. While the former correlated with inactivation of ESC enhancers and down-regulation of pluripotent genes, the latter might facilitate contacts between the putative new enhancers arising in 2CLCs and neighboring 2C genes. Importantly, disruption of chromatin organization by depleting CTCF or the cohesin complex promoted the ESC to 2CLC transition. Our results thus establish a critical role of 3D genome organization in totipotency acquisition.

Transposable element sequence fragments incorporated into coding and noncoding transcripts modulate the transcriptome of human pluripotent stem cells.

Transposable elements (TEs) occupy nearly 40% of mammalian genomes and, whilst most are fragmentary and no longer capable of transposition, they can nevertheless contribute to cell function. TEs within genes transcribed by RNA polymerase II can be copied as parts of primary transcripts; however, their full contribution to mature transcript sequences remains unresolved. Here, using long and short read (LR and SR) RNA sequencing data, we show that 26% of coding and 65% of noncoding transcripts in human pluripotent stem cells (hPSCs) contain TE-derived sequences. Different TE families are incorporated into RNAs in unique patterns, with consequences to transcript structure and function. The presence of TE sequences within a transcript is correlated with TE-type specific changes in its subcellular distribution, alterations in steady-state levels and half-life, and differential association with RNA Binding Proteins (RBPs). We identify hPSC-specific incorporation of endogenous retroviruses (ERVs) and LINE:L1 into protein-coding mRNAs, which generate TE sequence-derived peptides. Finally, single cell RNA-seq reveals that hPSCs express ERV-containing transcripts, whilst differentiating subpopulations lack ERVs and express SINE and LINE-containing transcripts. Overall, our comprehensive analysis demonstrates that the incorporation of TE sequences into the RNAs of hPSCs is more widespread and has a greater impact than previously appreciated.

scDIOR: single cell RNA-seq data IO software.

BACKGROUND: Single-cell RNA sequencing is becoming a powerful tool to identify cell states, reconstruct developmental trajectories, and deconvolute spatial expression. The rapid development of computational methods promotes the insight of heterogeneous single-cell data. An increasing number of tools have been provided for biological analysts, of which two programming languages- R and Python are widely used among researchers. R and Python are complementary, as many methods are implemented specifically in R or Python. However, the different platforms immediately caused the data sharing and transformation problem, especially for Scanpy, Seurat, and SingleCellExperiemnt. Currently, there is no efficient and user-friendly software to perform data transformation of single-cell omics between platforms, which makes users spend unbearable time on data Input and Output (IO), significantly reducing the efficiency of data analysis. RESULTS: We developed scDIOR for single-cell data transformation between platforms of R and Python based on Hierarchical Data Format Version 5 (HDF5). We have created a data IO ecosystem between three R packages (Seurat, SingleCellExperiment, Monocle) and a Python package (Scanpy). Importantly, scDIOR accommodates a variety of data types across programming languages and platforms in an ultrafast way, including single-cell RNA-seq and spatial resolved transcriptomics data, using only a few codes in IDE or command line interface. For large scale datasets, users can partially load the needed information, e.g., cell annotation without the gene expression matrices. scDIOR connects the analytical tasks of different platforms, which makes it easy to compare the performance of algorithms between them. CONCLUSIONS: scDIOR contains two modules, dior in R and diopy in Python. scDIOR is a versatile and user-friendly tool that implements single-cell data transformation between R and Python rapidly and stably. The software is freely accessible at https://github.com/JiekaiLab/scDIOR .

AP-1 activity is a major barrier of human somatic cell reprogramming.

Human induced pluripotent stem cells (iPSCs) technology has been widely applied to cell regeneration and disease modeling. However, most mechanism of somatic reprogramming is studied on mouse system, which is not always generic in human. Consequently, the generation of human iPSCs remains inefficient. Here, we map the chromatin accessibility dynamics during the induction of human iPSCs from urine cells. Comparing to the mouse system, we found that the closing of somatic loci is much slower in human. Moreover, a conserved AP-1 motif is highly enriched among the closed loci. The introduction of AP-1 repressor, JDP2, enhances human reprogramming and facilitates the reactivation of pluripotent genes. However, ESRRB, KDM2B and SALL4, several known pluripotent factors promoting mouse somatic reprogramming fail to enhance human iPSC generation. Mechanistically, we reveal that JDP2 promotes the closing of somatic loci enriching AP-1 motifs to enhance human reprogramming. Furthermore, JDP2 can rescue reprogramming deficiency without MYC or KLF4. These results indicate AP-1 activity is a major barrier to prevent chromatin remodeling during somatic cell reprogramming.

Concurrent binding to DNA and RNA facilitates the pluripotency reprogramming activity of Sox2.

Some transcription factors that specifically bind double-stranded DNA appear to also function as RNA-binding proteins. Here, we demonstrate that the transcription factor Sox2 is able to directly bind RNA in vitro as well as in mouse and human cells. Sox2 targets RNA via a 60-amino-acid RNA binding motif (RBM) positioned C-terminally of the DNA binding high mobility group (HMG) box. Sox2 can associate with RNA and DNA simultaneously to form ternary RNA/Sox2/DNA complexes. Deletion of the RBM does not affect selection of target genes but mitigates binding to pluripotency related transcripts, switches exon usage and impairs the reprogramming of somatic cells to a pluripotent state. Our findings designate Sox2 as a multi-functional factor that associates with RNA whilst binding to cognate DNA sequences, suggesting that it may co-transcriptionally regulate RNA metabolism during somatic cell reprogramming.

Epigenetic control of transposable elements and cell fate decision.

Transposable elements (TEs) are the most prevalent elements in mammalian genomes. Although potential risks for genome stability, they are a pool of potential regulatory sequences, chromatin control elements, protein-coding genes, and substrates for evolutionary processes. Consequently, a delicate balance is maintained between the potential benefits and deleterious aspects of TEs, and this balance is mediated by the epigenetic regulatory system. In this review, we introduce the role of heterochromatin associated epigentic modifications such as histone 3 lysine 9 trimethylation (H3K9me3) and DNA methylation in the silencing of TEs as well as epigenetic modifications such as histone 3 lysine 4 monomethylation (H3K4me1) and histone 3 lysine 27 acetylation (H3K27ac) in activation of TEs. Further, we elaborate the functions of TEs as binding sites of transcription factors and as anchors of chromosomal conformation in regulation of gene expression. We introduce the impact of TEs on the process of cell fate determination including natural embryonic development in vivo and artificial cell fate transition in vitro. We discuss the main challenges associated with computational TEs analysis and TEs functions exploration, as well as the different experimental and computational strategies in studying these processes. In all, this article provides a comprehensive review of the research advances and existing problems in study of transposable elements in epigenetic regulatory mechanisms, gene transcriptional regulation, and cell fate determination, thereby providing some references for researchers in the fields.

Vitamin C-dependent lysine demethylase 6 (KDM6)-mediated demethylation promotes a chromatin state that supports the endothelial-to-hematopoietic transition.

Hematopoietic stem cells (HSCs)/progenitor cells (HPCs) are generated from hemogenic endothelial cells (HECs) during the endothelial-to-hematopoietic transition (EHT); however, the underlying mechanism remains poorly understood. Here, using an array of approaches, including CRSPR/Cas9 gene knockouts, RNA-Seq, ChIP-Seq, ATAC-Seq etc., we report that vitamin C (Vc) is essential in HPC generation during human pluripotent stem cell (hPSC) differentiation in defined culture conditions. Mechanistically, we found that the endothelial cells generated in the absence of Vc fail to undergo the EHT because of an apparent failure in opening up genomic loci essential for hematopoiesis. Under Vc deficiency, these loci exhibited abnormal accumulation of histone H3 trimethylation at Lys-27 (H3K27me3), a repressive histone modification that arose because of lower activities of demethylases that target H3K27me3. Consistently, deletion of the two H3K27me3 demethylases, Jumonji domain-containing 3 (JMJD3 or KDM6B) and histone demethylase UTX (UTX or KDM6A), impaired HPC generation even in the presence of Vc. Furthermore, we noted that Vc and jmjd3 are also important for HSC generation during zebrafish development. Together, our findings reveal an essential role for Vc in the EHT for hematopoiesis, and identify KDM6-mediated chromatin demethylation as an important regulatory mechanism in hematopoietic cell differentiation.

Models of global gene expression define major domains of cell type and tissue identity.

The current classification of cells in an organism is largely based on their anatomic and developmental origin. Cells types and tissues are traditionally classified into those that arise from the three embryonic germ layers, the ectoderm, mesoderm and endoderm, but this model does not take into account the organization of cell type-specific patterns of gene expression. Here, we present computational models for cell type and tissue specification derived from a collection of 921 RNA-sequencing samples from 272 distinct mouse cell types or tissues. In an unbiased fashion, this analysis accurately predicts the three known germ layers. Unexpectedly, this analysis also suggests that in total there are eight major domains of cell type-specification, corresponding to the neurectoderm, neural crest, surface ectoderm, endoderm, mesoderm, blood mesoderm, germ cells and the embryonic domain. Further, we identify putative genes responsible for specifying the domain and the cell type. This model has implications for understanding trans-lineage differentiation for stem cells, developmental cell biology and regenerative medicine.

Passive DNA demethylation preferentially up-regulates pluripotency-related genes and facilitates the generation of induced pluripotent stem cells.

A high proliferation rate has been observed to facilitate somatic cell reprogramming, but the pathways that connect proliferation and reprogramming have not been reported. DNA methyltransferase 1 (DNMT1) methylates hemimethylated CpG sites produced during S phase and maintains stable inheritance of DNA methylation. Impairing this process results in passive DNA demethylation. In this study, we show that the cell proliferation rate positively correlated with the expression of Dnmt1 in G1 phase. In addition, as determined by whole-genome bisulfate sequencing and high-performance liquid chromatography, global DNA methylation of mouse embryonic fibroblasts was significantly higher in G1 phase than in G2/M phase. Thus, we suspected that high cellular proliferation requires more Dnmt1 expression in G1 phase to prevent passive DNA demethylation. The methylation differences of individual CpG sites between G1 and G2/M phase were related to the methylation status and the positions of their surrounding CpG sites. In addition, larger methylation differences were observed on the promoters of pluripotency-related genes; for example, Oct4, Nanog, Sox2, Esrrb, Cdh1, and Epcam When such methylation differences or passive DNA demethylation accumulated with Dnmt1 suppression and proliferation acceleration, DNA methylation on pluripotency-related genes was decreased, and their expression was up-regulated, which subsequently promoted pluripotency and mesenchymal-epithelial transition, a necessary step for reprogramming. We infer that high cellular proliferation rates promote generation of induced pluripotent stem cells at least partially by inducing passive DNA demethylation and up-regulating pluripotency-related genes. Therefore, these results uncover a connection between cell reprogramming and DNA methylation.

Chemical reprogramming of mouse embryonic and adult fibroblast into endoderm lineage.

We report here an approach to redirecting somatic cell fate under chemically defined conditions without transcription factors. We start by converting mouse embryonic fibroblasts to epithelial-like cells with chemicals and growth factors. Subsequent cell fate mapping reveals a robust induction of SOX17 in the resulting epithelial-like cells that can be further reprogrammed to endodermal progenitor cells. Interestingly, these cells can self-renew in vitro and further differentiate into albumin-producing hepatocytes that can rescue mice from acute liver injury. Our results demonstrate a rational approach to convert mouse embryonic fibroblasts to hepatocytes and suggest that this mechanism-driven approach may be generalized for other cells.

Gadd45a is a heterochromatin relaxer that enhances iPS cell generation.

Reprogramming of somatic cells to induced pluripotent stem cells rewrites the code of cell fate at the chromatin level. Yet, little is known about this process physically. Here, we describe a fluorescence recovery after photobleaching method to assess the dynamics of heterochromatin/euchromatin and show significant heterochromatin loosening at the initial stage of reprogramming. We identify growth arrest and DNA damage-inducible protein a (Gadd45a) as a chromatin relaxer in mouse embryonic fibroblasts, which also enhances somatic cell reprogramming efficiency. We show that residue glycine 39 (G39) in Gadd45a is essential for interacting with core histones, opening chromatin and enhancing reprogramming. We further demonstrate that Gadd45a destabilizes histone-DNA interactions and facilitates the binding of Yamanaka factors to their targets for activation. Our study provides a method to screen factors that impact on chromatin structure in live cells, and identifies Gadd45a as a chromatin relaxer.