HEARTSVG: a fast and accurate method for identifying spatially variable genes in large-scale spatial transcriptomics.

Identifying spatially variable genes (SVGs) is crucial for understanding the spatiotemporal characteristics of diseases and tissue structures, posing a distinctive challenge in spatial transcriptomics research. We propose HEARTSVG, a distribution-free, test-based method for fast and accurately identifying spatially variable genes in large-scale spatial transcriptomic data. Extensive simulations demonstrate that HEARTSVG outperforms state-of-the-art methods with higher F 1 scores (average F 1 Score=0.948), improved computational efficiency, scalability, and reduced false positives (FPs). Through analysis of twelve real datasets from various spatial transcriptomic technologies, HEARTSVG identifies a greater number of biologically significant SVGs (average AUC = 0.792) than other comparative methods without prespecifying spatial patterns. Furthermore, by clustering SVGs, we uncover two distinct tumor spatial domains characterized by unique spatial expression patterns, spatial-temporal locations, and biological functions in human colorectal cancer data, unraveling the complexity of tumors.

CTHRC1+ fibroblasts and SPP1+ macrophages synergistically contribute to pro-tumorigenic tumor microenvironment in pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is an extremely lethal cancer that accounts for over 90% of all pancreatic cancer cases. With a 5-year survival rate of only 13%, PDAC has proven to be extremely desmoplastic and immunosuppressive to most current therapies, including chemotherapy and surgical resection. In recent years, focus has shifted to understanding the tumor microenvironment (TME) around PDAC, enabling a greater understanding of biological pathways and intercellular interactions that can ultimately lead to potential for future drug targets. In this study, we leverage a combination of single-cell and spatial transcriptomics to further identify cellular populations and interactions within the highly heterogeneous TME. We demonstrate that SPP1+APOE+ tumor-associated macrophages (TAM) and CTHRC1+GREM1+ cancer-associated myofibroblasts (myCAF) not only act synergistically to promote an immune-suppressive TME through active extracellular matrix (ECM) deposition and epithelial mesenchymal transition (EMT), but are spatially colocalized and correlated, leading to worse prognosis. Our results highlight the crosstalk between stromal and myeloid cells as a significant area of study for future therapeutic targets to treat cancer.

A machine learning and drug repurposing approach to target ferroptosis in colorectal cancer stratified by sex and KRAS.

The landscape of sex differences in Colorectal Cancer (CRC) has not been well characterized with respect to the mechanisms of action for oncogenes such as KRAS. However, our recent study showed that tumors from male patients with KRAS mutations have decreased iron-dependent cell death called ferroptosis. Building on these findings, we further examined ferroptosis in CRC, considering both sex of the patient and KRAS mutations, using public databases and our in-house CRC tumor cohort. Through subsampling inference and variable importance analysis (VIMP), we identified significant differences in gene expression between KRAS mutant and wild type tumors from male patients. These genes suppress (e.g., SLC7A11 ) or drive (e.g., SLC1A5 ) ferroptosis, and these findings were further validated with Gaussian mixed models. Furthermore, we explored the prognostic value of ferroptosis regulating genes and discovered sex- and KRAS-specific differences at both the transcriptional and metabolic levels by random survival forest with backward elimination algorithm (RSF-BE). Of note, genes and metabolites involved in arginine synthesis and glutathione metabolism were uniquely associated with prognosis in tumors from males with KRAS mutations. Additionally, drug repurposing is becoming popular due to the high costs, attrition rates, and slow pace of new drug development, offering a way to treat common and rare diseases more efficiently. Furthermore, increasing evidence has shown that ferroptosis inhibition or induction can improve drug sensitivity or overcome chemotherapy drug resistance. Therefore, we investigated the correlation between gene expression, metabolite levels, and drug sensitivity across all CRC primary tumor cell lines using data from the Genomics of Drug Sensitivity in Cancer (GDSC) resource. We observed that ferroptosis suppressor genes such as DHODH , GCH1 , and AIFM2 in KRAS mutant CRC cell lines were resistant to cisplatin and paclitaxel, underscoring why these drugs are not effective for these patients. The comprehensive map generated here provides valuable biological insights for future investigations, and the findings are supported by rigorous analysis of large-scale publicly available data and our in-house cohort. The study also emphasizes the potential application of VIMP, Gaussian mixed models, and RSF-BE models in the multi-omics research community. In conclusion, this comprehensive approach opens doors for leveraging precision molecular feature analysis and drug repurposing possibilities in KRAS mutant CRC.

Incorporating prior information in gene expression network-based cancer heterogeneity analysis.

Cancer is molecularly heterogeneous, with seemingly similar patients having different molecular landscapes and accordingly different clinical behaviors. In recent studies, gene expression networks have been shown as more effective/informative for cancer heterogeneity analysis than some simpler measures. Gene interconnections can be classified as "direct" and "indirect," where the latter can be caused by shared genomic regulators (such as transcription factors, microRNAs, and other regulatory molecules) and other mechanisms. It has been suggested that incorporating the regulators of gene expressions in network analysis and focusing on the direct interconnections can lead to a deeper understanding of the more essential gene interconnections. Such analysis can be seriously challenged by the large number of parameters (jointly caused by network analysis, incorporation of regulators, and heterogeneity) and often weak signals. To effectively tackle this problem, we propose incorporating prior information contained in the published literature. A key challenge is that such prior information can be partial or even wrong. We develop a two-step procedure that can flexibly accommodate different levels of prior information quality. Simulation demonstrates the effectiveness of the proposed approach and its superiority over relevant competitors. In the analysis of a breast cancer dataset, findings different from the alternatives are made, and the identified sample subgroups have important clinical differences.

Partial Hepatectomy and Ablation for Survival of Early-Stage Hepatocellular Carcinoma Patients: A Bayesian Emulation Analysis.

Partial hepatectomy and ablation therapy are two widely used surgical procedures for localized early-stage hepatocellular carcinoma (HCC) patients. This article aimed to evaluate their relative effectiveness in terms of overall survival. An emulation analysis approach was first developed based on the Bayesian technique. We estimated propensity scores via Bayesian logistic regression and adopted a weighted Bayesian Weibull accelerated failure time (AFT) model incorporating prior information contained in the published literature. With the Surveillance, Epidemiology, and End Results (SEER)-Medicare data, an emulated target trial with rigorously defined inclusion/exclusion criteria and treatment regimens for early-stage HCC patients over 66 years old was developed. For the main cohort with tumor size less than or equal to 5 cm, a total of 1146 patients were enrolled in the emulated trial, with 301 and 845 in the partial hepatectomy and ablation arms, respectively. The analysis suggested ablation to be significantly associated with inferior overall survival (hazard ratio [HR] = 1.35; 95% credible interval [CrI]: 1.14, 1.60). For the subgroup with tumor size less than or equal to 3 cm, there was no significant difference in overall survival between the two arms (HR = 1.15; 95% CrI: 0.88, 1.52). Overall, the comparative treatment effect of ablation and partial hepatectomy on survival remains inconclusive. This finding may provide further insight into HCC clinical treatment.

Prediction Consistency Regularization for Learning with Noise Labels Based on Contrastive Clustering.

In the classification task, label noise has a significant impact on models' performance, primarily manifested in the disruption of prediction consistency, thereby reducing the classification accuracy. This work introduces a novel prediction consistency regularization that mitigates the impact of label noise on neural networks by imposing constraints on the prediction consistency of similar samples. However, determining which samples should be similar is a primary challenge. We formalize the similar sample identification as a clustering problem and employ twin contrastive clustering (TCC) to address this issue. To ensure similarity between samples within each cluster, we enhance TCC by adjusting clustering prior to distribution using label information. Based on the adjusted TCC's clustering results, we first construct the prototype for each cluster and then formulate a prototype-based regularization term to enhance prediction consistency for the prototype within each cluster and counteract the adverse effects of label noise. We conducted comprehensive experiments using benchmark datasets to evaluate the effectiveness of our method under various scenarios with different noise rates. The results explicitly demonstrate the enhancement in classification accuracy. Subsequent analytical experiments confirm that the proposed regularization term effectively mitigates noise and that the adjusted TCC enhances the quality of similar sample recognition.

Estimation of multiple networks with common structures in heterogeneous subgroups.

Network estimation has been a critical component of high-dimensional data analysis and can provide an understanding of the underlying complex dependence structures. Among the existing studies, Gaussian graphical models have been highly popular. However, they still have limitations due to the homogeneous distribution assumption and the fact that they are only applicable to small-scale data. For example, cancers have various levels of unknown heterogeneity, and biological networks, which include thousands of molecular components, often differ across subgroups while also sharing some commonalities. In this article, we propose a new joint estimation approach for multiple networks with unknown sample heterogeneity, by decomposing the Gaussian graphical model (GGM) into a collection of sparse regression problems. A reparameterization technique and a composite minimax concave penalty are introduced to effectively accommodate the specific and common information across the networks of multiple subgroups, making the proposed estimator significantly advancing from the existing heterogeneity network analysis based on the regularized likelihood of GGM directly and enjoying scale-invariant, tuning-insensitive, and optimization convexity properties. The proposed analysis can be effectively realized using parallel computing. The estimation and selection consistency properties are rigorously established. The proposed approach allows the theoretical studies to focus on independent network estimation only and has the significant advantage of being both theoretically and computationally applicable to large-scale data. Extensive numerical experiments with simulated data and the TCGA breast cancer data demonstrate the prominent performance of the proposed approach in both subgroup and network identifications.

Endocrine disrupting chemical mixture exposure and risk of papillary thyroid cancer in U.S. military personnel: A nested case-control study.

Single-pollutant methods to evaluate associations between endocrine disrupting chemicals (EDCs) and thyroid cancer risk may not reflect realistic human exposures. Therefore, we evaluated associations between exposure to a mixture of 18 EDCs, including polychlorinated biphenyls (PCBs), brominated flame retardants, and organochlorine pesticides, and risk of papillary thyroid cancer (PTC), the most common thyroid cancer histological subtype. We conducted a nested case-control study among U.S. military servicemembers of 652 histologically-confirmed PTC cases diagnosed between 2000 and 2013 and 652 controls, matched on birth year, sex, race/ethnicity, military component (active duty/reserve), and serum sample timing. We estimated mixture odds ratios (OR), 95% confidence intervals (95% CI), and standard errors (SE) for associations between pre-diagnostic serum EDC mixture concentrations, overall PTC risk, and risk of histological subtypes of PTC (classical, follicular), adjusted for body mass index and military branch, using quantile g-computation. Additionally, we identified relative contributions of individual mixture components to PTC risk, represented by positive and negative weights (w). A one-quartile increase in the serum mixture concentration was associated with a non-statistically significant increase in overall PTC risk (OR = 1.19; 95% CI = 0.91, 1.56; SE = 0.14). Stratified by histological subtype and race (White, Black), a one-quartile increase in the mixture was associated with increased classical PTC risk among those of White race (OR = 1.59; 95% CI = 1.06, 2.40; SE = 0.21), but not of Black race (OR = 0.95; 95% CI = 0.34, 2.68; SE = 0.53). PCBs 180, 199, and 118 had the greatest positive weights driving this association among those of White race (w = 0.312, 0.255, and 0.119, respectively). Findings suggest that exposure to an EDC mixture may be associated with increased classical PTC risk. These findings warrant further investigation in other study populations to better understand PTC risk by histological subtype and race.

Organochlorine pesticides and risk of papillary thyroid cancer in U.S. military personnel: a nested case-control study.

BACKGROUND: The effects of organochlorine pesticide (OCP) exposure on the development of human papillary thyroid cancer (PTC) are not well understood. A nested case-control study was conducted with data from the U.S. Department of Defense Serum Repository (DoDSR) cohort between 2000 and 2013 to assess associations of individual OCPs serum concentrations with PTC risk. METHODS: This study included 742 histologically confirmed PTC cases (341 females, 401 males) and 742 individually-matched controls with pre-diagnostic serum samples selected from the DoDSR. Associations between categories of lipid-corrected serum concentrations of seven OCPs and PTC risk were evaluated for classical PTC and follicular PTC using conditional logistic regression, adjusted for body mass index category and military branch to compute odds ratios (OR) and 95% confidence intervals (CIs). Effect modification by sex, birth cohort, and race was examined. RESULTS: There was no evidence of associations between most of the OCPs and PTC, overall or stratified by histological subtype. Overall, there was no evidence of an association between hexachlorobenzene (HCB) and PTC, but stratified by histological subtype HCB was associated with significantly increased risk of classical PTC (third tertile above the limit of detection (LOD) vs.

Information-incorporated sparse hierarchical cancer heterogeneity analysis.

Cancer heterogeneity analysis is essential for precision medicine. Most of the existing heterogeneity analyses only consider a single type of data and ignore the possible sparsity of important features. In cancer clinical practice, it has been suggested that two types of data, pathological imaging and omics data, are commonly collected and can produce hierarchical heterogeneous structures, in which the refined sub-subgroup structure determined by omics features can be nested in the rough subgroup structure determined by the imaging features. Moreover, sparsity pursuit has extraordinary significance and is more challenging for heterogeneity analysis, because the important features may not be the same in different subgroups, which is ignored by the existing heterogeneity analyses. Fortunately, rich information from previous literature (for example, those deposited in PubMed) can be used to assist feature selection in the present study. Advancing from the existing analyses, in this study, we propose a novel sparse hierarchical heterogeneity analysis framework, which can integrate two types of features and incorporate prior knowledge to improve feature selection. The proposed approach has satisfactory statistical properties and competitive numerical performance. A TCGA real data analysis demonstrates the practical value of our approach in analyzing data heterogeneity and sparsity.

Neurogenetic underpinnings of nicotine use severity: Integrating the brain transcriptomes and GWAS variants via network approaches.

Our study focused on human brain transcriptomes and the genetic risks of cigarettes per day (CPD) to investigate the neurogenetic mechanisms of individual variation in nicotine use severity. We constructed whole-brain and intramodular region-specific coexpression networks using BrainSpan's transcriptomes, and the genomewide association studies identified risk variants of CPD, confirmed the associations between CPD and each gene set in the region-specific subnetworks using an independent dataset, and conducted bioinformatic analyses. Eight brain-region-specific coexpression subnetworks were identified in association with CPD: amygdala, hippocampus, medial prefrontal cortex (MPFC), orbitofrontal cortex (OPFC), dorsolateral prefrontal cortex, striatum, mediodorsal nucleus of the thalamus (MDTHAL), and primary motor cortex (M1C). Each gene set in the eight subnetworks was associated with CPD. We also identified three hub proteins encoded by GRIN2A in the amygdala, PMCA2 in the hippocampus, MPFC, OPFC, striatum, and MDTHAL, and SV2B in M1C. Intriguingly, the pancreatic secretion pathway appeared in all the significant protein interaction subnetworks, suggesting pleiotropic effects between cigarette smoking and pancreatic diseases. The three hub proteins and genes are implicated in stress response, drug memory, calcium homeostasis, and inhibitory control. These findings provide novel evidence of the neurogenetic underpinnings of smoking severity.

Hierarchical False Discovery Rate Control for High-dimensional Survival Analysis with Interactions.

With the development of data collection techniques, analysis with a survival response and high-dimensional covariates has become routine. Here we consider an interaction model, which includes a set of low-dimensional covariates, a set of high-dimensional covariates, and their interactions. This model has been motivated by gene-environment (G-E) interaction analysis, where the E variables have a low dimension, and the G variables have a high dimension. For such a model, there has been extensive research on estimation and variable selection. Comparatively, inference studies with a valid false discovery rate (FDR) control have been very limited. The existing high-dimensional inference tools cannot be directly applied to interaction models, as interactions and main effects are not "equal". In this article, for high-dimensional survival analysis with interactions, we model survival using the Accelerated Failure Time (AFT) model and adopt a "weighted least squares + debiased Lasso" approach for estimation and selection. A hierarchical FDR control approach is developed for inference and respect of the "main effects, interactions" hierarchy. The asymptotic distribution properties of the debiased Lasso estimators are rigorously established. Simulation demonstrates the satisfactory performance of the proposed approach, and the analysis of a breast cancer dataset further establishes its practical utility.

The Bayesian Regularized Quantile Varying Coefficient Model.

The quantile varying coefficient (VC) model can flexibly capture dynamical patterns of regression coefficients. In addition, due to the quantile check loss function, it is robust against outliers and heavy-tailed distributions of the response variable, and can provide a more comprehensive picture of modeling via exploring the conditional quantiles of the response variable. Although extensive studies have been conducted to examine variable selection for the high-dimensional quantile varying coefficient models, the Bayesian analysis has been rarely developed. The Bayesian regularized quantile varying coefficient model has been proposed to incorporate robustness against data heterogeneity while accommodating the non-linear interactions between the effect modifier and predictors. Selecting important varying coefficients can be achieved through Bayesian variable selection. Incorporating the multivariate spike-and-slab priors further improves performance by inducing exact sparsity. The Gibbs sampler has been derived to conduct efficient posterior inference of the sparse Bayesian quantile VC model through Markov chain Monte Carlo (MCMC). The merit of the proposed model in selection and estimation accuracy over the alternatives has been systematically investigated in simulation under specific quantile levels and multiple heavy-tailed model errors. In the case study, the proposed model leads to identification of biologically sensible markers in a non-linear gene-environment interaction study using the NHS data.