系统生物学

Non-small cell lung cancer (NSCLC) is the most common type of lung cancer. According to 2020 reports, globally, 2.2 million cases are reported every year, with the mortality number being as high as 1.8 million patients. To study NSCLC, systems biology offers mathematical modeling as a tool to understand complex pathways and provide insights into the identification of biomarkers and potential therapeutic targets, which aids precision therapy. Mathematical modeling, specifically ordinary differential equations (ODEs), is used to better understand the dynamics of cancer growth and immunological interactions in the tumor microenvironment. This study highlighted the dual role of the cyclic GMP-AMP synthase–stimulator of interferon genes (cGAS/STING) pathway's classical involvement in regulating type 1 interferon (IFN I) and pro-inflammatory responses to promote tumor regression through senescence and apoptosis. Alternative signaling was induced by nuclear factor kappa B (NF-κB), mutated tumor protein p53 (p53), and programmed death-ligand1 (PD-L1), which lead to tumor growth. We identified key regulators in cancer progression by simulating the model and validating it with the following model estimation parameters: local sensitivity analysis, principal component analysis, rate of flow of metabolites, and model reduction. Integration of multiple signaling axes revealed that cGAS-STING, phosphoinositide 3-kinases (PI3K), and Ak strain transforming (AKT) may be potential targets that can be validated for cancer therapy.

Mitochondrial genomes (mitogenomes) display relatively rapid mutation rates, low sequence recombination, high copy numbers, and maternal inheritance patterns, rendering them valuable blueprints for mapping lineages, uncovering historical migration patterns, understanding intraspecific population dynamics, and investigating how environmental pressures shape traits underpinned by genetic variation. Here, we present the bioinformatic pipeline and code used to assemble and annotate the complete mitogenomes of five houndsharks (Chondrichthyes: Triakidae) and compare them to the mitogenomes of other closely related species. We demonstrate the value of a combined assembly approach for detecting deviations in mitogenome structure and describe how to select an assembly approach that best suits the sequencing data. The datasets required to run our analyses are available on the GitHub and Dryad repositories.

The limited standards for the rigorous and objective use of mitochondrial genomes (mitogenomes) can lead to uncertainties regarding the phylogenetic relationships of taxa under varying evolutionary constraints. The mitogenome exhibits heterogeneity in base composition, and evolutionary rates may vary across different regions, which can cause empirical data to violate assumptions of the applied evolutionary models. Consequently, the unique evolutionary signatures of the dataset must be carefully evaluated before selecting an appropriate approach for phylogenomic inference. Here, we present the bioinformatic pipeline and code used to expand the mitogenome phylogeny of the order Carcharhiniformes (groundsharks), with a focus on houndsharks (Chondrichthyes: Triakidae). We present a rigorous approach for addressing difficult-to-resolve phylogenies, incorporating multi-species coalescent modelling (MSCM) to address gene/species tree discordance. The protocol describes carefully designed approaches for preparing alignments, partitioning datasets, assigning models of evolution, inferring phylogenies based on traditional site-homogenous concatenation approaches as well as under multispecies coalescent and site heterogenous models, and generating statistical data for comparison of different topological outcomes. The datasets required to run our analyses are available on GitHub and Dryad repositories.

Capturing produced, consumed, or exchanged metabolites (metabolomics) and the result of gene expression (transcriptomics) require the extraction of metabolites and RNA. Multi-omics approaches and, notably, the combination of metabolomics and transcriptomic analyses are required for understanding the functional changes and adaptation of microorganisms to different physico-chemical and environmental conditions. A protocol was developed to extract total RNA and metabolites from less than 6 mg of a kind of phototrophic biofilm: oxygenic photogranules. These granules are aggregates of several hundred micrometers up to several millimeters. They harbor heterotrophic bacteria and phototrophs. After a common step for cell disruption by bead-beating, a part of the volume was recovered for RNA extraction, and the other half was used for the methanol- and dichloromethane-based extraction of metabolites. The solvents enabled the separation of two phases (aqueous and lipid) containing hydrophilic and lipophilic metabolites, respectively. The ¹H nuclear magnetic resonance (NMR) analysis of these extracts produced spectra that contained over a hundred signals with a signal-to-noise ratio higher than 10. The quality of the spectra enabled the identification of dozens of metabolites per sample. Total RNA was purified using a commercially available kit, yielding sufficient concentration and quality for metatranscriptomic analysis. This novel method enables the co-extraction of RNA and metabolites from the same sample, as opposed to the parallel extraction from two samples. Using the same sample for both extractions is particularly advantageous when working with inherently heterogeneous complex biofilm. In heterogeneous systems, differences between samples may be substantial. The co-extraction will enable a holistic analysis of the metabolomics and metatranscriptomics data generated, minimizing experimental biases, including technical variations and, notably, biological variability. As a result, it will ensure more robust multi-omics analyses, particularly by improving the correlation between metabolic changes and transcript modifications.

Many small molecules require derivatization to increase their volatility and to be amenable to gas chromatographic (GC) separation. Derivatization is usually time-consuming, and typical batch-wise procedures increase sample variability. Sequential automation of derivatization via robotic liquid handling enables the overlapping of sample preparation and analysis, maximizing time efficiency and minimizing variability. Herein, a protocol for the fully automated, two-stage derivatization of human blood–based samples in line with GC–[Orbitrap] mass spectrometry (MS)-based metabolomics is described. The protocol delivers a sample-to-sample runtime of 31 min, being suitable for better throughput routine metabolomic analysis.

Bone repair is a complex regenerative process relying on skeletal stem/progenitor cells (SSPCs) recruited predominantly from the periosteum. Activation and differentiation of periosteal SSPCs occur in a heterogeneous environment, raising the need for single cell/nucleus transcriptomics to decipher the response of the periosteum to injury. Enzymatic cell dissociation can induce a stress response affecting the transcriptome and lead to overrepresentation of certain cell types (i.e., immune and endothelial cells) and low coverage of other cell types of interest. To counteract these limitations, we optimized a protocol to isolate nuclei directly from the intact periosteum and from the fracture callus to perform single-nucleus RNA sequencing. This protocol is adapted for fresh murine periosteum, fracture callus, and frozen human periosteum. Nuclei are isolated using mechanical extraction combined with fluorescence-based nuclei sorting to obtain high-quality nucleus suspensions. This protocol allows the capture of the full diversity of cell types in the periosteum and fracture environment to better reflect the in vivo tissue composition.

Dual RNA-Seq technology has significantly advanced the study of biological interactions between two organisms by allowing parallel transcriptomic analysis. Existing analysis methods employ various combinations of open-source bioinformatics tools to process dual RNA-Seq data. Upon reviewing these methods, we intend to explore crucial criteria for selecting standard tools and methods, especially focusing on critical steps such as trimming and mapping reads to the reference genome. In order to validate the different combinatorial approaches, we performed benchmarking using top-ranking tools and a publicly available dual RNA-Seq Sequence Read Archive (SRA) dataset. An important observation while evaluating the mapping approach is that when the adapter trimmed reads are first mapped to the pathogen genome, more reads align to the pathogen genome than the unmapped reads derived from the traditional host-first mapping approach. This mapping method prevents the misalignment of pathogen reads to the host genome due to their shorter length. In this way, the pathogenic read information found at lesser proportions in a complex eukaryotic dataset is precisely obtained. This protocol presents a comprehensive comparison of these possible approaches, resulting in a robust unified standard methodology.

Glioblastoma (GBM) is the most aggressive brain tumor, and different efforts have been employed in the search for new drugs and therapeutic protocols for GBM. A label-free, mass spectrometry–based quantitative proteomics has been developed to identify and characterize proteins that are differentially expressed in GBM to gain a better understanding of the interactions and functions that lead to the pathological state focusing on the extracellular matrix (ECM). The main challenge in GBM research has been to identify novel molecular therapeutic targets and accurate diagnostic/prognostic biomarkers. To better investigate the GBM secretome upon in vitro treatment with histone deacetylase inhibitor (iHDAC), we employed a high-throughput label-free methodology of protein identification and quantification based on mass spectrometry followed by in silico studies. Our analysis revealed significant changes in the ECM protein profile, particularly those associated with the angiogenic matrisome. Proteins such as decorin, ADAM10, ADAM12, and ADAM15 were differentially regulated upon in silico analysis. In contrast, key angiogenesis markers such as VEGF and ECM proteins like fibronectin and integrins did not display significant changes. These results suggest that iHDAC inhibitors may modulate or suppress tumor behavior growth by targeting ECM proteins’ secretion rather than directly inhibiting angiogenesis.

Cellular communication relies on the intricate interplay of signaling molecules, which come together to form the cell–cell interaction (CCI) network that orchestrates tissue behavior. Researchers have shown that shallow neural networks can effectively reconstruct the CCI from the abundant molecular data captured in spatial transcriptomics (ST). However, in scenarios characterized by sparse connections and excessive noise within the CCI, shallow networks are often susceptible to inaccuracies, leading to suboptimal reconstruction outcomes. To achieve a more comprehensive and precise CCI reconstruction, we propose a novel method called triple-enhancement-based graph neural network (TENET). The TENET framework has been implemented and evaluated on both real and synthetic ST datasets. This protocol primarily introduces our network architecture and its implementation.

Stable-isotope resolved metabolomics (SIRM) is a powerful approach for characterizing metabolic states in cells and organisms. By incorporating isotopes, such as ¹³C, into substrates, researchers can trace reaction rates across specific metabolic pathways. Integrating metabolomics data with gene expression profiles further enriches the analysis, as we demonstrated in our prior study on glioblastoma metabolic symbiosis. However, the bioinformatics tools for analyzing tracer metabolomics data have been limited. In this protocol, we encourage the researchers to use SIRM and transcriptomics data and to perform the downstream analysis using our software tool DIMet. Indeed, DIMet is the first comprehensive tool designed for the differential analysis of tracer metabolomics data, alongside its integration with transcriptomics data. DIMet facilitates the analysis of stable-isotope labeling and metabolic abundances, offering a streamlined approach to infer metabolic changes without requiring complex flux analysis. Its pathway-based "metabologram" visualizations effectively integrate metabolomics and transcriptomics data, offering a versatile platform capable of analyzing corrected tracer datasets across diverse systems, organisms, and isotopes. We provide detailed steps for sample preparation and data analysis using DIMet through its intuitive, web-based Galaxy interface. To showcase DIMet's capabilities, we analyzed LDHA/B knockout glioblastoma cell lines compared to controls. Accessible to all researchers through Galaxy, DIMet is free, user-friendly, and open source, making it a valuable resource for advancing metabolic research.