SpaNCMG: improving spatial domains identification of spatial transcriptomics using neighborhood-complementary mixed-view graph convolutional network.

The advancement of spatial transcriptomics (ST) technology contributes to a more profound comprehension of the spatial properties of gene expression within tissues. However, due to challenges of high dimensionality, pronounced noise and dynamic limitations in ST data, the integration of gene expression and spatial information to accurately identify spatial domains remains challenging. This paper proposes a SpaNCMG algorithm for the purpose of achieving precise spatial domain description and localization based on a neighborhood-complementary mixed-view graph convolutional network. The algorithm enables better adaptation to ST data at different resolutions by integrating the local information from KNN and the global structure from r-radius into a complementary neighborhood graph. It also introduces an attention mechanism to achieve adaptive fusion of different reconstructed expressions, and utilizes KPCA method for dimensionality reduction. The application of SpaNCMG on five datasets from four sequencing platforms demonstrates superior performance to eight existing advanced methods. Specifically, the algorithm achieved highest ARI accuracies of 0.63 and 0.52 on the datasets of the human dorsolateral prefrontal cortex and mouse somatosensory cortex, respectively. It accurately identified the spatial locations of marker genes in the mouse olfactory bulb tissue and inferred the biological functions of different regions. When handling larger datasets such as mouse embryos, the SpaNCMG not only identified the main tissue structures but also explored unlabeled domains. Overall, the good generalization ability and scalability of SpaNCMG make it an outstanding tool for understanding tissue structure and disease mechanisms. Our codes are available at https://github.com/ZhihaoSi/SpaNCMG.

Integrating somatic mutation profiles with structural deep clustering network for metabolic stratification in pancreatic cancer: a comprehensive analysis of prognostic and genomic landscapes.

Pancreatic cancer is a globally recognized highly aggressive malignancy, posing a significant threat to human health and characterized by pronounced heterogeneity. In recent years, researchers have uncovered that the development and progression of cancer are often attributed to the accumulation of somatic mutations within cells. However, cancer somatic mutation data exhibit characteristics such as high dimensionality and sparsity, which pose new challenges in utilizing these data effectively. In this study, we propagated the discrete somatic mutation data of pancreatic cancer through a network propagation model based on protein-protein interaction networks. This resulted in smoothed somatic mutation profile data that incorporate protein network information. Based on this smoothed mutation profile data, we obtained the activity levels of different metabolic pathways in pancreatic cancer patients. Subsequently, using the activity levels of various metabolic pathways in cancer patients, we employed a deep clustering algorithm to establish biologically and clinically relevant metabolic subtypes of pancreatic cancer. Our study holds scientific significance in classifying pancreatic cancer based on somatic mutation data and may provide a crucial theoretical basis for the diagnosis and immunotherapy of pancreatic cancer patients.

A computational framework of routine test data for the cost-effective chronic disease prediction.

Chronic diseases, because of insidious onset and long latent period, have become the major global disease burden. However, the current chronic disease diagnosis methods based on genetic markers or imaging analysis are challenging to promote completely due to high costs and cannot reach universality and popularization. This study analyzed massive data from routine blood and biochemical test of 32 448 patients and developed a novel framework for cost-effective chronic disease prediction with high accuracy (AUC 87.32%). Based on the best-performing XGBoost algorithm, 20 classification models were further constructed for 17 types of chronic diseases, including 9 types of cancers, 5 types of cardiovascular diseases and 3 types of mental illness. The highest accuracy of the model was 90.13% for cardia cancer, and the lowest was 76.38% for rectal cancer. The model interpretation with the SHAP algorithm showed that CREA, R-CV, GLU and NEUT% might be important indices to identify the most chronic diseases. PDW and R-CV are also discovered to be crucial indices in classifying the three types of chronic diseases (cardiovascular disease, cancer and mental illness). In addition, R-CV has a higher specificity for cancer, ALP for cardiovascular disease and GLU for mental illness. The association between chronic diseases was further revealed. At last, we build a user-friendly explainable machine-learning-based clinical decision support system (DisPioneer: http://bioinfor.imu.edu.cn/dispioneer) to assist in predicting, classifying and treating chronic diseases. This cost-effective work with simple blood tests will benefit more people and motivate clinical implementation and further investigation of chronic diseases prevention and surveillance program.

EmAtlas: a comprehensive atlas for exploring spatiotemporal activation in mammalian embryogenesis.

The emerging importance of embryonic development research rapidly increases the volume for a professional resource related to multi-omics data. However, the lack of global embryogenesis repository and systematic analysis tools limits the preceding in stem cell research, human congenital diseases and assisted reproduction. Here, we developed the EmAtlas, which collects the most comprehensive multi-omics data and provides multi-scale tools to explore spatiotemporal activation during mammalian embryogenesis. EmAtlas contains data on multiple types of gene expression, chromatin accessibility, DNA methylation, nucleosome occupancy, histone modifications, and transcription factors, which displays the complete spatiotemporal landscape in mouse and human across several time points, involving gametogenesis, preimplantation, even fetus and neonate, and each tissue involves various cell types. To characterize signatures involved in the tissue, cell, genome, gene and protein levels during mammalian embryogenesis, analysis tools on these five scales were developed. Additionally, we proposed EmRanger to deliver extensive development-related biological background annotations. Users can utilize these tools to analyze, browse, visualize, and download data owing to the user-friendly interface. EmAtlas is freely accessible at http://bioinfor.imu.edu.cn/ematlas.

iProbiotics: a machine learning platform for rapid identification of probiotic properties from whole-genome primary sequences.

Lactic acid bacteria consortia are commonly present in food, and some of these bacteria possess probiotic properties. However, discovery and experimental validation of probiotics require extensive time and effort. Therefore, it is of great interest to develop effective screening methods for identifying probiotics. Advances in sequencing technology have generated massive genomic data, enabling us to create a machine learning-based platform for such purpose in this work. This study first selected a comprehensive probiotics genome dataset from the probiotic database (PROBIO) and literature surveys. Then, k-mer (from 2 to 8) compositional analysis was performed, revealing diverse oligonucleotide composition in strain genomes and apparently more probiotic (P-) features in probiotic genomes than non-probiotic genomes. To reduce noise and improve computational efficiency, 87 376 k-mers were refined by an incremental feature selection (IFS) method, and the model achieved the maximum accuracy level at 184 core features, with a high prediction accuracy (97.77%) and area under the curve (98.00%). Functional genomic analysis using annotations from gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Rapid Annotation using Subsystem Technology (RAST) databases, as well as analysis of genes associated with host gastrointestinal survival/settlement, carbohydrate utilization, drug resistance and virulence factors, revealed that the distribution of P-features was biased toward genes/pathways related to probiotic function. Our results suggest that the role of probiotics is not determined by a single gene, but by a combination of k-mer genomic components, providing new insights into the identification and underlying mechanisms of probiotics. This work created a novel and free online bioinformatic tool, iProbiotics, which would facilitate rapid screening for probiotics.

Designing a structure-function alphabet of helix based on reduced amino acid clusters.

Several simple secondary structures could form complex and diverse functional proteins, meaning that secondary structures may contain a lot of hidden information and are arranged according to certain principles, to carry enough information of functional specificity and diversity. However, these inner information and principles have not been understood systematically. In our study, we designed a structure-function alphabet of helix based on reduced amino acid clusters to describe the typical features of helices and delve into the information. Firstly, we selected 480 typical helices from membrane proteins, zymoproteins, transcription factors, and other proteins to define and calculate the interval range, and the helices are classified in terms of hydrophilicity, charge and length: (1) hydrophobic helix (≤43%), amphiphilic helix (43%∼71%), and hydrophilic helix (≥71%). (2) positive helix, negative helix, electrically neutral helix and uncharged helix. (3) short helix (≤8 aa), medium-length helix (9-28 aa), and long helix (≥29 aa). Then, we designed an alphabet containing 36 triplet codes according to the above classification, so that the main features of each helix can be represented by only three letters. This alphabet not only preliminarily defined the helix characteristics, but also greatly reduced the informational dimension of protein structure. Finally, we present an application example to demonstrate the value of the structure-function alphabet in protein functional determination and differentiation.

RaacFold: a webserver for 3D visualization and analysis of protein structure by using reduced amino acid alphabets.

Protein structure exhibits greater complexity and diversity than DNA structure, and usually affects the interpretation of the function, interactions and biological annotations. Reduced amino acid alphabets (Raaa) exhibit a powerful ability to decrease protein complexity and identify functional conserved regions, which motivated us to create RaacFold. The RaacFold provides 687 reduced amino acid clusters (Raac) based on 58 reduction methods and offers three analysis tools: Protein Analysis, Align Analysis, and Multi Analysis. The Protein Analysis and Align Analysis provide reduced representations of sequence-structure according to physicochemical similarities and computational biology strategies. With the simplified representations, the protein structure can be viewed more concise and clearer to capture biological insight than the unreduced structure. Thus, the design of artificial protein will be more convenient, and redundant interference is avoided. In addition, Multi Analysis allows users to explore biophysical variation and conservation in the evolution of protein structure and function. This supplies important information for the identification and exploration of the nonhomologous functions of paralogs. Simultaneously, RaacFold provides powerful 2D and 3D rendering performance with advanced parameters for sequences, structures, and related annotations. RaacFold is freely available at http://bioinfor.imu.edu.cn/raacfold.

Identification of Key lncRNAs Associated with Immune Infiltration and Prognosis in Gastric Cancer.

Long non-coding RNAs (lncRNAs), as promising novel biomarkers for cancer treatment and prognosis, can function as tumor suppressors and oncogenes in the occurrence and development of many types of cancer, including gastric cancer (GC). However, little is known about the complex regulatory system of lncRNAs in GC. In this study, we systematically analyzed lncRNA and miRNA transcriptomic profiles of GC based on bioinformatics methods and experimental validation. An lncRNA-miRNA interaction network related to GC was constructed, and the nine crucial lncRNAs were identified. These 9 lncRNAs were found to be associated with the prognosis of GC patients by Cox proportional hazards regression analysis. Among them, the expression of lncRNA SNHG14 can affect the survival of GC patients as a potential prognostic marker. Moreover, it was shown that SNHG14 was involved in immune-related pathways and significantly correlated with immune cell infiltration in GC. Meanwhile, we found that SNHG14 affected immune function in many cancers, such as breast cancer and esophageal carcinoma. Such information revealed that SNHG14 may serve as a potential target for cancer immunotherapy. As well, our study could provide practical and theoretical guiding significance for clinical application of non-coding RNAs.

Integrated analysis of the ubiquitination mechanism reveals the specific signatures of tissue and cancer.

BACKGROUND: Ubiquitination controls almost all cellular processes. The dysregulation of ubiquitination signals is closely associated with the initiation and progression of multiple diseases. However, there is little comprehensive research on the interaction and potential function of ubiquitination regulators (UBRs) in spermatogenesis and cancer. METHODS: We systematically characterized the mRNA and protein expression of UBRs across tissues and further evaluated their roles in testicular development and spermatogenesis. Subsequently, we explored the genetic alterations, expression perturbations, cancer hallmark-related pathways, and clinical relevance of UBRs in pan-cancer. RESULTS: This work reveals heterogeneity in the expression patterns of UBRs across tissues, and the expression pattern in testis is the most distinct. UBRs are dynamically expressed during testis development, which are critical for normal spermatogenesis. Furthermore, UBRs have widespread genetic alterations and expression perturbations in pan-cancer. The expression of 79 UBRs was identified to be closely correlated with the activity of 32 cancer hallmark-related pathways, and ten hub genes were screened for further clinical relevance analysis by a network-based method. More than 90% of UBRs can affect the survival of cancer patients, and hub genes have an excellent prognostic classification for specific cancer types. CONCLUSIONS: Our study provides a comprehensive analysis of UBRs in spermatogenesis and pan-cancer, which can build a foundation for understanding male infertility and developing cancer drugs in the aspect of ubiquitination.

Dppa2/4 as a trigger of signaling pathways to promote zygote genome activation by binding to CG-rich region.

Developmental pluripotency-associated 2 (Dppa2) and developmental pluripotency-associated 4 (Dppa4) as positive drivers were helpful for transcriptional regulation of zygotic genome activation (ZGA). Here, we systematically assessed the cooperative interplay of Dppa2 and Dppa4 in regulating cell pluripotency and found that simultaneous overexpression of Dppa2/4 can make induced pluripotent stem cells closer to embryonic stem cells (ESCs). Compared with other pluripotency transcription factors, Dppa2/4 can regulate majorities of signaling pathways by binding on CG-rich region of proximal promoter (0-500 bp), of which 85% and 77% signaling pathways were significantly activated by Dppa2 and Dppa4, respectively. Notably, Dppa2/4 also can dramatically trigger the decisive signaling pathways for facilitating ZGA, including Hippo, MAPK and TGF-beta signaling pathways and so on. At last, we found alkaline phosphatase, placental-like 2 (Alppl2) was completely silenced when Dppa2 and 4 single- or double-knockout in ESC, which is consistent with Dux. Moreover, Alppl2 was significantly activated in mouse 2-cell embryos and 4-8 cells stage of human embryos, further predicted that Alppl2 was directly regulated by Dppa2/4 as a ZGA candidate driver to facilitate pre-embryonic development.

RaacLogo: a new sequence logo generator by using reduced amino acid clusters.

Sequence logos give a fast and concise display in visualizing consensus sequence. Protein exhibits greater complexity and diversity than DNA, which usually affects the graphical representation of the logo. Reduced amino acids perform powerful ability for simplifying complexity of sequence alignment, which motivated us to establish RaacLogo. As a new sequence logo generator by using reduced amino acid alphabets, RaacLogo can easily generate many different simplified logos tailored to users by selecting various reduced amino acid alphabets that consisted of more than 40 clustering algorithms. This current web server provides 74 types of reduced amino acid alphabet, which were manually extracted to generate 673 reduced amino acid clusters (RAACs) for dealing with protein alignment. A two-dimensional selector was proposed for easily selecting desired RAACs with underlying biology knowledge. It is anticipated that the RaacLogo web server will play more high-potential roles for protein sequence alignment, topological estimation and protein design experiments. RaacLogo is freely available at http://bioinfor.imu.edu.cn/raaclogo.

HelPredictor models single-cell transcriptome to predict human embryo lineage allocation.

The in-depth understanding of cellular fate decision of human preimplantation embryos has prompted investigations on how changes in lineage allocation, which is far from trivial and remains a time-consuming task by experimental methods. It is desirable to develop a novel effective bioinformatics strategy to consider transitions of coordinated embryo lineage allocation and stage-specific patterns. There are rapidly growing applications of machine learning models to interpret complex datasets for identifying candidate development-related factors and lineage-determining molecular events. Here we developed the first machine learning platform, HelPredictor, that integrates three feature selection methods, namely, principal components analysis, F-score algorithm and squared coefficient of variation, and four classical machine learning classifiers that different combinations of methods and classifiers have independent outputs by increment feature selection method. With application to single-cell sequencing data of human embryo, HelPredictor not only achieved 94.9% and 90.9% respectively with cross-validation and independent test, but also fast classified different embryonic lineages and their development trajectories using less HelPredictor-predicted factors. The above-mentioned candidate lineage-specific genes were discussed in detail and were clustered for exploring transitions of embryonic heterogeneity. Our tool can fast and efficiently reveal potential lineage-specific and stage-specific biomarkers and provide insights into how advanced computational tools contribute to development research. The source code is available at https://github.com/liameihao/HelPredictor.

Immune cell infiltration-based signature for prognosis and immunogenomic analysis in breast cancer.

Breast cancer is one of the most human malignant diseases and the leading cause of cancer-related death in the world. However, the prognostic and therapeutic benefits of breast cancer patients cannot be predicted accurately by the current stratifying system. In this study, an immune-related prognostic score was established in 22 breast cancer cohorts with a total of 6415 samples. An extensive immunogenomic analysis was conducted to explore the relationships between immune score, prognostic significance, infiltrating immune cells, cancer genotypes and potential immune escape mechanisms. Our analysis revealed that this immune score was a promising biomarker for estimating overall survival in breast cancer. This immune score was associated with important immunophenotypic factors, such as immune escape and mutation load. Further analysis revealed that patients with high immune scores exhibited therapeutic benefits from chemotherapy and immunotherapy. Based on these results, we can conclude that this immune score may be a useful tool for overall survival prediction and treatment guidance for patients with breast cancer.

Modular arrangements of sequence motifs determine the functional diversity of KDM proteins.

Histone lysine demethylases (KDMs) play a vital role in regulating chromatin dynamics and transcription. KDM proteins are given modular activities by its sequence motifs with obvious roles division, which endow the complex and diverse functions. In our review, according to functional features, we classify sequence motifs into four classes: catalytic motifs, targeting motifs, regulatory motifs and potential motifs. JmjC, as the main catalytic motif, combines to Fe2+ and α-ketoglutarate by residues H-D/E-H and S-N-N/Y-K-N/Y-T/S. Targeting motifs make catalytic motifs recognize specific methylated lysines, such as PHD that helps KDM5 to demethylate H3K4me3. Regulatory motifs consist of a functional network. For example, NLS, Ser-rich, TPR and JmjN motifs regulate the nuclear localization. And interactions through the CW-type-C4H2C2-SWIRM are necessary to the demethylase activity of KDM1B. Additionally, many conservative domains that have potential functions but no deep exploration are reviewed for the first time. These conservative domains are usually amino acid-rich regions, which have great research value. The arrangements of four types of sequence motifs generate that KDM proteins diversify toward modular activities and biological functions. Finally, we draw a blueprint of functional mechanisms to discuss the modular activity of KDMs.

Clinical significance and immunogenomic landscape analyses of the immune cell signature based prognostic model for patients with breast cancer.

Breast cancer is one of the most common types of cancers and the leading cause of death from malignancy among women worldwide. Tumor-infiltrating lymphocytes are a source of important prognostic biomarkers for breast cancer patients. In this study, based on the tumor-infiltrating lymphocytes in the tumor immune microenvironment, a risk score prognostic model was developed in the training cohort for risk stratification and prognosis prediction in breast cancer patients. The prognostic value of this risk score prognostic model was also verified in the two testing cohorts and the TCGA pan cancer cohort. Nomograms were also established in the training and testing cohorts to validate the clinical use of this model. Relationships between the risk score, intrinsic molecular subtypes, immune checkpoints, tumor-infiltrating immune cell abundances and the response to chemotherapy and immunotherapy were also evaluated. Based on these results, we can conclude that this risk score model could serve as a robust prognostic biomarker, provide therapeutic benefits for the development of novel chemotherapy and immunotherapy, and may be helpful for clinical decision making in breast cancer patients.

Characterized the diversity of ABCB1 subtypes in immunogenomic landscape for predicting the drug response in breast cancer.

ABCB1 is an important gene that closely related to analgesic tolerance to opioids, and plays an important role in their postoperative treatment. Recent studies have demonstrated that ABCB1 genotype is significantly associated with the chemico-resistance and chemical sensitivity in breast cancer patients. So, it is become very important to investigate the important role of ABCB1 for predicting drug response in breast cancer patients. In this study, by conducting the Cox proportional hazards regression analysis in breast cancer patients, significant differences were found in prognosis between the ABCB1 high- and low-expression subtypes. Meanwhile, by using immune infiltration profiles as well as transcriptomics datasets, the ABCB1 high subtype was found to be significantly enriched in many immune-related KEGG pathways and biological processes, and was characterized by the high infiltration levels of immune cell types. Furthermore, bioinformatics inference revealed that the ABCB1 subtypes were associated with the therapeutic effect of immunotherapy, which would be important for patient prognosis. In conclusion, these findings may provide useful help for recognizing the diversity between ABCB1 subtypes in tumor immune microenvironment, and may unravel prognosis outcomes and immunotherapy utility for ABCB1 in breast cancer.

Reprogramming barriers in bovine cells nuclear transfer revealed by single-cell RNA-seq analysis.

Many progresses have recently been achieved in animal somatic cell nuclear transfer (SCNT). However, embryos derived from SCNT rarely result in live births. Single-cell RNA sequencing (scRNA-seq) can be used to investigate the development details of SCNT embryos. Here, bovine fibroblasts and three factors bovine iPSCs (3F biPSCs) were used as donors for bovine nuclear transfer, and the single blastomere transcriptome was analysed by scRNA-seq. Compared to in vitro fertilization (IVF) embryos, SCNT embryos exhibited many defects. Abnormally expressed genes were found at each stage of embryos, which enriched in metabolism, and epigenetic modification. The DEGs of the adjacent stage in SCNT embryos did not follow the temporal expression pattern similar to that of IVF embryos. Particularly, SCNT 8-cell stage embryos showed failures in some gene activation, including ZSCAN4, and defects in protein association networks which cored as POLR2K, GRO1, and ANKRD1. Some important signalling pathways also showed incomplete activation at SCNT zygote to morula stage. Interestingly, 3F biPSCNT embryos exhibited more dysregulated genes than SCNT embryos at zygote and 2-cell stage, including genes in KDM family. Pseudotime analysis of 3F biPSCNT embryos showed the different developmental fate from SCNT and IVF embryos. These findings suggested partial reprogrammed 3F biPS cells as donors for bovine nuclear transfer hindered the reprogramming of nuclear transfer embryos. Our studies revealed the abnormal gene expression and pathway activation of SCNT embryos, which could increase our understanding of the development of SCNT embryos and give hints to improve the efficiency of nuclear transfer.

eHSCPr discriminating the cell identity involved in endothelial to hematopoietic transition.

MOTIVATION: Hematopoietic stem cells (HSCs) give rise to all blood cells and play a vital role throughout the whole lifespan through their pluripotency and self-renewal properties. Accurately identifying the stages of early HSCs is extremely important, as it may open up new prospects for extracorporeal blood research. Existing experimental techniques for identifying the early stages of HSCs development are time-consuming and expensive. Machine learning has shown its excellence in massive single-cell data processing and it is desirable to develop related computational models as good complements to experimental techniques. RESULTS: In this study, we presented a novel predictor called eHSCPr specifically for predicting the early stages of HSCs development. To reveal the distinct genes at each developmental stage of HSCs, we compared F-score with three state-of-art differential gene selection methods (limma, DESeq2, edgeR) and evaluated their performance. F-score captured the more critical surface markers of endothelial cells and hematopoietic cells, and the area under receiver operating characteristic curve (ROC) value was 0.987. Based on SVM, the 10-fold cross-validation accuracy of eHSCpr in the independent dataset and the training dataset reached 94.84% and 94.19%, respectively. Importantly, we performed transcription analysis on the F-score gene set, which indeed further enriched the signal markers of HSCs development stages. eHSCPr can be a powerful tool for predicting early stages of HSCs development, facilitating hypothesis-driven experimental design and providing crucial clues for the in vitro blood regeneration studies. AVAILABILITY AND IMPLEMENTATION: http://bioinfor.imu.edu.cn/ehscpr. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Multi-substrate selectivity based on key loops and non-homologous domains: new insight into ALKBH family.

AlkB homologs (ALKBH) are a family of specific demethylases that depend on Fe2+ and α-ketoglutarate to catalyze demethylation on different substrates, including ssDNA, dsDNA, mRNA, tRNA, and proteins. Previous studies have made great progress in determining the sequence, structure, and molecular mechanism of the ALKBH family. Here, we first review the multi-substrate selectivity of the ALKBH demethylase family from the perspective of sequence and structural evolution. The construction of the phylogenetic tree and the comparison of key loops and non-homologous domains indicate that the paralogs with close evolutionary relationship have similar domain compositions. The structures show that the lack and variations of four key loops change the shape of clefts to cause the differences in substrate affinity, and non-homologous domains may be related to the compatibility of multiple substrates. We anticipate that the new insights into selectivity determinants of the ALKBH family are useful for understanding the demethylation mechanisms.

The Cumulative Formation of R-loop Interacts with Histone Modifications to Shape Cell Reprogramming.

R-loop, a three-stranded RNA/DNA structure, plays important roles in modulating genome stability and gene expression, but the molecular mechanism of R-loops in cell reprogramming remains elusive. Here, we comprehensively profiled the genome-wide landscape of R-loops during cell reprogramming. The results showed that the R-loop formation on most different types of repetitive elements is stage-specific in cell reprogramming. We unveiled that the cumulative deposition of an R-loop subset is positively correlated with gene expression during reprogramming. More importantly, the dynamic turnover of this R-loop subset is accompanied by the activation of the pluripotent transcriptional regulatory network (TRN). Moreover, the large accumulation of the active histone marker H3K4me3 and the reduction in H3K27me3 were also observed in these R-loop regions. Finally, we characterized the dynamic network of R-loops that facilitates cell fate transitions in reprogramming. Together, our study provides a new clue for deciphering the interplay mechanism between R-loops and HMs to control cell reprogramming.