Distinct fibroblast subpopulations associated with bone, brain or intrapulmonary metastasis in advanced non-small-cell lung cancer.

BACKGROUND: Bone or brain metastases may develop in 20-40% of individuals with late-stage non-small-cell lung cancer (NSCLC), resulting in a median overall survival of only 4-6 months. However, the primary lung cancer tissue's distinctions between bone, brain and intrapulmonary metastases of NSCLC at the single-cell level have not been underexplored. METHODS: We conducted a comprehensive analysis of 14 tissue biopsy samples obtained from treatment-naïve advanced NSCLC patients with bone (n = 4), brain (n = 6) or intrapulmonary (n = 4) metastasis using single-cell sequencing originating from the lungs. Following quality control and the removal of doublets, a total of 80 084 cells were successfully captured. RESULTS: The most significant inter-group differences were observed in the fraction and function of fibroblasts. We identified three distinct cancer-associated fibroblast (CAF) subpopulations: myofibroblastic CAF (myCAF), inflammatory CAF (iCAF) and antigen-presenting CAF (apCAF). Notably, apCAF was prevalent in NSCLC with bone metastasis, while iCAF dominated in NSCLC with brain metastasis. Intercellular signalling network analysis revealed that apCAF may play a role in bone metastasis by activating signalling pathways associated with cancer stemness, such as SPP1-CD44 and SPP1-PTGER4. Conversely, iCAF was found to promote brain metastasis by activating invasion and metastasis-related molecules, such as MET hepatocyte growth factor. Furthermore, the interaction between CAFs and tumour cells influenced T-cell exhaustion and signalling pathways within the tumour microenvironment. CONCLUSIONS: This study unveils the direct interplay between tumour cells and CAFs in NSCLC with bone or brain metastasis and identifies potential therapeutic targets for inhibiting metastasis by disrupting these critical cell-cell interactions.

eccDNA-pipe: an integrated pipeline for identification, analysis and visualization of extrachromosomal circular DNA from high-throughput sequencing data.

Extrachromosomal circular DNA (eccDNA) is currently attracting considerable attention from researchers due to its significant impact on tumor biogenesis. High-throughput sequencing (HTS) methods for eccDNA identification are continually evolving. However, an efficient pipeline for the integrative and comprehensive analysis of eccDNA obtained from HTS data is still lacking. Here, we introduce eccDNA-pipe, an accessible software package that offers a user-friendly pipeline for conducting eccDNA analysis starting from raw sequencing data. This dataset includes data from various sequencing techniques such as whole-genome sequencing (WGS), Circle-seq and Circulome-seq, obtained through short-read sequencing or long-read sequencing. eccDNA-pipe presents a comprehensive solution for both upstream and downstream analysis, encompassing quality control and eccDNA identification in upstream analysis and downstream tasks such as eccDNA length distribution analysis, differential analysis of genes enriched with eccDNA and visualization of eccDNA structures. Notably, eccDNA-pipe automatically generates high-quality publication-ready plots. In summary, eccDNA-pipe provides a comprehensive and user-friendly pipeline for customized analysis of eccDNA research.

RNAkines are secreted messengers shaping health and disease.

Extracellular noncoding RNAs (ncRNAs) have crucial roles in intercellular communications. The process of ncRNA secretion is highly regulated, with specific ncRNA profiles produced under different physiological and pathological circumstances. These ncRNAs are transported primarily via extracellular vesicles (EVs) from their origin cells to target cells, utilising both endocrine and paracrine pathways. The intercellular impacts of extracellular ncRNAs are essential for maintaining homeostasis and the pathogenesis of various diseases. Given the unique aspects of extracellular ncRNAs, here we propose the term 'RNAkine' to describe these recently identified secreted factors. We explore their roles as intercellular modulators, particularly in their ability to regulate metabolism and influence tumorigenesis, highlighting their definition and importance as a distinct class of secreted factors.

A single-cell landscape of pre- and post-menopausal high-grade serous ovarian cancer ascites.

High-grade serous ovarian cancer (HGSOC) is a hormone-related cancer with high mortality and poor prognosis. Based on the transcriptome of 57,444 cells in ascites from 10 patients with HGSOC (including 5 pre-menopausal and 5 post-menopausal patients), we identified 14 cell clusters which were further classified into 6 cell types, including T cells, B cells, NK cells, myeloid cells, epithelial cells, and stromal cells. We discovered an increased proportion of epithelial cells and a decreased proportion of T cells in pre-menopausal ascites compared with post-menopausal ascites. GO analysis revealed the pre-menopausal tumor microenvironments (TME) are closely associated with viral infection, while the post-menopausal TME are mostly related to the IL-17 immune pathway. SPP1/CD44-mediated crosstalk between myeloid cells and B cells, NK cells, and stromal cells mainly present in the pre-menopausal group, while SPP1/PTGER4 -mediated crosstalk between myeloid cells and epithelial cells mostly present in the post-menopausal group.

The kava chalcone flavokawain B exerts inhibitory activity and synergizes with BCL-2 inhibition in malignant B-cell lymphoma.

BACKGROUND: B-cell lymphoma, which originates from B cells at diverse differentiation stages, is the most common non-Hodgkin lymphoma with tremendous treatment challenges and unsatisfactory clinical outcomes. Flavokawain B (FKB), a naturally occurring chalcone extracted from kava, possesses promising anticancer properties. However, evidence on the effects of FKB on hematological malignancies, particularly lymphomas, remains scarce. PURPOSE: This study aimed to investigate the antilymphoma effect of FKB and its underlying mechanisms. STUDY DESIGN/METHODS: Proliferation assays, flow cytometry, and western blotting were employed to determine whether and how FKB affected B-cell lymphoma cell lines in vitro. Xenograft mouse models were established to evaluate the antilymphoma efficacy of FKB in vivo. RESULTS: FKB reduced the viability of a panel of B-cell lymphoma cell lines in a dose- and time-dependent manner. Mitochondrial apoptosis was markedly induced by FKB, as evidenced by an increased percentage of annexin V-positive cells, a loss of mitochondrial membrane potential, and cleavage of caspase-3 and PARP. Moreover, FKB inhibited BCL-XL expression and synergized with the BCL-2 inhibitor ABT-199. Mechanistically, FKB treatment decreased the phosphorylation of Akt, mammalian target of rapamycin (mTOR), glycogen synthase kinase-3ß (GSK3ß), and ribosomal protein S6 (RPS6). Pharmacological blockage of phosphoinositide 3-kinase (PI3K), Akt, or GSK3ß potentiated the activity of FKB, indicating the involvement of the PI3K/Akt cascade in FKB-mediated inhibitory effects. In mouse xenograft models, the intraperitoneal administration of FKB significantly decreased lymphoma growth, accompanied by diminished mitosis and Ki-67 staining of tumor tissues. CONCLUSION: Our data demonstrate the robust therapeutic potential of FKB in the treatment of B-cell lymphoma.

Dynamics and regulation of mitotic chromatin accessibility bookmarking at single-cell resolution.

Although mitotic chromosomes are highly compacted and transcriptionally inert, some active chromatin features are retained during mitosis to ensure the proper postmitotic reestablishment of maternal transcriptional programs, a phenomenon termed "mitotic bookmarking." However, the dynamics and regulation of mitotic bookmarking have not been systemically surveyed. Using single-cell transposase-accessible chromatin sequencing (scATAC-seq), we examined 6538 mitotic L02 human liver cells of variable stages and found that chromatin accessibility remained changing throughout cell division, with a constant decrease until metaphase and a gradual increase as chromosomes segregated. In particular, a subset of chromatin regions were identified to remain open throughout mitosis, and genes associated with these bookmarked regions are primarily linked to rapid reactivation upon mitotic exit. We also demonstrated that nuclear transcription factor Y subunit α (NF-YA) preferentially occupied bookmarked regions and contributed to transcriptional reactivation after mitosis. Our study uncovers the dynamic and regulatory blueprint of mitotic bookmarking.

The gut microbiome promotes arsenic metabolism and alleviates the metabolic disorder for their mammal host under arsenic exposure.

Gut microbiome can participate in arsenic metabolism. However, its efficacy in the host under arsenic stress is still controversial. To clarify their roles in fecal arsenic excretion, tissue arsenic accumulation, host physiological states and metabolism, in this study, ninety-six C57BL/6 male mice were randomly divided to four groups, groups A and B were given sterile water, and groups C and D were given the third generation of broad-spectrum antibiotic (ceftriaxone) to erase the background gut microbiome. Subsequently, groups B and D were subchronicly exposed to arsenic containing feed prepared by adding arsenical mixture (rice arsenic composition) into control feed. In group D, the fecal total arsenic (CtAs) decreased by 25.5 %, iAsIII composition increased by 46.9 %, unclarified As (uAs) composition decreased by 92.4 %, and the liver CtAs increased by 26.7 %; the fecal CtAs was positively correlated with microbial richness and some metabolites (organic acids, amino acids, carbohydrates, SCFAs, hydrophilic bile acids and their derivatives); and fecal DMA was positively correlated with microbial richness and some metabolites (ferulic acid, benzenepropanoic acid and pentanoic acid); network analysis showed that the numbers of modules, nodes, links were decreased and vulnerability was increased; some SCFAs and hydrophilic bile acid decreased, and hydrophobic bile acids increased (Ps < 0.05). In the tissue samples of group D, Il-18 and Ifn-γ gene expression increased and intestinal barrier-related genes Muc2, Occludin and Zo-1 expression decreased (Ps < 0.05); serum glutathione and urine malondialdehyde significantly increased (Ps < 0.05); urine metabolome significantly changed and the variation was correlated with six SCFAs-producing bacteria, and some SCFAs including isobutyric acid, valeric acid and heptanoic acid decreased (Ps < 0.05). Therefore, the normal gut microbiome increases fecal arsenic excretion and biotransformation, which can maintain a healthier microbiome and metabolic functions, and alleviate the metabolic disorder for their mammal host under arsenic exposure.

Single-cell transcriptome profiling and chromatin accessibility reveal an exhausted regulatory CD4+ T cell subset in systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease, and CD4+ T cells are known to promote SLE development. Here, we explore heterogeneities in the CD4+ T cell regulome and their associations with SLE pathogenesis by performing assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) and single-cell transcriptome sequencing (single-cell RNA sequencing [scRNA-seq]) of peripheral CD4+ T cells from 72 SLE patients and 30 healthy controls. Chromatin accessibility signatures of CD4+ T cells are correlated with disease severity. Further, we generate 34,176 single-cell transcriptomes of healthy and SLE CD4+ T cells and reveal transcriptional dysfunction of regulatory T (Treg) cells, identifying two Treg subpopulations, among which the CCR7lowCD74hi Treg subgroup features type I interferon-induced functional exhaustion in SLE patients. These transcriptome-level findings for SLE Tregs are mirrored in trends from the ATAC-seq data. Our study establishes a rich empirical foundation for understanding SLE and uncovers previously unknown contributions of Treg with exhaustion-like properties to SLE pathogenesis.

Single-cell transcriptomics reveal a unique memory-like NK cell subset that accumulates with ageing and correlates with disease severity in COVID-19.

BACKGROUND: Natural killer (NK) cells are innate lymphoid cells that mediate antitumour and antiviral responses. However, very little is known about how ageing influences human NK cells, especially at the single-cell level. METHODS: We applied single-cell sequencing (scRNA-seq) to human lymphocytes and NK cells from 4 young and 4 elderly individuals and then analysed the transcriptome data using Seurat. We detected the proportion and phenotype of NK cell subsets in peripheral blood samples from a total of 62 young and 52 elderly healthy donors by flow cytometry. We also used flow cytometry to examine the effector functions of NK cell subsets upon IFN-α/IL-12+IL-15/K562/IL-2 stimulation in vitro in peripheral blood samples from a total of 64 young and 63 elderly healthy donors. We finally studied and integrated single-cell transcriptomes of NK cells from 15 young and 41 elderly COVID-19 patients with those from 12 young and 6 elderly healthy control individuals to investigate the impacts of ageing on NK cell subsets in COVID-19 disease. RESULTS: We discovered a memory-like NK subpopulation (NK2) exhibiting the largest distribution change between elderly and young individuals among lymphocytes. Notably, we discovered a unique NK subset that was predominantly CD52+ NK2 cells (NK2.1). These memory-like NK2.1 cells accumulated with age, exhibited proinflammatory characteristics, and displayed a type I interferon response state. Integrative analyses of a large-cohort COVID-19 dataset and our datasets revealed that NK2.1 cells from elderly COVID-19 patients are enriched for type I interferon signalling, which is positively correlated with disease severity in COVID-19. CONCLUSIONS: We identified a unique memory-like NK cell subset that accumulates with ageing and correlates with disease severity in COVID-19. Our results identify memory-like NK2.1 cells as a potential target for developing immunotherapies for infectious diseases and for addressing age-related dysfunctions of the immune system.

miR-146a promotes M2 macrophage polarization and accelerates diabetic wound healing by inhibiting the TLR4/NF-κB axis.

We tried to unveil the clinical significance of miR-146a as a biomarker in M2 macrophage polarization in diabetic wound healing. Initially, we found reduced miR-146a in macrophages of diabetic patients. Next, dual-luciferase assay verified that toll-like receptor 4 (TLR4) was a target gene of miR-146 and was negatively regulated by miR-146. Moreover, after ectopic expression and depletion experiments of miR-146 and/or TLR4, lipopolysaccharide-induced inflammatory response of macrophages was detected. The results revealed that overexpression of miR-146a promoted the M2 macrophage polarization by suppressing the TLR4/nuclear factor-kappaB (NF-κB) axis, so as to enhance wound healing in diabetic ulcers. Further, mouse models with diabetic ulcers were established to investigate the effects of miR-146a on diabetic wound healing in vivo, which revealed that miR-146a promoted wound healing in diabetic ulcers by inhibiting the TLR4/NF-κB axis. In conclusion, we demonstrate that miR-146a can induce M2 macrophage polarization to enhance wound healing in diabetic ulcers by inhibiting the TLR4/NF-κB axis.

MicroRNA-10a/b inhibit TGF-ß/Smad-induced renal fibrosis by targeting TGF-ß receptor 1 in diabetic kidney disease.

TGF-ß/Smad signaling plays a vital role in the development of fibrosis in diabetic kidney disease (DKD). However, remedies targeting key elements in TGF-ß/Smad signaling are lacking. Here, we found that TGF-ß receptor 1 (TGFBR1), a key protein in TGF-ß/Smad signaling, was upregulated in kidney from diabetic mice and patients with DKD. Induction of TGFBR1 was regulated by microRNA-10a and -10b (miR-10a/b) by a post-transcriptional mechanism. Furthermore, the decreased XRN2, an exoribonuclease, was identified to contribute to affecting miR-10a/b maturation in vitro. In streptozotocin (STZ)-induced DKD mice, preventing the reduction of miR-10a/b in the kidney by an in situ lentivirus-injection method attenuated collagen deposition and foot process effacement, whereas deprivation of miR-10a/b aggravated renal fibrosis. Mechanistically, manipulating miR-10a/b in the kidney influenced TGFBR1 protein expression, TGF-ß/Smad signaling activation, and downstream pro-fibrotic genes expression including fibronectin (FN) and α-smooth muscle actin (α-SMA). In a cohort of patients diagnosed DKD, renal miR-10a/b expressions were downregulated, whereas both TGFBR1 and fibrosis were enhanced. Our finding suggests that overexpressing miR-10a/b in kidney may be a promising method for the treatment of fibrosis in DKD.

In Vivo Visualized Tracking of Tumor-Derived Extracellular Vesicles Using CRISPR-Cas9 System.

Introduction: Tumor extracellular vesicles (EVs) and their relevance to various processes of tumor growth have been vigorously investigated over the past decade. However, obtaining direct evidence of spontaneous EV transfer in vivo remains challenging. In our previous study, a single-guide RNA (sgRNA): Cas9 ribonucleoprotein complex, which can efficiently delete target genes, was delivered into recipient cells using an engineered EV. Aim: Applying this newly discovered exosomal bio-cargo to track the uptake and distribution of tumor EVs. Methods: Tumor cells of interest were engineered to express and release the sgRNA:Cas9 complex, and a reporter cell/system containing STOP-fluorescent protein (FP) elements was also generated. EV-delivered Cas9 proteins from donor cells were programmed by a pair of sgRNAs to completely delete a blockade sequence and, in turn, recuperated the expression of FP in recipient reporter cells. Thus, fluorescently illuminated cells indicate the uptake of EVs. To improve the efficiency and sensitivity of this tracking system in vivo, we optimized the sgRNA design, which could more efficiently trigger the expression of reporter proteins. Results: We demonstrated the EV-mediated crosstalk between tumor cells, and between tumor cells and normal cells in vitro. In vivo, we showed that intravenously administered EVs can be taken up by the liver. Moreover, we showed that EVs derived from melanoma xenografts in vivo preferentially target the brain and liver. This distribution resembles the manifestation of organotrophic metastasis of melanoma. Conclusion: This study provides an alternative tool to study the distribution and uptake of tumor EVs.

Single-cell analysis of COVID-19, sepsis, and HIV infection reveals hyperinflammatory and immunosuppressive signatures in monocytes.

The mortality risk of coronavirus disease 2019 (COVID-19) patients has been linked to the cytokine storm caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Understanding the inflammatory responses shared between COVID-19 and other infectious diseases that feature cytokine storms may therefore help in developing improved therapeutic strategies. Here, we use integrative analysis of single-cell transcriptomes to characterize the inflammatory signatures of peripheral blood mononuclear cells from patients with COVID-19, sepsis, and HIV infection. We identify ten hyperinflammatory cell subtypes in which monocytes are the main contributors to the transcriptional differences in these infections. Monocytes from COVID-19 patients share hyperinflammatory signatures with HIV infection and immunosuppressive signatures with sepsis. Finally, we construct a "three-stage" model of heterogeneity among COVID-19 patients, related to the hyperinflammatory and immunosuppressive signatures in monocytes. Our study thus reveals cellular and molecular insights about inflammatory responses to SARS-CoV-2 infection and provides therapeutic guidance to improve treatments for subsets of COVID-19 patients.

Profiling of the immune repertoire in COVID-19 patients with mild, severe, convalescent, or retesting-positive status.

Forty-seven samples of peripheral blood mononuclear cells from four groups of coronavirus disease (COVID)-19 patients (mild, severe, convalescent, retesting-positive) and healthy controls were applied to profile the immune repertoire of COVID-19 patients in acute infection or convalescence by transcriptome sequencing and immune-receptor repertoire (IRR) sequencing. Transcriptome analyses showed that genes within principal component group 1 (PC1) were associated with infection and disease severity whereas genes within PC2 were associated with recovery from COVID-19. A "dual-injury mechanism" of COVID-19 severity was related to an increased number of proinflammatory pathways and activated hypercoagulable pathways. A machine-learning model based on the genes associated with inflammatory and hypercoagulable pathways had the potential to be employed to monitor COVID-19 severity. Signature analyses of B-cell receptors (BCRs) and T-cell receptors (TCRs) revealed the dominant selection of longer V-J pairs (e.g., IGHV3-9-IGHJ6 and IGHV3-23-IGHJ6) and continuous tyrosine motifs in BCRs and lower diversity of TCRs. These findings provide potential predictors for COVID-19 outcomes, and new potential targets for COVID-19 treatment.

Single-cell profiling of the human decidual immune microenvironment in patients with recurrent pregnancy loss.

Maintaining homeostasis of the decidual immune microenvironment at the maternal-fetal interface is essential for placentation and reproductive success. Although distinct decidual immune cell subpopulations have been identified under normal conditions, systematic understanding of the spectrum and heterogeneity of leukocytes under recurrent miscarriage in human deciduas remains unclear. To address this, we profiled the respective transcriptomes of 18,646 primary human decidual immune cells isolated from patients with recurrent pregnancy loss (RPL) and healthy controls at single-cell resolution. We discovered dramatic differential distributions of immune cell subsets in RPL patients compared with the normal decidual immune microenvironment. Furthermore, we found a subset of decidual natural killer (NK) cells that support embryo growth were diminished in proportion due to abnormal NK cell development in RPL patients. We also elucidated the altered cellular interactions between the decidual immune cell subsets in the microenvironment and those of the immune cells with stromal cells and extravillous trophoblast under disease state. These results provided deeper insights into the RPL decidual immune microenvironment disorder that are potentially applicable to improve the diagnosis and therapeutics of this disease.

Salvianolic acid A inhibits the growth of diffuse large B-cell lymphoma through MAPK pathways.

Treatment options are limited in patients with diffuse large B-cell lymphoma (DLBCL). Salvianolic acid A (SAA) is a water-soluble phenolic acid extracted from Salvia miltiorrhiza (Danshen) with anti-tumor properties. The anti-leukemic activity of SAA found in our recent research prompted us to investigate the therapeutic effect and mechanism of action of SAA in DLBCL. In the work described here, we found that SAA inhibited the viability of DLBCL cells by inducing cellular apoptosis, which was accompanied by upregulation of Bax and cleavage of PARP. Pre-incubation of SAA increased the phosphorylation of JNK, while it decreased the phosphorylation of p38 and ERK in DLBCL cells. Importantly, pharmacologic JNK inhibition partially mitigated the anti-survival effect of SAA, and inhibition of p38 and ERK synergized with SAA. Furthermore, SAA suppressed DLBCL tumor growth in a xenograft mouse model in vivo. Therefore, our data suggest the therapeutic utility of SAA in the management of DLBCL.

Diagnostic value of optic disc retinal nerve fiber layer thickness for diabetic peripheral neuropathy.

OBJECTIVE: To investigate the value of optic disc retinal nerve fiber layer (RNFL) thickness in the diagnosis of diabetic peripheral neuropathy (DPN). METHODS: Ninety patients with type 2 diabetes, including 60 patients without DPN (NDPN group) and 30 patients with DPN (DPN group), and 30 healthy participants (normal group) were enrolled. Optical coherence tomography (OCT) was used to measure the four quadrants and the overall average RNFL thickness of the optic disc. The receiver operator characteristic curve was drawn and the area under the curve (AUC) was calculated to evaluate the diagnostic value of RNFL thickness in the optic disc area for DPN. RESULTS: The RNFL thickness of the DPN group was thinner than those of the normal and NDPN groups in the overall average ((101.07± 12.40) µm vs. (111.07±6.99) µm and (109.25±6.90) µm), superior quadrant ((123.00±19.04) µm vs. (138.93±14.16) µm and (134.47±14.34) µm), and inferior quadrant ((129.37±17.50) µm vs. (143.60±12.22) µm and (144.48±14.10) µm), and the differences were statistically significant. The diagnostic efficiencies of the overall average, superior quadrant, and inferior quadrant RNFL thicknesses, and a combined index of superior and inferior quadrant RNFL thicknesses were similar, and the AUCs were 0.739 (95% confidence interval (CI) 0.635-0.826), 0.683 (95% CI 0.576-0.778), 0.755 (95% CI 0.652-0.840), and 0.773 (95% CI 0.672-0.854), respectively. The diagnostic sensitivity of RNFL thickness in the superior quadrant reached 93.33%. CONCLUSIONS: The thickness of the RNFL in the optic disc can be used as a diagnostic method for DPN.

Single Cell Sequencing: A New Dimension in Cancer Diagnosis and Treatment.

Cancer is one of the leading causes of death worldwide and well known for its complexity. Cancer cells within the same tumor or from different tumors are highly heterogeneous. Furthermore, stromal and immune cells within tumor microenvironment interact with cancer cells to play important roles in how tumors progress and respond to different treatments. Recent advances in single cell technologies, especially massively parallel single cell sequencing, have made it possible to analyze cancer cells and cells in its tumor microenvironment in parallel with unprecedented high resolution. In this chapter, we will review recent developments in single cell sequencing technologies and their applications in cancer research. We will also explain how insights generated from single cell sequencing can be used to develop novel diagnostic and therapeutic approaches to conquer cancer.

Single-cell analysis of two severe COVID-19 patients reveals a monocyte-associated and tocilizumab-responding cytokine storm.

Several studies show that the immunosuppressive drugs targeting the interleukin-6 (IL-6) receptor, including tocilizumab, ameliorate lethal inflammatory responses in COVID-19 patients infected with SARS-CoV-2. Here, by employing single-cell analysis of the immune cell composition of two severe-stage COVID-19 patients prior to and following tocilizumab-induced remission, we identify a monocyte subpopulation that contributes to the inflammatory cytokine storms. Furthermore, although tocilizumab treatment attenuates the inflammation, immune cells, including plasma B cells and CD8+ T cells, still exhibit robust humoral and cellular antiviral immune responses. Thus, in addition to providing a high-dimensional dataset on the immune cell distribution at multiple stages of the COVID-19, our work also provides insights into the therapeutic effects of tocilizumab, and identifies potential target cell populations for treating COVID-19-related cytokine storms.

APEC: an accesson-based method for single-cell chromatin accessibility analysis.

The development of sequencing technologies has promoted the survey of genome-wide chromatin accessibility at single-cell resolution. However, comprehensive analysis of single-cell epigenomic profiles remains a challenge. Here, we introduce an accessibility pattern-based epigenomic clustering (APEC) method, which classifies each cell by groups of accessible regions with synergistic signal patterns termed "accessons". This python-based package greatly improves the accuracy of unsupervised single-cell clustering for many public datasets. It also predicts gene expression, identifies enriched motifs, discovers super-enhancers, and projects pseudotime trajectories. APEC is available at https://github.com/QuKunLab/APEC.