The transcription factor TCF-1 enforces commitment to the innate lymphoid cell lineage.

Innate lymphoid cells (ILCs) play important functions in immunity and tissue homeostasis, but their development is poorly understood. Through the use of single-cell approaches, we examined the transcriptional and functional heterogeneity of ILC progenitors, and studied the precursor-product relationships that link the subsets identified. This analysis identified two successive stages of ILC development within T cell factor 1-positive (TCF-1+) early innate lymphoid progenitors (EILPs), which we named 'specified EILPs' and 'committed EILPs'. Specified EILPs generated dendritic cells, whereas this potential was greatly decreased in committed EILPs. TCF-1 was dispensable for the generation of specified EILPs, but required for the generation of committed EILPs. TCF-1 used a pre-existing regulatory landscape established in upstream lymphoid precursors to bind chromatin in EILPs. Our results provide insight into the mechanisms by which TCF-1 promotes developmental progression of ILC precursors, while constraining their dendritic cell lineage potential and enforcing commitment to ILC fate.

Transcription factors TCF-1 and GATA3 are key factors for the epigenetic priming of early innate lymphoid progenitors toward distinct cell fates.

The differentiation of innate lymphoid cells (ILCs) from hematopoietic stem cells needs to go through several multipotent progenitor stages. However, it remains unclear whether the fates of multipotent progenitors are predefined by epigenetic states. Here, we report the identification of distinct accessible chromatin regions in all lymphoid progenitors (ALPs), EILPs, and ILC precursors (ILCPs). Single-cell MNase-seq analyses revealed that EILPs contained distinct subpopulations epigenetically primed toward either dendritic cell lineages or ILC lineages. We found that TCF-1 and GATA3 co-bound to the lineage-defining sites for ILCs (LDS-Is), whereas PU.1 binding was enriched in the LDSs for alternative dendritic cells (LDS-As). TCF-1 and GATA3 were indispensable for the epigenetic priming of LDSs at the EILP stage. Our results suggest that the multipotency of progenitor cells is defined by the existence of a heterogeneous population of cells epigenetically primed for distinct downstream lineages, which are regulated by key transcription factors.

GLYR1 transcriptionally regulates PER3 expression to promote the proliferation and migration of multiple myeloma.

Period circadian regulator 3 (PER3) functions as a tumor suppressor in various cancers. However, the role of PER3 in multiple myeloma (MM) has not been reported yet. Through this study, we aimed to investigate the potential role of PER3 in MM and the underlying mechanisms. RT-qPCR and western blotting were used to determine the mRNA and protein expression levels of PER3. Glyoxylate reductase 1 homolog (GLYR1) was predicted to be a transcription factor of PER3. The binding sites of GLYR1 on the promoter region of PER3 were analyzed using UCSC and confirmed using luciferase and chromatin immunoprecipitation assays. Viability, apoptosis, and metathesis were determined using CCK-8, colony formation, TUNEL, and transwell assays. We found that PER3 expression decreased in MM. Low PER3 levels may predict poor survival rates; PER3 overexpression suppresses the viability and migration of MM cells and promotes apoptosis. Moreover, GLYR1 transcriptionally activates PER3, and the knockdown of PER3 alleviates the effects of GLYR1 and induces its malignant behavior in MM cells. To conclude, GLYR1 upregulates PER3 and suppresses the aggressive behavior of MM cells, suggesting that GLYR1/PER3 signaling may be a potential therapeutic target for MM.

Principles of nucleosome organization revealed by single-cell micrococcal nuclease sequencing.

Nucleosome positioning is critical to chromatin accessibility and is associated with gene expression programs in cells1-3. Previous nucleosome mapping methods assemble profiles from cell populations and reveal a cell-averaged pattern: nucleosomes are positioned and form a phased array that surrounds the transcription start sites of active genes3-6 and DNase I hypersensitive sites7. However, even in a homogenous population of cells, cells exhibit heterogeneity in expression in response to active signalling8,9 that may be related to heterogeneity in chromatin accessibility10-12. Here we report a technique, termed single-cell micrococcal nuclease sequencing (scMNase-seq), that can be used to simultaneously measure genome-wide nucleosome positioning and chromatin accessibility in single cells. Application of scMNase-seq to NIH3T3 cells, mouse primary naive CD4 T cells and mouse embryonic stem cells reveals two principles of nucleosome organization: first, nucleosomes in heterochromatin regions, or that surround the transcription start sites of silent genes, show large variation in positioning across different cells but are highly uniformly spaced along the nucleosome array; and second, nucleosomes that surround the transcription start sites of active genes and DNase I hypersensitive sites show little variation in positioning across different cells but are relatively heterogeneously spaced along the nucleosome array. We found a bimodal distribution of nucleosome spacing at DNase I hypersensitive sites, which corresponds to inaccessible and accessible states and is associated with nucleosome variation and variation in accessibility across cells. Nucleosome variation is smaller within single cells than across cells, and smaller within the same cell type than across cell types. A large fraction of naive CD4 T cells and mouse embryonic stem cells shows depleted nucleosome occupancy at the de novo enhancers detected in their respective differentiated lineages, revealing the existence of cells primed for differentiation to specific lineages in undifferentiated cell populations.

Publisher Correction: Principles of nucleosome organization revealed by single-cell micrococcal nuclease sequencing.

Change history: In Fig. 1c of this Letter, the two graphs were duplicates. The right panel of Fig. 1c has been corrected online.

Survival prediction in acute myeloid leukemia using gene expression profiling.

BACKGROUND: Acute myeloid leukemia (AML) is a genetically heterogeneous blood disorder. AML patients are associated with a relatively poor overall survival. The objective of this study was to establish a machine learning model to accurately perform the prognosis prediction in AML patients. METHODS: We first screened for prognosis-related genes using Kaplan-Meier survival analysis in The Cancer Genome Atlas dataset and validated the results in the Oregon Health & Science University dataset. With a random forest model, we built a prognostic risk score using patient's age, TP53 mutation, ELN classification and normalized 197 gene expression as predictor variable. Gene set enrichment analysis was implemented to determine the dysregulated gene sets between the high-risk and low-risk groups. Similarity Network Fusion (SNF)-based integrative clustering was performed to identify subgroups of AML patients with different clinical features. RESULTS: The random forest model was deemed the best model (area under curve value, 0.75). The random forest-derived risk score exhibited significant association with shorter overall survival in AML patients. The gene sets of pantothenate and coa biosynthesis, glycerolipid metabolism, biosynthesis of unsaturated fatty acids were significantly enriched in phenotype high risk score. SNF-based integrative clustering indicated three distinct subsets of AML patients in the TCGA cohort. The cluster3 AML patients were characterized by older age, higher risk score, more frequent TP53 mutations, higher cytogenetics risk, shorter overall survival. CONCLUSIONS: The random forest-based risk score offers an effective method to perform prognosis prediction for AML patients.

Trac-looping measures genome structure and chromatin accessibility.

Long-range chromatin interactions play critical roles in genome organization and regulation of transcription. We now report transposase-mediated analysis of chromatin looping (Trac-looping) for simultaneous detection of multiscale genome-wide chromatin interactions among regulatory elements and chromatin accessibility. With this technique, a bivalent oligonucleotide linker is inserted between two interacting regions such that the chromatin interactions are captured without prior chromatin fragmentation and proximity-based ligation. Application of Trac-looping to human CD4+ T cells revealed substantial reorganization of enhancer-promoter interactions associated with changes in gene expression after T cell receptor stimulation.

First case of bloodstream infection of Trichosporon loubieri in a patient with B-cell lymphocytic leukemia in China.

Trichosporon loubieri is an opportunistic pathogenic fungus that could causes invasive infections for human, which rarely been reported. In the present study, we reported a case of bloodstream infection from a patient with Ph-positive B-cell acute lymphocytic leukemia (B-ALL) due to Trichosporon loubieri. Trichosporon loubieri was identified by Internal Transcribed Spacer (ITS) gene sequencing. The patient was treated by intravenous voriconazole (VCZ) and amphotericin B (AmB) according to antifungal susceptibility testing and he was still alive so far. To the best of our knowledge, this is the fourth report of human bloodstream infection due to Trichosporon loubieri and the first survival case of its kind in an immunocompromised patient.

MLL3/MLL4 are required for CBP/p300 binding on enhancers and super-enhancer formation in brown adipogenesis.

Histone H3K4me1/2 methyltransferases MLL3/MLL4 and H3K27 acetyltransferases CBP/p300 are major enhancer epigenomic writers. To understand how these epigenomic writers orchestrate enhancer landscapes in cell differentiation, we have profiled genomic binding of MLL4, CBP, lineage-determining transcription factors (EBF2, C/EBPß, C/EBPα, PPARγ), coactivator MED1, RNA polymerase II, as well as epigenome (H3K4me1/2/3, H3K9me2, H3K27me3, H3K36me3, H3K27ac), transcriptome and chromatin opening during adipogenesis of immortalized preadipocytes derived from mouse brown adipose tissue (BAT). We show that MLL4 and CBP drive the dynamic enhancer epigenome, which correlates with the dynamic transcriptome. MLL3/MLL4 are required for CBP/p300 binding on enhancers activated during adipogenesis. Further, MLL4 and CBP identify super-enhancers (SEs) of adipogenesis and that MLL3/MLL4 are required for SE formation. Finally, in brown adipocytes differentiated in culture, MLL4 identifies primed SEs of genes fully activated in BAT such as Ucp1. Comparison of MLL4-defined SEs in brown and white adipogenesis identifies brown-specific SE-associated genes that could be involved in BAT functions. These results establish MLL3/MLL4 and CBP/p300 as master enhancer epigenomic writers and suggest that enhancer-priming by MLL3/MLL4 followed by enhancer-activation by CBP/p300 sequentially shape dynamic enhancer landscapes during cell differentiation. Our data also provide a rich resource for understanding epigenomic regulation of brown adipogenesis.

Enhancer priming by H3K4 methyltransferase MLL4 controls cell fate transition.

Transcriptional enhancers control cell-type-specific gene expression. Primed enhancers are marked by histone H3 lysine 4 (H3K4) mono/di-methylation (H3K4me1/2). Active enhancers are further marked by H3K27 acetylation (H3K27ac). Mixed-lineage leukemia 4 (MLL4/KMT2D) is a major enhancer H3K4me1/2 methyltransferase with functional redundancy with MLL3 (KMT2C). However, its role in cell fate maintenance and transition is poorly understood. Here, we show in mouse embryonic stem cells (ESCs) that MLL4 associates with, but is surprisingly dispensable for the maintenance of, active enhancers of cell-identity genes. As a result, MLL4 is dispensable for cell-identity gene expression and self-renewal in ESCs. In contrast, MLL4 is required for enhancer-binding of H3K27 acetyltransferase p300, enhancer activation, and induction of cell-identity genes during ESC differentiation. MLL4 protein, rather than MLL4-mediated H3K4 methylation, controls p300 recruitment to enhancers. We also show that, in somatic cells, MLL4 is dispensable for maintaining cell identity but essential for reprogramming into induced pluripotent stem cells. These results indicate that, although enhancer priming by MLL4 is dispensable for cell-identity maintenance, it controls cell fate transition by orchestrating p300-mediated enhancer activation.

[Analysis of mutations of 14 genes among 87 patients with myelodysplastic syndrome].

OBJECTIVE: To explore the correlation of genetic mutations and clinical features of myelodysplastic syndromes (MDS) with scores of Revised International Prognostic Scoring System (IPSS-R). METHODS: Eighty-seven patients with de novo MDS were enrolled. Mutations of MDS-related genes and clinical features were used to determine the incidence and subtype of mutations. Clinical features and IPSS-R scores of the patients with high frequency mutations involving TET2, TP53, ASXL1, RUNX1 and SF3B1 genes were compared. RESULTS: Fifty-four patients (62.1%) harbored at least one point mutation. The incidences of various mutations were significantly different, with the incidence of MDS-EB-2 being 100% and MDS-SLD being only 38.9%. Compared with the wild types, patients harboring mutations had higher lactate dehydrogenase, higher ß2 microglobulin, higher percentage of bone marrow blast cells and lower hemoglobin levels (P=0.027, <0.01, <0.01, 0.046, respectively). The IPSS-R scores of MDS patients with mutations were significantly higher than the wild types (P<0.01). The IPSS-R scores of the TP53 mutation groups were 7.82±1.83, which was significantly higher than the control group (3.77±1.66, P<0.01). No difference was found between the IPSS-R between patients carrying TET2, ASXL1, RUNX1, and SF3B1 mutations or the wild types (P>0.05). CONCLUSION: Genetic mutations are commonly found in MDS. MDS patients with mutations have unique clinical laboratory characteristics. Although the prognostic value of most genes is controversial, TP53 is an definite indicator of poor prognosis.

Gcn5 and PCAF negatively regulate interferon-ß production through HAT-independent inhibition of TBK1.

Viral infection triggers innate immune signaling, which in turn induces interferon-ß (IFN-ß) production to establish innate antiviral immunity. Previous studies showed that Gcn5 (Kat2a), a histone acetyltransferase (HAT) with partial functional redundancy with PCAF (Kat2b), and Gcn5/PCAF-mediated histone H3K9 acetylation (H3K9ac) are enriched on the active IFNB gene promoter. However, whether Gcn5/PCAF and H3K9ac regulate IFN-ß production is unknown. Here, we show that Gcn5/PCAF-mediated H3K9ac correlates well with, but is surprisingly dispensable for, the expression of endogenous IFNB and the vast majority of active genes in fibroblasts. Instead, Gcn5/PCAF repress IFN-ß production and innate antiviral immunity in several cell types in a HAT-independent and non-transcriptional manner: by inhibiting the innate immune signaling kinase TBK1 in the cytoplasm. Our results thus identify Gcn5 and PCAF as negative regulators of IFN-ß production and innate immune signaling.

InteMAP: Integrated metagenomic assembly pipeline for NGS short reads.

BACKGROUND: Next-generation sequencing (NGS) has greatly facilitated metagenomic analysis but also raised new challenges for metagenomic DNA sequence assembly, owing to its high-throughput nature and extremely short reads generated by sequencers such as Illumina. To date, how to generate a high-quality draft assembly for metagenomic sequencing projects has not been fully addressed. RESULTS: We conducted a comprehensive assessment on state-of-the-art de novo assemblers and revealed that the performance of each assembler depends critically on the sequencing depth. To address this problem, we developed a pipeline named InteMAP to integrate three assemblers, ABySS, IDBA-UD and CABOG, which were found to complement each other in assembling metagenomic sequences. Making a decision of which assembling approaches to use according to the sequencing coverage estimation algorithm for each short read, the pipeline presents an automatic platform suitable to assemble real metagenomic NGS data with uneven coverage distribution of sequencing depth. By comparing the performance of InteMAP with current assemblers on both synthetic and real NGS metagenomic data, we demonstrated that InteMAP achieves better performance with a longer total contig length and higher contiguity, and contains more genes than others. CONCLUSIONS: We developed a de novo pipeline, named InteMAP, that integrates existing tools for metagenomics assembly. The pipeline outperforms previous assembly methods on metagenomic assembly by providing a longer total contig length, a higher contiguity and covering more genes. InteMAP, therefore, could potentially be a useful tool for the research community of metagenomics.

Transcriptional profile of human thymus reveals IGFBP5 is correlated with age-related thymic involution.

Thymus is the main immune organ which is responsible for the production of self-tolerant and functional T cells, but it shrinks rapidly with age after birth. Although studies have researched thymus development and involution in mouse, the critical regulators that arise with age in human thymus remain unclear. We collected public human single-cell transcriptomic sequencing (scRNA-seq) datasets containing 350,678 cells from 36 samples, integrated them as a cell atlas of human thymus. Clinical samples were collected and experiments were performed for validation. We found early thymocyte-specific signaling and regulons which played roles in thymocyte migration, proliferation, apoptosis and differentiation. Nevertheless, signaling patterns including number, strength and path completely changed during aging, Transcription factors (FOXC1, MXI1, KLF9, NFIL3) and their target gene, IGFBP5, were resolved and up-regulated in aging thymus and involved in promoting epithelial-mesenchymal transition (EMT), responding to steroid and adipogenesis process of thymic epithelial cell (TECs). Furthermore, we validated that IGFBP5 protein increased at TECs and Hassall's corpuscle in both human and mouse aging thymus and knockdown of IGFBP5 significantly increased the expression of proliferation-related genes in thymocytes. Collectively, we systematically explored cell-cell communications and regulons of early thymocytes as well as age-related differences in human thymus by using both bioinformatic and experimental verification, indicating IGFBP5 as a functional marker of thymic involution and providing new insights into the mechanisms of thymus involution.

A de novo metagenomic assembly program for shotgun DNA reads.

MOTIVATION: A high-quality assembly of reads generated from shotgun sequencing is a substantial step in metagenome projects. Although traditional assemblers have been employed in initial analysis of metagenomes, they cannot surmount the challenges created by the features of metagenomic data. RESULT: We present a de novo assembly approach and its implementation named MAP (metagenomic assembly program). Based on an improved overlap/layout/consensus (OLC) strategy incorporated with several special algorithms, MAP uses the mate pair information, resulting in being more applicable to shotgun DNA reads (recommended as >200 bp) currently widely used in metagenome projects. Results of extensive tests on simulated data show that MAP can be superior to both Celera and Phrap for typical longer reads by Sanger sequencing, as well as has an evident advantage over Celera, Newbler and the newest Genovo, for typical shorter reads by 454 sequencing. AVAILABILITY AND IMPLEMENTATION: The source code of MAP is distributed as open source under the GNU GPL license, the MAP program and all simulated datasets can be freely available at http://bioinfo.ctb.pku.edu.cn/MAP/

Single-Cell RNA Sequencing Reveals Cellular Heterogeneity in an Acral Amelanotic Melanoma After Immunotherapy Treatment.

Background: Anti-programmed cell death ligand-1 (anti-PD-L1) immunotherapy is often used for advanced urothelial carcinoma and melanoma, including amelanotic melanoma, a relatively rare subtype with little to no pigment in the tumor cells. However, cellular heterogeneity of amelanotic melanoma during or after anti-PD-L1 immunotherapy treatments has not been described. Purpose: To investigate cellular heterogeneity in acral amelanotic melanoma after immunotherapy exposure. Methods: We evaluated subtle visual changes of the melanoma by dermoscopy and performed a pathological examination to analyze the heterogeneity of microscopic morphological and immunohistochemistry changes. The cellular transcriptional heterogeneity and corresponding biological function profiles of the melanoma were determined by single-cell RNA sequencing (scRNA-seq). Results: The dermoscopic examination revealed black globules and scar-like depigmentation areas against a homogeneous red background. Pigmented and amelanotic melanoma cells were observed microscopically. The pigmented cells were large and contained melanin granules expressing Melan-A and HMB45; the amelanotic cells were small and did not express HMB45. Ki-67 immunohistochemical staining revealed that the pigmented melanoma cells had a higher proliferative ability than the amelanotic cells. scRNA-seq identified three cell clusters: amelanotic cell cluster 1, amelanotic cell cluster 2, and pigmented cell cluster. Furthermore, a pseudo-time trajectory analysis showed that amelanotic cell cluster 2 originated from amelanotic cell cluster 1 and transformed into the pigmented melanoma cell cluster. The expression pattern of melanin synthesis-related and lysosome-endosome-related genes in different cell clusters supported the cell cluster transformation results. Also, upregulated expression of cell cycle genes indicated that the pigmented melanoma cells had a high proliferative ability. Conclusion: Coexisting amelanotic and pigmented melanoma cells indicated cellular heterogeneity in an acral amelanotic melanoma from a patient who underwent immunotherapy treatment. Additionally, the pigmented melanoma cells acquired a higher proliferative ability than the amelanotic melanoma cells.

Assessing Biomaterial-Induced Stem Cell Lineage Fate by Machine Learning-Based Artificial Intelligence.

Current functional assessment of biomaterial-induced stem cell lineage fate in vitro mainly relies on biomarker-dependent methods with limited accuracy and efficiency. Here a "Mesenchymal stem cell Differentiation Prediction (MeD-P)" framework for biomaterial-induced cell lineage fate prediction is reported. MeD-P contains a cell-type-specific gene expression profile as a reference by integrating public RNA-seq data related to tri-lineage differentiation (osteogenesis, chondrogenesis, and adipogenesis) of human mesenchymal stem cells (hMSCs) and a predictive model for classifying hMSCs differentiation lineages using the k-nearest neighbors (kNN) strategy. It is shown that MeD-P exhibits an overall accuracy of 90.63% on testing datasets, which is significantly higher than the model constructed based on canonical marker genes (80.21%). Moreover, evaluations of multiple biomaterials show that MeD-P provides accurate prediction of lineage fate on different types of biomaterials as early as the first week of hMSCs culture. In summary, it is demonstrated that MeD-P is an efficient and accurate strategy for stem cell lineage fate prediction and preliminary biomaterial functional evaluation.

Adoption value of support vector machine algorithm-based computed tomography imaging in the diagnosis of secondary pulmonary fungal infections in patients with malignant hematological disorders.

This study aimed to assess the feasibility of diagnosing secondary pulmonary fungal infections (PFIs) in patients with hematological malignancies (HM) using computerized tomography (CT) imaging and a support vector machine (SVM) algorithm. A total of 100 patients with HM complicated by secondary PFI underwent CT scans, and they were included in the training group. Concurrently, 80 patients with the same underlying disease who were treated at our institution were included in the test group. The types of pathogens among different PFI patients and the CT imaging features were compared. Radiomic features were extracted from the CT imaging data of patients, and a diagnostic SVM model was constructed by integrating these features with clinical characteristics. Aspergillus was the most common pathogen responsible for PFIs, followed by Candida, Pneumocystis jirovecii, Mucor, and Cryptococcus, in descending order of occurrence. Patients typically exhibited bilateral diffuse lung lesions. Within the SVM algorithm model, six radiomic features, namely the square root of the inverse covariance of the gray-level co-occurrence matrix (square root IV), the square root of the inverse covariance of the gray-level co-occurrence matrix, and small dependency low gray-level emphasis, significantly influenced the diagnosis of secondary PFIs in patients with HM. The area under the curve values for the training and test sets were 0.902 and 0.891, respectively. Therefore, CT images based on the SVM algorithm demonstrated robust predictive capability in diagnosing secondary PFIs in conjunction with HM.

Transcriptomic Heterogeneity of Skin Across Different Anatomic Sites.

Multiomic studies, including RNA sequencing, single-cell RNA sequencing, and epigenomics, can provide insight into the connection between anatomically heterogeneous gene expression profile of the skin and dermatoses-predisposed sites, in which RNA sequencing is essential. Therefore, in this study, 159 skin samples collected mainly from discarded normal skin tissue during surgical treatment for benign skin tumors were used for RNA sequencing. On the basis of cluster analysis, the skin was divided into four regions, with each region showing specific physiological characteristics through differentially expressed gene analysis. The results showed that the head and neck region, perineum, and palmoplantar area were closely associated with lipid metabolism, hormone metabolism, blood circulation, and related neural regulation, respectively. Transcription factor enrichment indicated that different regions were associated with the development of adjacent tissues. Specifically, the head and neck region, trunk and extremities, perineum, and palmoplantar area were associated with the central nervous, axial, urogenital, and vascular systems, respectively. The results were imported into an open website (https://dermvis.github.io/) for retrieval. Our transcriptomic data elucidated that human skin exhibits transcriptomic heterogeneity reflecting physiological and developmental variation at different anatomic sites and provided guidance for further studies on skin development and dermatoses predisposed sites.

The gut microbiota correlate with the disease characteristics and immune status of patients with untreated diffuse large B-cell lymphoma.

Background: Gut microbiota characteristics in patients with diffuse large B-cell lymphoma (DLBCL) are reportedly different when compared with the healthy population and it remains unclear if the gut microbiota affects host immunity and clinical disease features. This research investigated the gut microbiota in patients with untreated DLBCL and analyzed its correlation with patient clinical characteristics, humoral, and cell immune status. Methods: Thirty-five patients with untreated DLBCL and 20 healthy controls (HCs) were recruited to this study and microbiota differences in stool samples were analyzed by 16S rDNA sequencing. Absolute ratios of immune cell subset counts in peripheral blood were detected by flow cytometry and peripheral blood cytokine levels were detected by enzyme-linked immunosorbent assay. Relationships between changes in patient microbiomes and clinical characteristics, such as clinical stage, international prognostic index (IPI) risk stratification, cell origin, organ involved and treatment responses were investigated and correlations between differential microbiota and host immune indices were analyzed. Results: The alpha-diversity index of intestinal microecology in DLBCL patients was not significantly different when compared with HCs (P>0.05), nonetheless beta-diversity was significantly decreased (P=0.001). p_Proteobacteria were dominant in DLBCL, while p_Bacteroidetes abundance was significantly decreased when compared with HCs (P<0.05). Gut microbiota characteristics were identified that were associated with clinical features, such as tumor load, risk stratification and cell origin, and correlation analyses were performed between differential flora abundance associated with these clinical features and host immune status. The p_Firmicutes was positively correlated with absolute lymphocyte values, g_Prevotella_2 and s_un_g_Prevotella_2 were negatively correlated with absolute lymphocyte values, T cell counts and CD4 cell counts, while g_Pyramidobacter, s_un_g_Pyramidobacter, and f_Peptostreptococcaceae were negatively correlated with IgA. Conclusions: Dominant gut microbiota, abundance, diversity, and structure in DLBCL were influenced by the disease, correlated with patient immune status and this suggested that the microecology-immune axis may be involved in regulating lymphoma development. In the future, it may be possible to improve immune function in patients with DLBCL by regulating the gut microbiota, improve treatment response rates and increase patient survival rates.