PancanQTLv2.0: a comprehensive resource for expression quantitative trait loci across human cancers.

Expression quantitative trait locus (eQTL) analysis is a powerful tool used to investigate genetic variations in complex diseases, including cancer. We previously developed a comprehensive database, PancanQTL, to characterize cancer eQTLs using The Cancer Genome Atlas (TCGA) dataset, and linked eQTLs with patient survival and GWAS risk variants. Here, we present an updated version, PancanQTLv2.0 (https://hanlaboratory.com/PancanQTLv2/), with advancements in fine-mapping causal variants for eQTLs, updating eQTLs overlapping with GWAS linkage disequilibrium regions and identifying eQTLs associated with drug response and immune infiltration. Through fine-mapping analysis, we identified 58 747 fine-mapped eQTLs credible sets, providing mechanic insights of gene regulation in cancer. We further integrated the latest GWAS Catalog and identified a total of 84 592 135 linkage associations between eQTLs and the existing GWAS loci, which represents a remarkable ∼50-fold increase compared to the previous version. Additionally, PancanQTLv2.0 uncovered 659516 associations between eQTLs and drug response and identified 146948 associations between eQTLs and immune cell abundance, providing potentially clinical utility of eQTLs in cancer therapy. PancanQTLv2.0 expanded the resources available for investigating gene expression regulation in human cancers, leading to advancements in cancer research and precision oncology.

FLUXestimator: a webserver for predicting metabolic flux and variations using transcriptomics data.

Quantitative assessment of single cell fluxome is critical for understanding the metabolic heterogeneity in diseases. Unfortunately, laboratory-based single cell fluxomics is currently impractical, and the current computational tools for flux estimation are not designed for single cell-level prediction. Given the well-established link between transcriptomic and metabolomic profiles, leveraging single cell transcriptomics data to predict single cell fluxome is not only feasible but also an urgent task. In this study, we present FLUXestimator, an online platform for predicting metabolic fluxome and variations using single cell or general transcriptomics data of large sample-size. The FLUXestimator webserver implements a recently developed unsupervised approach called single cell flux estimation analysis (scFEA), which uses a new neural network architecture to estimate reaction rates from transcriptomics data. To the best of our knowledge, FLUXestimator is the first web-based tool dedicated to predicting cell-/sample-wise metabolic flux and metabolite variations using transcriptomics data of human, mouse and 15 other common experimental organisms. The FLUXestimator webserver is available at http://scFLUX.org/, and stand-alone tools for local use are available at https://github.com/changwn/scFEA. Our tool provides a new avenue for studying metabolic heterogeneity in diseases and has the potential to facilitate the development of new therapeutic strategies.

A graph neural network model to estimate cell-wise metabolic flux using single-cell RNA-seq data.

The metabolic heterogeneity and metabolic interplay between cells are known as significant contributors to disease treatment resistance. However, with the lack of a mature high-throughput single-cell metabolomics technology, we are yet to establish systematic understanding of the intra-tissue metabolic heterogeneity and cooperative mechanisms. To mitigate this knowledge gap, we developed a novel computational method, namely, single-cell flux estimation analysis (scFEA), to infer the cell-wise fluxome from single-cell RNA-sequencing (scRNA-seq) data. scFEA is empowered by a systematically reconstructed human metabolic map as a factor graph, a novel probabilistic model to leverage the flux balance constraints on scRNA-seq data, and a novel graph neural network-based optimization solver. The intricate information cascade from transcriptome to metabolome was captured using multilayer neural networks to capitulate the nonlinear dependency between enzymatic gene expressions and reaction rates. We experimentally validated scFEA by generating an scRNA-seq data set with matched metabolomics data on cells of perturbed oxygen and genetic conditions. Application of scFEA on this data set showed the consistency between predicted flux and the observed variation of metabolite abundance in the matched metabolomics data. We also applied scFEA on five publicly available scRNA-seq and spatial transcriptomics data sets and identified context- and cell group-specific metabolic variations. The cell-wise fluxome predicted by scFEA empowers a series of downstream analyses including identification of metabolic modules or cell groups that share common metabolic variations, sensitivity evaluation of enzymes with regards to their impact on the whole metabolic flux, and inference of cell-tissue and cell-cell metabolic communications.

A generalized phase 1-2-3 design integrating dose optimization with confirmatory treatment comparison.

A generalized phase 1-2-3 design, Gen 1-2-3, that includes all phases of clinical treatment evaluation is proposed. The design extends and modifies the design of Chapple and Thall (2019), denoted by CT. Both designs begin with a phase 1-2 trial including dose acceptability and optimality criteria, and both select an optimal dose for phase 3. The Gen 1-2-3 design has the following key differences. In stage 1, it uses phase 1-2 criteria to identify a set of candidate doses rather than 1 dose. In stage 2, which is intermediate between phase 1-2 and phase 3, it randomizes additional patients fairly among the candidate doses and an active control treatment arm and uses survival time data from both stage 1 and stage 2 patients to select an optimal dose. It then makes a Go/No Go decision of whether or not to conduct phase 3 based on the predictive probability that the selected optimal dose will provide a specified substantive improvement in survival time over the control. A simulation study shows that the Gen 1-2-3 design has desirable operating characteristics compared to the CT design and 2 conventional designs.

Adaptive phase I-II clinical trial designs identifying optimal biological doses for targeted agents and immunotherapies.

Targeted agents and immunotherapies have revolutionized cancer treatment, offering promising options for various cancer types. Unlike traditional therapies the principle of "more is better" is not always applicable to these new therapies due to their unique biomedical mechanisms. As a result, various phase I-II clinical trial designs have been proposed to identify the optimal biological dose that maximizes the therapeutic effect of targeted therapies and immunotherapies by jointly monitoring both efficacy and toxicity outcomes. This review article examines several innovative phase I-II clinical trial designs that utilize accumulated efficacy and toxicity outcomes to adaptively determine doses for subsequent patients and identify the optimal biological dose, maximizing the overall therapeutic effect. Specifically, we highlight three categories of phase I-II designs: efficacy-driven, utility-based, and designs incorporating multiple efficacy endpoints. For each design, we review the dose-outcome model, the definition of the optimal biological dose, the dose-finding algorithm, and the software for trial implementation. To illustrate the concepts, we also present two real phase I-II trial examples utilizing the EffTox and ISO designs. Finally, we provide a classification tree to summarize the designs discussed in this article.

T3 + 3: 3 + 3 Design With Delayed Outcomes.

Delayed outcome is common in phase I oncology clinical trials. It causes logistic difficulty, wastes resources, and prolongs the trial duration. This article investigates this issue and proposes the time-to-event 3 + 3 (T3 + 3) design, which utilizes the actual follow-up time for at-risk patients with pending toxicity outcomes. The T3 + 3 design allows continuous accrual without unnecessary trial suspension and is costless and implementable with pretabulated dose decision rules. Besides, the T3 + 3 design uses the isotonic regression to estimate the toxicity rates across dose levels and therefore can accommodate for any targeted toxicity rate for maximum tolerated dose (MTD). It dramatically facilitates the trial preparation and conduct without intensive computation and statistical consultation. The extension to other algorithm-based phase I dose-finding designs (e.g., i3 + 3 design) is also studied. Comprehensive computer simulation studies are conducted to investigate the performance of the T3 + 3 design under various dose-toxicity scenarios. The results confirm that the T3 + 3 design substantially shortens the trial duration compared with the conventional 3 + 3 design and yields much higher accuracy in MTD identification than the rolling six design. In summary, the T3 + 3 design addresses the delayed outcome issue while keeping the desirable features of the 3 + 3 design, such as simplicity, transparency, and costless implementation. It has great potential to accelerate early-phase drug development.

Characterization of Proteoform Post-Translational Modifications by Top-Down and Bottom-Up Mass Spectrometry in Conjunction with Annotations.

Many proteoforms can be produced from a gene due to genetic mutations, alternative splicing, post-translational modifications (PTMs), and other variations. PTMs in proteoforms play critical roles in cell signaling, protein degradation, and other biological processes. Mass spectrometry (MS) is the primary technique for investigating PTMs in proteoforms, and two alternative MS approaches, top-down and bottom-up, have complementary strengths. The combination of the two approaches has the potential to increase the sensitivity and accuracy in PTM identification and characterization. In addition, protein and PTM knowledge bases, such as UniProt, provide valuable information for PTM characterization and verification. Here, we present a software pipeline PTM-TBA (PTM characterization by Top-down and Bottom-up MS and Annotations) for identifying and localizing PTMs in proteoforms by integrating top-down and bottom-up MS as well as PTM annotations. We assessed PTM-TBA using a technical triplicate of bottom-up and top-down MS data of SW480 cells. On average, database search of the top-down MS data identified 2000 mass shifts, 814.5 (40.7%) of which were matched to 11 common PTMs and 423 of which were localized. Of the mass shifts identified by top-down MS, PTM-TBA verified 435 mass shifts using the bottom-up MS data and UniProt annotations.

TopFD: A Proteoform Feature Detection Tool for Top-Down Proteomics.

Top-down liquid chromatography-mass spectrometry (LC-MS) analyzes intact proteoforms and generates mass spectra containing peaks of proteoforms with various isotopic compositions, charge states, and retention times. An essential step in top-down MS data analysis is proteoform feature detection, which aims to group these peaks into peak sets (features), each containing all peaks of a proteoform. Accurate protein feature detection enhances the accuracy in MS-based proteoform identification and quantification. Here, we present TopFD, a software tool for top-down MS feature detection that integrates algorithms for proteoform feature detection, feature boundary refinement, and machine learning models for proteoform feature evaluation. We performed extensive benchmarking of TopFD, ProMex, FlashDeconv, and Xtract using seven top-down MS data sets and demonstrated that TopFD outperforms other tools in feature accuracy, reproducibility, and feature abundance reproducibility.

Supervised clustering of high-dimensional data using regularized mixture modeling.

Identifying relationships between genetic variations and their clinical presentations has been challenged by the heterogeneous causes of a disease. It is imperative to unveil the relationship between the high-dimensional genetic manifestations and the clinical presentations, while taking into account the possible heterogeneity of the study subjects.We proposed a novel supervised clustering algorithm using penalized mixture regression model, called component-wise sparse mixture regression (CSMR), to deal with the challenges in studying the heterogeneous relationships between high-dimensional genetic features and a phenotype. The algorithm was adapted from the classification expectation maximization algorithm, which offers a novel supervised solution to the clustering problem, with substantial improvement on both the computational efficiency and biological interpretability. Experimental evaluation on simulated benchmark datasets demonstrated that the CSMR can accurately identify the subspaces on which subset of features are explanatory to the response variables, and it outperformed the baseline methods. Application of CSMR on a drug sensitivity dataset again demonstrated the superior performance of CSMR over the others, where CSMR is powerful in recapitulating the distinct subgroups hidden in the pool of cell lines with regards to their coping mechanisms to different drugs. CSMR represents a big data analysis tool with the potential to resolve the complexity of translating the clinical representations of the disease to the real causes underpinning it. We believe that it will bring new understanding to the molecular basis of a disease and could be of special relevance in the growing field of personalized medicine.

SSMD: a semi-supervised approach for a robust cell type identification and deconvolution of mouse transcriptomics data.

Deconvolution of mouse transcriptomic data is challenged by the fact that mouse models carry various genetic and physiological perturbations, making it questionable to assume fixed cell types and cell type marker genes for different data set scenarios. We developed a Semi-Supervised Mouse data Deconvolution (SSMD) method to study the mouse tissue microenvironment. SSMD is featured by (i) a novel nonparametric method to discover data set-specific cell type signature genes; (ii) a community detection approach for fixing cell types and their marker genes; (iii) a constrained matrix decomposition method to solve cell type relative proportions that is robust to diverse experimental platforms. In summary, SSMD addressed several key challenges in the deconvolution of mouse tissue data, including: (i) varied cell types and marker genes caused by highly divergent genotypic and phenotypic conditions of mouse experiment; (ii) diverse experimental platforms of mouse transcriptomics data; (iii) small sample size and limited training data source and (iv) capable to estimate the proportion of 35 cell types in blood, inflammatory, central nervous or hematopoietic systems. In silico and experimental validation of SSMD demonstrated its high sensitivity and accuracy in identifying (sub) cell types and predicting cell proportions comparing with state-of-the-arts methods. A user-friendly R package and a web server of SSMD are released via https://github.com/xiaoyulu95/SSMD.

The INGENIOUS trial: Impact of pharmacogenetic testing on adverse events in a pragmatic clinical trial.

Adverse drug events (ADEs) account for a significant mortality, morbidity, and cost burden. Pharmacogenetic testing has the potential to reduce ADEs and inefficacy. The objective of this INGENIOUS trial (NCT02297126) analysis was to determine whether conducting and reporting pharmacogenetic panel testing impacts ADE frequency. The trial was a pragmatic, randomized controlled clinical trial, adapted as a propensity matched analysis in individuals (N = 2612) receiving a new prescription for one or more of 26 pharmacogenetic-actionable drugs across a community safety-net and academic health system. The intervention was a pharmacogenetic testing panel for 26 drugs with dosage and selection recommendations returned to the health record. The primary outcome was occurrence of ADEs within 1 year, according to modified Common Terminology Criteria for Adverse Events (CTCAE). In the propensity-matched analysis, 16.1% of individuals experienced any ADE within 1-year. Serious ADEs (CTCAE level ≥ 3) occurred in 3.2% of individuals. When combining all 26 drugs, no significant difference was observed between the pharmacogenetic testing and control arms for any ADE (Odds ratio 0.96, 95% CI: 0.78-1.18), serious ADEs (OR: 0.91, 95% CI: 0.58-1.40), or mortality (OR: 0.60, 95% CI: 0.28-1.21). However, sub-group analyses revealed a reduction in serious ADEs and death in individuals who underwent pharmacogenotyping for aripiprazole and serotonin or serotonin-norepinephrine reuptake inhibitors (OR 0.34, 95% CI: 0.12-0.85). In conclusion, no change in overall ADEs was observed after pharmacogenetic testing. However, limitations incurred during INGENIOUS likely affected the results. Future studies may consider preemptive, rather than reactive, pharmacogenetic panel testing.

Generalized phase I-II designs to increase long term therapeutic success rate.

Designs for early phase dose finding clinical trials typically are either phase I based on toxicity, or phase I-II based on toxicity and efficacy. These designs rely on the implicit assumption that the dose of an experimental agent chosen using these short-term outcomes will maximize the agent's long-term therapeutic success rate. In many clinical settings, this assumption is not true. A dose selected in an early phase oncology trial may give suboptimal progression-free survival or overall survival time, often due to a high rate of relapse following response. To address this problem, a new family of Bayesian generalized phase I-II designs is proposed. First, a conventional phase I-II design based on short-term outcomes is used to identify a set of candidate doses, rather than selecting one dose. Additional patients then are randomized among the candidates, patients are followed for a predefined longer time period, and a final dose is selected to maximize the long-term therapeutic success rate, defined in terms of duration of response. Dose-specific sample sizes in the randomization are determined adaptively to obtain a desired level of selection reliability. The design was motivated by a phase I-II trial to find an optimal dose of natural killer cells as targeted immunotherapy for recurrent or treatment-resistant B-cell hematologic malignancies. A simulation study shows that, under a range of scenarios in the context of this trial, the proposed design has much better performance than two conventional phase I-II designs.

A Bayesian phase I/II design to determine subgroup-specific optimal dose for immunotherapy sequentially combined with radiotherapy.

Sequential administration of immunotherapy following radiotherapy (immunoRT) has attracted much attention in cancer research. Due to its unique feature that radiotherapy upregulates the expression of a predictive biomarker for immunotherapy, novel clinical trial designs are needed for immunoRT to identify patient subgroups and the optimal dose for each subgroup. In this article, we propose a Bayesian phase I/II design for immunotherapy administered after standard-dose radiotherapy for this purpose. We construct a latent subgroup membership variable and model it as a function of the baseline and pre-post radiotherapy change in the predictive biomarker measurements. Conditional on the latent subgroup membership of each patient, we jointly model the continuous immune response and the binary efficacy outcome using plateau models, and model toxicity using the equivalent toxicity score approach to account for toxicity grades. During the trial, based on accumulating data, we continuously update model estimates and adaptively randomize patients to admissible doses. Simulation studies and an illustrative trial application show that our design has good operating characteristics in terms of identifying both patient subgroups and the optimal dose for each subgroup.

Classification of Fluorescently Labelled Maize Kernels Using Convolutional Neural Networks.

Accurate real-time classification of fluorescently labelled maize kernels is important for the industrial application of its advanced breeding techniques. Therefore, it is necessary to develop a real-time classification device and recognition algorithm for fluorescently labelled maize kernels. In this study, a machine vision (MV) system capable of identifying fluorescent maize kernels in real time was designed using a fluorescent protein excitation light source and a filter to achieve optimal detection. A high-precision method for identifying fluorescent maize kernels based on a YOLOv5s convolutional neural network (CNN) was developed. The kernel sorting effects of the improved YOLOv5s model, as well as other YOLO models, were analysed and compared. The results show that using a yellow LED light as an excitation light source combined with an industrial camera filter with a central wavelength of 645 nm achieves the best recognition effect for fluorescent maize kernels. Using the improved YOLOv5s algorithm can increase the recognition accuracy of fluorescent maize kernels to 96%. This study provides a feasible technical solution for the high-precision, real-time classification of fluorescent maize kernels and has universal technical value for the efficient identification and classification of various fluorescently labelled plant seeds.

Evaluation of Machine Learning Models for Proteoform Retention and Migration Time Prediction in Top-Down Mass Spectrometry.

Reversed-phase liquid chromatography (RPLC) and capillary zone electrophoresis (CZE) are two primary proteoform separation methods in mass spectrometry (MS)-based top-down proteomics. Proteoform retention time (RT) prediction in RPLC and migration time (MT) prediction in CZE provide additional information for accurate proteoform identification and quantification. While existing methods are mainly focused on peptide RT and MT prediction in bottom-up MS, there is still a lack of methods for proteoform RT and MT prediction in top-down MS. We systematically evaluated eight machine learning models and a transfer learning method for proteoform RT prediction and five models and the transfer learning method for proteoform MT prediction. Experimental results showed that a gated recurrent unit (GRU)-based model with transfer learning achieved a high accuracy (R = 0.978) for proteoform RT prediction and that the GRU-based model and a fully connected neural network model obtained a high accuracy of R = 0.982 and 0.981 for proteoform MT prediction, respectively.

Local and distant brain control in melanoma and NSCLC brain metastases with concurrent radiosurgery and immune checkpoint inhibition.

INTRODUCTION: The treatment of brain metastases with stereotactic radiosurgery (SRS) in combination with immune checkpoint inhibitors (ICI) has become more common in recent years, but there is a lack of prospective data on cancer control outcomes when these therapies are administered concurrently. METHODS: Data were retrospectively reviewed for patients with non-small cell lung cancer (NSCLC) and melanoma brain metastases treated with SRS at a single institution from May 2008 to January 2017. A parametric proportional hazard model is used to detect the effect of concurrent ICI within 30, 60, or 90 days of ICI administration on local control and distant in-brain control. Other patient and lesion characteristics are treated as covariates and adjusted in the regression. A frailty term is added in the baseline hazard to capture the within-patient correlation. RESULTS: We identified 144 patients with 477 total lesions, including 95 NSCLC patients (66.0%), and 49 (34.0%) melanoma patients. On multivariate analysis, concurrent SRS and ICI (SRS within 30 days of ICI administration) was not associated with local control but was associated with distant brain control. When controlling for prior treatment to lesion, number of lesions, and presence of extracranial metastases, patients receiving this combination had a statistically significant decrease in distant brain failure compared to patients that received non-concurrent ICI or no ICI (HR 0.15; 95% CI 0.05-0.47, p = 0.0011). CONCLUSION: Concurrent ICI can enhance the efficacy of SRS. Prospective studies would allow for stronger evidence to support the impact of concurrent SRS and ICI on disease outcomes.

Evaluating risk for alcohol use disorder: Polygenic risk scores and family history.

BACKGROUND: Early identification of individuals at high risk for alcohol use disorder (AUD) coupled with prompt interventions could reduce the incidence of AUD. In this study, we investigated whether Polygenic Risk Scores (PRS) can be used to evaluate the risk for AUD and AUD severity (as measured by the number of DSM-5 AUD diagnostic criteria met) and compared their performance with a measure of family history of AUD. METHODS: We studied individuals of European ancestry from the Collaborative Study on the Genetics of Alcoholism (COGA). DSM-5 diagnostic criteria were available for 7203 individuals, of whom 3451 met criteria for DSM-IV alcohol dependence or DSM-5 AUD and 1616 were alcohol-exposed controls aged ≥21 years with no history of AUD or drug dependence. Further, 4842 individuals had a positive first-degree family history of AUD (FH+), 2722 had an unknown family history (FH?), and 336 had a negative family history (FH-). PRS were derived from a meta-analysis of a genome-wide association study of AUD from the Million Veteran Program and scores from the problem subscale of the Alcohol Use Disorders Identification Test in the UK Biobank. We used mixed models to test the association between PRS and risk for AUD and AUD severity. RESULTS: AUD cases had higher PRS than controls with PRS increasing as the number of DSM-5 diagnostic criteria increased (p-values ≤ 1.85E-05 ) in the full COGA sample, the FH+ subsample, and the FH? subsample. Individuals in the top decile of PRS had odds ratios (OR) for developing AUD of 1.96 (95% CI: 1.54 to 2.51, p-value = 7.57E-08 ) and 1.86 (95% CI: 1.35 to 2.56, p-value = 1.32E-04 ) in the full sample and the FH+ subsample, respectively. These values are comparable to previously reported ORs for a first-degree family history (1.91 to 2.38) estimated from national surveys. PRS were also significantly associated with the DSM-5 AUD diagnostic criterion count in the full sample, the FH+ subsample, and the FH? subsample (p-values ≤6.7E-11 ). PRS remained significantly associated with AUD and AUD severity after accounting for a family history of AUD (p-values ≤6.8E-10 ). CONCLUSIONS: Both PRS and family history were associated with AUD and AUD severity, indicating that these risk measures assess distinct aspects of liability to AUD traits.

BIPSE: A biomarker-based phase I/II design for immunotherapy trials with progression-free survival endpoint.

A Bayesian biomarker-based phase I/II design (BIPSE) is presented for immunotherapy trials with a progression-free survival (PFS) endpoint. The objective is to identify the subgroup-specific optimal dose, defined as the dose with the best risk-benefit tradeoff in each biomarker subgroup. We jointly model the immune response, toxicity outcome, and PFS with information borrowing across subgroups. A plateau model is used to describe the marginal distribution of the immune response. Conditional on the immune response, we model toxicity using probit regression and model PFS using the mixture cure rate model. During the trial, based on the accumulating data, we continuously update model estimates and adaptively randomize patients to doses with high desirability within each subgroup. Simulation studies show that the BIPSE design has desirable operating characteristics in selecting the subgroup-specific optimal doses and allocating patients to those optimal doses, and outperforms conventional designs.

Diversity synthesis of indole derivatives via catalyst control cyclization reaction of 2-indolylmethanols and azonaphthalene.

A series of indole-fused scaffolds and derivatives was synthesized via the cyclization reaction of 2-indolylmethanols with azonaphthalene. These reactions were realized under mild reaction conditions through catalyst control, providing structurally diverse indole derivatives with moderate to excellent yields. This protocol also shows good substrate adaptability, especially in six-membered ring products.

SCI: A Bayesian adaptive phase I/II dose-finding design accounting for semi-competing risks outcomes for immunotherapy trials.

An immunotherapy trial often uses the phase I/II design to identify the optimal biological dose, which monitors the efficacy and toxicity outcomes simultaneously in a single trial. The progression-free survival rate is often used as the efficacy outcome in phase I/II immunotherapy trials. As a result, patients developing disease progression in phase I/II immunotherapy trials are generally seriously ill and are often treated off the trial for ethical consideration. Consequently, the happening of disease progression will terminate the toxicity event but not vice versa, so the issue of the semi-competing risks arises. Moreover, this issue can become more intractable with the late-onset outcomes, which happens when a relatively long follow-up time is required to ascertain progression-free survival. This paper proposes a novel Bayesian adaptive phase I/II design accounting for semi-competing risks outcomes for immunotherapy trials, referred to as the dose-finding design accounting for semi-competing risks outcomes for immunotherapy trials (SCI) design. To tackle the issue of the semi-competing risks in the presence of late-onset outcomes, we re-construct the likelihood function based on each patient's actual follow-up time and develop a data augmentation method to efficiently draw posterior samples from a series of Beta-binomial distributions. We propose a concise curve-free dose-finding algorithm to adaptively identify the optimal biological dose using accumulated data without making any parametric dose-response assumptions. Numerical studies show that the proposed SCI design yields good operating characteristics in dose selection, patient allocation, and trial duration.