Identification of the metaphyseal skeletal stem cell building trabecular bone.

Skeletal stem cells (SSCs) that are capable of self-renewal and multipotent differentiation contribute to bone development and homeostasis. Several populations of SSCs at different skeletal sites have been reported. Here, we identify a metaphyseal SSC (mpSSC) population whose transcriptional landscape is distinct from other bone mesenchymal stromal cells (BMSCs). These mpSSCs are marked by Sstr2 or Pdgfrb+Kitl-, located just underneath the growth plate, and exclusively derived from hypertrophic chondrocytes (HCs). These HC-derived mpSSCs have properties of self-renewal and multipotency in vitro and in vivo, producing most HC offspring postnatally. HC-specific deletion of Hgs, a component of the endosomal sorting complex required for transport, impairs the HC-to-mpSSC conversion and compromises trabecular bone formation. Thus, mpSSC is the major source of BMSCs and osteoblasts in bone marrow, supporting the postnatal trabecular bone formation.

Transformer for one stop interpretable cell type annotation.

Consistent annotation transfer from reference dataset to query dataset is fundamental to the development and reproducibility of single-cell research. Compared with traditional annotation methods, deep learning based methods are faster and more automated. A series of useful single cell analysis tools based on autoencoder architecture have been developed but these struggle to strike a balance between depth and interpretability. Here, we present TOSICA, a multi-head self-attention deep learning model based on Transformer that enables interpretable cell type annotation using biologically understandable entities, such as pathways or regulons. We show that TOSICA achieves fast and accurate one-stop annotation and batch-insensitive integration while providing biologically interpretable insights for understanding cellular behavior during development and disease progressions. We demonstrate TOSICA's advantages by applying it to scRNA-seq data of tumor-infiltrating immune cells, and CD14+ monocytes in COVID-19 to reveal rare cell types, heterogeneity and dynamic trajectories associated with disease progression and severity.

De novo generation of macrophage from placenta-derived hemogenic endothelium.

Macrophages play pivotal roles in immunity, hematopoiesis, and tissue homeostasis. In mammals, macrophages have been shown to originate from yolk-sac-derived erythro-myeloid progenitors and aorta-gonad-mesonephros (AGM)-derived hematopoietic stem cells. However, whether macrophages can arise from other embryonic sites remains unclear. Here, using single-cell RNA sequencing, we profile the transcriptional landscape of mouse fetal placental hematopoiesis. We uncover and experimentally validate that a CD44+ subpopulation of placental endothelial cells (ECs) exhibits hemogenic potential. Importantly, lineage tracing using the newly generated Hoxa13 reporter line shows that Hoxa13-labeled ECs can produce placental macrophages, named Hofbauer cell (HBC)-like cells. Furthermore, we identify two subtypes of HBC-like cells, and cell-cell interaction analysis identifies their potential roles in angiogenesis and antigen presentation, separately. Our study provides a comprehensive understanding of placental hematopoiesis and highlights the placenta as a source of macrophages, which has important implications for both basic and translational research.

A single-cell resolution developmental atlas of hematopoietic stem and progenitor cell expansion in zebrafish.

During vertebrate embryogenesis, fetal hematopoietic stem and progenitor cells (HSPCs) exhibit expansion and differentiation properties in a supportive hematopoietic niche. To profile the developmental landscape of fetal HSPCs and their local niche, here, using single-cell RNA-sequencing, we deciphered a dynamic atlas covering 28,777 cells and 9 major cell types (23 clusters) of zebrafish caudal hematopoietic tissue (CHT). We characterized four heterogeneous HSPCs with distinct lineage priming and metabolic gene signatures. Furthermore, we investigated the regulatory mechanism of CHT niche components for HSPC development, with a focus on the transcription factors and ligand-receptor networks involved in HSPC expansion. Importantly, we identified an endothelial cell-specific G protein-coupled receptor 182, followed by in vivo and in vitro functional validation of its evolutionally conserved role in supporting HSPC expansion in zebrafish and mice. Finally, comparison between zebrafish CHT and human fetal liver highlighted the conservation and divergence across evolution. These findings enhance our understanding of the regulatory mechanism underlying hematopoietic niche for HSPC expansion in vivo and provide insights into improving protocols for HSPC expansion in vitro.

The Effects of Graded Levels of Calorie Restriction: XVI. Metabolomic Changes in the Cerebellum Indicate Activation of Hypothalamocerebellar Connections Driven by Hunger Responses.

Calorie restriction (CR) remains the most robust intervention to extend life span and improve healthspan. Though the cerebellum is more commonly associated with motor control, it has strong links with the hypothalamus and is thought to be associated with nutritional regulation and adiposity. Using a global mass spectrometry-based metabolomics approach, we identified 756 metabolites that were significantly differentially expressed in the cerebellar region of the brain of C57BL/6J mice, fed graded levels of CR (10, 20, 30, and 40 CR) compared to mice fed ad libitum for 12 hours a day. Pathway enrichment indicated changes in the pathways of adenosine and guanine (which are precursors of DNA production), aromatic amino acids (tyrosine, phenylalanine, and tryptophan) and the sulfur-containing amino acid methionine. We also saw increases in the tricarboxylic acid cycle (TCA) cycle, electron donor, and dopamine and histamine pathways. In particular, changes in l-histidine and homocarnosine correlated positively with the level of CR and food anticipatory activity and negatively with insulin and body temperature. Several metabolic and pathway changes acted against changes seen in age-associated neurodegenerative disorders, including increases in the TCA cycle and reduced l-proline. Carnitine metabolites contributed to discrimination between CR groups, which corroborates previous work in the liver and plasma. These results indicate the conservation of certain aspects of metabolism across tissues with CR. Moreover, this is the first study to indicate CR alters the cerebellar metabolome, and does so in a graded fashion, after only a short period of restriction.

Ribosome 18S m6A Methyltransferase METTL5 Promotes Translation Initiation and Breast Cancer Cell Growth.

N6 methylation at adenosine 1832 (m6A1832) of mammalian 18S rRNA, occupying a critical position within the decoding center, is modified by a conserved methyltransferase, METTL5. Here, we find that METTL5 shows strong substrate preference toward the 18S A1832 motif but not the other reported m6A motifs. Comparison with a yeast ribosome structural model unmodified at this site indicates that the modification may facilitate mRNA binding by inducing conformation changes in the mammalian ribosomal decoding center. METTL5 promotes p70-S6K activation and proper translation initiation, and the loss of METTL5 significantly reduces the abundance of polysome. METTL5 expression is elevated in breast cancer patient samples and is required for growth of several breast cancer cell lines. We further find that Caenorhabditis elegans lacking the homolog metl-5 develop phenotypes known to be associated with impaired translation. Altogether, our findings uncover critical and conserved roles of METTL5 in the regulation of translation.

The precursor of PI(3,4,5)P3 alleviates aging by activating daf-18(Pten) and independent of daf-16.

Aging is characterized by the loss of homeostasis and the general decline of physiological functions, accompanied by various degenerative diseases and increased rates of mortality. Aging targeting small molecule screens have been performed many times, however, few have focused on endogenous metabolic intermediates-metabolites. Here, using C. elegans lifespan assays, we conducted a worm metabolite screen and identified an eukaryotes conserved metabolite, myo-inositol (MI), to extend lifespan, increase mobility and reduce fat content. Genetic analysis of enzymes in MI metabolic pathway suggest that MI alleviates aging through its derivative PI(4,5)P2. MI and PI(4,5)P2 are precursors of PI(3,4,5)P3, which is negatively related to longevity. The longevity effect of MI is dependent on the tumor suppressor gene, daf-18 (homologous to mouse Pten), independent of its classical pathway downstream genes, akt or daf-16. Furthermore, we found MI effects on aging and lifespan act through mitophagy regulator PTEN induced kinase-1 (pink-1) and mitophagy. MI's anti-aging effect is also conserved in mouse, indicating a conserved mechanism in mammals.

Individual variation of the SARS-CoV-2 receptor ACE2 gene expression and regulation.

The COVID-19 coronavirus is now spreading worldwide. Its pathogen, SARS-CoV-2, has been shown to use angiotensin-converting enzyme 2 (ACE2) as its host cell receptor, same as the severe acute respiratory syndrome coronavirus (SARS-CoV) in 2003. Epidemiology studies found males although only slightly more likely to be infected than females account for the majority of the severely ill and fatality, which also bias for people older than 60 years or with metabolic and cardiovascular diseases. Here by analyzing GTEx and other public data in 30 tissues across thousands of individuals, we found a significantly higher level in Asian females, an age-dependent decrease in all ethnic groups, and a highly significant decrease in type II diabetic patients of ACE2 expression. Consistently, the most significant expression quantitative loci (eQTLs) contributing to high ACE2 expression are close to 100% in East Asians, >30% higher than other ethnic groups. A shockingly common enrichment of viral infection pathways was found among ACE2 anti-expressed genes, and multiple binding sites of virus infection related transcription factors and sex hormone receptors locate at ACE2 regulatory regions. Human and mice data analysis further revealed ACE2 expression is reduced in T2D patients and with inflammatory cytokine treatment and upregulated by estrogen and androgen (both decrease with age). Our findings revealed a negative correlation between ACE2 expression and COVID-19 fatality at both population and molecular levels. These results will be instrumental when designing potential prevention and treatment strategies for ACE2 binding coronaviruses in general.

Characterization of functional transposable element enhancers in acute myeloid leukemia.

Transposable elements (TEs) have been shown to have important gene regulatory functions and their alteration could lead to disease phenotypes. Acute myeloid leukemia (AML) develops as a consequence of a series of genetic changes in hematopoietic precursor cells, including mutations in epigenetic factors. Here, we set out to study the gene regulatory role of TEs in AML. We first explored the epigenetic landscape of TEs in AML patients using ATAC-seq data. We show that a large number of TEs in general, and more specifically mammalian-wide interspersed repeats (MIRs), are more enriched in AML cells than in normal blood cells. We obtained a similar finding when analyzing histone modification data in AML patients. Gene Ontology enrichment analysis showed that genes near MIRs in open chromatin regions are involved in leukemogenesis. To functionally validate their regulatory role, we selected 19 MIR regions in AML cells, and tested them for enhancer activity in an AML cell line (Kasumi-1) and a chronic myeloid leukemia (CML) cell line (K562); the results revealed several MIRs to be functional enhancers. Taken together, our results suggest that TEs are potentially involved in myeloid leukemogenesis and highlight these sequences as potential candidates harboring AML-associated variation.

The Effects of Graded Levels of Calorie Restriction: XIV. Global Metabolomics Screen Reveals Brown Adipose Tissue Changes in Amino Acids, Catecholamines, and Antioxidants After Short-Term Restriction in C57BL/6 Mice.

Animals undergoing calorie restriction (CR) often lower their body temperature to conserve energy. Brown adipose tissue (BAT) is stimulated through norepinephrine when rapid heat production is needed, as it is highly metabolically active due to the uncoupling of the electron transport chain from ATP synthesis. To better understand how BAT metabolism changes with CR, we used metabolomics to identify 883 metabolites that were significantly differentially expressed in the BAT of C57BL/6 mice, fed graded CR (10%, 20%, 30%, and 40% CR relative to their individual baseline intake), compared with mice fed ad libitum (AL) for 12 hours a day. Pathway analysis revealed that graded CR had an impact on the TCA cycle and fatty acid degradation. In addition, an increase in nucleic acids and catecholamine pathways was seen with graded CR in the BAT metabolome. We saw increases in antioxidants with CR, suggesting a beneficial effect of mitochondrial uncoupling. Importantly, the instigator of BAT activation, norepinephrine, was increased with CR, whereas its precursors l-tyrosine and dopamine were decreased, indicating a shift of metabolites through the activation pathway. Several of these key changes were correlated with food anticipatory activity and body temperature, indicating BAT activation may be driven by responses to hunger.

A 3D Atlas of Hematopoietic Stem and Progenitor Cell Expansion by Multi-dimensional RNA-Seq Analysis.

In vertebrates, hematopoiesis occurring in different niches is orchestrated by intrinsic and extrinsic regulators. Previous studies have revealed numerous linear and planar regulatory mechanisms. However, a multi-dimensional transcriptomic atlas of any given hematopoietic organ has not yet been established. Here, we use multiple RNA sequencing (RNA-seq) approaches, including cell type-specific, temporal bulk RNA-seq, in vivo GEO-seq, and single-cell RNA-seq (scRNA-seq), to characterize the detailed spatiotemporal transcriptome during hematopoietic stem and progenitor cell (HSPC) expansion in the caudal hematopoietic tissue (CHT) of zebrafish. Combinatorial expression profiling reveals that, in the CHT niche, HSPCs and their neighboring supporting cells are co-regulated by shared signaling pathways and intrinsic factors, such as integrin signaling and Smchd1. Moreover, scRNA-seq analysis unveils the strong association between cell cycle status and HSPC differentiation. Taken together, we report a global transcriptome landscape that provides valuable insights and a rich resource to understand HSPC expansion in an intact vertebrate hematopoietic organ.

Immune Cell Types and Secreted Factors Contributing to Inflammation-to-Cancer Transition and Immune Therapy Response.

Although chronic inflammation increases many cancers' risk, how inflammation facilitates cancer development is still not well studied. Recognizing whether and when inflamed tissues transition to cancerous tissues is of utmost importance. To unbiasedly infer molecular events, immune cell types, and secreted factors contributing to the inflammation-to-cancer (I2C) transition, we develop a computational package called "SwitchDetector" based on liver, gastric, and colon cancer I2C data. Using it, we identify angiogenesis associated with a common critical transition stage for multiple I2C events. Furthermore, we infer infiltrated immune cell type composition and their secreted or suppressed extracellular proteins to predict expression of important transition stage genes. This identifies extracellular proteins that may serve as early-detection biomarkers for pre-cancer and early-cancer stages. They alone or together with I2C hallmark angiogenesis genes are significantly related to cancer prognosis and can predict immune therapy response. The SwitchDetector and I2C database are publicly available at www.inflammation2cancer.org.

HER2/EGFR-AKT Signaling Switches TGFß from Inhibiting Cell Proliferation to Promoting Cell Migration in Breast Cancer.

TGFß signaling inhibits cell proliferation to block cancer initiation, yet it also enhances metastasis to promote malignancy during breast cancer development. The mechanisms underlying these differential effects are still unclear. Here, we report that HER2/EGFR signaling switches TGFß function in breast cancer cells from antiproliferation to cancer promotion. Inhibition of HER2/EGFR activity attenuated TGFß-induced epithelial-mesenchymal transition and migration but enhanced the antiproliferative activity of TGFß. Activation of HER2/EGFR induced phosphorylation of Smad3 at Ser208 of the linker region through AKT, which promoted the nuclear accumulation of Smad3 and subsequent expression of the genes related to EMT and cell migration. In contrast, HER2/EGFR signaling had no effects on the nuclear localization of Smad2. Knockdown of Smad3, but not Smad2, blocked TGFß-induced breast cancer cell migration. We observed a positive correlation between the nuclear localization of Smad3 and HER2 activation in advanced human breast cancers. Our results demonstrate a key role for HER2/EGFR in differential regulation of Smad3 activity to shift TGFß function from antitumorigenic to protumorigenic during breast cancer development.Significance: TGFß signaling can shift from inhibiting to promoting breast cancer development via HER2/EGFR AKT-mediated phosphorylation of Smad3 at S208, enhancing its nuclear accumulation and upregulation of EMT-related genes.Graphical Abstract: http://cancerres.aacrjournals.org/content/canres/78/21/6073/F1.large.jpg Cancer Res; 78(21); 6073-85. ©2018 AACR.

The Effects of Graded Levels of Calorie Restriction: X. Transcriptomic Responses of Epididymal Adipose Tissue.

Calorie restriction (CR) leads to a remarkable decrease in adipose tissue mass and increases longevity in many taxa. Since the discovery of leptin, the secretory abilities of adipose tissue have gained prominence in the responses to CR. We quantified transcripts of epididymal white adipose tissue of male C57BL/6 mice exposed to graded levels of CR (0-40% CR) for 3 months. The numbers of differentially expressed genes (DEGs) involved in NF-κB, HIF1-α, and p53 signaling increased with increasing levels of CR. These pathways were all significantly downregulated at 40% CR relative to 12 h ad libitum feeding. In addition, graded CR had a substantial impact on DEGs associated with pathways involved in angiogenesis. Of the 497 genes differentially expressed with graded CR, 155 of these genes included a signal peptide motif. These putative signaling proteins were involved in the response to ketones, TGF-ß signaling, negative regulation of insulin secretion, and inflammation. This accords with the previously established effects of graded CR on glucose homeostasis in the same mice. Overall these data suggest reduced levels of adipose tissue under CR may contribute to the protective impact of CR in multiple ways linked to changes in a large population of secreted proteins.

CXXC5 suppresses hepatocellular carcinoma by promoting TGF-ß-induced cell cycle arrest and apoptosis.

Evading TGF-ß-mediated growth inhibition is often associated with tumorigenesis in liver, including hepatocellular carcinoma (HCC). To better understand the functions and the underlying molecular mechanisms of TGF-ß in HCC initiation and progression, we carried out transcriptome sequencing (RNA-Seq) to identify the target genes of TGF-ß. CXXC5, a member of the CXXC-type zinc finger domain-containing protein family, was identified as a novel TGF-ß target gene in Hep3B HCC cells. Knockdown of CXXC5 attenuated the expression of a substantial portion of TGF-ß target genes and ameliorated TGF-ß-induced growth inhibition or apoptosis of Hep3B cells, suggesting that CXXC5 is required for TGF-ß-mediated inhibition of HCC progression. Analysis of the TCGA database indicated that CXXC5 expression is reduced in the majority of HCC tissue samples in comparison to that in normal tissues. Furthermore, CXXC5 associates with the histone deacetylase HDAC1 and competes its interaction with Smad2/3, thereby abolishing the inhibitory effect of HDAC1 on TGF-ß signaling. These observations together suggest that CXXC5 may act as a tumor suppressor by promoting TGF-ß signaling via a positive feedback loop, and reveal a strategy for HCC to bypass TGF-ß-mediated cytostasis by disrupting the positive feedback regulation. Our findings shed new light on TGF-ß signaling regulation and demonstrate the function of CXXC5 in HCC development.

Inference of differentiation time for single cell transcriptomes using cell population reference data.

Single-cell RNA sequencing (scRNA-seq) is a powerful method for dissecting intercellular heterogeneity during development. Conventional trajectory analysis provides only a pseudotime of development, and often discards cell-cycle events as confounding factors. Here using matched cell population RNA-seq (cpRNA-seq) as a reference, we developed an "iCpSc" package for integrative analysis of cpRNA-seq and scRNA-seq data. By generating a computational model for reference "biological differentiation time" using cell population data and applying it to single-cell data, we unbiasedly associated cell-cycle checkpoints to the internal molecular timer of single cells. Through inferring a network flow from cpRNA-seq to scRNA-seq data, we predicted a role of M phase in controlling the speed of neural differentiation of mouse embryonic stem cells, and validated it through gene knockout (KO) experiments. By linking temporally matched cpRNA-seq and scRNA-seq data, our approach provides an effective and unbiased approach for identifying developmental trajectory and timing-related regulatory events.

Transcriptome analysis reveals determinant stages controlling human embryonic stem cell commitment to neuronal cells.

Proper neural commitment is essential for ensuring the appropriate development of the human brain and for preventing neurodevelopmental diseases such as autism spectrum disorders, schizophrenia, and intellectual disorders. However, the molecular mechanisms underlying the neural commitment in humans remain elusive. Here, we report the establishment of a neural differentiation system based on human embryonic stem cells (hESCs) and on comprehensive RNA sequencing analysis of transcriptome dynamics during early hESC differentiation. Using weighted gene co-expression network analysis, we reveal that the hESC neurodevelopmental trajectory has five stages: pluripotency (day 0); differentiation initiation (days 2, 4, and 6); neural commitment (days 8-10); neural progenitor cell proliferation (days 12, 14, and 16); and neuronal differentiation (days 18, 20, and 22). These stages were characterized by unique module genes, which may recapitulate the early human cortical development. Moreover, a comparison of our RNA-sequencing data with several other transcriptome profiling datasets from mice and humans indicated that Module 3 associated with the day 8-10 stage is a critical window of fate switch from the pluripotency to the neural lineage. Interestingly, at this stage, no key extrinsic signals were activated. In contrast, using CRISPR/Cas9-mediated gene knockouts, we also found that intrinsic hub transcription factors, including the schizophrenia-associated SIX3 gene and septo-optic dysplasia-related HESX1 gene, are required to program hESC neural determination. Our results improve the understanding of the mechanism of neural commitment in the human brain and may help elucidate the etiology of human mental disorders and advance therapies for managing these conditions.

The effects of graded levels of calorie restriction: XI. Evaluation of the main hypotheses underpinning the life extension effects of CR using the hepatic transcriptome.

Calorie restriction (CR) may extend longevity by modulating the mechanisms involved in aging. Different hypotheses have been proposed for its main mode of action. We quantified hepatic transcripts of male C57BL/6 mice exposed to graded levels of CR (0% to 40% CR) for three months, and evaluated the responses relative to these various hypotheses. Of the four main signaling pathways implied to be linked to the impact of CR on lifespan (insulin/insulin like growth factor 1 (IGF-1), nuclear factor-kappa beta (NF-ĸB), mechanistic target of rapamycin (mTOR) and sirtuins (SIRTs)), all the pathways except SIRT were altered in a manner consistent with increased lifespan. However, the expression levels of SIRT4 and SIRT7 were decreased with increasing levels of CR. Changes consistent with altered fuel utilization under CR may reduce reactive oxygen species production, which was paralleled by reduced protection. Downregulated major urinary protein (MUP) transcription suggested reduced reproductive investment. Graded CR had a positive effect on autophagy and xenobiotic metabolism, and was protective with respect to cancer signaling. CR had no significant effect on fibroblast growth factor-21 (FGF21) transcription but affected transcription in the hydrogen sulfide production pathway. Responses to CR were consistent with several different hypotheses, and the benefits of CR on lifespan likely reflect the combined impact on multiple aging related processes.

The metabolic ER stress sensor IRE1α suppresses alternative activation of macrophages and impairs energy expenditure in obesity.

Obesity is associated with metabolic inflammation and endoplasmic reticulum (ER) stress, both of which promote metabolic disease progression. Adipose tissue macrophages (ATMs) are key players orchestrating metabolic inflammation, and ER stress enhances macrophage activation. However, whether ER stress pathways underlie ATM regulation of energy homeostasis remains unclear. Here, we identified inositol-requiring enzyme 1α (IRE1α) as a critical switch governing M1-M2 macrophage polarization and energy balance. Myeloid-specific IRE1α abrogation in Ern1f/f; Lyz2-Cre mice largely reversed high-fat diet (HFD)-induced M1-M2 imbalance in white adipose tissue (WAT) and blocked HFD-induced obesity, insulin resistance, hyperlipidemia and hepatic steatosis. Brown adipose tissue (BAT) activity, WAT browning and energy expenditure were significantly higher in Ern1f/f; Lyz2-Cre mice. Furthermore, IRE1α ablation augmented M2 polarization of macrophages in a cell-autonomous manner. Thus, IRE1α senses protein unfolding and metabolic and immunological states, and consequently guides ATM polarization. The macrophage IRE1α pathway drives obesity and metabolic syndrome through impairing BAT activity and WAT browning.

LIN-28 balances longevity and germline stem cell number in Caenorhabditis elegans through let-7/AKT/DAF-16 axis.

The RNA-binding protein LIN-28 was first found to control developmental timing in Caenorhabditis elegans. Later, it was found to play important roles in pluripotency, metabolism, and cancer in mammals. Here we report that a low dosage of lin-28 enhanced stress tolerance and longevity, and reduced germline stem/progenitor cell number in C. elegans. The germline LIN-28-regulated microRNA let-7 was required for these effects by targeting akt-1/2 and decreasing their protein levels. AKT-1/2 and the downstream DAF-16 transcription factor were both required for the lifespan and germline stem cell effects of lin-28. The pathway also mediated dietary restriction induced lifespan extension and reduction in germline stem cell number. Thus, the LIN-28/let-7/AKT/DAF-16 axis we delineated here is a program that plays an important role in balancing reproduction and somatic maintenance and their response to the environmental energy level-a central dogma of the 'evolutionary optimization' of resource allocation that modulates aging.