COVID-19 immune features revealed by a large-scale single-cell transcriptome atlas.

A dysfunctional immune response in coronavirus disease 2019 (COVID-19) patients is a recurrent theme impacting symptoms and mortality, yet a detailed understanding of pertinent immune cells is not complete. We applied single-cell RNA sequencing to 284 samples from 196 COVID-19 patients and controls and created a comprehensive immune landscape with 1.46 million cells. The large dataset enabled us to identify that different peripheral immune subtype changes are associated with distinct clinical features, including age, sex, severity, and disease stages of COVID-19. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA was found in diverse epithelial and immune cell types, accompanied by dramatic transcriptomic changes within virus-positive cells. Systemic upregulation of S100A8/A9, mainly by megakaryocytes and monocytes in the peripheral blood, may contribute to the cytokine storms frequently observed in severe patients. Our data provide a rich resource for understanding the pathogenesis of and developing effective therapeutic strategies for COVID-19.

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.

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.

Regeneration of the human segmentation clock in somitoids in vitro.

Each vertebrate species appears to have a unique timing mechanism for forming somites along the vertebral column, and the process in human remains poorly understood at the molecular level due to technical and ethical limitations. Here, we report the reconstitution of human segmentation clock by direct reprogramming. We first reprogrammed human urine epithelial cells to a presomitic mesoderm (PSM) state capable of long-term self-renewal and formation of somitoids with an anterior-to-posterior axis. By inserting the RNA reporter Pepper into HES7 and MESP2 loci of these iPSM cells, we show that both transcripts oscillate in the resulting somitoids at ~5 h/cycle. GFP-tagged endogenous HES7 protein moves along the anterior-to-posterior axis during somitoid formation. The geo-sequencing analysis further confirmed anterior-to-posterior polarity and revealed the localized expression of WNT, BMP, FGF, and RA signaling molecules and HOXA-D family members. Our study demonstrates the direct reconstitution of human segmentation clock from somatic cells, which may allow future dissection of the mechanism and components of such a clock and aid regenerative medicine.

Sirtuin 6 ameliorates bleomycin-induced pulmonary fibrosis via activation of lipid catabolism.

Pulmonary fibrosis is a chronic and serious interstitial lung disease with little effective therapies currently. Our incomplete understanding of its pathogenesis remains obstacles in therapeutic developments. Sirtuin 6 (SIRT6) has been shown to mitigate multiple organic fibrosis. However, the involvement of SIRT6-mediated metabolic regulation in pulmonary fibrosis remains unclear. Here, we demonstrated that SIRT6 was predominantly expressed in alveolar epithelial cells in human lung tissues by using a single-cell sequencing database. We showed that SIRT6 protected against bleomycin-induced injury of alveolar epithelial cells in vitro and pulmonary fibrosis of mice in vivo. High-throughput sequencing revealed enriched lipid catabolism in Sirt6 overexpressed lung tissues. Mechanismly, SIRT6 ameliorates bleomycin-induced ectopic lipotoxicity by enhancing lipid degradation, thereby increasing the energy supply and reducing the levels of lipid peroxides. Furthermore, we found that peroxisome proliferator-activated receptor α (PPARα) was essential for SIRT6-mediated lipid catabolism, anti-inflammatory responses, and antifibrotic signaling. Our data suggest that targeting SIRT6-PPARα-mediated lipid catabolism could be a potential therapeutic strategy for diseases complicated with pulmonary fibrosis.

Single-cell landscape of peripheral immune response in patients with anti-melanoma differentiation-associated gene 5 dermatomyositis.

OBJECTIVE: Anti-melanoma differentiation-associated gene 5 (MDA5)-positive dermatomyositis (DM) is a rare but life-threatening autoimmune disorder with a high risk to develop rapidly progressive interstitial lung disease. Current empirical therapies have limited improvement on patients' survival, as little is known about the aetiology of MDA5 DM. To best understand its immune landscape, we applied single-cell RNA sequencing (scRNA-seq) to peripheral blood samples from DM patients and healthy controls. METHODS: Peripheral blood mononuclear cells (PBMCs) from eight DM patients, comprising three distinct subtypes, as well as two healthy donors, were sequenced by 10X Genomics platform. Additional scRNA-seq data of four healthy donors were incorporated for further bioinformatic analysis. RESULTS: Aberrant increased proportions of CD14+ monocyte and plasma cells were observed in MDA5 DM samples. Moreover, we found overactivated type I interferon response and antiviral immunity in both innate and adaptive immune cells derived from MDA5 DM patients, which was positively correlated with disease severity. Importantly, a unique subset of CD14+ monocyte that highly expressed interferon alpha-inducible protein 27 (IFI27, a biomarker for viral infection) and interferon induced with helicase C domain 1 (IFIH1, encodes MDA5) was specifically identified in MDA5 DM samples for the first time. CONCLUSION: Our study demonstrates the peripheral immune cell atlas of different DM subtypes, provides compelling evidence for viral infection-derived origin of MDA5 DM, and offers potential targets for innovative therapeutic interventions.

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.

Long-term Remission of Acute Intermittent Porphyria Treated with Gonadotropin-Releasing Hormone Analogues and Estrogen: a Case Report.

BACKGROUND: Acute intermittent porphyria (AIP) is an autosomal dominant hepatic porphyria characterized by a partial deficiency of hydroxymethylbilane synthase involved in heme biosynthesis. It is difficult for all patients to achieve complete control of AIP episodes. METHOD: We report on a 20-year-old female woman who suffered from recurrent abdominal pain and was diagnosed as "acute intermittent porphyria". She failed to respond to conventional symptomatic treatment and subsequently was treated with gonadotropin-releasing hormone analogues (GnRH) combined with estrogen for one year. RESULT: The case did not experience acute attacks and obtained long-term clinical remission to date. CONCLUSIONS: GnRH combined with estrogen, one of the treatment options for menstrual-associated AIP, might induce long-term remission.

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.

Homocystinuria in a Family with Novel Cystathionine Beta Synthase Gene Mutations.

BACKGROUND: Classic homocystinuria is caused by cystathionine beta synthase deficiency owing to genetic mutations. The most common symptoms are ectopia lentis, osteoporosis, thrombosis, and mental retardation. This disease is prone to misdiagnosis and delayed diagnosis. METHODS: Here, we report a 19-year-old woman with Marfan's morphotype, high blood homocysteine, and a history of ectopia lentis. Total homocysteine levels became normal following treatment with vitamin therapy. RESULTS: Genetic analysis revealed two heterozygous nucleotide mutations in the parents. The mutation from the patient's father had not been described previously. CONCLUSIONS: Screening for blood homocysteine should be performed early. Early diagnosis and treatment can prevent related symptoms.

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.

The Relative Importance of Clinical and Socio-demographic Variables in Prognostic Prediction in Non-Small Cell Lung Cancer: A Variable Importance Approach.

BACKGROUND: Prognostic modeling in health care has been predominantly statistical, despite a rapid growth of literature on machine-learning approaches in biological data analysis. We aim to assess the relative importance of variables in predicting overall survival among patients with non-small cell lung cancer using a Variable Importance (VIMP) approach in a machine-learning Random Survival Forest (RSF) model for posttreatment planning and follow-up. METHODS: A total of 935 non-small cell lung cancer patients were randomly and equally divided into 2 training and testing cohorts in an RFS model. The prognostic variables included age, sex, race, the TNM Classification of Malignant Tumors (TNM) stage, smoking history, Eastern Cooperative Oncology Group performance status, histologic type, treatment category, maximum standard uptake value of whole-body tumor (SUVmaxWB), whole-body metabolic tumor volume (MTVwb), and Charlson Comorbidity Index. The VIMP was calculated using a permutation method in the RSF model. We further compared the VIMP of the RSF model to that of the standard Cox survival model. We examined the order of VIMP with the differential functional forms of the variables. RESULTS: In both the RSF and the standard Cox models, the most important variables are treatment category, TNM stage, and MTVwb. The order of VIMP is more robust in RSF model than in Cox model regarding the differential functional forms of the variables. CONCLUSIONS: The RSF VIMP approach can be applied alongside with the Cox model to further advance the understanding of the roles of prognostic factors, and improve prognostic precision and care efficiency.

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.

Osteoclast stimulatory transmembrane protein induces a phenotypic switch in macrophage polarization suppressing an M1 pro-inflammatory state.

Macrophages are the key cells in metabolic syndrome and are also a risk factor for metabolic disease. Macrophages have different functions and transcriptional profiles, but all are required for maintaining homeostasis. It is well known that macrophages play a key role in inflammation and early atherogenesis, and are present in two phenotypes: pro-inflammatory (M1) and anti-inflammatory (M2). Osteoclast stimulatory transmembrane protein (oc-stamp) is a multiple-pass transmembrane protein; however, its function remains unclear. In this study, we explored the role of oc-stamp in macrophages physiology. The results showed that oc-stamp was notably decreased under LPS and IFN-γ stimulation, while it was increased with IL-4 treatment. Furthermore, oc-stamp induced a phenotypic switch in macrophage polarization, suppressing the M1 pro-inflammatory state in the overexpression group, and promoting the M1 pro-inflammatory state in the knockdown group. Further study revealed that oc-stamp regulated macrophage polarization possibly via STAT6. Taken together, our results are the first to demonstrate that oc-stamp may play an important role in macrophage polarization and inhibit the M1 pro-inflammatory state.

TFE3 Alleviates Hepatic Steatosis through Autophagy-Induced Lipophagy and PGC1α-Mediated Fatty Acid ß-Oxidation.

Autophagy flux deficiency is closely related to the development of hepatic steatosis. Transcription factor E3 (TFE3) is reported to be a crucial gene that regulates autophagy flux and lysosome function. Therefore, we investigated the role of TFE3 in a cell model of hepatic steatosis. We constructed L02 hepatocyte lines that stably over-expressed or knocked down the expression of TFE3. Subsequently, the effects of TFE3 on hepatocellular lipid metabolism were determined by autophagy flux assay, lipid oil red O (ORO) staining, immunofluorescence staining, and mitochondrial ß-oxidation assessment. Finally, we analyzed whether peroxisome proliferative activated receptor gamma coactivator 1α (PGC1α) was the potential target gene of TFE3 in the regulation of hepatic steatosis using a chromatin immunoprecipitation (CHIP) assay and a luciferase reporter system. We found that overexpression of TFE3 markedly alleviated hepatocellular steatosis. On the contrary, downregulation of TFE3 resulted in an aggravated steatosis. The mechanistic studies revealed that the TFE3-manipulated regulatory effects on hepatocellular steatosis are dependent on autophagy-induced lipophagy and PGC1α-mediated fatty acid ß-oxidation because blocking these pathways with an Atg5 small interfering RNA (siRNA) or PGC1α siRNA dramatically blunted the TFE3-mediated regulation of steatosis. In conclusion, TFE3 gene provides a novel insight into the treatment of hepatic steatosis and other metabolic disease.

[Three-dimensional CT liver image segmentation based on hierarchical contextual active contour].

In this paper, we propose a new active contour algorithm, i. e. hierarchical contextual active contour (HCAC), and apply it to automatic liver segmentation from three-dimensional CT (3D-CT) images. HCAC is a learning-based method and can be divided into two stages. At the first stage, i.e. the training stage, given a set of abdominal 3D-CT training images and the corresponding manual liver labels, we tried to establish a mapping between automatic segmentations (in each round) and manual reference segmentations via context features, and obtained a series of self-correcting classifiers. At the second stage, i.e. the segmentation stage, we firstly used the basic active contour to segment the image and subsequently used the contextual active contour (CAC) iteratively, which combines the image information and the current shape model, to improve the segmentation result. The current shape model is produced by the corresponding self-correcting classifier (the input is the previous automatic segmentation result). The proposed method was evaluated on the datasets of MICCAI 2007 liver segmentation challenge. The experimental results showed that we would get more and more accurate segmentation results by the iterative steps and the satisfied results would be obtained after about six rounds of iterations.

Protocol for siRNA-mediated TET1 knockdown during differentiation of human embryonic stem cells into definitive endoderm cells.

TET1-mediated active DNA demethylation is required for endogenous retrovirus (ERV) enhancer activation during human ES differentiation into definitive endoderm (DE) cells. Here we present a protocol for siRNA-mediated TET1 knockdown during this process to decipher TET1's role in ERV activation and DE differentiation. We describe steps for inducing ES into DE cells. We then detail steps for knocking down TET1 during differentiation and for examining the effects of TET1 knockdown on LTR6B methylation, cell morphology, and gene expression. For complete details on the use and execution of this protocol, please refer to Wu et al. (2022).1.

Exploring the correlation between genetic transcription and multi-temporal developmental autism spectrum disorder using resting-state functional magnetic resonance imaging.

Introduction: The present investigation aimed to explore the neurodevelopmental trajectory of autism spectrum disorder (ASD) by identifying the changes in brain function and gene expression associated with the disorder. Previous studies have indicated that ASD is a highly inherited neurodevelopmental disorder of the brain that displays symptom heterogeneity across different developmental periods. However, the transcriptomic changes underlying these developmental differences remain largely unknown. Methods: To address this gap in knowledge, our study employed resting-state functional magnetic resonance imaging (rs-fMRI) data from a large sample of male participants across four representative age groups to stratify the abnormal changes in brain function associated with ASD. Partial least square regression (PLSr) was utilized to identify unique changes in gene expression in brain regions characterized by aberrant functioning in ASD. Results: Our results revealed that ASD exhibits distinctive developmental trajectories in crucial brain regions such as the default mode network (DMN), temporal lobe, and prefrontal lobes during critical periods of neurodevelopment when compared to the control group. These changes were also associated with genes primarily located in synaptic tissues. Discussion: The findings of this study suggest that the neurobiology of ASD is uniquely heterogeneous across different ages and may be accompanied by distinct molecular mechanisms related to gene expression.

Transposable elements activation triggers necroptosis in mouse embryonic stem cells.

Deficiency of the histone H3K9 methyltransferase SETDB1 induces RIPK3-dependent necroptosis in mouse embryonic stem cells (mESCs). However, how necroptosis pathway is activated in this process remains elusive. Here we report that the reactivation of transposable elements (TEs) upon SETDB1 knockout is responsible for the RIPK3 regulation through both cis and trans mechanisms. IAPLTR2_Mm and MMERVK10c-int, both of which are suppressed by SETDB1-dependent H3K9me3, act as enhancer-like cis-regulatory elements and their RIPK3 nearby members enhance RIPK3 expression when SETDB1 is knockout. Moreover, reactivated endogenous retroviruses generate excessive viral mimicry, which promotes necroptosis mainly through Z-DNA-binding protein 1 (ZBP1). These results indicate TEs play an important role in regulating necroptosis.