Clinicopathological role of Cyclin A2 in uterine corpus endometrial carcinoma: Integration of tissue microarrays and ScRNA-Seq

BACKGROUND

The comprehensive expression level and potential molecular role of Cyclin A2 (CCNA2) in uterine corpus endometrial carcinoma (UCEC) remains undiscovered.

METHODS

UCEC and normal endometrium tissues from in-house and public databases were collected for investigating protein and messenger RNA expression of CCNA2. The transcription factors of CCNA2 were identified by the Cistrome database. The prognostic significance of CCNA2 in UCEC was evaluated through univariate and multivariate Cox regression as well as Kaplan-Meier curve analysis. Single-cell RNA-sequencing (scRNA-seq) analysis was performed to explore cell types in UCEC, and the AUCell algorithm was used to investigate the activity of CCNA2 in different cell types.

RESULTS

A total of 32 in-house UCEC and 30 normal endometrial tissues as well as 720 UCEC and 165 control samples from public databases were eligible and collected. Integrated calculation showed that the CCNA2 expression was up-regulated in the UCEC tissues (SMD = 2.43, 95% confidence interval 2.23∼2.64). E2F1 and FOXM1 were identified as transcription factors due to the presence of binding peaks on transcription site of CCNA2. CCNA2 predicted worse prognosis in UCEC. However, CCNA2 was not an independent prognostic factor in UCEC. The scRNA-seq analysis disclosed five cell types: B cells, T cells, monocytes, natural killer cells, and epithelial cells in UCEC. The expression of CCNA2 was mainly located in B cells and T cells. Moreover, CCNA2 was active in T cells and B cells using the AUCell algorithm.

CONCLUSION

CCNA2 was up-regulated and mainly located in T cells and B cells in UCEC. Overexpression of CCNA2 predicted unfavorable prognosis of UCEC.

Dual topologies of myotomal collagen XV and Tenascin C act in concert to guide and shape developing motor axons

During development, motor axons are guided toward muscle target by various extrinsic cues including extracellular matrix (ECM) proteins whose identities and cellular source remain poorly characterized. Here, using single-cell RNAseq of sorted GFP+ cells from smyhc1:gfp-injected zebrafish embryos, we unravel the slow muscle progenitors (SMP) pseudotemporal trajectory at the single-cell level and show that differentiating SMPs are a major source of ECM proteins. The SMP core-matrisome was characterized and computationally predicted to form a basement membrane-like structure tailored for motor axon guidance, including basement membrane-associated ECM proteins, as collagen XV-B, one of the earliest core-matrisome gene transcribed in differentiating SMPs and the glycoprotein Tenascin C. To investigate how contact-mediated guidance cues are organized along the motor path to exert their function in vivo, we used microscopy-based methods to analyze and quantify motor axon navigation in tnc and col15a1b knock-out fish. We show that motor axon shape and growth rely on the timely expression of the attractive cue Collagen XV-B that locally provides axons with a permissive soft microenvironment and separately organizes the repulsive cue Tenascin C into a unique functional dual topology. Importantly, bioprinted micropatterns that mimic this in vivo ECM topology were sufficient to drive directional motor axon growth. Our study offers evidence that not only the composition of ECM cues but their topology critically influences motor axon navigation in vertebrates with potential applications in regenerative medicine for peripheral nerve injury as regenerating nerves follow their original path.

Lineage tracing of T cell differentiation from T-iPSC by 2D feeder-free culture and 3D organoid culture

Introduction

T cells induced from induced pluripotent stem cells(iPSCs) derived from antigen-specific T cells (T-iPS-T cells) are an attractive tool for T cell immunotherapy. The induction of cytotoxic T-iPS-T cells is well established in feeder-free condition for the aim of off-the-shelf production, however, the induction of helper T-iPS-T cells remains challenging.

Methods

We analyzed T-iPS-T cells matured in 3D organoid culture at different steps in the culture process at the single-cell level. T-iPS-T cell datasets were merged with an available human thymocyte dataset based in single-cell RNA sequencing (scRNA-seq). Particularly, we searched for genes crucial for generation CD4+ T-iPS-T cells by comparing T-iPS-T cells established in 2D feeder-free or 3D organoid culture.

Results

The scRNA-seq data indicated that T-iPS-T cells are similar to T cells transitioning to human thymocytes, with SELENOW, GIMAP4, 7, SATB1, SALMF1, IL7R, SYTL2, S100A11, STAT1, IFITM1, LZTFL1 and SOX4 identified as candidate genes for the 2D feeder-free induction of CD4+ T-iPS-T cells.

Discussion

This study provides single cell transcriptome datasets of iPS-T cells and leads to further analysis for CD4+ T cell generation from T-iPSCs.

Single-cell RNA seq analysis of erythroid cells reveals a specific sub-population of stress erythroid progenitors

BACKGROUND

: Erythroid cells play important roles in hemostasis and disease. However, there is still significant knowledge gap regarding stress erythropoiesis.

METHODS

: Two single-cell RNAseq datasets of erythroid cells on GEO with accession numbers GSE149938 and GSE184916 were obtained. The datasets from two sources, bone marrow and peripheral blood were analyzed using Seurat v4.1.1, and other tools in R. QC metrics were performed, data were normalized and scaled. Principal components that capture the variation of the data were determined. In clustering the cells, KNN graph was constructed and Louvain algorithm was applied to optimize the standard modularity function. Clusters were defined via differential expression of features.

RESULTS

We identified 9 different cell types, with a particular cluster representing the stress erythroids. The clusters showed differentially expressed genes as observed from the gene signature plot. The stress erythroid cluster differentially expressed some genes including ALAS2, HEMGN, and GUK1.

CONCLUSION

The erythroid population was found to be heterogeneous, with a distinct sub-cell type constituting the stress erythroids; this may have important implications for our knowledge of steady-state and stress erythropoiesis, and the markers found in this cluster may prove useful for future research into the dynamics of stress erythroid progenitor cell differentiation.

The Simple prEservatioN of Single cElls method for cryopreservation enables the generation of single-cell immune profiles from whole blood

Introduction

Current multistep methods utilized for preparing and cryopreserving single-cell suspensions from blood samples for single-cell RNA sequencing (scRNA-seq) are time-consuming, requiring trained personnel and special equipment, so limiting their clinical adoption. We developed a method, Simple prEservatioN of Single cElls (SENSE), for single-step cryopreservation of whole blood (WB) along with granulocyte depletion during single-cell assay, to generate high quality single-cell profiles (SCP).

Methods

WB was cryopreserved using the SENSE method and peripheral blood mononuclear cells (PBMCs) were isolated and cryopreserved using the traditional density-gradient method (PBMC method) from the same blood sample (n=6). The SCPs obtained from both methods were processed using a similar pipeline and quality control parameters. Further, entropy calculation, differential gene expression, and cellular communication analysis were performed to compare cell types and subtypes from both methods.

Results

Highly viable (86.3 ± 1.51%) single-cell suspensions (22,353 cells) were obtained from the six WB samples cryopreserved using the SENSE method. In-depth characterization of the scRNA-seq datasets from the samples processed with the SENSE method yielded high-quality profiles of lymphoid and myeloid cell types which were in concordance with the profiles obtained with classical multistep PBMC method processed samples. Additionally, the SENSE method cryopreserved samples exhibited significantly higher T-cell enrichment, enabling deeper characterization of T-cell subtypes. Overall, the SENSE and PBMC methods processed samples exhibited transcriptional, and cellular communication network level similarities across cell types with no batch effect except in myeloid lineage cells.

Discussion

Comparative analysis of scRNA-seq datasets obtained with the two cryopreservation methods i.e., SENSE and PBMC methods, yielded similar cellular and molecular profiles, confirming the suitability of the former method's incorporation in clinics/labs for cryopreserving and obtaining high-quality single-cells for conducting critical translational research.

Acute irradiation causes a long-term disturbance in the heterogeneity and gene expression profile of medullary thymic epithelial cells

The thymus has the ability to regenerate from acute injury caused by radiation, infection, and stressors. In addition to thymocytes, thymic epithelial cells in the medulla (mTECs), which are crucial for T cell self-tolerance by ectopically expressing and presenting thousands of tissue-specific antigens (TSAs), are damaged by these insults and recover thereafter. However, given recent discoveries on the high heterogeneity of mTECs, it remains to be determined whether the frequency and properties of mTEC subsets are restored during thymic recovery from radiation damage. Here we demonstrate that acute total body irradiation with a sublethal dose induces aftereffects on heterogeneity and gene expression of mTECs. Single-cell RNA-sequencing (scRNA-seq) analysis showed that irradiation reduces the frequency of mTECs expressing AIRE, which is a critical regulator of TSA expression, 15 days after irradiation. In contrast, transit-amplifying mTECs (TA-mTECs), which are progenitors of AIRE-expressing mTECs, and Ccl21a-expressing mTECs, were less affected. Interestingly, a detailed analysis of scRNA-seq data suggested that the proportion of a unique mTEC cluster expressing Ccl25 and a high level of TSAs was severely decreased by irradiation. In sum, we propose that the effects of acute irradiation disrupt the heterogeneity and properties of mTECs over an extended period, which potentially leads to an impairment of thymic T cell selection.

Using combined single-cell gene expression, TCR sequencing and cell surface protein barcoding to characterize and track CD4+ T cell clones from murine tissues

Single-cell gene expression analysis using sequencing (scRNA-seq) has gained increased attention in the past decades for studying cellular transcriptional programs and their heterogeneity in an unbiased manner, and novel protocols allow the simultaneous measurement of gene expression, T-cell receptor clonality and cell surface protein expression. In this article, we describe the methods to isolate scRNA/TCR-seq-compatible CD4+ T cells from murine tissues, such as skin, spleen, and lymph nodes. We describe the processing of cells and quality control parameters during library preparation, protocols for multiplexing of samples, and strategies for sequencing. Moreover, we describe a step-by-step bioinformatic analysis pipeline from sequencing data generated using these protocols. This includes quality control, preprocessing of sequencing data and demultiplexing of individual samples. We perform quantification of gene expression and extraction of T-cell receptor alpha and beta chain sequences, followed by quality control and doublet detection, and methods for harmonization and integration of datasets. Next, we describe the identification of highly variable genes and dimensionality reduction, clustering and pseudotemporal ordering of data, and we demonstrate how to visualize the results with interactive and reproducible dashboards. We will combine different analytic R-based frameworks such as Bioconductor and Seurat, illustrating how these can be interoperable to optimally analyze scRNA/TCR-seq data of CD4+ T cells from murine tissues.

Orthogonal outlier detection and dimension estimation for improved MDS embedding of biological datasets

Conventional dimensionality reduction methods like Multidimensional Scaling (MDS) are sensitive to the presence of orthogonal outliers, leading to significant defects in the embedding. We introduce a robust MDS method, called DeCOr-MDS (Detection and Correction of Orthogonal outliers using MDS), based on the geometry and statistics of simplices formed by data points, that allows to detect orthogonal outliers and subsequently reduce dimensionality. We validate our methods using synthetic datasets, and further show how it can be applied to a variety of large real biological datasets, including cancer image cell data, human microbiome project data and single cell RNA sequencing data, to address the task of data cleaning and visualization.

150 risk variants for diverticular disease of intestine prioritize cell types and enable polygenic prediction of disease susceptibility

We conducted a genome-wide association study (GWAS) analysis of diverticular disease (DivD) of intestine within 724,372 individuals and identified 150 independent genome-wide significant DNA variants. Integration of the GWAS results with human gut single-cell RNA sequencing data implicated gut myocyte, mesothelial and stromal cells, and enteric neurons and glia in DivD development. Ninety-five genes were prioritized based on multiple lines of evidence, including SLC9A3, a drug target gene of tenapanor used for the treatment of the constipation subtype of irritable bowel syndrome. A DivD polygenic score (PGS) enables effective risk prediction (area under the curve [AUC], 0.688; 95% confidence interval [CI], 0.645-0.732) and the top 20% PGS was associated with ∼3.6-fold increased DivD risk relative to the remaining population. Our statistical and bioinformatic analyses suggest that the mechanism of DivD is through colon structure, gut motility, gastrointestinal mucus, and ionic homeostasis. Our analyses reinforce the link between gastrointestinal disorders and the enteric nervous system through genetics.

Shared Differential Expression-Based Distance Reflects Global Cell Type Relationships in Single-Cell RNA Sequencing Data

Unsupervised cell clustering on the basis of meaningful biological variation in single-cell RNA sequencing (scRNA seq) data has received significant attention, as it assists with ontological subpopulation identification among the data. A key step in the clustering process is to compute distances between the cells under a specified distance measure. Although particular distance measures may successfully separate cells into biologically relevant clusters, they may fail to retain global structure of the data, such as relative similarity between the cell clusters. In this article, we modify a biologically motivated distance measure, SIDEseq, for use of aggregate comparisons of cell types in large single-cell assays, and demonstrate that, across simulated and real scRNA seq data, the distance matrix more consistently retains global cell type relationships than commonly used distance measures for scRNA seq clustering. We call the modified distance measure "SIDEREF." We explore spectral dimension reduction of the SIDEREF distance matrix as a means of noise filtering, similar to principal components analysis applied directly to expression data. We utilize a summary measure of relative cell type distances to better display the cell group relationships. SIDEREF visualizations more consistently reflect global structures in the data than other commonly considered distance measures. We utilize relative cell type distances and the SIDEREF distance measure to uncover compositional differences between annotated leukocyte cell groups in a compendium of Mus musculus scRNA seq assays comprising 12 tissues. SIDEREF and associated analysis is openly available on GitHub.

Integration of scRNA-Seq and TCGA RNA-Seq to Analyze the Heterogeneity of HPV+ and HPV- Cervical Cancer Immune Cells and Establish Molecular Risk Models

Background

Numerous studies support that Human papillomavirus (HPV) can cause cervical cancer. However, few studies have surveyed the heterogeneity of HPV infected or uninfected (HPV+ and HPV-) cervical cancer (CESC) patients. Integration of scRNA-seq and TCGA data to analyze the heterogeneity of HPV+ and HPV- cervical cancer patients on a single-cell level could improve understanding of the cellular mechanisms during HPV-induced cervical cancer.

Methods

CESC scRNA-seq data obtained from the Gene Expression Omnibus (GEO) database and the Seurat, Monocle3 package were used for scRNA-seq data analysis. The ESTIMATE package was used for single-sample gene immune score, CIBERSORT package was used to identify immune scores of cells, and the “WGCNA” package for the weighted correlation network analysis. Univariate Cox and LASSO regression were performed to establish survival and relapse signatures. KEGG and GO analyses were performed for the signature gene. Gene Expression Profiling Interactive Analysis was used for Pan-cancer analysis.

Results

In the HPV+ CESC group, CD8+ T cells and B cells were down-regulated, whereas T reg cells, CD4+ T cells, and epithelial cells were up-regulated according to scRNA-seq data. Survival analysis of TCGA-CESC revealed that increased expression of naive B cells or CD8+ T cells favors the survival probability of CESC patients. WGCNA, univariate Cox, and LASSO Cox regression established a 9-genes survival signature and a 7-gene relapse model. Pan-cancer analysis identified IKZF3, FOXP3, and JAK3 had a similar distribution and effects in HPV-associated HNSC.

Conclusion

Analysis of scRNA-seq and bulk RNA-seq of HPV+ and HPV- CESC samples revealed heterogeneity from transcriptional state to immune infiltration. Survival and relapse models were adjusted according to the heterogeneity of HPV+ and HPV- CESC immune cells to assess the prognostic risk accurately. Hub genes represent similar protection in HPV- associated HNSC while showing irrelevant to other potential HPV-related cancers.

AhR Activation Leads to Attenuation of Murine Autoimmune Hepatitis: Single-Cell RNA-Seq Analysis Reveals Unique Immune Cell Phenotypes and Gene Expression Changes in the Liver

The aryl hydrocarbon receptor (AhR) is a ubiquitously expressed ligand-activated transcription factor. While initially identified as an environmental sensor, this receptor has been shown more recently to regulate a variety of immune functions. AhR ligands vary in structure and source from environmental chemicals such as 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) and indoles found in cruciferous vegetables to endogenous ligands derived from tryptophan metabolism. In the current study, we used TCDD, a high affinity AhR ligand to study the impact of AhR activation in the murine model of autoimmune hepatitis (AIH). Primarily, we used single-cell RNA-sequencing (scRNA-seq) technology to study the nature of changes occurring in the immune cells in the liver at the cellular and molecular level. We found that AhR activation attenuated concanavalin A (ConA)-induced AIH by limiting chemotaxis of pro-inflammatory immune cell subsets, promoting anti-inflammatory cytokine production, and suppressing pro-inflammatory cytokine production. scRNA-seq analysis showed some unusual events upon ConA injection such as increased presence of mature B cells, natural killer (NK) T cells, CD4+ or CD8+ T cells, Kupffer cells, memory CD8+ T cells, and activated T cells while TCDD treatment led to the reversal of most of these events. Additionally, the immune cells showed significant alterations in the gene expression profiles. Specifically, we observed downregulation of inflammation-associated genes including Ptma, Hspe1, and CD52 in TCDD-treated AIH mice as well as alterations in the expression of migratory markers such as CXCR2. Together, the current study characterizes the nature of inflammatory changes occurring in the liver during AIH, and sheds light on how AhR activation during AIH attenuates liver inflammation by inducing phenotypic and genotypic changes in immune cells found in the liver.

Postnatal Smad3 Inactivation in Murine Smooth Muscle Cells Elicits a Temporally and Regionally Distinct Transcriptional Response

Heterozygous, loss of function mutations in positive regulators of the Transforming Growth Factor-β (TGF-β) pathway cause hereditary forms of thoracic aortic aneurysm. It is unclear whether and how the initial signaling deficiency triggers secondary signaling upregulation in the remaining functional branches of the pathway, and if this contributes to maladaptive vascular remodeling. To examine this process in a mouse model in which time-controlled, partial interference with postnatal TGF-β signaling in vascular smooth muscle cells (VSMCs) could be assessed, we used a VSMC-specific tamoxifen-inducible system, and a conditional allele, to inactivate Smad3 at 6 weeks of age, after completion of perinatal aortic development. This intervention induced dilation and histological abnormalities in the aortic root, with minor involvement of the ascending aorta. To analyze early and late events associated with disease progression, we performed a comparative single cell transcriptomic analysis at 10- and 18-weeks post-deletion, when aortic dilation is undetectable and moderate, respectively. At the early time-point, Smad3-inactivation resulted in a broad reduction in the expression of extracellular matrix components and critical components of focal adhesions, including integrins and anchoring proteins, which was reflected histologically by loss of connections between VSMCs and elastic lamellae. At the later time point, however, expression of several transcripts belonging to the same functional categories was normalized or even upregulated; this occurred in association with upregulation of transcripts coding for TGF-β ligands, and persistent downregulation of negative regulators of the pathway. To interrogate how VSMC heterogeneity may influence this transition, we examined transcriptional changes in each of the four VSMC subclusters identified, regardless of genotype, as partly reflecting the proximal-to-distal anatomic location based on in situ RNA hybridization. The response to Smad3-deficiency varied depending on subset, and VSMC subsets over-represented in the aortic root, the site most vulnerable to dilation, most prominently upregulated TGF-β ligands and pro-pathogenic factors such as thrombospondin-1, angiotensin converting enzyme, and pro-inflammatory mediators. These data suggest that Smad3 is required for maintenance of focal adhesions, and that loss of contacts with the extracellular matrix has consequences specific to each VSMC subset, possibly contributing to the regional susceptibility to dilation in the aorta.

Loss of Function of the Neural Cell Adhesion Molecule NrCAM Regulates Differentiation, Proliferation and Neurogenesis in Early Postnatal Hypothalamic Tanycytes

Hypothalamic tanycytes are neural stem and progenitor cells, but little is known of how they are regulated. Here we provide evidence that the cell adhesion molecule, NrCAM, regulates tanycytes in the adult niche. NrCAM is strongly expressed in adult mouse tanycytes. Immunohistochemical and in situ hybridization analysis revealed that NrCAM loss of function leads to both a reduced number of tanycytes and reduced expression of tanycyte-specific cell markers, along with a small reduction in tyrosine hydroxylase-positive arcuate neurons. Similar analyses of NrCAM mutants at E16 identify few changes in gene expression or cell composition, indicating that NrCAM regulates tanycytes, rather than early embryonic hypothalamic development. Neurosphere and organotypic assays support the idea that NrCAM governs cellular homeostasis. Single-cell RNA sequencing (scRNA-Seq) shows that tanycyte-specific genes, including a number that are implicated in thyroid hormone metabolism, show reduced expression in the mutant mouse. However, the mild tanycyte depletion and loss of markers observed in NrCAM-deficient mice were associated with only a subtle metabolic phenotype.

Ductal Macrophages Predominate in the Immune Landscape of the Lactating Mammary Gland

The mammary gland is unique in female mammals. Mammary tissue undergoes development and remodeling during lactation, a stage associated with high susceptibility to bacterial infections, inducing an inflammatory condition called mastitis. Although the immune response of the mammary gland has been the subject of intense research to improve prevention and treatment efficacy, the precise definition of its immune composition at this particular physiological stage is still missing. We combined single-cell RNA-Seq, flow cytometry, and three-dimensional confocal microscopy techniques to characterize the immune landscape of lactating murine mammary tissue. Macrophages dominated the immune cell repertoire and could be subdivided into at least two subsets: ductal and stromal macrophages. Ductal macrophages represented approximately 80% of the total CD45^pos immune cells and co-expressed F4/80 and CD11c, with high levels of MHC class II molecules. They were strategically poised below the alveolar basal cells in contact with the myoepithelial cell network. Adaptive T and B lymphocytes were remarkably less numerous at this stage, which could explain the limited efficacy of vaccination against mastitis. These results support the view that new strategies to increase mammary immunity and prevent mastitis should be devised.

Transcriptomic Signatures of Airway Epithelium Infected With SARS-CoV-2: A Balance Between Anti-infection and Virus Load

COVID-19 pneumonia requires effective medical therapies. However, it is a challenge to find therapeutic drugs that not only inhibit viral replication, but also inhibit the accompanying cytokine storm and maintain an appropriate immune response. In this study, the effects of SARS-CoV-2 on gene expression in lung epithelial cells from patients with COVID-19 were systematically evaluated with bioinformatics analysis methods. Transcriptome expression specific to bystander (exposed but uninfected) and infected cells were found, and the vital pathways were identified by conducting differentially expressed gene analysis regarding the relationship between gene signatures of COVID-19 infection and disease severity. We found that a high viral load did not necessarily imply a low response of epithelial cells or a poor disease convalescence. The ability to distinguish the role of virus-correlated genes facilitates the development of potential new medicines and therapies for COVID-19 infection.

Single-Cell Transcriptome Profiling Reveals the Suppressive Role of Retinal Neurons in Microglia Activation Under Diabetes Mellitus

Diabetic retinopathy, as one of the common complications of diabetes mellitus, is the leading cause of blindness in the working-age population worldwide. The disease is characterized by damage to retinal vasculature, which is associated with the activation of retina microglial and induces chronic neurodegeneration. Previous studies have identified the effects of activated microglial on the retinal neurons, but the cellular and molecular mechanisms underlying microglial activation is largely unknown. Here, we performed scRNA-seq on the retina of non-human primates with diabetes mellitus, and identified cell-type-specific molecular changes of the six major cell types. By identifying the ligand-receptor expression patterns among different cells, we established the interactome of the whole retina. The data showed that TNF-α signal mediated the activation of microglia through an autocrine manner. And we found TGFβ2, which was upregulated in cone dramatically by hyperglycemia, inhibited microglia activation at the early stage of diabetic retinopathy. In summary, our study is the first to profile cell-specific molecular changes and the cell-cell interactome of retina under diabetes mellitus, paving a way to dissect the cellular and molecular mechanisms underlying early-stage diabetic retinopathy.

Transcriptional, Epigenetic, and Functional Reprogramming of Monocytes From Non-Human Primates Following Chronic Alcohol Drinking

Chronic heavy drinking (CHD) of alcohol is a known risk factor for increased susceptibility to bacterial and viral infection as well as impaired wound healing. Evidence suggests that these defects are mediated by a dysregulated inflammatory response originating from myeloid cells, notably monocytes and macrophages, but the mechanisms remain poorly understood. Our ability to study CHD is impacted by the complexities of human drinking patterns and behavior as well as comorbidities and confounding risk factors for patients with alcohol use disorders. To overcome these challenges, we utilized a translational rhesus macaque model of voluntary ethanol self-administration that closely recapitulates human drinking patterns and chronicity. In this study, we examined the effects of CHD on blood monocytes in control and CHD female macaques after 12 months of daily ethanol consumption. While monocytes from CHD female macaques generated a hyper-inflammatory response to ex vivo LPS stimulation, their response to E. coli was dampened. In depth scRNA-Seq analysis of purified monocytes revealed significant shifts in classical monocyte subsets with accumulation of cells expressing markers of hypoxia (HIF1A) and inflammation (NFkB signaling pathway) in CHD macaques. The increased presence of monocyte subsets skewed towards inflammatory phenotypes was complemented by epigenetic analysis, which revealed higher accessibility of promoter regions that regulate genes involved in cytokine signaling pathways. Collectively, data presented in this manuscript demonstrate that CHD shifts classical monocyte subset composition and primes the monocytes towards a more hyper-inflammatory response to LPS, but compromised pathogen response.

S100A9 Derived From Myeloma Associated Myeloid Cells Promotes TNFSF13B/TNFRSF13B-Dependent Proliferation and Survival of Myeloma Cells

Multiple myeloma (MM) is a lethal hematological malignancy characterized by abundant myeloid cells in the microenvironment that fuel tumor progression. But the mechanism by which myeloid cells support myeloma cells has not been fully explored. We aimed to examine their effect on bone marrow cells of MM patients by scRNA-seq transcriptome analysis and reveal a high-resolution gene profile of myeloma cells and myeloma-associated myeloid cells. Based on correlation analysis of integrated scRNA-seq and bulk RNA-seq datasets from patients, we confirmed that myeloid-derived S100A9 was involved in TNFSF13B-dependent myeloma cell proliferation and survival. In the animal experiments, S100A9 was found to be critical for MM cell proliferation and survival via TNFSF13B production by myeloid cells, neutrophils, and macrophages. In-vitro analysis of patient primary myeloma cells further demonstrated that enhanced TNFSF13B signaling triggered the canonical NF-κB pathway to boost tumor cell proliferation. All these results suggest that myeloid-derived S100A9 is required for TNFSF13B/TNFRSF13B-dependent cell-fate specification, which provides fresh insights into MM progression.

Novel Targets of SARS-CoV-2 Spike Protein in Human Fetal Brain Development Suggest Early Pregnancy Vulnerability

Pregnant women are at greater risk of infection by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), because of their altered immunity and strained cardiovascular system. Emerging studies of placenta, embryos, and cerebral organoids suggest that fetal organs including brain could also be vulnerable to coronavirus disease 2019 (COVID-19). Additionally, a case study from Paris has reported transient neurological complications in neonates born to pregnant mothers. However, it remains poorly understood whether the fetal brain expresses cellular components that interact with Spike protein (S) of coronaviruses, which facilitates fusion of virus and host cell membrane and is the primary protein in viral entry. To address this question, we analyzed the expression of known (ACE2, TMPRSS2, and FURIN) and novel (ZDHHC5, GOLGA7, and ATP1A1) S protein interactors in publicly available fetal brain bulk and single cell RNA sequencing datasets. Bulk RNA sequencing analysis across multiple regions of fetal brain spanning 8 weeks post conception (wpc)-37wpc indicates that two of the known S protein interactors are expressed at low levels with median normalized gene expression values ranging from 0.08 to 0.06 (ACE2) and 0.01-0.02 (TMPRSS2). However, the third known S protein interactor FURIN is highly expressed (11.1-44.09) in fetal brain. Interestingly, all three novel S protein interactors are abundantly expressed throughout fetal brain development with median normalized gene expression values ranging from 20.38-21.60 (ZDHHC5), 92.47-68.35 (GOLGA7), and 65.45-194.5 (ATP1A1). Moreover, the peaks of expression of novel interactors is around 12-26wpc. Using publicly available single cell RNA sequencing datasets, we further show that novel S protein interactors show higher co-expression with neurons than with neural progenitors and astrocytes. These results suggest that even though two of the known S protein interactors are present at low levels in fetal brain, novel S protein interactors are abundantly present and could play a direct or indirect role in SARS-CoV-2 fetal brain pathogenesis, especially during the 2nd and 3rd trimesters of pregnancy.