Search

Summary: 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.

Summary: 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.

Summary: 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.

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. • Outer membrane protein K (OmpK) demonstrates immunogenic potential in silico. • OmpK antigen for immunization is produced in cost-effective plant bio-factories. • OmpK produced in plants induced a strong antibody response in mice and fish. • Vibriosis in aquaculture may be eradicated by cost-effective subunit vaccination. [ABSTRACT FROM AUTHOR]

Summary: 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 (PBH = 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.

Summary: 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.

Summary: 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.

Summary: 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/.

Summary: 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.

Summary: 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.

Summary: 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.

Summary: 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.

Summary: 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.

Summary: 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.

Summary: 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.

Summary: 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.

Summary: 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.

Summary: Small nucleolar RNA host genes (SNHGs) are a special family of long non-coding RNAs (lncRNAs), which not only function in a way similar to other lncRNAs but also influence the intracellular level of small nucleolar RNAs to modulate cancers. However, the features of SNHGs and their role in the prognosis and immunotherapeutic response of human cancer have not been explored. We found that SNHGs were commonly deregulated and correlated with patient survival in various cancers. The critical role of DNA methylation and somatic alterations on deregulation was also identified. SNHG family score was significantly associated with survival, multiple tumor characteristics, and tumor microenvironment. SNHG-related risk score could serve as a prognostic and immunotherapeutic response biomarker based on multiple databases. This study emphasizes the potential of SNHGs as biomarkers for prognosis and immunotherapeutic response, enabling further research into the immune regulatory mechanism and therapeutic potentials of SNHGs in cancer.

Summary: Extensive metabolic heterogeneity in breast cancers has limited the deployment of metabolic therapies. To enable patient stratification, we studied the metabolic landscape in breast cancers (∼3000 patients combined) and identified three subtypes with increasing degrees of metabolic deregulation. Subtype M1 was found to be dependent on bile-acid biosynthesis, whereas M2 showed reliance on methionine pathway, and M3 engaged fatty-acid, nucleotide, and glucose metabolism. The extent of metabolic alterations correlated strongly with tumor aggressiveness and patient outcome. This pattern was reproducible in independent datasets and using in vivo tumor metabolite data. Using machine-learning, we identified robust and generalizable signatures of metabolic subtypes in tumors and cell lines. Experimental inhibition of metabolic pathways in cell lines representing metabolic subtypes revealed subtype-specific sensitivity, therapeutically relevant drugs, and promising combination therapies. Taken together, metabolic stratification of breast cancers can thus aid in predicting patient outcome and designing precision therapies.

Summary: Machine learning (ML) has the potential to identify subsets of patients with distinct phenotypes from gene expression data. However, phenotype prediction using ML has often relied on identifying important genes without a systems biology context. To address this, we created an interpretable ML approach based on blood transcriptomics to predict phenotype in systemic lupus erythematosus (SLE), a heterogeneous autoimmune disease. We employed a sequential grouped feature importance algorithm to assess the performance of gene sets, including immune and metabolic pathways and cell types, known to be abnormal in SLE in predicting disease activity and organ involvement. Gene sets related to interferon, tumor necrosis factor, the mitoribosome, and T cell activation were the best predictors of phenotype with excellent performance. These results suggest potential relationships between the molecular pathways identified in each model and manifestations of SLE. This ML approach to phenotype prediction can be applied to other diseases and tissues.

Summary: Bile acid (BA) metabolism is a complex system that includes a wide variety of primary and secondary, as well as conjugated and unconjugated BAs that undergo continuous enterohepatic circulation (EHC). Alterations in both composition and dynamics of BAs have been associated with various diseases. However, a mechanistic understanding of the relationship between altered BA metabolism and related diseases is lacking. Computational modeling may support functional analyses of the physiological processes involved in the EHC of BAs along the gut-liver axis. In this study, we developed a physiologically based model of murine BA metabolism describing synthesis, hepatic and microbial transformations, systemic distribution, excretion, and EHC of BAs at the whole-body level. For model development, BA metabolism of specific pathogen-free (SPF) mice was characterized in vivo by measuring BA levels and composition in various organs, expression of transporters along the gut, and cecal microbiota composition. We found significantly different BA levels between male and female mice that could only be explained by adjusted expression of the hepatic enzymes and transporters in the model. Of note, this finding was in agreement with experimental observations. The model for SPF mice could also describe equivalent experimental data in germ-free mice by specifically switching off microbial activity in the intestine. The here presented model can therefore facilitate and guide functional analyses of BA metabolism in mice, e.g., the effect of pathophysiological alterations on BA metabolism and translation of results from mouse studies to a clinically relevant context through cross-species extrapolation.

Summary: Early onset of type 2 diabetes and cardiovascular disease are common complications for women diagnosed with gestational diabetes. Prediabetes refers to a condition in which blood glucose levels are higher than normal, but not yet high enough to be diagnosed as type 2 diabetes. Currently, there is no accurate way of knowing which women with gestational diabetes are likely to develop postpartum prediabetes. This study aims to predict the risk of postpartum prediabetes in women diagnosed with gestational diabetes. Our sparse logistic regression approach selects only two variables – antenatal fasting glucose at OGTT and HbA1c soon after the diagnosis of GDM – as relevant, but gives an area under the receiver operating characteristic curve of 0.72, outperforming all other methods. We envision this to be a practical solution, which coupled with a targeted follow-up of high-risk women, could yield better cardiometabolic outcomes in women with a history of GDM.

Summary: While lipid traits are known essential mediators of cardiovascular disease, few approaches have taken advantage of their shared genetic effects. We apply a Bayesian multivariate size estimator, mash, to GWAS of four lipid traits in the Million Veterans Program (MVP) and provide posterior mean and local false sign rates for all effects. These estimates borrow information across traits to improve effect size accuracy. We show that controlling local false sign rates accurately and powerfully identifies replicable genetic associations and that multivariate control furthers the ability to explain complex diseases. Our application yields high concordance between independent datasets, more accurately prioritizes causal genes, and significantly improves polygenic prediction beyond state-of-the-art methods by up to 59% for lipid traits. The use of Bayesian multivariate genetic shrinkage has yet to be applied to human quantitative trait GWAS results, and we present a staged approach to prediction on a polygenic scale.

Summary: Presence/absence variation (PAV) is a well-known phenomenon in prokaryotes that was described for the first time in bivalves in 2020 in Mytilus galloprovincialis. The objective of the present study was to further our understanding of the PAV phenomenon in mussel biology. The distribution of PAV was studied in a mussel chromosome-level genome assembly, revealing a widespread distribution but with hotspots of dispensability. Special attention was given to the effect of PAV in gene expression, since dispensable genes were found to be inherently subject to distortions due to their sparse distribution among individuals. Furthermore, the high expression and strong tissue specificity of some dispensable genes, such as myticins, strongly supported their biological relevance. The significant differences in the repertoire of dispensable genes associated with two geographically distinct populations suggest that PAV is involved in local adaptation. Overall, the PAV phenomenon would provide a key selective advantage at the population level.

Summary: The increased amount of tertiary lymphoid structures (TLSs) is associated with a favorable prognosis in patients with lung adenocarcinoma (LUAD). However, evaluating TLSs manually is an experience-dependent and time-consuming process, which limits its clinical application. In this multi-center study, we developed an automated computational workflow for quantifying the TLS density in the tumor region of routine hematoxylin and eosin (H&E)-stained whole-slide images (WSIs). The association between the computerized TLS density and disease-free survival (DFS) was further explored in 802 patients with resectable LUAD of three cohorts. Additionally, a Cox proportional hazard regression model, incorporating clinicopathological variables and the TLS density, was established to assess its prognostic ability. The computerized TLS density was an independent prognostic biomarker in patients with resectable LUAD. The integration of the TLS density with clinicopathological variables could support individualized clinical decision-making by improving prognostic stratification.

Summary: Adenosquamous carcinoma (ASC) is frequently misdiagnosed or overlooked in clinical practice due to its dual histological components and potential transformation from either adenocarcinoma (ADC) or squamous cell carcinoma (SCC). Our study aimed to differentiate ASC from ADC and SCC by incorporating features of enhanced CTs and clinical characteristics to build radiomics and deep learning models. The classification models were trained in Xiangya Hospital and validated in two other independent hospitals. The areas under the receiver operating characteristic curves (AUC), accuracy, sensitivity, specificity, positive predictive value, and negative predictive value were used to estimate the performance. The optimal three-class classification model achieved a maximum AUC of 0.89 and accuracy of 0.81 in external validation sets, AUC of 0.99 and accuracy of 0.99 in the internal test set. These findings highlight the efficacy of our models in differentiating ASC, providing a non-invasive, timely, and accurate diagnostic approach before and during the treatment.

Summary: Personalized prediction is ideal in chronic lymphocytic leukemia (CLL). Although refined models have been developed, stratifying patients in risk groups, it is required to accommodate time-dependent information of patients, to address the clinical heterogeneity observed within these groups. In this direction, this study proposes a personalized stepwise dynamic predictive algorithm (PSDPA) for the time-to-first-treatment of the individual patient. The PSDPA introduces a personalized Score, reflecting the evolution in the patient’s follow-up, employed to develop a reference pool of patients. Score evolution’s similarity is used to predict, at a selected time point, the time-to-first-treatment for a new patient. Additional patient’s biological information may be utilized. The algorithm was applied to 20 CLL patients, indicating that stricter assessment criteria for the Score evolution’s similarity, and biological similarity exploitation, may improve prediction. The PSDPA capitalizes on both the follow-up and the biological background of the individual patient, dynamically promoting personalized prediction in CLL.

Summary: Acute myocardial infarction dominates coronary artery disease mortality. Identifying bio-signatures for plaque destabilization and rupture is important for preventing the transition from coronary stability to instability and the occurrence of thrombosis events. This computational systems biology study enrolled 2,235 samples from 22 independent bulks cohorts and 14 samples from two single-cell cohorts. A machine-learning integrative program containing nine learners was developed to generate a warning classifier linked to atherosclerotic plaque vulnerability signature (APVS). The classifier displays the reliable performance and robustness for distinguishing ST-elevation myocardial infarction from chronic coronary syndrome at presentation, and revealed higher accuracy to 33 pathogenic biomarkers. We also developed an APVS-based quantification system (APVSLevel) for comprehensively quantifying atherosclerotic plaque vulnerability, empowering early-warning capabilities, and accurate assessment of atherosclerosis severity. It unraveled the multidimensional dysregulated mechanisms at high resolution. This study provides a potential tool for macro-level differential diagnosis and evaluation of subtle genetic pathological changes in atherosclerosis.

Summary: Trilobites were one of the first animals on Earth to leave their imprints on the seafloor. Such imprints represent behavioral traces related to feeding or protection, in both cases implying different types of locomotion. Modeling how trilobites moved is essential to understand their evolutionary history and ecological impact on marine substrates. Herein, locomotion in trilobites is approached by means of three-dimensional models, which yielded two main gait types. These two gaits reflect basic behaviors: burrowing and walking. This model reveals that trilobites could change their gait and consequently increase rapidly their speed varying the amplitude of the metachronal wave, a change independent from their biological structure. Fast increases in speed enhanced the protection of trilobites against predators and sudden environmental crises. The trilobite body pattern constrained their gaits, controlled by the distance between the pair of legs and between legs in a same segment.

Summary: Here, we exploit a deep serological profiling strategy coupled with an integrated, computational framework for the analysis of SARS-CoV-2 humoral immune responses. Applying a high-density peptide array (HDPA) spanning the entire proteomes of SARS-CoV-2 and endemic human coronaviruses allowed identification of B cell epitopes and relate them to their evolutionary and structural properties. We identify hotspots of pre-existing immunity and identify cross-reactive epitopes that contribute to increasing the overall humoral immune response to SARS-CoV-2. Using a public dataset of over 38,000 viral genomes from the early phase of the pandemic, capturing both inter- and within-host genetic viral diversity, we determined the evolutionary profile of epitopes and the differences across proteins, waves, and SARS-CoV-2 variants. Lastly, we show that mutations in spike and nucleocapsid epitopes are under stronger selection between than within patients, suggesting that most of the selective pressure for immune evasion occurs upon transmission between hosts.

Summary: This paper describes the development of a deep learning model for prediction of hip fractures on pelvic radiographs (X-rays). Developed using over 40,000 pelvic radiographs from a single institution, the model demonstrated high sensitivity and specificity when applied to a test set of emergency department radiographs. This study approximates the real-world application of a deep learning fracture detection model by including radiographs with sub-optimal image quality, other non-hip fractures, and metallic implants, which were excluded from prior published work. The study also explores the effect of ethnicity on model performance, as well as the accuracy of visualization algorithm for fracture localization.

Summary: Metabarcoding revolutionized our understanding of diversity and ecology of microorganisms in different habitats. However, it is also associated with several inherent biases, one of which is associated with intragenomic diversity of a molecular barcode. Here, we compare intragenomic variability of the V9 region of the 18S rRNA gene in 19 eukaryotic phyla abundant in marine plankton. The level of intragenomic variability is comparable across all the phyla, and in most genomes and transcriptomes one V9 sequence and one OTU is predominant. However, most of the variability observed at the barcode level is probably caused by sequencing errors and is mitigated by using a denoising tool, DADA2. The SWARM algorithm commonly used in metabarcoding studies is not optimal for collapsing genuine and erroneous sequences into a single OTU, leading to an overestimation of diversity in metabarcoding data. For an unknown reason, SWARM inflates diversity of eupelagonemids more than that of other eukaryotes.

Summary: To understand the clinical significance of the tumor microenvironment (TME), it is essential to study the interactions between malignant and non-malignant cells in clinical specimens. Here, we established a computational framework for a multiplex imaging system to comprehensively characterize spatial contexts of the TME at multiple scales, including close and long-distance spatial interactions between cell type pairs. We applied this framework to a total of 1,393 multiplex imaging data newly generated from 88 primary central nervous system lymphomas with complete follow-up data and identified significant prognostic subgroups mainly shaped by the spatial context. A supervised analysis confirmed a significant contribution of spatial context in predicting patient survival. In particular, we found an opposite prognostic value of macrophage infiltration depending on its proximity to specific cell types. Altogether, we provide a comprehensive framework to analyze spatial cellular interaction that can be broadly applied to other technologies and tumor contexts.

Summary: Droplet-based single-cell RNA-seq (scRNA-seq) data are plagued by ambient contaminations caused by nucleic acid material released by dead and dying cells. This material is mixed into the buffer and is co-encapsulated with cells, leading to a lower signal-to-noise ratio. Although there exist computational methods to remove ambient contaminations post-hoc, the reliability of algorithms in generating high-quality data from low-quality sources remains uncertain. Here, we assess data quality before data filtering by a set of quantitative, contamination-based metrics that assess data quality more effectively than standard metrics. Through a series of controlled experiments, we report improvements that can minimize ambient contamination outside of tissue dissociation, via cell fixation, improved cell loading, microfluidic dilution, and nuclei versus cell preparation; many of these parameters are inaccessible on commercial platforms. We provide end-users with insights on factors that can guide their decision-making regarding optimizations that minimize ambient contamination, and metrics to assess data quality.

Summary: Exploring early detection methods through comprehensive evaluation of DNA methylation for lung squamous cell carcinoma (LUSC) patients is of great significance. By using different machine learning algorithms for feature selection and model construction based on The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, five methylation biomarkers in LUSC (along with mapped genes) were identified including cg14823851 (TBX4), cg02772121 (TRIM15), cg10424681 (C6orf201), cg12910906 (ARHGEF4), and cg20181079 (OR4D11), achieving extremely high sensitivity and specificity in distinguishing LUSC from normal samples in independent cohorts. Pyrosequencing assay verified DNA methylation levels, meanwhile qRT-PCR and immunohistochemistry results presented their accordant methylation-related gene expression statuses in paired LUSC and normal lung tissues. The five methylation-based biomarkers proposed in this study have great potential for the diagnosis of LUSC and could guide studies in methylation-regulated tumor development and progression.

Odor perception in humans is initiated by activation of odorant receptors (ORs) in the nose. However, the ORs linked to specific olfactory percepts are unknown, unlike in vision or taste where receptors are linked to perception of different colors and tastes. The large family of ORs (~400) and multiple receptors activated by an odorant present serious challenges. Here, we first use machine learning to screen ~0.5 million compounds for new ligands and identify enriched structural motifs for ligands of 34 human ORs. We next demonstrate that the activity of ORs successfully predicts many of the 146 different perceptual qualities of chemicals. Although chemical features have been used to model odor percepts, we show that biologically relevant OR activity is often superior. Interestingly, each odor percept could be predicted with very few ORs, implying they contribute more to each olfactory percept. A similar model is observed in Drosophila where comprehensive OR-neuron data are available.

Summary: Lineage plasticity, especially transdifferentiation between neural/neuroendocrine (NE) and non-NE lineage, has been observed in multiple cancer types and linked to increased tumor aggressiveness. However, existing NE/non-NE subtype classifications in various cancer types were established through ad hoc approaches in different studies, making it difficult to align findings across cancer types and extend investigations to new datasets. To address this issue, we developed a generalized strategy to generate quantitative NE scores and a web application to facilitate its implementation. We applied this method to nine datasets covering seven cancer types, including two neural cancers, two neuroendocrine cancers, and three non-NE cancers. Our analysis revealed significant NE inter-tumoral heterogeneity and identified strong associations between NE scores and molecular, histological, and clinical features, including prognosis in different cancer types. These results support the translational utility of NE scores. Overall, our work demonstrated a broadly applicable strategy for determining the NE properties of tumors.

Summary: Monocyte-derived macrophages help maintain tissue homeostasis and defend the organism against pathogens. In tumors, recent studies have uncovered complex macrophage populations, including tumor-associated macrophages, which support tumorigenesis through cancer hallmarks such as immunosuppression, angiogenesis, or matrix remodeling. In the case of chronic lymphocytic leukemia, these macrophages are known as nurse-like cells (NLCs) and they protect leukemic cells from spontaneous apoptosis, contributing to their chemoresistance. We propose an agent-based model of monocyte differentiation into NLCs upon contact with leukemic B cells in vitro. We performed patient-specific model optimization using cultures of peripheral blood mononuclear cells from patients. Using our model, we were able to reproduce the temporal survival dynamics of cancer cells in a patient-specific manner and to identify patient groups related to distinct macrophage phenotypes. Our results show a potentially important role of phagocytosis in the polarization process of NLCs and in promoting cancer cells’ enhanced survival.

Summary: The last decade has witnessed massive advancements in high-throughput techniques capable of producing increasingly complex gene expression datasets across time and space and at the resolution of single cells. Yet, the large volume of big data available and the complexity of experimental designs hamper an easy understanding and effective communication of the results.We present expressyouRcell, an easy-to-use R package to map the multi-dimensional variations of transcript and protein levels in dynamic cell pictographs. expressyouRcell visualizes gene expression variations as pictographic representations of cell-type thematic maps. expressyouRcell visually reduces the complexity of displaying gene expression and protein level changes across multiple measurements (time points or single-cell trajectories) by generating dynamic representations of cellular pictographs.We applied expressyouRcell to single cell, bulk RNA sequencing (RNA-seq), and proteomics datasets, demonstrating its flexibility and usability in the visualization of complex variations in gene expression. Our approach improves the standard quantitative interpretation and communication of relevant results.

Summary: Estrogen receptor positive (ER+) breast cancer is responsive to a number of targeted therapies used clinically. Unfortunately, the continuous application of targeted therapy often results in resistance, driving the consideration of combination and alternating therapies. Toward this end, we developed a mathematical model that can simulate various mono, combination, and alternating therapies for ER + breast cancer cells at different doses over long time scales. The model is used to look for optimal drug combinations and predicts a significant synergism between Cdk4/6 inhibitors in combination with the anti-estrogen fulvestrant, which may help explain the clinical success of adding Cdk4/6 inhibitors to anti-estrogen therapy. Furthermore, the model is used to optimize an alternating treatment protocol so it works as well as monotherapy while using less total drug dose.

Summary: The domains of contemporary medicine and biology have generated substantial high-dimensional genetic data. Identifying representative genes and decreasing the dimensionality of the data can be challenging. The goal of gene selection is to minimize computing costs and enhance classification precision. Therefore, this article designs a new wrapper gene selection algorithm named artificial bee bare-bone hunger games search (ABHGS), which is the hunger games search (HGS) integrated with an artificial bee strategy and a Gaussian bare-bone structure to address this issue. To evaluate and validate the performance of our proposed method, ABHGS is compared to HGS and a single strategy embedded in HGS, six classic algorithms, and ten advanced algorithms on the CEC 2017 functions. The experimental results demonstrate that the bABHGS outperforms the original HGS. Compared to peers, it increases classification accuracy and decreases the number of selected features, indicating its actual engineering utility in spatial search and feature selection.

Summary: With the rapid expansion of transcriptome studies in many fishes, a great number of RNA-seq data have been published, allowing for a more systematic understanding of the general profiles and details of gene expression in fish. FishGET is dedicated to gathering and curating fish RNA-seq data to discover more new RNAs, including mRNA and lncRNA, thereby getting a more complete reference transcriptome and providing more comprehensive and accurate transcriptome annotations. We obtained a total of 1362 RNA-seq paired-end data of 8 fishes from 97 different studies, and then we performed transcript assembly, meta-assembly, weighted gene co-expression network analysis (WGCNA), functional annotations, neighbor location annotation, lncRNA type annotation, homology annotation. To promote research into fish genes at the transcriptional level, we developed a user-friendly web interface that allows users to view all information and makes use of multiple types of dynamic interactive visualization services.

Summary: In the absence of targetable hormonal axes, chemoresistance for triple-negative breast cancer (TNBC) often compromises patient outcomes. To investigate the underlying tumor dynamics, we performed trajectory analysis on the single-nuclei RNA-seq (snRNA-seq) of chemoresistant tumor clones during neoadjuvant chemotherapy (NAC). It revealed a common tumor trajectory across multiple patients with HER2-like expansions during NAC. Genome-wide CRISPR-Cas9 knock-out on mammary epithelial cells revealed chemosensitivity-promoting knock-outs were up-regulated along the tumor trajectory. Furthermore, we derived a consensus gene signature of TNBC chemoresistance by comparing the trajectory transcriptome with chemoresistant transcriptomes from TNBC cell lines and poor prognosis patient samples to predict FDA-approved drugs, including afatinib (pan-HER inhibitor), targeting the consensus signature. We validated the synergistic efficacy of afatinib and paclitaxel in chemoresistant TNBC cells and confirmed pharmacological suppression of the consensus signature. The study provides a dynamic model of chemoresistant tumor transcriptome, and computational framework for pharmacological intervention.

Summary: The major histocompatibility complex class I (MHC-I) genes are highly polymorphic. MHC-I genotyping is required for determining the peptide epitopes available to an individual’s T-cell repertoire. Current genotyping software tools do not work for the dog, due to very limited known canine alleles. To address this, we developed a Kmer-based paired-end read (KPR) de novo assembler and genotyper, which assemble paired-end RNA-seq reads from MHC-I regions into contigs, and then genotype each contig and estimate its expression level. KPR tools outperform other popular software examined in typing new alleles. We used KPR tools to successfully genotype152 dogs from a published dataset. The study discovers 33 putative new alleles, finds dominant alleles in 4 dog breeds, and builds allele diversity and expression landscapes among the 152 dogs. Our software meets a significant need in biomedical research.

Summary: The metabolic activity of all the micro-organism composing the human microbiome interacts with the host metabolism contributing to human health and disease in a way that is not fully understood. Here, we introduce STELLA, a computational method to derive the spectrum of metabolites associated with the microbiome of an individual. STELLA integrates known information on metabolic pathways associated with each bacterial species and extracts from these the list of metabolic products of each singular reaction by means of automatic text analysis. By comparing the result obtained on a single subject with the metabolic profile data of a control set of healthy subjects, we are able to identify individual metabolic alterations. To illustrate the method, we present applications to autism spectrum disorder and multiple sclerosis.

Summary: The chromatin loop plays a critical role in the study of gene expression and disease. Supervised learning-based algorithms to predict the chromatin loops require large priori information to satisfy the model construction, while the prediction sensitivity of unsupervised learning-based algorithms is still unsatisfactory. Therefore, we propose an unsupervised algorithm, Ecomap-loop. It takes advantage of extrusion complex-associated patterns, including CTCF, RAD21, and SMC enrichments, as well as the orientation distribution of CTCF motif of loops to build feature matrices; then the eigen decomposition model is employed to obtain the cell type-specific loops. We compare the performance of Ecomap-loop with the state-of-the-art unsupervised algorithm using Hi-C, ChIA-PET, expression quantitative trait locus (eQTL), and CRISPR interference (CRISPRi) screen data; the results show that Ecomap-loop achieves the best in four cell types. In addition, the functional analysis reveals the ability of Ecomap-loop to predict active functionality-related and cell type-specific loops.

Summary: Graph and image are two common representations of Hi-C cis-contact maps. Existing computational tools have only adopted Hi-C data modeled as unitary data structures but neglected the potential advantages of synergizing the information of different views. Here we propose GILoop, a dual-branch neural network that learns from both representations to identify genome-wide CTCF-mediated loops. With GILoop, we explore the combined strength of integrating the two view representations of Hi-C data and corroborate the complementary relationship between the views. In particular, the model outperforms the state-of-the-art loop calling framework and is also more robust against low-quality Hi-C libraries. We also uncover distinct preferences for matrix density by graph-based and image-based models, revealing interesting insights into Hi-C data elucidation. Finally, along with multiple transfer-learning case studies, we demonstrate that GILoop can accurately model the organizational and functional patterns of CTCF-mediated looping across different cell lines.

Summary: Drug resistance and tumor relapse in patients with melanoma is attributed to a combination of genetic and non-genetic mechanisms. Dedifferentiation, a common mechanism of non-genetic resistance in melanoma is characterized by the loss of melanocytic markers. While various molecular attributes of de-differentiation have been identified, the transition dynamics remain poorly understood. Here, we construct cell-state transition landscapes, to quantify the stochastic dynamics driving phenotypic switching in melanoma based on its underlying regulatory network. These landscapes reveal the existence of multiple alternative paths to resistance—de-differentiation and transition to a hyper-pigmented phenotype. Finally, by visualizing the changes in the landscape during in silico molecular perturbations, we identify combinatorial strategies that can lead to the most optimal outcome—a landscape with the minimum occupancy of the two drug-resistant states. Therefore, we present these landscapes as platforms to screen possible therapeutic interventions in terms of their ability to lead to the most favorable patient outcomes.

Summary: Predicting associations between microRNAs (miRNAs) and diseases from the viewpoint of function modules has become increasingly popular. However, existing methods obtained the relations between diseases and miRNAs only through the construction of similarity networks and neglected the complex network characteristic. In this paper, a new method named combining miRNA function similarities and network topology similarities based on module identification in networks (ComSim-MINE) was developed. Combined similarity is calculated from the harmonic mean between miRNA function similarities and network topology similarities. Experimental results showed that ComSim-MINE can compete with several state-of-the-art weighted function module algorithms, such as ClusterONE, MCODE, NEMO, and SPICi, and achieved the satisfactory results in terms of the composite score of F-measure, sensitivity, and accuracy based on the generated miRNA function interaction network. From the analysis of case studies, some new findings obtained from our proposed method provide clinicians new clues for epidemic diseases, such as COVID-19.

Summary: CRISPR effector Cas13 recognizes and degrades RNA molecules that are complementary to its guide RNA (gRNA) and possesses potential as an antiviral biotechnology because it can degrade viral mRNA and RNA genomes. Because multiplexed targeting is a critical strategy to improve viral suppression, we developed a strategy to design of gRNAs where individual gRNAs have maximized activity at multiple viral targets, simultaneously, by exploiting the molecular biophysics of promiscuous target recognition by Cas13. These “polyvalent” gRNA sequences (“pgRNAs”) provide superior antiviral elimination across tissue/organ scales in a higher organism (Nicotiana benthamiana) compared to conventionally-designed gRNAs—reducing detectable viral RNA by >30-fold, despite lacking perfect complementarity with either of their targets and, when multiplexed, reducing viral RNA by >99.5%. Pairs of pgRNA-targetable sequences are abundant in the genomes of RNA viruses, and this work highlights the need for specific approaches to the challenges of targeting viruses in eukaryotes using CRISPR.