Citations

Summary Long-read nanopore sequencing enables simultaneous detection of germline variation and native DNA base modifications on individual DNA molecules, providing a unique opportunity to investigate allele-specific epigenetic regulation. Here, we performed whole-genome nanopore sequencing on normal and tumor prostate tissues to characterize differential methylation, methylation entropy, and allele-specific methylation (ASM) associated with noncoding genetic variants. Genome-wide analysis identified extensive cancer-associated differentially methylated regions (DMRs), with hypermethylated DMRs significantly enriched near transcription start sites and transcriptional regulatory regions. Integration with transcriptomic datasets revealed strong inverse relationships between promoter methylation and gene expression, while 5-hydroxymethylcytosine (5hmC) levels positively correlated with transcriptional activity across gene bodies. Using fragment-level methylation patterns enabled by long-read sequencing, we further quantified methylation entropy incorporating both 5mCG and 5hmCG states. Cancer-hypermethylated DMRs exhibited markedly reduced entropy, consistent with clonal fixation of methylation states during tumor progression. Entropy profiling across chromatin annotations demonstrated maximal epigenetic heterogeneity at partially modified enhancer-associated regions. To investigate cis-regulatory genetic effects, we developed a simple ASM framework (nanoASM) that can partition sequencing reads by allelic state and identifies allele-specific DMRs directly from long-read data. Compared with conventional population-level mQTL analysis, ASM demonstrated substantially improved statistical efficiency by leveraging within-individual contrasts and reducing sample-level heterogeneity. Although germline single nucleotide polymorphisms (SNPs) were largely shared between normal and tumor tissues, ASM patterns differed substantially, with tumor-associated ASM regions displaying significantly larger genomic span and stronger allelic methylation differences. Comparative analysis with TCGA prostate mQTL and GTEx prostate eQTL datasets demonstrated substantial concordance between ASM directionality and downstream transcriptional effects, particularly for variants located within DMRs and near transcription start sites. At the IRX4 prostate cancer risk locus, ASM identified an androgen-responsive regulatory domain overlapping AR ChIP-seq and H3K27ac peaks, nominating rs6885084 as a candidate functional variant. At the PSCA locus, ASM anchored by rs4736369 was associated with allele-specific methylation, chromatin activation, transcript abundance, and isoform usage. Together, these findings establish nanopore-based ASM analysis as a powerful approach for resolving functional noncoding variants and their regulatory domains they control in prostate cancer.

Abstract Background Lumbar puncture is a routinely performed procedure, yet it presents significant challenges in patients with abnormal spinal anatomy, such as severe scoliosis or obesity. Conventional ultrasound-assisted neuraxial techniques, while superior to landmark-based methods, are constrained by two fundamental limitations: inaccurate projection of the probe’s geometric center onto the skin when the probe is tilted or rocked, and the difficulty of consistently replicating the exact probe angle during needle insertion. These limitations contribute to multiple needle passes, prolonged procedure time, and increased patient discomfort and complication risk. To address these issues, we developed a novel angle-free ultrasound-assisted technique designed to eliminate the need for replicating probe angulation during needle insertion. Case presentation: This angle-free approach was utilized for serial lumbar punctures in two patients with spinal muscular atrophy requiring repeated intrathecal nusinersen administration. Both patients presented with extreme anatomical complexity, including high body mass index (30 and 35.1 kg/m²) and severe scoliosis (Cobb angles 90° and 83°). The technique follows a three-step protocol: (1) Preprocedural imaging to select the optimal interlaminar space; (2) Ultrasound-guided postural adjustment until the probe could be held perpendicular to the bed edge and parallel to the floor while maintaining target visualization; (3) Needle insertion along this orthogonal trajectory. All 22 procedures (11 per patient) were completed successfully without persistent neurological complications. Procedural outcomes demonstrated a first-attempt success rate of 81.8% (Case 1) and 100% (Case 2), with first-pass success rates of 45.5% and 54.5%, respectively. Longitudinal analysis across sequential procedures revealed a consistent reduction in both the number of needle passes and needling time, reflecting a favorable learning curve. Conclusions The angle-free ultrasound-assisted neuraxial technique offers a valuable alternative for patients with complex spinal anatomy. By using real-time postural adjustments to achieve a needle trajectory, it enhances procedural precision and eliminates the technical hurdle of replicating probe angulation. The structured "4P" framework (Planning, Positioning, Perpendicular, Parallel) offers a practical guide for clinical adoption. The observed learning curve suggests that clinicians can readily acquire and apply this technique to manage such challenging cases.

Abstract Network-based methods utilize protein–protein interaction information to identify significantly perturbed subnetworks in cancer and to propose key molecular pathways. Numerous methods have been developed, but to date, a rigorous benchmark analysis to compare the performance of existing approaches is lacking. In this paper, we proposed a novel benchmarking framework using synthetic data and conducted a comprehensive analysis to investigate the ability of existing methods to detect target genes and subnetworks and to control false positives, and how they perform in the presence of topological biases at both gene and subnetwork levels. Our analysis revealed insights into algorithmic performance that were previously unattainable. Based on the results of the benchmark study, we presented a practical guide for users on how to select appropriate detection methods and protein–protein interaction networks for cancer pathway identification, and provided suggestions for future algorithm development.

Abstract Prostate cancer is a significant health concern and the most commonly diagnosed cancer in men worldwide. Understanding the complex process of prostate tumor evolution and progression is crucial for improved diagnosis, treatments, and patient outcomes. Previous studies have focused on unraveling the dynamics of prostate cancer evolution using phylogenetic or lineage analysis approaches. However, those approaches have limitations in capturing the complete disease process or incorporating genomic and transcriptomic variations comprehensively. In this study, we applied a novel computational approach to derive a prostate cancer progression model using multidimensional data from 497 prostate tumor samples and 52 tumor-adjacent normal samples obtained from The Cancer Genome Atlas study. The model was validated using data from an independent cohort of 545 primary tumor samples. By integrating transcriptomic and genomic data, our model provides a comprehensive view of prostate tumor progression, identifies crucial signaling pathways and genetic events, and uncovers distinct transcription signatures associated with disease progression. Our findings have significant implications for cancer research and hold promise for guiding personalized treatment strategies in prostate cancer. Significance: We developed and validated a progression model of prostate cancer using >1,000 tumor and normal tissue samples. The model provided a comprehensive view of prostate tumor evolution and progression.

Background: The identification of cancer driver genes and key molecular pathways has been the focus of large-scale cancer genome studies. Network-based methods detect significantly perturbed subnetworks as putative cancer pathways by incorporating genomics data with the topological information of PPI networks. However, commonly used PPI networks have distinct topological structures, making the results of the same method vary widely when applied to different networks. Furthermore, emerging context-specific PPI networks often have incomplete topological structures, which pose serious challenges for existing subnetwork detection algorithms. Methods: In this paper, we propose a novel method, referred to as MultiFDRnet, to address the above issues. The basic idea is to model a set of PPI networks as a multiplex network to preserve the topological structure of individual networks, while introducing dependencies among them, and, then, to detect significantly perturbed subnetworks on the modeled multiplex network using all the structural information simultaneously. Results: To illustrate the effectiveness of the proposed approach, an extensive benchmark analysis was conducted on both simulated and real cancer data. The experimental results showed that the proposed method is able to detect significantly perturbed subnetworks jointly supported by multiple PPI networks and to identify novel modular structures in context-specific PPI networks.

Cross-talk between Rho- and Arf-family guanosine triphosphatases (GTPases) plays an important role in linking the actin cytoskeleton to membrane protrusions, organelle morphology, and vesicle trafficking. The central actin regulator, WAVE regulatory complex (WRC), integrates Rac1 (a Rho-family GTPase) and Arf signaling to promote Arp2/3-mediated actin polymerization in many processes, but how WRC senses Arf signaling is unknown. Here, we have reconstituted a direct interaction between Arf and WRC. This interaction is greatly enhanced by Rac1 binding to the D site of WRC. Arf1 binds to a previously unidentified, conserved surface on the Sra1 subunit of WRC, which, in turn, drives WRC activation using a mechanism distinct from that of Rac1. Mutating the Arf binding site abolishes Arf1-WRC interaction, disrupts Arf1-mediated WRC activation, and impairs lamellipodia formation and cell migration. This work uncovers a new mechanism underlying WRC activation and provides a mechanistic foundation for studying how WRC-mediated actin polymerization links Arf and Rac signaling in cells.

Abstract Summary Quantifying pairwise sequence similarities is a key step in metagenomics studies. Alignment-free methods provide a computationally efficient alternative to alignment-based methods for large-scale sequence analysis. Several neural network-based methods have recently been developed for this purpose. However, existing methods do not perform well on sequences of varying lengths and are sensitive to the presence of insertions and deletions. In this article, we describe the development of a new method, referred to as AsMac that addresses the aforementioned issues. We proposed a novel neural network structure for approximate string matching for the extraction of pertinent information from biological sequences and developed an efficient gradient computation algorithm for training the constructed neural network. We performed a large-scale benchmark study using real-world data that demonstrated the effectiveness and potential utility of the proposed method. Availability and implementation The open-source software for the proposed method and trained neural-network models for some commonly used metagenomics marker genes were developed and are freely available at www.acsu.buffalo.edu/~yijunsun/lab/AsMac.html. Supplementary information Supplementary data are available at Bioinformatics online.

Bladder cancer is a biologically heterogeneous disease with variable clinical presentations, outcomes and responses to therapy. Thus, the clinical utility of single biomarkers for the detection and prediction of biological behavior of bladder cancer is limited. We have previously identified and validated a bladder cancer diagnostic signature composed of 10 biomarkers, which has been incorporated into a multiplex immunoassay bladder cancer test, Oncuria™. In this study, we evaluate whether these 10 biomarkers can assist in the prediction of bladder cancer clinical outcomes. Tumor gene expression and patient survival data from bladder cancer cases from The Cancer Genome Atlas (TCGA) were analyzed. Alignment between the mRNA expression of 10 biomarkers and the TCGA 2017 subtype classification was assessed. Kaplan–Meier analysis of multiple gene expression datasets indicated that high expression of the combined 10 biomarkers correlated with a significant reduction in overall survival. The analysis of three independent, publicly available gene expression datasets confirmed that multiplex prognostic models outperformed single biomarkers. In total, 8 of the 10 biomarkers from the Oncuria™ test were significantly associated with either luminal or basal molecular subtypes, and thus, the test has the potential to assist in the prediction of clinical outcome.

BACKGROUND: Intravesical Bacillus Calmette-Guerin (BCG), a live attenuated tuberculosis vaccine that acts as a non-specific immune system stimulant, is the most effective adjuvant treatment for patients with intermediate or high-risk non-muscle-invasive bladder cancer (NMIBC). However, to date, there are no reliable tests that are predictive of BCG treatment response. In this study, we evaluated the performance of OncuriaTM, a bladder cancer detection test, to predict response to intravesical BCG. METHODS: OncuriaTM data was evaluated in voided urine samples obtained from a prospectively collected cohort of 64 subjects with intermediate or high risk NMIBC prior to treatment with intravesical BCG. The OncuriaTM test, which measures 10 cancer-associated biomarkers was performed in an independent clinical laboratory. The ability of the test to identify those patients in whom BCG is ineffective against tumor recurrence was tested. Predictive models were derived using supervised learning and cross-validation analyses. Model performance was assessed using ROC curves. RESULTS: Pre-treatment urinary concentrations of MMP9, VEGFA, CA9, SDC1, PAI1, APOE, A1AT, ANG and MMP10 were increased in patients who developed disease recurrence. A combinatorial predictive model of treatment outcome achieved an AUROC 0.89 [95% CI: 0.80–0.99], outperforming any single biomarker, with a test sensitivity of 81.8% and a specificity of 84.9%. Hazard ratio analysis revealed that patients with higher urinary levels of ANG, CA9 and MMP10 had a significantly higher risk of disease recurrence. CONCLUSIONS: Monitoring the urinary levels of a cancer-associated biomarker panel enabled the discrimination of patients who did not respond to intravesical BCG therapy. With further study, the multiplex OncuriaTM test may be applicable for the clinical evaluation of bladder cancer patients considering intravesical BCG treatment.

Background: The Oncuria™ urine test for the detection of bladder cancer measures a multiplex protein signature. In this study, we investigated the influence of urinary cellularity, protein, and hematuria on the performance of the Oncuria™ test in an ex vivo experimental model. Materials and Methods: Pooled urine from healthy subjects was spiked with cultured benign (UROtsa) or malignant cells (T24), cellular proteins, or whole blood. The resulting samples were analyzed using the Oncuria™ test following the manufacturer’s instructions. Results: Urine samples obtained from healthy subjects were negative for bladder cancer by Oncuria™ test criteria. The majority of the manipulated conditions did not result in a false-positive test. The addition of whole blood (high concentration) did result in a false-positive result, but this was abrogated by sample centrifugation prior to analysis. The addition of cellular proteins (high concentration) resulted in a positive Oncuria™ test, and this was unaffected by pre-analysis sample centrifugation. Conclusions: The Oncuria™ multiplex test performed well in the ex vivo experimental model and shows promise for clinical application. The identification of patients who require additional clinical evaluation could reduce the need to subject patients who do not have bladder cancer to frequent, uncomfortable and expensive cystoscopic examinations, thus benefiting both patients and the healthcare system.