Approximated gene expression trajectories for gene regulatory network inference on cell tracks.

The study of pattern formation has benefited from our ability to reverse-engineer gene regulatory network (GRN) structure from spatiotemporal quantitative gene expression data. Traditional approaches have focused on systems where the timescales of pattern formation and morphogenesis can be separated. Unfortunately, this is not the case in most animal patterning systems, where pattern formation and morphogenesis are co-occurring and tightly linked. To elucidate patterning mechanisms in such systems we need to adapt our GRN inference methodologies to include cell movements. In this work, we fill this gap by integrating quantitative data from live and fixed embryos to approximate gene expression trajectories (AGETs) in single cells and use these to reverse-engineer GRNs. This framework generates candidate GRNs that recapitulate pattern at the tissue level, gene expression dynamics at the single cell level, recover known genetic interactions and recapitulate experimental perturbations while incorporating cell movements explicitly for the first time.

Single-cell transcriptome and crosstalk analysis reveals immune alterations and key pathways in the bone marrow of knee OA patients.

Knee osteoarthritis (OA) is a significant medical and economic burden. To understand systemic immune effects, we performed deep exploration of bone marrow aspirate concentrates (BMACs) from knee-OA patients via single-cell RNA sequencing and proteomic analyses from a randomized clinical trial (MILES: NCT03818737). We found significant cellular and immune alterations in the bone marrow, specifically in MSCs, T cells and NK cells, along with changes in intra-tissue cellular crosstalk during OA progression. Unlike previous studies focusing on injury sites or peripheral blood, our probe into the bone marrow-an inflammation and immune regulation hub-highlights remote organ impact of OA, identifying cell types and pathways for potential therapeutic targeting. Our findings highlight increased cellular senescence and inflammatory pathways, revealing key upstream genes, transcription factors, and ligands. Additionally, we identified significant enrichment in key biological pathways like PI3-AKT-mTOR signaling and IFN responses, showing their potentially crucial role in OA onset and progression.

Panera: An innovative framework for surmounting uncertainty in microbial community modeling using pan-genera metabolic models.

Utilization of 16S rRNA data in constraint-based modeling to characterize microbial communities confronts a major hurdle of lack of species-level resolution, impeding the construction of community models. We introduce "Panera," an innovative framework designed to model communities under this uncertainty and yet perform metabolic inferences using pan-genus metabolic models (PGMMs). We demonstrated PGMMs' utility for comprehending the metabolic capabilities of a genus and in characterizing community models using amplicon data. The unique, adaptable nature of PGMMs unlocks their potential in building hybrid communities, combining genome-scale metabolic models (GSMMs) and PGMMs. Notably, these models provide predictions comparable to the standard GSMM-based community models, while achieving a nearly 46% reduction in error compared to the genus model-based communities. In essence, "Panera" presents a potent and effective approach to aid in metabolic modeling by enabling robust predictions of community metabolic potential when dealing with amplicon data, and offers insights into genus-level metabolic landscapes.

Module control of network analysis in psychopathology.

The network approach to characterizing psychopathology departs from traditional latent categorical and dimensional approaches. Causal interplay among symptoms contributed to dynamic psychopathology system. Therefore, analyzing the symptom clusters is critical for understanding mental disorders. Furthermore, despite extensive research studying the topological features of symptom networks, the control relationships between symptoms remain largely unclear. Here, we present a novel systematizing concept, module control, to analyze the control principle of the symptom network at a module level. We introduce Module Control Network (MCN) to identify key modules that regulate the network's behavior. By applying our approach to a multivariate psychological dataset, we discover that non-emotional modules, such as sleep-related and stress-related modules, are the primary controlling modules in the symptom network. Our findings indicate that module control can expose central symptom cluster governing psychopathology network, offering novel insights into the underlying mechanisms of mental disorders and individualized approach to psychological interventions.

Assessment of circulating proteins in thyroid cancer: Proteome-wide Mendelian randomization and colocalization analysis.

The causality between circulating proteins and thyroid cancer (TC) remains unclear. We employed five large-scale circulating proteomic genome-wide association studies (GWASs) with up to 100,000 participants and a TC meta-GWAS (nCase = 3,418, nControl = 292,703) to conduct proteome-wide Mendelian randomization (MR) and Bayesian colocalization analysis. Protein and gene expressions were validated in thyroid tissue. Through MR analysis, we identified 26 circulating proteins with a putative causal relationship with TCs, among which NANS protein passed multiple corrections (P BH = 3.28e-5, 0.05/1,525). These proteins were involved in amino acids and organic acid synthesis pathways. Colocalization analysis further identified six proteins associated with TCs (VCAM1, LGMN, NPTX1, PLEKHA7, TNFAIP3, and BMP1). Tissue validation confirmed BMP1, LGMN, and PLEKHA7's differential expression between normal and TC tissues. We found limited evidence for linking circulating proteins and the risk of TCs. Our study highlighted the contribution of proteins, particularly those involved in amino acid metabolism, to TCs.

Decoupling the correlation between cytotoxic and exhausted T lymphocyte states enhances melanoma immunotherapy response prediction.

Cytotoxic T lymphocyte (CTL) and terminal exhausted T lymphocyte (ETL) activities crucially influence immune checkpoint inhibitor (ICI) response. Despite this, the efficacy of ETL and CTL transcriptomic signatures for response prediction remains limited. Investigating this across the TCGA and publicly available single-cell cohorts, we find a strong positive correlation between ETL and CTL expression signatures in most cancers. We hence posited that their limited predictability arises due to their mutually canceling effects on ICI response. Thus, we developed DETACH, a computational method to identify a gene set whose expression pinpoints to a subset of melanoma patients where the CTL and ETL correlation is low. DETACH enhances CTL's prediction accuracy, outperforming existing signatures. DETACH signature genes activity also demonstrates a positive correlation with lymphocyte infiltration and the prevalence of reactive T cells in the tumor microenvironment (TME), advancing our understanding of the CTL cell state within the TME.

A statistical approach for identifying single nucleotide variants that affect transcription factor binding.

Non-coding variants located within regulatory elements may alter gene expression by modifying transcription factor (TF) binding sites, thereby leading to functional consequences. Different TF models are being used to assess the effect of DNA sequence variants, such as single nucleotide variants (SNVs). Often existing methods are slow and do not assess statistical significance of results. We investigated the distribution of absolute maximal differential TF binding scores for general computational models that affect TF binding. We find that a modified Laplace distribution can adequately approximate the empirical distributions. A benchmark on in vitro and in vivo datasets showed that our approach improves upon an existing method in terms of performance and speed. Applications on eQTLs and on a genome-wide association study illustrate the usefulness of our statistics by highlighting cell type-specific regulators and target genes. An implementation of our approach is freely available on GitHub and as bioconda package.

ProteoMixture: A cell type deconvolution tool for bulk tissue proteomic data.

Numerous multi-omic investigations of cancer tissue have documented varying and poor pairwise transcript:protein quantitative correlations, and most deconvolution tools aiming to predict cell type proportions (cell admixture) have been developed and credentialed using transcript-level data alone. To estimate cell admixture using protein abundance data, we analyzed proteome and transcriptome data generated from contrived admixtures of tumor, stroma, and immune cell models or those selectively harvested from the tissue microenvironment by laser microdissection from high grade serous ovarian cancer (HGSOC) tumors. Co-quantified transcripts and proteins performed similarly to estimate stroma and immune cell admixture (r ≥ 0.63) in two commonly used deconvolution algorithms, ESTIMATE or ConsensusTME. We further developed and optimized protein-based signatures estimating cell admixture proportions and benchmarked these using bulk tumor proteomic data from over 150 patients with HGSOC. The optimized protein signatures supporting cell type proportion estimates from bulk tissue proteomic data are available at https://lmdomics.org/ProteoMixture/.

scGREAT: Transformer-based deep-language model for gene regulatory network inference from single-cell transcriptomics.

Gene regulatory networks (GRNs) involve complex and multi-layer regulatory interactions between regulators and their target genes. Precise knowledge of GRNs is important in understanding cellular processes and molecular functions. Recent breakthroughs in single-cell sequencing technology made it possible to infer GRNs at single-cell level. Existing methods, however, are limited by expensive computations, and sometimes simplistic assumptions. To overcome these obstacles, we propose scGREAT, a framework to infer GRN using gene embeddings and transformer from single-cell transcriptomics. scGREAT starts by constructing gene expression and gene biotext dictionaries from scRNA-seq data and gene text information. The representation of TF gene pairs is learned through optimizing embedding space by transformer-based engine. Results illustrated scGREAT outperformed other contemporary methods on benchmarks. Besides, gene representations from scGREAT provide valuable gene regulation insights, and external validation on spatial transcriptomics illuminated the mechanism behind scGREAT annotation. Moreover, scGREAT identified several TF target regulations corroborated in studies.

Safeguarding authenticity for mitigating the harms of generative AI: Issues, research agenda, and policies for detection, fact-checking, and ethical AI.

As the influence of transformer-based approaches in general and generative artificial intelligence (AI) in particular continues to expand across various domains, concerns regarding authenticity and explainability are on the rise. Here, we share our perspective on the necessity of implementing effective detection, verification, and explainability mechanisms to counteract the potential harms arising from the proliferation of AI-generated inauthentic content and science. We recognize the transformative potential of generative AI, exemplified by ChatGPT, in the scientific landscape. However, we also emphasize the urgency of addressing associated challenges, particularly in light of the risks posed by disinformation, misinformation, and unreproducible science. This perspective serves as a response to the call for concerted efforts to safeguard the authenticity of information in the age of AI. By prioritizing detection, fact-checking, and explainability policies, we aim to foster a climate of trust, uphold ethical standards, and harness the full potential of AI for the betterment of science and society.

FOTF-CPI: A compound-protein interaction prediction transformer based on the fusion of optimal transport fragments.

Compound-protein interaction (CPI) affinity prediction plays an important role in reducing the cost and time of drug discovery. However, the interpretability of how fragments function in CPI is impacted by the fact that current methods ignore the affinity relationships between fragments of compounds and fragments of proteins in CPI modeling. This article introduces an improved Transformer called FOTF-CPI (a Fusion of Optimal Transport Fragments compound-protein interaction prediction model). We use an optimal transport-based fragmentation approach to improve the model's understanding of compound and protein sequences. Additionally, a fused attention mechanism is employed, which combines the features of fragments to capture full affinity information. This fused attention redistributes higher attention scores to fragments with higher affinity. Experimental results show FOTF-CPI achieves an average 2% higher performance than other models on all three datasets. Furthermore, the visualization confirms the potential of FOTF-CPI for drug discovery applications.

HLAEquity: Examining biases in pan-allele peptide-HLA binding predictors.

Peptide-HLA (pHLA) binding prediction is essential in screening peptide candidates for personalized peptide vaccines. Machine learning (ML) pHLA binding prediction tools are trained on vast amounts of data and are effective in screening peptide candidates. Most ML models report the ability to generalize to HLA alleles unseen during training ("pan-allele" models). However, the use of datasets with imbalanced allele content raises concerns about biased model performance. First, we examine the data bias of two ML-based pan-allele pHLA binding predictors. We find that the pHLA datasets overrepresent alleles from geographic populations of high-income countries. Second, we show that the identified data bias is perpetuated within ML models, leading to algorithmic bias and subpar performance for alleles expressed in low-income geographic populations. We draw attention to the potential therapeutic consequences of this bias, and we challenge the use of the term "pan-allele" to describe models trained with currently available public datasets.

In silico designing and characterization of outer membrane protein K (OmpK) from Vibrio anguillarum and its expression in Nicotiana tabacum for the development of a plant-based vaccine against fish vibriosis.

Vibriosis is caused by Vibrio anguillarum in various species of aquaculture. A novel, secure, and stable vaccine is needed to eradicate vibriosis. Here, for reverse vaccinology and plant-based expression, the outer membrane protein K (OmpK) of V. anguillarum was chosen due to its conserved nature in all Vibrio species. OmpK, an ideal vaccine candidate against vibriosis, demonstrated immunogenic, non-allergic, and non-toxic behavior by using various bioinformatics tools. Docking showed the interaction of the OmpK model with TLR-5. In comparison to costly platforms, plants can be used as alternative and economic bio-factories to produce vaccine antigens. We expressed OmpK antigen in Nicotiana tabacum using Agrobacterium-mediated transformation. The expression vector was constructed using Gateway® cloning. Transgene integration was verified by polymerase chain reaction (PCR), and the copy number via qRT-PCR, which showed two copies of transgenes. Western blotting detected monomeric form of OmpK protein. The total soluble protein (TSP) fraction of OmpK was equivalent to 0.38% as detected by ELISA. Mice and fish were immunized with plant-derived OmpK antigen, which showed a significantly high level of anti-OmpK antibodies. The present study is the first report of OmpK antigen expression in higher plants for the potential use as vaccine in aquaculture against vibriosis, which could provide protection against multiple Vibrio species due to the conserved nature OmpK antigen.

Epitranscriptome analysis of NAD-capped RNA by spike-in-based normalization and prediction of chronological age.

Nicotinamide adenine dinucleotide (NAD) can be used as an initiating nucleotide in RNA transcription to produce NAD-capped RNA (NAD-RNA). RNA modification by NAD that links metabolite with expressed transcript is a poorly studied epitranscriptomic modification. Current NAD-RNA profiling methods involve multi-steps of chemo-enzymatic labeling and affinity-based enrichment, thus presenting a critical analytical challenge to remove unwanted variations, particularly batch effects. Here, we propose a computational framework, enONE, to remove unwanted variations. We demonstrate that designed spike-in RNA, together with modular normalization procedures and evaluation metrics, can mitigate technical noise, empowering quantitative and comparative assessment of NAD-RNA across different datasets. Using enONE and a human aging cohort, we reveal age-associated features of NAD-capping and further develop an accurate RNA-based aging clock that combines signatures from both transcriptome and NAD-modified epitranscriptome. enONE facilitates the discovery of NAD-RNA responsive to physiological changes, laying an important foundation for functional investigations into this modification.

Gene-metabolite annotation with shortest reactional distance enhances metabolite genome-wide association studies results.

Metabolite genome-wide association studies (mGWAS) have advanced our understanding of the genetic control of metabolite levels. However, interpreting these associations remains challenging due to a lack of tools to annotate gene-metabolite pairs beyond the use of conservative statistical significance threshold. Here, we introduce the shortest reactional distance (SRD) metric, drawing from the comprehensive KEGG database, to enhance the biological interpretation of mGWAS results. We applied this approach to three independent mGWAS, including a case study on sickle cell disease patients. Our analysis reveals an enrichment of small SRD values in reported mGWAS pairs, with SRD values significantly correlating with mGWAS p values, even beyond the standard conservative thresholds. We demonstrate the utility of SRD annotation in identifying potential false negatives and inaccuracies within current metabolic pathway databases. Our findings highlight the SRD metric as an objective, quantitative and easy-to-compute annotation for gene-metabolite pairs, suitable to integrate statistical evidence to biological networks.

Review: Computational analysis of human skeletal remains in ancient DNA and forensic genetics.

Degraded DNA is used to answer questions in the fields of ancient DNA (aDNA) and forensic genetics. While aDNA studies typically center around human evolution and past history, and forensic genetics is often more concerned with identifying a specific individual, scientists in both fields face similar challenges. The overlap in source material has prompted periodic discussions and studies on the advantages of collaboration between fields toward mutually beneficial methodological advancements. However, most have been centered around wet laboratory methods (sampling, DNA extraction, library preparation, etc.). In this review, we focus on the computational side of the analytical workflow. We discuss limitations and considerations to consider when working with degraded DNA. We hope this review provides a framework to researchers new to computational workflows for how to think about analyzing highly degraded DNA and prompts an increase of collaboration between the forensic genetics and aDNA fields.

A data-driven computational model for obesity-driven diabetes onset and remission through weight loss.

Obesity is a major risk factor for the development of type 2 diabetes (T2D), where a sustained weight loss may result in T2D remission in individuals with obesity. To design effective and feasible intervention strategies to prevent or reverse T2D, it is imperative to study the progression of T2D and remission together. Unfortunately, this is not possible through experimental and observational studies. To address this issue, we introduce a data-driven computational model and use human data to investigate the progression of T2D with obesity and remission through weight loss on the same timeline. We identify thresholds for the emergence of T2D and necessary conditions for remission. We explain why remission is only possible within a window of opportunity and the way that window depends on the progression history of T2D, individual's metabolic state, and calorie restrictions. These findings can help to optimize therapeutic intervention strategies for T2D prevention or treatment.

Genetic overlap for ten cardiovascular diseases: A comprehensive gene-centric pleiotropic association analysis and Mendelian randomization study.

Recent studies suggest that pleiotropic effects may explain the genetic architecture of cardiovascular diseases (CVDs). We conducted a comprehensive gene-centric pleiotropic association analysis for ten CVDs using genome-wide association study (GWAS) summary statistics to identify pleiotropic genes and pathways that may underlie multiple CVDs. We found shared genetic mechanisms underlying the pathophysiology of CVDs, with over two-thirds of the diseases exhibiting common genes and single-nucleotide polymorphisms (SNPs). Significant positive genetic correlations were observed in more than half of paired CVDs. Additionally, we investigated the pleiotropic genes shared between different CVDs, as well as their functional pathways and distribution in different tissues. Moreover, six hub genes, including ALDH2, XPO1, HSPA1L, ESR2, WDR12, and RAB1A, as well as 26 targeted potential drugs, were identified. Our study provides further evidence for the pleiotropic effects of genetic variants on CVDs and highlights the importance of considering pleiotropy in genetic association studies.

A multi-center performance assessment for automated histopathological classification and grading of glioma using whole slide images.

Accurate pathological classification and grading of gliomas is crucial in clinical diagnosis and treatment. The application of deep learning techniques holds promise for automated histological pathology diagnosis. In this study, we collected 733 whole slide images from four medical centers, of which 456 were used for model training, 150 for internal validation, and 127 for multi-center testing. The study includes 5 types of common gliomas. A subtask-guided multi-instance learning image-to-label training pipeline was employed. The pipeline leveraged "patch prompting" for the model to converge with reasonable computational cost. Experiments showed that an overall accuracy of 0.79 in the internal validation dataset. The performance on the multi-center testing dataset showed an overall accuracy to 0.73. The findings suggest a minor yet acceptable performance decrease in multi-center data, demonstrating the model's strong generalizability and establishing a robust foundation for future clinical applications.

Preliminary study of all-solid-state batteries: Evaluation of blast formation during the thermal runaway.

All-solid-state batteries have been developed to increase energy density by replacing the lithiated graphite negative electrode by a lithium metal foil and to increase safety by removing the organic compounds. However, the safety issues of these batteries have received little attention up to now. The behavior of a reassembled all-solid-state battery under thermal stress was recorded by X-ray radiography and a high-speed camera. The thermal runaway (TR) lasted about 5 ms, thus extremely fast reaction kinetics. In comparison, the TR of a lithium-ion battery is about 500 ms. Furthermore, a 188-mbar aerial overpressure was measured using a piezoelectric sensor. Although this cell is not an explosive, 2.7 g TNT equivalent was calculated for it. This atypical behavior could have an impact on the casing or the battery pack. Therefore, it must be studied in greater detail.