Systems Biology Perspective Research Articles

Laparoscopic sigmoid colectomy with primary anastomosis for experimental modeling in the nonhuman primate.

Laparoscopic colon surgery is performed frequently in the clinical setting for a multitude of reasons including cancer, infection, and autoimmune disease. As a result, extensive research has been conducted in relation to clinical outcomes after surgery, but more recently, in relation to the impact of surgery and other patient factors on physiologic homeostasis including the host microbiome. Despite this, experimental surgical models for laparoscopic colon surgery are scarce in the literature with most studies utilizing rodents. While rodent studies provide valuable insights into basic mechanistic processes, the translation of novel therapeutic approaches to clinical practice often requires the use of large animal models. In exploring the intricate systems biology linking surgery and medicine, sophisticated models such as nonhuman primates (NHPs) play a pivotal role. By closely resembling human anatomical, physiological, and behavioral characteristics, NHPs facilitate the development and refinement of complex surgical techniques and peri-operative practices. Furthermore, they enable longitudinal studies that comprehensively assess both immediate and long-term outcomes. The availability and utilization of multiple robust models enhance the validity of surgical research, leading to more successful translation to human clinical practice. Here we describe our technique for performing a laparoscopic sigmoid colectomy with a primary anastomosis in an NHP. The entire procedure was well tolerated without significant ventilation or hemodynamic issue. To our knowledge, this represents the first laparoscopic sigmoid colectomy with primary anastomosis performed in an NHP. Furthermore, this demonstrates the feasibility of the technique and provides a relevant, preclinical model for the study of surgical colon disease. Although the surgical colectomy model in NHPs closely resembles the clinical scenario, it is crucial to recognize that a 'model' inherently comes with limitations. The intended use of any model should be carefully evaluated concerning the target patient population with the consideration of potential disparities in anatomy, physiology, environmental factors, and disease to properly interpret results. This model provides an opportunity to study mechanisms, from a systems biology perspective, underlying both innovative surgical treatments and their effects on diseases such as colon cancer, as well as benign conditions like inflammatory bowel disease, diverticulitis, and anastomotic leak, offering high predictive value.

Research Progress on Genomic Regions and Candidate Genes Related to Milk Composition Traits of Dairy Goats Based on Functional Genomics: A Narrative Review.

Goat milk has gained global attention for its unique nutritional properties and potential health benefits. Advancements in functional genomic technologies have significantly progressed genetic research on milk composition traits in dairy goats. This review summarizes various research methodologies applied in this field. Genome-wide association studies (GWAS) have identified genomic regions associated with major milk components, with the diacylglycerol acyltransferase 1 (DGAT1) gene and casein gene cluster consistently linked to milk composition traits. Transcriptomics has revealed gene expression patterns in mammary tissue across lactation stages, while the role of non-coding RNAs (such as miRNAs and circRNAs) in regulating milk composition has been confirmed. Proteomic and metabolomic studies have not only helped us gain a more comprehensive understanding of goat milk composition characteristics but have also provided crucial support for the functional validation of genes related to milk components. The integration of multi-omics data has emerged as an effective strategy for elucidating complex regulatory networks from a systems biology perspective. Despite progress, challenges remain, including refining reference genomes, collecting large-scale phenotypic data, and conducting functional validations. Future research should focus on improving reference genomes, expanding study populations, investigating functional milk components, exploring epigenetic regulation and non-coding RNAs, and studying microbiome-host genome interactions. These efforts will inform more precise genomic and marker-assisted selection strategies, advancing genetic improvements in milk composition traits in dairy goats.

Network modeling links kidney developmental programs and the cancer type-specificity of VHL mutations

Elucidating the molecular dependencies behind the cancer-type specificity of driver mutations may reveal new therapeutic opportunities. We hypothesized that developmental programs would impact the transduction of oncogenic signaling activated by a driver mutation and shape its cancer-type specificity. Therefore, we designed a computational analysis framework by combining single-cell gene expression profiles during fetal organ development, latent factor discovery, and information theory-based differential network analysis to systematically identify transcription factors that selectively respond to driver mutations under the influence of organ-specific developmental programs. After applying this approach to VHL mutations, which are highly specific to clear cell renal cell carcinoma (ccRCC), we revealed important regulators downstream of VHL mutations in ccRCC and used their activities to cluster patients with ccRCC into three subtypes. This classification revealed a more significant difference in prognosis than the previous mRNA profile-based method and was validated in an independent cohort. Moreover, we found that EP300, a key epigenetic factor maintaining the regulatory network of the subtype with the worst prognosis, can be targeted by a small inhibitor, suggesting a potential treatment option for a subset of patients with ccRCC. This work demonstrated an intimate relationship between organ development and oncogenesis from the perspective of systems biology, and the method can be generalized to study the influence of other biological processes on cancer driver mutations.

Integrating “fast” and “slow” systems therapeutics for physiological renormalization in neurodegenerative disorders

In our previous paper, we provided evidence that simultaneously activating multiple cannabinergic pathways, CB1 and CB2, a systems therapeutic approach, was harnessed to provide “fast” protection from neurodegenerative diseases by renormalizing the physiology of the retina. We had used an animal model of glaucoma, using glutamate-induced neurodegeneration of the retina. Other forms of chronic or acute neurotoxicity throughout the CNS may likewise benefit from a systems approach to cannabinergic mediated rescue of neural tissue. However, from a systems biology perspective, cannabinergic pathways act in many ways, not fully explained in our previous paper, and must act in synergy with many other pathways, including the molecules (e.g. HSP) released by adjacent stem cells, to optimally yield neuroprotection. In addition to the “fast” pathways previously described, we’ll address “slow” acting, lifestyle and nutritional pathways too. Building on our previous paper, we highlight some of the important “fast and “slow” neuroprotective pathways that should be activated together through a systems therapeutic approach.

Severe Allergy as a Chronic Inflammatory Condition From a Systems Biology Perspective.

Persistent and unresolved inflammation is a common underlying factor observed in several and seemingly unrelated human diseases, including cardiovascular and neurodegenerative diseases. Particularly, in atopic conditions, acute inflammatory responses such as those triggered by insect venom, food or drug allergies possess also a life-threatening potential. However, respiratory allergies predominantly exhibit late immune responses associated with chronic inflammation, that can eventually progress into a severe phenotype displaying similar features as those observed in other chronic inflammatory diseases, as is the case of uncontrolled severe asthma. This review aims to explore the different facets and systems involved in chronic allergic inflammation, including processes such as tissue remodelling and immune cell dysregulation, as well as genetic, metabolic and microbiota alterations, which are common to other inflammatory conditions. Our goal here was to deepen on the understanding of an entangled disease as is chronic allergic inflammation and expose potential avenues for the development of better diagnostic and intervention strategies.

Phenotypic Plasticity and Cancer: A System Biology Perspective.

Epithelial-to-mesenchymal transition (EMT) is a major axis of phenotypic plasticity not only in diseased conditions such as cancer metastasis and fibrosis but also during normal development and wound healing. Yet-another important axis of plasticity with metastatic implications includes the cancer stem cell (CSCs) and non-CSC transitions. However, in both processes, epithelial (E) and mesenchymal (M) phenotypes are not merely binary states. Cancer cells acquire a spectrum of phenotypes with traits, properties, and markers of both E and M phenotypes, giving rise to intermediary hybrid (E/M) phenotypes. E/M cells play an important role in tumor initiation, metastasis, and disease progression in multiple cancers. Furthermore, the hybrid phenotypes also play a major role in causing therapeutic resistance in cancer. Here, we discuss how a systems biology perspective on the problem, which is implicit in the 'Team Medicine' approach outlined in the theme of this Special Issue of The Journal of Clinical Medicine and includes an interdisciplinary team of experts, is more likely to shed new light on EMT in cancer and help us to identify novel therapeutics and strategies to target phenotypic plasticity in cancer.

LipidSig 2.0: integrating lipid characteristic insights into advanced lipidomics data analysis.

In the field of lipidomics, where the complexity of lipid structures and functions presents significant analytical challenges, LipidSig stands out as the first web-based platform providing integrated, comprehensive analysis for efficient data mining of lipidomic datasets. The upgraded LipidSig 2.0 (https://lipidsig.bioinfomics.org/) simplifies the process and empowers researchers to decipher the complex nature of lipids and link lipidomic data to specific characteristics and biological contexts. This tool markedly enhances the efficiency and depth of lipidomic research by autonomously identifying lipid species and assigning 29 comprehensive characteristics upon data entry. LipidSig 2.0 accommodates 24 data processing methods, streamlining diverse lipidomic datasets. The tool's expertise in automating intricate analytical processes, including data preprocessing, lipid ID annotation, differential expression, enrichment analysis, and network analysis, allows researchers to profoundly investigate lipid properties and their biological implications. Additional innovative features, such as the 'Network' function, offer a system biology perspective on lipid interactions, and the 'Multiple Group' analysis aids in examining complex experimental designs. With its comprehensive suite of features for analyzing and visualizing lipid properties, LipidSig 2.0 positions itself as an indispensable tool for advanced lipidomics research, paving the way for new insights into the role of lipids in cellular processes and disease development.

A Network Toxicology Approach for Mechanistic Modelling of Nanomaterial Hazard and Adverse Outcomes.

Hazard assessment is the first step in evaluating the potential adverse effects of chemicals. Traditionally, toxicological assessment has focused on the exposure, overlooking the impact of the exposed system on the observed toxicity. However, systems toxicology emphasizes how system properties significantly contribute to the observed response. Hence, systems theory states that interactions store more information than individual elements, leading to the adoption of network based models to represent complex systems in many fields of life sciences. Here, they develop a network-based approach to characterize toxicological responses in the context of a biological system, inferring biological system specific networks. They directly link molecular alterations to the adverse outcome pathway (AOP) framework, establishing direct connections between omics data and toxicologically relevant phenotypic events. They apply this framework to a dataset including 31 engineered nanomaterials with different physicochemical properties in two different in vitro and one in vivo models and demonstrate how the biological system is the driving force of the observed response. This work highlights the potential of network-based methods to significantly improve their understanding of toxicological mechanisms from a systems biology perspective and provides relevant considerations and future data-driven approaches for the hazard assessment of nanomaterials and other advanced materials.

Computational single-cell methods for predicting cancer risk.

Despite recent biotechnological breakthroughs, cancer risk prediction remains a formidable computational and experimental challenge. Addressing it is critical in order to improve prevention, early detection and survival rates. Here, I briefly summarize some key emerging theoretical and computational challenges as well as recent computational advances that promise to help realize the goals of cancer-risk prediction. The focus is on computational strategies based on single-cell data, in particular on bottom-up network modeling approaches that aim to estimate cancer stemness and dedifferentiation at single-cell resolution from a systems-biological perspective. I will describe two promising methods, a tissue and cell-lineage independent one based on the concept of diffusion network entropy, and a tissue and cell-lineage specific one that uses transcription factor regulons. Application of these tools to single-cell and single-nucleus RNA-seq data from stages prior to invasive cancer reveal that they can successfully delineate the heterogeneous inter-cellular cancer-risk landscape, identifying those cells that are more likely to turn cancerous. Bottom-up systems biological modeling of single-cell omic data is a novel computational analysis paradigm that promises to facilitate the development of preventive, early detection and cancer-risk prediction strategies.

CSF amino acid profiles in ICV-streptozotocin-induced sporadic Alzheimer's disease in male Wistar rat: a metabolomics and systems biology perspective.

Alzheimer's disease (AD) is an increasingly important public health concern due to the increasing proportion of older individuals within the general population. The impairment of processes responsible for adequate brain energy supply primarily determines the early features of the aging process. Restricting brain energy supply results in brain hypometabolism prior to clinical symptoms and is anatomically and functionally associated with cognitive impairment. The present study investigated changes in metabolic profiles induced by intracerebroventricular-streptozotocin (ICV-STZ) in an AD-like animal model. To this end, male Wistar rats received a single injection of STZ (3 mg·kg-1) by ICV (2.5 μL into each ventricle for 5 min on each side). In the second week after receiving ICV-STZ, rats were tested for cognitive performance using the Morris Water Maze test and subsequently prepared for positron emission tomography (PET) to confirm AD-like symptoms. Tandem Mass Spectrometry (MS/MS) analysis was used to detect amino acid changes in cerebrospinal fluid (CFS) samples. Our metabolomics study revealed a reduction in the concentrations of various amino acids (alanine, arginine, aspartic acid, glutamic acid, glycine, isoleucine, methionine, phenylalanine, proline, serine, threonine, tryptophane, tyrosine, and valine) in CSF of ICV-STZ-treated animals as compared to controls rats. The results of the current study indicate amino acid levels could potentially be considered targets of nutritional and/or pharmacological interventions to interfere with AD progression.

Strategic pathway analysis for dual management of epilepsy and comorbid depression: a systems biology perspective

Strategic pathway analysis for dual management of epilepsy and comorbid depression: a systems biology perspective

A review of mathematical modeling of bone remodeling from a systems biology perspective

Bone remodeling is an essential, delicately balanced physiological process of coordinated activity of bone cells that remove and deposit new bone tissue in the adult skeleton. Due to the complex nature of this process, many mathematical models of bone remodeling have been developed. Each of these models has unique features, but they have underlying patterns. In this review, the authors highlight the important aspects frequently found in mathematical models for bone remodeling and discuss how and why these aspects are included when considering the physiology of the bone basic multicellular unit, which is the term used for the collection of cells responsible for bone remodeling. The review also emphasizes the view of bone remodeling from a systems biology perspective. Understanding the systemic mechanisms involved in remodeling will help provide information on bone pathology associated with aging, endocrine disorders, cancers, and inflammatory conditions and enhance systems pharmacology. Furthermore, some features of the bone remodeling cycle and interactions with other organ systems that have not yet been modeled mathematically are discussed as promising future directions in the field.

Genetic factors underlying host resistance to Rhipicephalus microplus tick infestation in Braford cattle: a systems biology perspective.

Approximately 80% of the world's cattle are raised in regions with a high risk of tick-borne diseases, resulting in significant economic losses due to parasitism by Rhipicephalus (Boophilus) microplus. However, the lack of a systemic biology approach hampers a comprehensive understanding of tick-host interactions that mediate tick resistance phenotypes. Here, we conducted a genome-wide association study (GWAS) of 2933 Braford cattle and found 340 single-nucleotide polymorphisms (SNPs) associated with tick counts. Gene expression analyses were performed on skin samples obtained from previously tick-exposed heifers with extremely high or low estimated breeding values for R. microplus counts. Evaluations were performed both before and after artificial infestation with ticks. Differentially expressed genes were found within 1-Mb windows centered at significant SNPs from GWAS. A total of 330 genes were related to the breakdown of homeostasis that was induced by larval attachment to bovine skin. Enrichment analysis pointed to a key role of proteolysis and signal transduction via JAK/STAT, NFKB and WNT/beta catenin signaling pathways. Integrative analysis on matrixEQTL revealed two cis-eQTLs and four significant SNPs in the genes peptidyl arginine deiminase type IV (PADI4) and LOC11449251. The integration of genomic data from QTL maps and transcriptome analyses has identified a set of twelve key genes that show significant associations with tick loads. These genes could be key candidates to improve the accuracy of genomic predictions for tick resistance in Braford cattle.

Nitrogen sensing and regulatory networks: it's about time and space.

A plant's response to external and internal nitrogen signals/status relies on sensing and signaling mechanisms that operate across spatial and temporal dimensions. From a comprehensive systems biology perspective, this involves integrating nitrogen responses in different cell types and over long distances to ensure organ coordination in real time and yield practical applications. In this prospective review, we focus on novel aspects of nitrogen (N) sensing/signaling uncovered using temporal and spatial systems biology approaches, largely in the model Arabidopsis. The temporal aspects span: transcriptional responses to N-dose mediated by Michaelis-Menten kinetics, the role of the master NLP7 transcription factor as a nitrate sensor, its nitrate-dependent TF nuclear retention, its "hit-and-run" mode of target gene regulation, and temporal transcriptional cascade identified by "network walking." Spatial aspects of N-sensing/signaling have been uncovered in cell type-specific studies in roots and in root-to-shoot communication. We explore new approaches using single-cell sequencing data, trajectory inference, and pseudotime analysis as well as machine learning and artificial intelligence approaches. Finally, unveiling the mechanisms underlying the spatial dynamics of nitrogen sensing/signaling networks across species from model to crop could pave the way for translational studies to improve nitrogen-use efficiency in crops. Such outcomes could potentially reduce the detrimental effects of excessive fertilizer usage on groundwater pollution and greenhouse gas emissions.

The Dual Role of Sulforaphane-Induced Cellular Stress-A Systems Biological Study.

The endoplasmic reticulum (ER) plays a crucial role in cellular homeostasis. When ER stress is generated, an autophagic self-digestive process is activated to promote cell survival; however, cell death is induced in the case of excessive levels of ER stress. The aim of the present study was to investigate the effect of a natural compound called sulforaphane (SFN) upon ER stress. Our goal was to investigate how SFN-dependent autophagy activation affects different stages of ER stress induction. We approached our scientific analysis from a systems biological perspective using both theoretical and molecular biological techniques. We found that SFN induced the various cell-death mechanisms in a concentration- and time-dependent manner. The short SFN treatment at low concentrations promoted autophagy, whereas the longer treatment at higher concentrations activated cell death. We proved that SFN activated autophagy in a mTORC1-dependent manner and that the presence of ULK1 was required for its function. A low concentration of SFN pre- or co-treatment combined with short and long ER stress was able to promote cell survival via autophagy induction in each treatment, suggesting the potential medical importance of SFN in ER stress-related diseases.

The MSC-EV-microRNAome: A Perspective on Therapeutic Mechanisms of Action in Sepsis and ARDS.

Mesenchymal stromal cells (MSCs) and MSC-derived extracellular vesicles (EVs) have emerged as innovative therapeutic agents for the treatment of sepsis and acute respiratory distress syndrome (ARDS). Although their potential remains undisputed in pre-clinical models, this has yet to be translated to the clinic. In this review, we focused on the role of microRNAs contained in MSC-derived EVs, the EV microRNAome, and their potential contribution to therapeutic mechanisms of action. The evidence that miRNA transfer in MSC-derived EVs has a role in the overall therapeutic effects is compelling. However, several questions remain regarding how to reconcile the stochiometric issue of the low copy numbers of the miRNAs present in the EV particles, how different miRNAs delivered simultaneously interact with their targets within recipient cells, and the best miRNA or combination of miRNAs to use as therapy, potency markers, and biomarkers of efficacy in the clinic. Here, we offer a molecular genetics and systems biology perspective on the function of EV microRNAs, their contribution to mechanisms of action, and their therapeutic potential.

Exploration of the link between COVID-19 and gastric cancer from the perspective of bioinformatics and systems biology.

Coronavirus disease 2019 (COVID-19), an infectious disease caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), has caused a global pandemic. Gastric cancer (GC) poses a great threat to people's health, which is a high-risk factor for COVID-19. Previous studies have found some associations between GC and COVID-19, whereas the underlying molecular mechanisms are not well understood. We employed bioinformatics and systems biology to explore these links between GC and COVID-19. Gene expression profiles of COVID-19 (GSE196822) and GC (GSE179252) were obtained from the Gene Expression Omnibus (GEO) database. After identifying the shared differentially expressed genes (DEGs) for GC and COVID-19, functional annotation, protein-protein interaction (PPI) network, hub genes, transcriptional regulatory networks and candidate drugs were analyzed. We identified 209 shared DEGs between COVID-19 and GC. Functional analyses highlighted immune-related pathways as key players in both diseases. Ten hub genes (CDK1, KIF20A, TPX2, UBE2C, HJURP, CENPA, PLK1, MKI67, IFI6, IFIT2) were identified. The transcription factor/gene and miRNA/gene interaction networks identified 38 transcription factors (TFs) and 234 miRNAs. More importantly, we identified ten potential therapeutic agents, including ciclopirox, resveratrol, etoposide, methotrexate, trifluridine, enterolactone, troglitazone, calcitriol, dasatinib and deferoxamine, some of which have been reported to improve and treat GC and COVID-19. This research offer valuable insights into the molecular interplay between COVID-19 and GC, potentially guiding future therapeutic strategies.

Cellular control of protein levels: A systems biology perspective.

How cells regulate protein levels is a central question of biology. Over the past decades, molecular biology research has provided profound insights into the mechanisms and the molecular machinery governing each step of the gene expression process, from transcription to protein degradation. Recent advances in transcriptomics and proteomics have complemented our understanding of these fundamental cellular processes with a quantitative, systems-level perspective. Multi-omic studies revealed significant quantitative, kinetic and functional differences between the genome, transcriptome and proteome. While protein levels often correlate with mRNA levels, quantitative investigations have demonstrated a substantial impact of translation and protein degradation on protein expression control. In addition, protein-level regulation appears to play a crucial role in buffering protein abundances against undesirable mRNA expression variation. These findings have practical implications for many fields, including gene function prediction and precision medicine.

Untargeted Lipidomic Approach for Studying Different Nervous System Tissues of the Murine Model of Krabbe Disease.

Krabbe disease is a rare neurodegenerative disease with an autosomal recessive character caused by a mutation in the GALC gene. The mutation leads to an accumulation of psychosine and a subsequent degeneration of oligodendrocytes and Schwann cells. Psychosine is the main biomarker of the disease. The Twitcher mouse is the most commonly used animal model to study Krabbe disease. Although there are many references to this model in the literature, the lipidomic study of nervous system tissues in the Twitcher model has received little attention. This study focuses on the comparison of the lipid profiles of four nervous system tissues (brain, cerebellum, spinal cord, and sciatic nerve) in the Twitcher mouse compared to the wild-type mouse. Altogether, approximately 230 molecular species belonging to 19 lipid classes were annotated and quantified. A comparison at the levels of class, molecular species, and lipid building blocks showed significant differences between the two groups, particularly in the sciatic nerve. The in-depth study of the lipid phenotype made it possible to hypothesize the genes and enzymes involved in the changes. The integration of metabolic data with genetic data may be useful from a systems biology perspective to gain a better understanding of the molecular basis of the disease.

Brain-wide genome-wide colocalization study for integrating genetics, transcriptomics and brain morphometry in Alzheimer's disease

Alzheimer's disease (AD) is one of the most common neurodegenerative diseases. However, the AD mechanism has not yet been fully elucidated to date, hindering the development of effective therapies. In our work, we perform a brain imaging genomics study to link genetics, single-cell gene expression data, tissue-specific gene expression data, brain imaging-derived volumetric endophenotypes, and disease diagnosis to discover potential underlying neurobiological pathways for AD. To do so, we perform brain-wide genome-wide colocalization analyses to integrate multidimensional imaging genomic biobank data. Specifically, we use (1) the individual-level imputed genotyping data and magnetic resonance imaging (MRI) data from the UK Biobank, (2) the summary statistics of the genome-wide association study (GWAS) from multiple European ancestry cohorts, and (3) the tissue-specific cis-expression quantitative trait loci (cis-eQTL) summary statistics from the GTEx project. We apply a Bayes factor colocalization framework and mediation analysis to these multi-modal imaging genomic data. As a result, we derive the brain regional level GWAS summary statistics for 145 brain regions with 482,831 single nucleotide polymorphisms (SNPs) followed by posthoc functional annotations. Our analysis yields the discovery of a potential AD causal pathway from a systems biology perspective: the SNP chr10:124165615:G>A (rs6585827) mutation upregulates the expression of BTBD16 gene in oligodendrocytes, a specialized glial cells, in the brain cortex, leading to a reduced risk of volumetric loss in the entorhinal cortex, resulting in the protective effect on AD. We substantiate our findings with multiple evidence from existing imaging, genetic and genomic studies in AD literature. Our study connects genetics, molecular and cellular signatures, regional brain morphologic endophenotypes, and AD diagnosis, providing new insights into the mechanistic understanding of the disease. Our findings can provide valuable guidance for subsequent therapeutic target identification and drug discovery in AD.