Spatial transcriptomics combined with single-nucleus RNA sequencing reveals glial cell heterogeneity in the human spinal cord

Glial cells play crucial roles in regulating physiological and pathological functions,including sensation,the response to infection and acute injury,and chronic neurodegenerative disorders.Glial cells include astrocytes,microglia,and oligodendrocytes in the central nervous system,and satellite glial cells and Schwann cells in the peripheral nervous system.Despite the greater understanding of glial cell types and functional heterogeneity achieved through single-cell and single-nucleus RNA sequencing in animal models,few studies have investigated the transcriptomic profiles of glial cells in the human spinal cord.Here,we used high-throughput single-nucleus RNA sequencing and spatial transcriptomics to map the cellular and molecular heterogeneity of astrocytes,microglia,and oligodendrocytes in the human spinal cord.To explore the conservation and divergence across species,we compared these findings with those from mice.In the human spinal cord,astrocytes,microglia,and oligodendrocytes were each divided into six distinct transcriptomic subclusters.In the mouse spinal cord,astrocytes,microglia,and oligodendrocytes were divided into five,four,and five distinct transcriptomic subclusters,respectively.The comparative results revealed substantial heterogeneity in all glial cell types between humans and mice.Additionally,we detected sex differences in gene expression in human spinal cord glial cells.Specifically,in all astrocyte subtypes,the levels of NEAT1 and CHI3L1 were higher in males than in females,whereas the levels of CST3 were lower in males than in females.In all microglial subtypes,all differentially expressed genes were located on the sex chromosomes.In addition to sex-specific gene differences,the levels of MT-ND4,MT2A,MT-ATP6,MT-CO3,MT-ND2,MT-ND3,and MT-CO_(2) in all spinal cord oligodendrocyte subtypes were higher in females than in males.Collectively,the present dataset extensively characterizes glial cell heterogeneity and offers a valuable resource for exploring the cellular basis of spinal cordrelated illnesses,including chronic pain,amyotrophic lateral sclerosis,and multiple sclerosis.

Mitochondrial targeting sequence of magnetoreceptor MagR:More than just targeting

Iron-sulfur clusters(ISC)are essential cofactors for proteins involved in various biological processes,such as electron transport,biosynthetic reactions,DNA repair,and gene expression regulation.ISC assembly protein IscA1(or MagR)is found within the mitochondria of most eukaryotes.Magnetoreceptor(MagR)is a highly conserved A-type iron and iron-sulfur cluster-binding protein,characterized by two distinct types of iron-sulfur clusters,[2Fe-2S]and[3Fe-4S],each conferring unique magnetic properties.MagR forms a rod-like polymer structure in complex with photoreceptive cryptochrome(Cry)and serves as a putative magnetoreceptor for retrieving geomagnetic information in animal navigation.Although the N-terminal sequences of MagR vary among species,their specific function remains unknown.In the present study,we found that the N-terminal sequences of pigeon MagR,previously thought to serve as a mitochondrial targeting signal(MTS),were not cleaved following mitochondrial entry but instead modulated the efficiency with which iron-sulfur clusters and irons are bound.Moreover,the N-terminal region of MagR was required for the formation of a stable MagR/Cry complex.Thus,the N-terminal sequences in pigeon MagR fulfil more important functional roles than just mitochondrial targeting.These results further extend our understanding of the function of MagR and provide new insights into the origin of magnetoreception from an evolutionary perspective.

Identification and validation of a pyroptosis-related prognostic model for colorectal cancer based on bulk and single-cell RNA sequencing data

BACKGROUND Pyroptosis impacts the development of malignant tumors,yet its role in colorectal cancer(CRC)prognosis remains uncertain.AIM To assess the prognostic significance of pyroptosis-related genes and their association with CRC immune infiltration.METHODS Gene expression data were obtained from The Cancer Genome Atlas(TCGA)and single-cell RNA sequencing dataset GSE178341 from the Gene Expression Omnibus(GEO).Pyroptosis-related gene expression in cell clusters was analyzed,and enrichment analysis was conducted.A pyroptosis-related risk model was developed using the LASSO regression algorithm,with prediction accuracy assessed through K-M and receiver operating characteristic analyses.A nomo-gram predicting survival was created,and the correlation between the risk model and immune infiltration was analyzed using CIBERSORTx calculations.Finally,the differential expression of the 8 prognostic genes between CRC and normal samples was verified by analyzing TCGA-COADREAD data from the UCSC database.RESULTS An effective pyroptosis-related risk model was constructed using 8 genes-CHMP2B,SDHB,BST2,UBE2D2,GJA1,AIM2,PDCD6IP,and SEZ6L2(P<0.05).Seven of these genes exhibited differential expression between CRC and normal samples based on TCGA database analysis(P<0.05).Patients with higher risk scores demonstrated increased death risk and reduced overall survival(P<0.05).Significant differences in immune infiltration were observed between low-and high-risk groups,correlating with pyroptosis-related gene expression.CONCLUSION We developed a pyroptosis-related prognostic model for CRC,affirming its correlation with immune infiltration.This model may prove useful for CRC prognostic evaluation.

The Application of Nicotiana benthamiana as a Transient Expression Host to Clone the Coding Sequences of Plant Genes

Coding sequences (CDS) are commonly used for transient gene expression, in yeast two-hybrid screening, to verify protein interactions and in prokaryotic gene expression studies. CDS are most commonly obtained using complementary DNA (cDNA) derived from messenger RNA (mRNA) extracted from plant tissues and generated by reverse transcription. However, some CDS are difficult to acquire through this process as they are expressed at extremely low levels or have specific spatial and/or temporal expression patterns in vivo. These challenges require the development of alternative CDS cloning technologies. In this study, we found that the genomic intron-containing gene coding sequences (gDNA) from Arabidopsis thaliana, Oryza sativa, Brassica napus, and Glycine max can be correctly transcribed and spliced into mRNA in Nicotiana benthamiana. In contrast, gDNAs from Triticum aestivum and Sorghum bicolor did not function correctly. In transient expression experiments, the target DNA sequence is driven by a constitutive promoter. Theoretically, a sufficient amount of mRNA can be extracted from the N. benthamiana leaves, making it conducive to the cloning of CDS target genes. Our data demonstrate that N. benthamiana can be used as an effective host for the cloning CDS of plant genes.

Single-cell transcriptome sequencing reveals the mechanism regulating rice plumule development

Seed plumules comprise multiple developing tissues and are key sites for above-ground plant organ morphogenesis.Here,the spatial expression of genes in developing rice seed plumules was characterized by single-cell transcriptome sequencing in Zhongjiazao 17,a popular Chinese indica rice cultivar.Of 15 cell clusters,13 were assigned to cell types using marker genes and cluster-specific genes.Marker genes of multiple cell types were expressed in several clusters,suggesting a complex developmental system.Some genes for signaling by phytohormones such as abscisic acid were highly expressed in specific clusters.Various cis-elements in the promoters of genes specifically expressed in cell clusters were calculated,and some key hormone-related motifs were frequent in certain clusters.Spatial expression patterns of genes involved in rapid seed germination,seedling growth,and development were identified.These findings enhanced our understanding of cellular diversity and specialization within plumules of rice,a monocotyledonous model crop.

Possible effects of selecting different station distributions in the optimal sequence estimation method

Since the inception of the optimal sequence estimation (OSE) method,various research teams have substantiated its efficacy as the optimal stacking technique for handling array data,leading to its successful application in numerous geoscience studies.Nevertheless,concerns persist regarding the potential impact of aliasing resulting from the choice of distinct station distributions on the outcomes derived from OSE.In this investigation,I employ theoretical deduction and experimental analysis to elucidate the reasons behind the immunity of the Y_(l'm')-related common signal obtained through OSE to variations in station distribution selection.The primary objective of OSE is also underscored,i.e.,to restore/strip a Y_(l'm')-related common periodic signal from various stations.Furthermore,I provide additional clarification that the‘Y_(l'm')-related common signal’and the‘Y_(l'm')-related equivalent excitation sequence’are distinct concepts.These analyses will facilitate the utilization of the OSE technique by other researchers in investigating intriguing geophysical phenomena and attaining sound explanations.

Pig pangenome graph reveals functional features of non‑reference sequences

Background The reliance on a solitary linear reference genome has imposed a significant constraint on our compre-hensive understanding of genetic variation in animals.This constraint is particularly pronounced for non-reference sequences(NRSs),which have not been extensively studied.Results In this study,we constructed a pig pangenome graph using 21 pig assemblies and identified 23,831 NRSs with a total length of 105 Mb.Our findings revealed that NRSs were more prevalent in breeds exhibiting greater genetic divergence from the reference genome.Furthermore,we observed that NRSs were rarely found within coding sequences,while NRS insertions were enriched in immune-related Gene Ontology terms.Notably,our investigation also unveiled a close association between novel genes and the immune capacity of pigs.We observed substantial differences in terms of frequencies of NRSs between Eastern and Western pigs,and the heat-resistant pigs exhibited a substantial number of NRS insertions in an 11.6 Mb interval on chromosome X.Additionally,we discovered a 665 bp insertion in the fourth intron of the TNFRSF19 gene that may be associated with the ability of heat tolerance in South-ern Chinese pigs.Conclusions Our findings demonstrate the potential of a graph genome approach to reveal important functional features of NRSs in pig populations.

Phylogenetic, phylogeographic and divergence time analysis of Anopheles subpictus species complex using ITS2 and COI sequences

Objective:To address the phylogenetic and phylogeographic relationship between different lineages of Anopheles(An.)subpictus species complex in most parts of the Asian continent by maximum utilization of Internal Transcriber Spacer 2(ITS2)and cytochrome C oxidase I(COI)sequences deposited at the GenBank.Methods:Seventy-five ITS2,210 COI and 26 concatenated sequences available in the NCBI database were used.Phylogenetic analysis was performed using Bayesian likelihood trees,whereas median-joining haplotype networks and time-scale divergence trees were generated for phylogeographic analysis.Genetic diversity indices and genetic differentiation were also calculated.Results:Two genetically divergent molecular forms of An.subpictus species complex corresponding to sibling species A and B are established.Species A evolved around 37-82 million years ago in Sri Lanka,India,and the Netherlands,and species B evolved around 22-79 million years ago in Sri Lanka,India,and Myanmar.Vietnam,Thailand,and Cambodia have two molecular forms:one is phylogenetically similar to species B.Other forms differ from species A and B and evolved recently in the above mentioned countries,Indonesia and the Philippines.Genetic subdivision among Sri Lanka,India,and the Netherlands is almost absent.A substantial genetic differentiation was obtained for some populations due to isolation by large geographical distances.Genetic diversity indices reveal the presence of a long-established stable mosquito population,at mutation-drift equilibrium,regardless of population fluctuations.Conclusions:An.subpictus species complex consists of more than two genetically divergent molecular forms.Species A is highly divergent from the rest.Sri Lanka and India contain only species A and B.

On power series statistical convergence and new uniform integrability of double sequences

In the present paper,we mostly focus on P_(p)^(2)-statistical convergence.We will look into the uniform integrability via the power series method and its characterizations for double sequences.Also,the notions of P_(p)^(2)-statistically Cauchy sequence,P_(p)^(2)-statistical boundedness and core for double sequences will be described in addition to these findings.

Incorporating genomic annotation into single-step genomic prediction with imputed whole-genome sequence data

Single-step genomic best linear unbiased prediction(ss GBLUP) is now intensively investigated and widely used in livestock breeding due to its beneficial feature of combining information from both genotyped and ungenotyped individuals in the single model. With the increasing accessibility of whole-genome sequence(WGS) data at the population level, more attention is being paid to the usage of WGS data in ss GBLUP. The predictive ability of ss GBLUP using WGS data might be improved by incorporating biological knowledge from public databases. Thus, we extended ss GBLUP, incorporated genomic annotation information into the model, and evaluated them using a yellow-feathered chicken population as the examples. The chicken population consisted of 1 338 birds with 23 traits, where imputed WGS data including 5 127 612 single nucleotide polymorphisms(SNPs) are available for 895 birds. Considering different combinations of annotation information and models, original ss GBLUP, haplotype-based ss GHBLUP, and four extended ss GBLUP incorporating genomic annotation models were evaluated. Based on the genomic annotation(GRCg6a) of chickens, 3 155 524 and 94 837 SNPs were mapped to genic and exonic regions, respectively. Extended ss GBLUP using genic/exonic SNPs outperformed other models with respect to predictive ability in 15 out of 23 traits, and their advantages ranged from 2.5 to 6.1% compared with original ss GBLUP. In addition, to further enhance the performance of genomic prediction with imputed WGS data, we investigated the genotyping strategies of reference population on ss GBLUP in the chicken population. Comparing two strategies of individual selection for genotyping in the reference population, the strategy of evenly selection by family(SBF) performed slightly better than random selection in most situations. Overall, we extended genomic prediction models that can comprehensively utilize WGS data and genomic annotation information in the framework of ss GBLUP, and validated the idea that properly handling the genomic annotation information and WGS data increased the predictive ability of ss GBLUP. Moreover, while using WGS data, the genotyping strategy of maximizing the expected genetic relationship between the reference and candidate population could further improve the predictive ability of ss GBLUP. The results from this study shed light on the comprehensive usage of genomic annotation information in WGS-based single-step genomic prediction.

Applications of single-cell RNA sequencing in spermatogenesis and molecular evolution

Spermatogenic cell heterogeneity is determined by the complex process of spermatogenesis differentiation.However,effectively revealing the regulatory mechanisms underlying mammalian spermatogenic cell development and differentiation via traditional methods is difficult.Advances in technology have led to the emergence of many single-cell transcriptome sequencing protocols,which have partially addressed these challenges.In this review,we detail the principles of 10x Genomics technology and summarize the methods for downstream analysis of single-cell transcriptome sequencing data.Furthermore,we explore the role of single-cell transcriptome sequencing in revealing the heterogeneity of testicular ecological niche cells,delineating the establishment and disruption of testicular immune homeostasis during human spermatogenesis,investigating abnormal spermatogenesis in humans,and,ultimately,elucidating the molecular evolution of mammalian spermatogenesis.

Revolutionizing stem cell research:unbiased insights through single-cell sequencing

Stem cells have shown great application potential in wound repair,tissue regeneration,and disease treatment.Therefore,a full understanding of stem cells and their related regulatory mechanisms in disease treatment is conducive to improving the therapeutic effect of stem cells.However,thus far,there are still many unsolved mysteries in thefield of stem cells due to technical limitations,which hinder the in-depth exploration of stem cells and their wide clinical application.Single-cell sequencing(SCS)has provided very powerful and unbiased insights into cell gene expression profiles at the single-cell level,bringing exciting results to the stem cellfield.At present,SCS has been widely applied in thefield of stem cells,covering various aspects,including lineage tracing the development of stem cells,identifying new stem cell types,exploring cellular heterogeneity,and identifying internal functional subpopulations.In this paper,we focus on the latest research progress and discuss the application of SCS technology in stem cells.

A cloud model target damage effectiveness assessment algorithm based on spatio-temporal sequence finite multilayer fragments dispersion

To solve the problem of target damage assessment when fragments attack target under uncertain projectile and target intersection in an air defense intercept,this paper proposes a method for calculating target damage probability leveraging spatio-temporal finite multilayer fragments distribution and the target damage assessment algorithm based on cloud model theory.Drawing on the spatial dispersion characteristics of fragments of projectile proximity explosion,we divide into a finite number of fragments distribution planes based on the time series in space,set up a fragment layer dispersion model grounded in the time series and intersection criterion for determining the effective penetration of each layer of fragments into the target.Building on the precondition that the multilayer fragments of the time series effectively assail the target,we also establish the damage criterion of the perforation and penetration damage and deduce the damage probability calculation model.Taking the damage probability of the fragment layer in the spatio-temporal sequence to the target as the input state variable,we introduce cloud model theory to research the target damage assessment method.Combining the equivalent simulation experiment,the scientific and rational nature of the proposed method were validated through quantitative calculations and comparative analysis.

Dissection of triple-negative breast cancer microenvironment and identification of potential therapeutic drugs using single-cell RNA sequencing analysis

Breast cancer remains a leading cause of mortality in women worldwide.Triple-negative breast cancer(TNBC)is a particularly aggressive subtype characterized by rapid progression,poor prognosis,and lack of clear therapeutic targets.In the clinic,delineation of tumor heterogeneity and development of effective drugs continue to pose considerable challenges.Within the scope of our study,high heterogeneity inherent to breast cancer was uncovered based on the landscape constructed from both tumor and healthy breast tissue samples.Notably,TNBC exhibited significant specificity regarding cell proliferation,differentiation,and disease progression.Significant associations between tumor grade,prognosis,and TNBC oncogenes were established via pseudotime trajectory analysis.Consequently,we further performed comprehensive characterization of the TNBC microenvironment.A crucial epithelial subcluster,E8,was identified as highly malignant and strongly associated with tumor cell proliferation in TNBC.Additionally,epithelial-mesenchymal transition(EMT)-associated fibroblast and M2 macrophage subclusters exerted an influence on E8 through cellular interactions,contributing to tumor growth.Characteristic genes in these three cluster cells could therefore serve as potential therapeutic targets for TNBC.The collective findings provided valuable insights that assisted in the screening of a series of therapeutic drugs,such as pelitinib.We further confirmed the anti-cancer effect of pelitinib in an orthotopic 4T1 tumor-bearing mouse model.Overall,our study sheds light on the unique characteristics of TNBC at single-cell resolution and the crucial cell types associated with tumor cell proliferation that may serve as potent tools in the development of effective anti-cancer drugs.

The Upper Jurassic sediments,Marib-Shabwa Basin,Yemen:Lithofacies aspects and sequence stratigraphic analysis

The present study is devoted to understanding the evolution of the Upper Jurassic Sab'atayn Formation in the Marib-Shabwa Basin,Yemen,through a sequence stratigraphic analysis based on integrating datasets of sedimentology,seismic sections,and well logs.The Sab'atayn Formation(Tithonian age)is represented by a series of clastic and evaporites that were deposited under fluvio-deltaic to prodeltaic settings.It is divided into four members including Yah(at the base),upwards to Seen,Alif,and Safir at the top.Two third-order depositional sequences were determined for the Tithonian succession which were separated by three sequence boundaries.These sequences were classified into their systems tracts signifying several sedimentation patterns of progradational,aggradational,and retrogradational parasequence sets.The first depositional sequence corresponds to the early-middle Tithonian Yah and Seen units that can be classified into lowstand,transgressive,and highstand systems tracts.The second sequence comprises the late Tithonian Alif unit that can be subdivided into transgressive and highstand systems tracts.The sandy deposits of the Alif Member(highstand deposits)represent the most productive hydrocarbon reservoir in the basin.The Upper Jurassic sediments in the study area were resulted from a combination of eustatic and tectonic effects.

Single-cell sequencing technology in diabetic wound healing:New insights into the progenitors-based repair strategies

Diabetes mellitus(DM),an increasingly prevalent chronic metabolic disease,is characterised by prolonged hyperglycaemia,which leads to long-term health consequences.Although much effort has been put into understanding the pathogenesis of diabetic wounds,the underlying mechanisms remain unclear.The advent of single-cell RNA sequencing(scRNAseq)has revolutionised biological research by enabling the identification of novel cell types,the discovery of cellular markers,the analysis of gene expression patterns and the prediction of develop-mental trajectories.This powerful tool allows for an in-depth exploration of pathogenesis at the cellular and molecular levels.In this editorial,we focus on progenitor-based repair strategies for diabetic wound healing as revealed by scRNAseq and highlight the biological behaviour of various healing-related cells and the alteration of signalling pathways in the process of diabetic wound healing.ScRNAseq could not only deepen our understanding of the complex biology of diabetic wounds but also identify and validate new targets for inter-vention,offering hope for improved patient outcomes in the management of this challenging complication of DM.

Food nutrition and toxicology targeting on specific organs in the era of single-cell sequencing

Due to the complex natures of dietary food components,it is difficult to elucidate how the compounds affect host health.Dietary food often selectively presents its mechanism of action on different cell types,and participates in the modulation of targeted cells and their microenvironments within organs.However,the limitations of traditional in vitro assays or in vivo animal experiments cannot comprehensively examine cellular heterogeneity and the tissue-biased influences.Single-cell RNA sequencing(sc RNA-seq)has emerged as an indispensable methodology to decompose tissues into different cell types for the demonstration of transcriptional profiles of individual cells.Sc RNA-seq applications has been summarized on three typical organs(brain,liver,kidney),and two representative immune-and tumor related health problems.The everincreasing role of sc RNA-seq in dietary food research with further improvement can provide sub-cellular information and the coupling between other cellular modalities.In this review,we propose utilizing sc RNAseq to more effectively capture the subtle and complex effects of food chemicals,and how they may lead to health problems at single-cell resolution.This novel technique will be valuable to elucidate the underlying mechanism of both the health benefits of food nutrients and the detrimental consequences food toxicants at the cellular level.

Landscape of Sequence Variations in Homologous Copies of FAD2 and FAD3 in Rapeseed(Brassica napus L.)Germplasm with High/Low Linolenic Acid Trait

Genetic manipulation(either restraint or enhancement)of the biosynthesis pathway ofα-linolenic acid(ALA)in seed oil is an important goal in Brassica napus breeding.B.napus is a tetraploid plant whose genome often har-bors four and six homologous copies,respectively,of the two fatty acid desaturases FAD2 and FAD3,which con-trol the last two steps of ALA biosynthesis during seed oil accumulation.In this study,we compared their promoters,coding sequences,and expression levels in three high-ALA inbred lines 2006L,R8Q10,and YH25005,a low-ALA line A28,a low-ALA/high-oleic-acid accession SW,and the wildtype ZS11.The expression levels of most FAD2 and FAD3 homologs in the three high-ALA accessions were higher than those in ZS11 and much higher than those in A28 and SW.The three high-ALA accessions shared similar sequences with the pro-moters and CDSs of BnFAD3.C4 and BnFAD3.A3.In A28 and SW,substitution of three amino acid residues in BnFAD2.A5 and BnFAD2.C5,an absence of BnFAD2.C1 locus,and a 549 bp long deletion on the BnFAD3.A3 promoter were detected.The profile of BnFAD2 mutation in the two low-ALA accessions A28 and SW is different from that reported in previous studies.The mutations in BnFAD3 in the high-ALA accessions are reported for thefirst time.In identifying the sites of these mutations,we provide detailed information to aid the design of mole-cular markers for accelerated breeding schemes.

Adaptive Successive POI Recommendation via Trajectory Sequences Processing and Long Short-Term Preference Learning

Point-of-interest(POI)recommendations in location-based social networks(LBSNs)have developed rapidly by incorporating feature information and deep learning methods.However,most studies have failed to accurately reflect different users’preferences,in particular,the short-term preferences of inactive users.To better learn user preferences,in this study,we propose a long-short-term-preference-based adaptive successive POI recommendation(LSTP-ASR)method by combining trajectory sequence processing,long short-term preference learning,and spatiotemporal context.First,the check-in trajectory sequences are adaptively divided into recent and historical sequences according to a dynamic time window.Subsequently,an adaptive filling strategy is used to expand the recent check-in sequences of users with inactive check-in behavior using those of similar active users.We further propose an adaptive learning model to accurately extract long short-term preferences of users to establish an efficient successive POI recommendation system.A spatiotemporal-context-based recurrent neural network and temporal-context-based long short-term memory network are used to model the users’recent and historical checkin trajectory sequences,respectively.Extensive experiments on the Foursquare and Gowalla datasets reveal that the proposed method outperforms several other baseline methods in terms of three evaluation metrics.More specifically,LSTP-ASR outperforms the previously best baseline method(RTPM)with a 17.15%and 20.62%average improvement on the Foursquare and Gowalla datasets in terms of the Fβmetric,respectively.

Action of circulating and infiltrating B cells in the immune microenvironment of colorectal cancer by single-cell sequencing analysis

BACKGROUND The complexity of the immune microenvironment has an impact on the treatment of colorectal cancer(CRC),one of the most prevalent malignancies worldwide.In this study,multi-omics and single-cell sequencing techniques were used to investigate the mechanism of action of circulating and infiltrating B cells in CRC.By revealing the heterogeneity and functional differences of B cells in cancer immunity,we aim to deepen our understanding of immune regulation and provide a scientific basis for the development of more effective cancer treatment strategies.AIM To explore the role of circulating and infiltrating B cell subsets in the immune microenvironment of CRC,explore the potential driving mechanism of B cell development,analyze the interaction between B cells and other immune cells in the immune microenvironment and the functions of communication molecules,and search for possible regulatory pathways to promote the anti-tumor effects of B cells.METHODS A total of 69 paracancer(normal),tumor and peripheral blood samples were collected from 23 patients with CRC from The Cancer Genome Atlas database(https://portal.gdc.cancer.gov/).After the immune cells were sorted by multicolor flow cytometry,the single cell transcriptome and B cell receptor group library were sequenced using the 10X Genomics platform,and the data were analyzed using bioinformatics tools such as Seurat.The differences in the number and function of B cell infiltration between tumor and normal tissue,the interaction between B cell subsets and T cells and myeloid cell subsets,and the transcription factor regulatory network of B cell subsets were explored and analyzed.RESULTS Compared with normal tissue,the infiltrating number of CD20+B cell subsets in tumor tissue increased significantly.Among them,germinal center B cells(GCB)played the most prominent role,with positive clone expansion and heavy chain mutation level increasing,and the trend of differentiation into memory B cells increased.However,the number of plasma cells in the tumor microenvironment decreased significantly,and the plasma cells secreting IgA antibodies decreased most obviously.In addition,compared with the immune microenvironment of normal tissues,GCB cells in tumor tissues became more closely connected with other immune cells such as T cells,and communication molecules that positively regulate immune function were significantly enriched.CONCLUSION The role of GCB in CRC tumor microenvironment is greatly enhanced,and its affinity to tumor antigen is enhanced by its significantly increased heavy chain mutation level.Meanwhile,GCB has enhanced its association with immune cells in the microenvironment,which plays a positive anti-tumor effect.