I am Peisen Sun, a Ph.D. candidate in the Department of Control Science and Engineering at Xi'an Jiaotong University. My research focuses on bioinformatics, with an emphasis on quantifying the functional impact of somatic variants at single-cell resolution and dissecting tumor tissues through spatial multi-omics.
I employ machine learning and artificial intelligence approaches to address complex genomics problems, leveraging state-of-the-art single-cell and spatial omics technologies. My work aims to elucidate the molecular mechanisms underlying biological systems.
My research contributions have been published as first author (including co-first author) in venues such as Cell Genomics, Nature Communications, and bioRxiv.
") does not match the recommended repository name for your site ("").
", so that your site can be accessed directly at "http://".
However, if the current repository name is intended, you can ignore this message by removing "{% include widgets/debug_repo_name.html %}" in index.html.
",
which does not match the baseurl ("") configured in _config.yml.
baseurl in _config.yml to "".

Peisen Sun (孙培森), Peng Jia (贾鹏), Chuyuan Lou, Xiaona Li, Miaojing Li, Ruifen Zhang, Hui Lin, Xiujuan Li, Yunsong Zhou, Shuzheng Cai, Zexin Zhang, Pengyu Zhang (张鹏宇), Tun Xu (徐暾), Songbo Wang (王松渤), Yudong Wang, Mingxuan Li (黎明轩), Yuyan Gong, Guoqin Sang, Jie Yin, Jizhe Li, Qi Xie, Zongsheng Yue, Dan Xu, Xiaofei Yang (杨晓飞), Stephen J. Bush, Songmei Geng (耿松梅), Deyu Meng (孟德宇), Xuexin Li (李薛鑫)#, Kai Ye (叶凯)# (# corresponding author)
bioRxiv 2026
Somatic SNVs are commonly prioritized by recurrence or sequence annotation, leaving the molecular consequences of private variants unresolved. We developed PrismSNV and Ligo-Seq, a single-cell full-length transcriptome sequencing platform, to quantify transcriptomic displacement associated with within-cell SNV-state toggling. Across 26 samples from three cancer types, more than 60% of analyzed SNVs showed significant cell-type-specific displacement with over 30% variation in magnitude. Sequence-based deleteriousness was poorly aligned with cellular displacement, and fewer than 3% of top-displacement SNVs across three tumor types were catalogued in ClinVar. Displacement was strongest in transcriptionally plastic, less differentiated states. Analysis of 655,349 CD4+ T-cell transcriptomes further showed that SNV-gene associations changed in magnitude and sometimes direction during activation. Lack of recurrence or clinical annotation does not imply transcriptomic inertness; variant impact is shaped by cellular context rather than sequence alone, establishing the variant-cell pair as a critical unit for interpreting somatic variation.
Peisen Sun (孙培森), Peng Jia (贾鹏), Chuyuan Lou, Xiaona Li, Miaojing Li, Ruifen Zhang, Hui Lin, Xiujuan Li, Yunsong Zhou, Shuzheng Cai, Zexin Zhang, Pengyu Zhang (张鹏宇), Tun Xu (徐暾), Songbo Wang (王松渤), Yudong Wang, Mingxuan Li (黎明轩), Yuyan Gong, Guoqin Sang, Jie Yin, Jizhe Li, Qi Xie, Zongsheng Yue, Dan Xu, Xiaofei Yang (杨晓飞), Stephen J. Bush, Songmei Geng (耿松梅), Deyu Meng (孟德宇), Xuexin Li (李薛鑫)#, Kai Ye (叶凯)# (# corresponding author)
bioRxiv 2026
Somatic SNVs are commonly prioritized by recurrence or sequence annotation, leaving the molecular consequences of private variants unresolved. We developed PrismSNV and Ligo-Seq, a single-cell full-length transcriptome sequencing platform, to quantify transcriptomic displacement associated with within-cell SNV-state toggling. Across 26 samples from three cancer types, more than 60% of analyzed SNVs showed significant cell-type-specific displacement with over 30% variation in magnitude. Sequence-based deleteriousness was poorly aligned with cellular displacement, and fewer than 3% of top-displacement SNVs across three tumor types were catalogued in ClinVar. Displacement was strongest in transcriptionally plastic, less differentiated states. Analysis of 655,349 CD4+ T-cell transcriptomes further showed that SNV-gene associations changed in magnitude and sometimes direction during activation. Lack of recurrence or clinical annotation does not imply transcriptomic inertness; variant impact is shaped by cellular context rather than sequence alone, establishing the variant-cell pair as a critical unit for interpreting somatic variation.

Mingxuan Li (黎明轩)*, Peisen Sun (孙培森)*, Yisi Luo (罗倚斯), Guancheng Zhou (周冠程), Xiaofei Yang (杨晓飞), Deyu Meng (孟德宇), Kai Ye (叶凯) (* equal contribution)
Nature Communications 2026
Integrating spatial multi-omics data presents significant challenges, particularly in uncovering the spatial patterns of cells and deciphering the real regulatory mechanisms among various omics. These insights are critical for harnessing the full potential of each modality while minimizing the impact of biotechnological biases that will lead to unstable results. Here, we introduce SpatialCOC, a framework that treats spatial information as prior knowledge to learn omics-specific spatial distributions, then discovering nonlinear correlations among modalities. The effectiveness and robustness of SpatialCOC are validated using real-world datasets, encompassing diverse tissue sections analyzed with multiple experimental techniques. Compared to existing methods, SpatialCOC excels in identifying region-specific continuous spatial domains and maintains batch-consistency across trajectory inferences. By providing a novel perspective on the interplay between spatial information and multi-omics modalities, SpatialCOC offers a flexible approach that can accommodate modality data of arbitrary dimensions.
Mingxuan Li (黎明轩)*, Peisen Sun (孙培森)*, Yisi Luo (罗倚斯), Guancheng Zhou (周冠程), Xiaofei Yang (杨晓飞), Deyu Meng (孟德宇), Kai Ye (叶凯) (* equal contribution)
Nature Communications 2026
Integrating spatial multi-omics data presents significant challenges, particularly in uncovering the spatial patterns of cells and deciphering the real regulatory mechanisms among various omics. These insights are critical for harnessing the full potential of each modality while minimizing the impact of biotechnological biases that will lead to unstable results. Here, we introduce SpatialCOC, a framework that treats spatial information as prior knowledge to learn omics-specific spatial distributions, then discovering nonlinear correlations among modalities. The effectiveness and robustness of SpatialCOC are validated using real-world datasets, encompassing diverse tissue sections analyzed with multiple experimental techniques. Compared to existing methods, SpatialCOC excels in identifying region-specific continuous spatial domains and maintains batch-consistency across trajectory inferences. By providing a novel perspective on the interplay between spatial information and multi-omics modalities, SpatialCOC offers a flexible approach that can accommodate modality data of arbitrary dimensions.

Peisen Sun (孙培森), Stephen J. Bush, Songbo Wang (王松渤), Peng Jia (贾鹏), Mingxuan Li (黎明轩), Tun Xu (徐暾), Pengyu Zhang (张鹏宇), Xiaofei Yang (杨晓飞), Chengyao Wang (王澄瑶), Linfeng Xu (许林峰), Tingjie Wang (王庭杰), Kai Ye (叶凯)
Cell Genomics 2025 Featured content
Analyzing spatial transcriptomics data from tumor tissues poses several challenges beyond those of healthy samples, including unclear boundaries between different regions, uneven cell densities, and relatively higher cellular heterogeneity. Collectively, these bias the background against which spatially variable genes are identified, which can result in misidentification of spatial structures and hinder potential insight into complex pathologies. To overcome this problem, STMiner leverages 2D Gaussian mixture models and optimal transport theory to directly characterize the spatial distribution of genes rather than the capture locations of the cells expressing them (spots). By effectively mitigating the impacts of both background bias and data sparsity, STMiner reveals key gene sets and spatial structures overlooked by spot-based analytic tools, facilitating novel biological discoveries. The core concept of directly analyzing overall gene expression patterns also allows for a broader application beyond spatial transcriptomics, positioning STMiner for continuous expansion as spatial omics technologies evolve.
Peisen Sun (孙培森), Stephen J. Bush, Songbo Wang (王松渤), Peng Jia (贾鹏), Mingxuan Li (黎明轩), Tun Xu (徐暾), Pengyu Zhang (张鹏宇), Xiaofei Yang (杨晓飞), Chengyao Wang (王澄瑶), Linfeng Xu (许林峰), Tingjie Wang (王庭杰), Kai Ye (叶凯)
Cell Genomics 2025 Featured content
Analyzing spatial transcriptomics data from tumor tissues poses several challenges beyond those of healthy samples, including unclear boundaries between different regions, uneven cell densities, and relatively higher cellular heterogeneity. Collectively, these bias the background against which spatially variable genes are identified, which can result in misidentification of spatial structures and hinder potential insight into complex pathologies. To overcome this problem, STMiner leverages 2D Gaussian mixture models and optimal transport theory to directly characterize the spatial distribution of genes rather than the capture locations of the cells expressing them (spots). By effectively mitigating the impacts of both background bias and data sparsity, STMiner reveals key gene sets and spatial structures overlooked by spot-based analytic tools, facilitating novel biological discoveries. The core concept of directly analyzing overall gene expression patterns also allows for a broader application beyond spatial transcriptomics, positioning STMiner for continuous expansion as spatial omics technologies evolve.

Peisen Sun, Guanglin Li
Frontiers in Genetics 2019
Circular RNAs (circRNAs), which play vital roles in many regulatory pathways, are widespread in many species. Although many circRNAs have been discovered in plants and animals, the functions of these RNAs have not been fully investigated. In addition to the function of circRNAs as microRNA (miRNA) decoys, the translation potential of circRNAs is important for the study of their functions; yet, few tools are available to identify their translation potential. With the development of high-throughput sequencing technology and the emergence of ribosome profiling technology, it is possible to identify the coding ability of circRNAs with high sensitivity. To evaluate the coding ability of circRNAs, we first developed the CircCode tool and then used CircCode to investigate the translation potential of circRNAs from humans and Arabidopsis thaliana. Based on the ribosome profile databases downloaded from NCBI, we found 3,610 and 1,569 translated circRNAs in humans and A. thaliana, respectively. Finally, we tested the performance of CircCode and found a low false discovery rate and high sensitivity for identifying circRNA coding ability. CircCode, a Python 3–based framework for identifying the coding ability of circRNAs, is also a simple and powerful command line-based tool. To investigate the translation potential of circRNAs, the user can simply fill in the given configuration file and run the Python 3 scripts.
Peisen Sun, Guanglin Li
Frontiers in Genetics 2019
Circular RNAs (circRNAs), which play vital roles in many regulatory pathways, are widespread in many species. Although many circRNAs have been discovered in plants and animals, the functions of these RNAs have not been fully investigated. In addition to the function of circRNAs as microRNA (miRNA) decoys, the translation potential of circRNAs is important for the study of their functions; yet, few tools are available to identify their translation potential. With the development of high-throughput sequencing technology and the emergence of ribosome profiling technology, it is possible to identify the coding ability of circRNAs with high sensitivity. To evaluate the coding ability of circRNAs, we first developed the CircCode tool and then used CircCode to investigate the translation potential of circRNAs from humans and Arabidopsis thaliana. Based on the ribosome profile databases downloaded from NCBI, we found 3,610 and 1,569 translated circRNAs in humans and A. thaliana, respectively. Finally, we tested the performance of CircCode and found a low false discovery rate and high sensitivity for identifying circRNA coding ability. CircCode, a Python 3–based framework for identifying the coding ability of circRNAs, is also a simple and powerful command line-based tool. To investigate the translation potential of circRNAs, the user can simply fill in the given configuration file and run the Python 3 scripts.