Chun-Yu Lin

Background: With the advance in single-cell RNA sequencing (scRNA-seq) technology, deriving inherent biological system information from expression profiles at a single-cell resolution has become possible. It has been known that network modeling by estimating the associations between genes could better reveal dynamic changes in biological systems. However, accurately constructing a single-cell network (SCN) to capture the network architecture of each cell and further explore cell-to-cell heterogeneity remains challenging. Results: We introduce SINUM, a method for constructing the SIngle-cell Network Using Mutual information, which estimates mutual information between any two genes from scRNA-seq data to determine whether they are dependent or independent in a specific cell. Experiments on various scRNA-seq datasets with different cell numbers based on eight performance indexes (e.g., adjusted rand index and F-measure index) validated the accuracy and robustness of SINUM in cell type identification, superior to the state-of-the-art SCN inference method. Additionally, the SINUM SCNs exhibit high overlap with the human interactome and possess the scale-free property. Conclusions: SINUM presents a view of biological systems at the network level to detect cell-type marker genes/gene pairs and investigate time-dependent changes in gene associations during embryo development. Codes for SINUM are freely available at https://github.com/SysMednet/SINUM.

The management of chronic myelogenous leukemia (CML) has seen significant progress with the introduction of tyrosine kinase inhibitors (TKIs), particularly Imatinib. However, a notable proportion of CML patients develop resistance to Imatinib, often due to the persistence of leukemia stem cells and resistance mechanisms independent of BCR::ABL1 This study investigates the roles of IL6R, IL7R, and MYC in Imatinib resistance by employing CRISPR/Cas9 for gene editing and the Non-Invasive Apoptosis Detection Sensor version 2 (NIADS v2) for apoptosis assessment. The results indicate that Imatinib-resistant K562 cells (K562-IR) predominantly express IL6R, IL7R, and MYC, with IL6R and MYC playing crucial roles in cell survival and sensitivity to Imatinib. Conversely, IL7R does not significantly impact cytotoxicity, either alone or in combination with Imatinib. Further genetic editing experiments confirm the protective functions of IL6R and MYC in K562-IR cells, suggesting their potential as therapeutic targets for overcoming Imatinib resistance in CML. This study contributes to understanding the mechanisms of Imatinib resistance in CML, proposing IL6R and MYC as pivotal targets for therapeutic strategies. Moreover, the utilization of NIADS v2 enhances our capability to analyze apoptosis and drug responses, contributing to a deeper understanding of CML pathogenesis and treatment options.

Pathway analysis, including nontopology-based (non-TB) and topology-based (TB) methods, is widely used to interpret the biological phenomena underlying differences in expression data between two phenotypes. By considering dependencies and interactions between genes, TB methods usually perform better than non-TB methods in identifying pathways that include closely relevant or directly causative genes for a given phenotype. However, most TB methods may be limited by incomplete pathway data used as the reference network or by difficulties in selecting appropriate reference networks for different research topics. Here, we propose a gene set correlation enrichment analysis method, Gscore, based on an expression dataset-derived coexpression network to examine whether a differentially expressed gene (DEG) list (or each of its DEGs) is associated with a known gene set. Gscore is better able to identify target pathways in 89 human disease expression datasets than eight other state-of-the-art methods and offers insight into how disease-wide and pathway-wide associations reflect clinical outcomes. When applied to RNA-seq data from COVID-19-related cells and patient samples, Gscore provided a means for studying how DEGs are implicated in COVID-19-related pathways. In summary, Gscore offers a powerful analytical approach for annotating individual DEGs, DEG lists, and genome-wide expression profiles based on existing biological knowledge.

Head and neck squamous cell carcinoma (HNSC) exhibits genetic heterogeneity in etiologies, tumor sites, and biological processes, which significantly impact therapeutic strategies and prognosis. While the influence of human papillomavirus on clinical outcomes is established, the molecular subtypes determining additional treatment options for HNSC remain unclear and inconsistent. This study aims to identify distinct HNSC molecular subtypes to enhance diagnosis and prognosis accuracy. In this study, we collected three HNSC microarrays (n = 306) from the Gene Expression Omnibus (GEO), and HNSC RNA-Seq data (n = 566) from The Cancer Genome Atlas (TCGA) to identify differentially expressed genes (DEGs) and validate our results. Two scoring methods, representative score (RS) and perturbative score (PS), were developed for DEGs to summarize their possible activation functions and influence in tumorigenesis. Based on the RS and PS scoring, we selected candidate genes to cluster TCGA samples for the identification of molecular subtypes in HNSC. We have identified 289 up-regulated DEGs and selected 88 genes (called HNSC88) using the RS and PS scoring methods. Based on HNSC88 and TCGA samples, we determined three HNSC subtypes, including one HPV-associated subtype, and two HPV-negative subtypes. One of the HPV-negative subtypes showed a relationship to smoking behavior, while the other exhibited high expression in tumor immune response. The Kaplan–Meier method was used to compare overall survival among the three subtypes. The HPV-associated subtype showed a better prognosis compared to the other two HPV-negative subtypes (log rank, p = 0.0092 and 0.0001; hazard ratio, 1.36 and 1.39). Additionally, within the HPV-negative group, the smoking-related subgroup exhibited worse prognosis compared to the subgroup with high expression in immune response (log rank, p = 0.039; hazard ratio, 1.53). The HNSC88 not only enables the identification of HPV-associated subtypes, but also proposes two potential HPV-negative subtypes with distinct prognoses and molecular signatures. This study provides valuable strategies for summarizing the roles and influences of genes in tumorigenesis for identifying molecular signatures and subtypes of HNSC.

Analyzing multiple networks is important to understand relevant features among different networks. Although many studies have been conducted for that purpose, not much attention has been paid to the analysis of attractors (i.e., steady states) in multiple networks. Therefore, we study common attractors and similar attractors in multiple networks to uncover hidden similarities and differences among networks using Boolean networks (BNs), where BNs have been used as a mathematical model of genetic networks and neural networks. We define three problems on detecting common attractors and similar attractors, and theoretically analyze the expected number of such objects for random BNs, where we assume that given networks have the same set of nodes (i.e., genes). We also present four methods for solving these problems. Computational experiments on randomly generated BNs are performed to demonstrate the efficiency of our proposed methods. In addition, experiments on a practical biological system, a BN model of the TGF- $\beta$ signaling pathway, are performed. The result suggests that common attractors and similar attractors are useful for exploring tumor heterogeneity and homogeneity in eight cancers.

Recently, extracting inherent biological system information (e.g. cellular networks) from genome-wide expression profiles for developing personalized diagnostic and therapeutic strategies has become increasingly important. However, accurately constructing single-sample networks (SINs) to capture individual characteristics and heterogeneity in disease remains challenging. Here, we propose a sample-specific-weighted correlation network (SWEET) method to model SINs by integrating the genome-wide sample-to-sample correlation (i.e. sample weights) with the differential network between perturbed and aggregate networks. For a group of samples, the genome-wide sample weights can be assessed without prior knowledge of intrinsic subpopulations to address the network edge number bias caused by sample size differences. Compared with the state-of-the-art SIN inference methods, the SWEET SINs in 16 cancers more likely fit the scale-free property, display higher overlap with the human interactomes and perform better in identifying three types of cancer-related genes. Moreover, integrating SWEET SINs with a network proximity measure facilitates characterizing individual features and therapy in diseases, such as somatic mutation, mut-driver and essential genes. Biological experiments further validated two candidate repurposable drugs, albendazole for head and neck squamous cell carcinoma (HNSCC) and lung adenocarcinoma (LUAD) and encorafenib for HNSCC. By applying SWEET, we also identified two possible LUAD subtypes that exhibit distinct clinical features and molecular mechanisms. Overall, the SWEET method complements current SIN inference and analysis methods and presents a view of biological systems at the network level to offer numerous clues for further investigation and clinical translation in network medicine and precision medicine.

Chronic myelogenous leukemia (CML) is the most common type of leukemia in adults, and more than 90% of CML patients harbor the abnormal Philadelphia chromosome (Ph) that encodes the BCR-ABL oncoprotein. Although the ABL kinase inhibitor (imatinib) has proven to be very effective in achieving high remission rates and improving prognosis, up to 33% of CML patients still cannot achieve an optimal response. Here, we used CRISPR/Cas9 to specifically target the BCR-ABL junction region in K562 cells, resulting in the inhibition of cancer cell growth and oncogenesis. Due to the variety of BCR-ABL junctions in CML patients, we utilized gene editing of the human ABL gene for clinical applications. Using the ABL gene-edited virus in K562 cells, we detected 41.2% indels in ABL sgRNA_2-infected cells. The ABL-edited cells reveled significant suppression of BCR-ABL protein expression and downstream signals, inhibiting cell growth and increasing cell apoptosis. Next, we introduced the ABL gene-edited virus into a systemic K562 leukemia xenograft mouse model, and bioluminescence imaging of the mice showed a significant reduction in the leukemia cell population in ABL-targeted mice, compared to the scramble sgRNA virus-injected mice. In CML cells from clinical samples, infection with the ABL gene-edited virus resulted in more than 30.9% indels and significant cancer cell death. Notably, no off-target effects or bone marrow cell suppression was found using the ABL gene-edited virus, ensuring both user safety and treatment efficacy. This study demonstrated the critical role of the ABL gene in maintaining CML cell survival and tumorigenicity in vitro and in vivo. ABL gene editing-based therapy might provide a potential strategy for imatinib-insensitive or resistant CML patients.

Although many studies have shown the association between smoking and the increased incidence and adverse prognosis of head and neck squamous cell carcinoma (HNSCC), the mechanisms and pharmaceutical targets involved remain unclear. Here, we integrated gene expression signatures, genetic alterations, and survival analyses to identify prognostic indicators and therapeutic targets for smoking HNSCC patients, and we discovered that the FDA-approved drug varenicline inhibits the target for cancer cell migration/invasion. We first identified 18 smoking-related and prognostic genes for HNSCC by using RNA-Seq and clinical follow-up data. One of these genes, CHRNB4 (neuronal acetylcholine receptor subunit beta-4), increased the risk of death by approximately threefold in CHRNB4-high expression smokers compared to CHRNB4-low expression smokers (log rank, p = 0.00042; hazard ratio, 2.82; 95% CI, 1.55–5.14), former smokers, and non-smokers. Furthermore, we examined the functional enrichment of co-regulated genes of CHRNB4 and its 246 frequently occurring copy number alterations (CNAs). We found that these genes were involved in promoting angiogenesis, resisting cell death, and sustaining proliferation, and contributed to much worse outcomes for CHRNB4-high patients. Finally, we performed CHRNB4 gene editing and drug inhibition assays, and the results validate these observations. In summary, our study suggests that CHRNB4 is a prognostic indicator for smoking HNSCC patients and provides a potential new therapeutic drug to prevent recurrence or distant metastasis.

Alterations in membrane proteins (MPs) and their regulated pathways have been established as cancer hallmarks and extensively targeted in clinical applications. However, the analysis of MP-interacting proteins and downstream pathways across human malignancies remains challenging. Here, we present a systematically integrated method to generate a resource of cancer membrane protein-regulated networks (CaMPNets), containing 63,746 high-confidence protein–protein interactions (PPIs) for 1962 MPs, using expression profiles from 5922 tumors with overall survival outcomes across 15 human cancers. Comprehensive analysis of CaMPNets links MP partner communities and regulated pathways to provide MP-based gene sets for identifying prognostic biomarkers and druggable targets. For example, we identify CHRNA9 with 12 PPIs (e.g., ERBB2) can be a therapeutic target and find its anti-metastasis agent, bupropion, for treatment in nicotine-induced breast cancer. This resource is a study to systematically integrate MP interactions, genomics, and clinical outcomes for helping illuminate cancer-wide atlas and prognostic landscapes in tumor homo/heterogeneity.

Prostate cancer (PCa) is the second leading cause of cancer death among men worldwide. About 70% of PCa patients were diagnosed at later stage, and metastasis has been observed. Additionally, the cure rate of PCa closely relies on the early diagnosis with biomarkers. The identification of biomarkers for diagnosis and prognosis is an urgent clinical issue for PCa. Here, we developed a novel scoring strategy, including cluster score (CS) and predicting score (PS), to identify 29 PCa genes (called PCa29) for early diagnostic biomarkers from two datasets in Gene Expression Omnibus. The result indicates that PCa29 can discriminate between normal and tumor tissues and are specific for prostate cancer. To validate PCa29, we found that 97% of PCa29 were consistently significant with these gene expressions in The Cancer Genome Atlas; furthermore, [ ]70% of PCa29 are consensus to the protein expression in The Human Protein Atlas. Finally, we examined 10 genes in PCa29 on three PCa cell lines by real-time quantitative polymerase chain reaction. The experimental results show that the trend of the differential PCa29 expression is consistent with the analyzed results from our novel scoring method. We believe that our method is useful and PCa29 are potential biomarkers that provide the clues to develop targeting therapy for PCa.