Transposable elements (TEs) make up a majority of a typical eukaryote’s genome, and contribute to cell heterogeneity in unclear ways. Single-cell sequencing technologies are powerful tools to explore cells, however analysis is typically gene-centric and TE expression has not been addressed. Here, we develop a single-cell TE processing pipeline, scTE, and report the expression of TEs in single cells in a range of biological contexts. Specific TE types are expressed in subpopulations of embryonic stem cells and are dynamically regulated during pluripotency reprogramming, differentiation, and embryogenesis. Unexpectedly, TEs are expressed in somatic cells, including human disease-specific TEs that are undetectable in bulk analyses. Finally, we apply scTE to single-cell ATAC-seq data, and demonstrate that scTE can discriminate cell type using chromatin accessibly of TEs alone. Overall, our results classify the dynamic patterns of TEs in single cells and their contributions to cell heterogeneity.

The majority of mammalian genomes are devoted to transposable elements (TEs). Whilst TEs are increasingly recognized for their important biological functions, they are a potential danger to genomic stability and are carefully regulated by the epigenetic system. However, the full complexity of this regulatory system is not understood. Here, using mouse embryonic stem cells, we show that TEs are suppressed by heterochromatic marks like H3K9me3, and are also labelled by all major types of chromatin modification in complex patterns, including bivalent activatory and repressive marks. We identified 29 epigenetic modifiers that significantly deregulated at least one type of TE. The loss of Setdb1, Ncor2, Rnf2, Kat5, Prmt5, Uhrf1, and Rrp8 caused widespread changes in TE expression and chromatin accessibility. These effects were context-specific, with different chromatin modifiers regulating the expression and chromatin accessibility of specific subsets of TEs. Our work reveals the complex patterns of epigenetic regulation of TEs.

The coronavirus disease 2019 (COVID-19) pandemic is caused by infection with the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which is spread primary via respiratory droplets and infects the lungs. Currently widely used cell lines and animals are unable to accurately mimic human physiological conditions because of the abnormal status of cell lines (transformed or cancer cells) and species differences between animals and humans. Organoids are stem cell-derived self-organized three-dimensional culture in vitro and model the physiological conditions of natural organs. Here we showed that SARS-CoV-2 infected and extensively replicated in human embryonic stem cells (hESCs)-derived lung organoids, including airway and alveolar organoids which covered the complete infection and spread route for SARS-CoV-2 within lungs. The infected cells were ciliated, club, and alveolar type 2 (AT2) cells, which were sequentially located from the proximal to the distal airway and terminal alveoli, respectively. Additionally, RNA-seq revealed early cell response to virus infection including an unexpected downregulation of the metabolic processes, especially lipid metabolism, in addition to the well-known upregulation of immune response. Further, Remdesivir and a human neutralizing antibody potently inhibited SARS-CoV-2 replication in lung organoids. Therefore, human lung organoids can serve as a pathophysiological model to investigate the underlying mechanism of SARS-CoV-2 infection and to discover and test therapeutic drugs for COVID-19.