Symmetry-breaking instabilities play an important role in understanding the mechanisms underlying the diversity of patterns observed in nature, such as in Turing's reaction-diffusion theory, which connects cellular signalling and transport with the development of growth and form. Extensive literature focuses on the linear stability analysis of homogeneous equilibria in these systems, culminating in a set of conditions for transport-driven instabilities that are commonly presumed to initiate self-organisation. We demonstrate that a selection of simple, canonical transport models with only mild multistable non-linearities can satisfy the Turing instability conditions while also robustly exhibiting only transient patterns. Hence, a Turing-like instability is insufficient for the existence of a patterned state. While it is known that linear theory can fail to predict the formation of patterns, we demonstrate that such failures can appear robustly in systems with multiple stable homogeneous equilibria. Given that biological systems such as gene regulatory networks and spatially distributed ecosystems often exhibit a high degree of multistability and nonlinearity, this raises important questions of how to analyse prospective mechanisms for self-organisation.
- Klíčová slova
- Multistability, Pattern formation, Turing instabilities,
- MeSH
- biologické modely MeSH
- difuze MeSH
- ekosystém * MeSH
- genové regulační sítě MeSH
- matematické pojmy * MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
Pre-mRNA splicing is a highly coordinated process. While its dysregulation has been linked to neurological deficits, our understanding of the underlying molecular and cellular mechanisms remains limited. We implicated pathogenic variants in U2AF2 and PRPF19, encoding spliceosome subunits in neurodevelopmental disorders (NDDs), by identifying 46 unrelated individuals with 23 de novo U2AF2 missense variants (including 7 recurrent variants in 30 individuals) and 6 individuals with de novo PRPF19 variants. Eight U2AF2 variants dysregulated splicing of a model substrate. Neuritogenesis was reduced in human neurons differentiated from human pluripotent stem cells carrying two U2AF2 hyper-recurrent variants. Neural loss of function (LoF) of the Drosophila orthologs U2af50 and Prp19 led to lethality, abnormal mushroom body (MB) patterning, and social deficits, which were differentially rescued by wild-type and mutant U2AF2 or PRPF19. Transcriptome profiling revealed splicing substrates or effectors (including Rbfox1, a third splicing factor), which rescued MB defects in U2af50-deficient flies. Upon reanalysis of negative clinical exomes followed by data sharing, we further identified 6 patients with NDD who carried RBFOX1 missense variants which, by in vitro testing, showed LoF. Our study implicates 3 splicing factors as NDD-causative genes and establishes a genetic network with hierarchy underlying human brain development and function.
- Klíčová slova
- Development, Genetic diseases, Genetics, Neurodevelopment, iPS cells,
- MeSH
- enzymy opravy DNA genetika MeSH
- genové regulační sítě MeSH
- jaderné proteiny genetika MeSH
- lidé MeSH
- missense mutace MeSH
- neurovývojové poruchy * genetika MeSH
- sestřih RNA MeSH
- sestřihové faktory genetika MeSH
- spliceozomy * genetika MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- Názvy látek
- enzymy opravy DNA MeSH
- jaderné proteiny MeSH
- PRPF19 protein, human MeSH Prohlížeč
- sestřihové faktory MeSH
Lithium (Li) is one of the most effective drugs for treating bipolar disorder (BD), however, there is presently no way to predict response to guide treatment. The aim of this study is to identify functional genes and pathways that distinguish BD Li responders (LR) from BD Li non-responders (NR). An initial Pharmacogenomics of Bipolar Disorder study (PGBD) GWAS of lithium response did not provide any significant results. As a result, we then employed network-based integrative analysis of transcriptomic and genomic data. In transcriptomic study of iPSC-derived neurons, 41 significantly differentially expressed (DE) genes were identified in LR vs NR regardless of lithium exposure. In the PGBD, post-GWAS gene prioritization using the GWA-boosting (GWAB) approach identified 1119 candidate genes. Following DE-derived network propagation, there was a highly significant overlap of genes between the top 500- and top 2000-proximal gene networks and the GWAB gene list (Phypergeometric = 1.28E-09 and 4.10E-18, respectively). Functional enrichment analyses of the top 500 proximal network genes identified focal adhesion and the extracellular matrix (ECM) as the most significant functions. Our findings suggest that the difference between LR and NR was a much greater effect than that of lithium. The direct impact of dysregulation of focal adhesion on axon guidance and neuronal circuits could underpin mechanisms of response to lithium, as well as underlying BD. It also highlights the power of integrative multi-omics analysis of transcriptomic and genomic profiling to gain molecular insights into lithium response in BD.
- MeSH
- antimanika farmakologie terapeutické užití MeSH
- bipolární porucha * farmakoterapie genetika MeSH
- celogenomová asociační studie * metody MeSH
- farmakogenetika metody MeSH
- fokální adheze * účinky léků genetika MeSH
- genomika metody MeSH
- genové regulační sítě * účinky léků genetika MeSH
- indukované pluripotentní kmenové buňky účinky léků metabolismus MeSH
- lidé MeSH
- lithium * farmakologie terapeutické užití MeSH
- multiomika MeSH
- neurony metabolismus účinky léků MeSH
- sloučeniny lithia farmakologie terapeutické užití MeSH
- stanovení celkové genové exprese metody MeSH
- transkriptom * genetika účinky léků MeSH
- Check Tag
- lidé MeSH
- mužské pohlaví MeSH
- ženské pohlaví MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
- Research Support, N.I.H., Extramural MeSH
- Názvy látek
- antimanika MeSH
- lithium * MeSH
- sloučeniny lithia MeSH
MOTIVATION: The problem of model inference is of fundamental importance to systems biology. Logical models (e.g. Boolean networks; BNs) represent a computationally attractive approach capable of handling large biological networks. The models are typically inferred from experimental data. However, even with a substantial amount of experimental data supported by some prior knowledge, existing inference methods often focus on a small sample of admissible candidate models only. RESULTS: We propose Boolean network sketches as a new formal instrument for the inference of Boolean networks. A sketch integrates (typically partial) knowledge about the network's topology and the update logic (obtained through, e.g. a biological knowledge base or a literature search), as well as further assumptions about the properties of the network's transitions (e.g. the form of its attractor landscape), and additional restrictions on the model dynamics given by the measured experimental data. Our new BNs inference algorithm starts with an 'initial' sketch, which is extended by adding restrictions representing experimental data to a 'data-informed' sketch and subsequently computes all BNs consistent with the data-informed sketch. Our algorithm is based on a symbolic representation and coloured model-checking. Our approach is unique in its ability to cover a broad spectrum of knowledge and efficiently produce a compact representation of all inferred BNs. We evaluate the method on a non-trivial collection of real-world and simulated data. AVAILABILITY AND IMPLEMENTATION: All software and data are freely available as a reproducible artefact at https://doi.org/10.5281/zenodo.7688740.
Boolean networks (BNs) are a well-accepted modelling formalism in computational systems biology. Nevertheless, modellers often cannot identify only a single BN that matches the biological reality. The typical reasons for this is insufficient knowledge or a lack of experimental data. Formally, this uncertainty can be expressed using partially specified Boolean networks (PSBNs), which encode the wide range of network candidates into a single structure. In this paper, we target the control of PSBNs. The goal of BN control is to find perturbations which guarantee stabilisation of the system in the desired state. Specifically, we consider variable perturbations (gene knock-out and over-expression) with three types of application time-window: one-step, temporary, and permanent. While the control of fully specified BNs is a thoroughly explored topic, control of PSBNs introduces additional challenges that we address in this paper. In particular, the unspecified components of the model cause a significant amount of additional state space explosion. To address this issue, we propose a fully symbolic methodology that can represent the numerous system variants in a compact form. In fully specified models, the efficiency of a perturbation is characterised by the count of perturbed variables (the perturbation size). However, in the case of a PSBN, a perturbation might work only for a subset of concrete BN models. To that end, we introduce and quantify perturbation robustness. This metric characterises how efficient the given perturbation is with respect to the model uncertainty. Finally, we evaluate the novel control methods using non-trivial real-world PSBN models. We inspect the method's scalability and efficiency with respect to the size of the state space and the number of unspecified components. We also compare the robustness metrics for all three perturbation types. Our experiments support the hypothesis that one-step perturbations are significantly less robust than temporary and permanent ones.
- Klíčová slova
- Boolean network, Permanent control, Perturbation, Symbolic algorithm, Temporary control,
- MeSH
- algoritmy MeSH
- genové regulační sítě * MeSH
- systémová biologie * MeSH
- Publikační typ
- časopisecké články MeSH
OBJECTIVES: Clinical phenotyping and predicting treatment responses in SLE patients is challenging. Extensive blood transcriptional profiling has identified various gene modules that are promising for stratification of SLE patients. We aimed to translate existing transcriptomic data into simpler gene signatures suitable for daily clinical practice. METHODS: Real-time PCR of multiple genes from the IFN M1.2, IFN M5.12, neutrophil (NPh) and plasma cell (PLC) modules, followed by a principle component analysis, was used to identify indicator genes per gene signature. Gene signatures were measured in longitudinal samples from two childhood-onset SLE cohorts (n = 101 and n = 34, respectively), and associations with clinical features were assessed. Disease activity was measured using Safety of Estrogen in Lupus National Assessment (SELENA)-SLEDAI. Cluster analysis subdivided patients into three mutually exclusive fingerprint-groups termed (1) all-signatures-low, (2) only IFN high (M1.2 and/or M5.12) and (3) high NPh and/or PLC. RESULTS: All gene signatures were significantly associated with disease activity in cross-sectionally collected samples. The PLC-signature showed the highest association with disease activity. Interestingly, in longitudinally collected samples, the PLC-signature was associated with disease activity and showed a decrease over time. When patients were divided into fingerprints, the highest disease activity was observed in the high NPh and/or PLC group. The lowest disease activity was observed in the all-signatures-low group. The same distribution was reproduced in samples from an independent SLE cohort. CONCLUSIONS: The identified gene signatures were associated with disease activity and were indicated to be suitable tools for stratifying SLE patients into groups with similar activated immune pathways that may guide future treatment choices.
- Klíčová slova
- biomarkers, childhood-onset SLE, clustering analysis, disease activity, gene signatures, interferon, neutrophils, plasma cells,
- MeSH
- dítě MeSH
- genové regulační sítě MeSH
- lidé MeSH
- longitudinální studie MeSH
- shluková analýza MeSH
- systémový lupus erythematodes * MeSH
- transkriptom * MeSH
- Check Tag
- dítě MeSH
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- multicentrická studie MeSH
- práce podpořená grantem MeSH
Recent research has already shown that circular RNAs (circRNAs) are functional in gene expression regulation and potentially related to diseases. Due to their stability, circRNAs can also be used as biomarkers for diagnosis. However, the function of most circRNAs remains unknown, and it is expensive and time-consuming to discover it through biological experiments. In this paper, we predict circRNA annotations from the knowledge of their interaction with miRNAs and subsequent miRNA-mRNA interactions. First, we construct an interaction network for a target circRNA and secondly spread the information from the network nodes with the known function to the root circRNA node. This idea itself is not new; our main contribution lies in proposing an efficient and exact deterministic procedure based on the principle of probability-generating functions to calculate the p-value of association test between a circRNA and an annotation term. We show that our publicly available algorithm is both more effective and efficient than the commonly used Monte-Carlo sampling approach that may suffer from difficult quantification of sampling convergence and subsequent sampling inefficiency. We experimentally demonstrate that the new approach is two orders of magnitude faster than the Monte-Carlo sampling, which makes summary annotation of large circRNA files feasible; this includes their reannotation after periodical interaction network updates, for example. We provide a summary annotation of a current circRNA database as one of our outputs. The proposed algorithm could be generalized towards other types of RNA in way that is straightforward.
- Klíčová slova
- Annotation term, Circular RNA, Interaction network,
- MeSH
- biologické markery MeSH
- genové regulační sítě MeSH
- kruhová RNA * MeSH
- messenger RNA genetika metabolismus MeSH
- mikro RNA * genetika metabolismus MeSH
- pravděpodobnost MeSH
- stanovení celkové genové exprese metody MeSH
- Publikační typ
- časopisecké články MeSH
- Názvy látek
- biologické markery MeSH
- kruhová RNA * MeSH
- messenger RNA MeSH
- mikro RNA * MeSH
BACKGROUND: Oral squamous cell carcinoma (OSCC) severely affects the quality of life and the 5-year survival rate is low. Exploring the potential miRNA-mRNA regulatory network and analyzing hub genes and clinical data can provide a theoretical basis for further elucidating the pathogenesis of OSCC. METHODS: The miRNA expression datasets of GSE113956 and GSE124566 and mRNA expression datasets of GSE31056, GSE37991 and GSE13601 were obtained from the Gene Expression Omnibus databases. The differentially expressed miRNAs (DEMs) and mRNAs (DEGs) were screened using GEO2R. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed by DAVID database. The PPI network was established through STRING database and the hub genes were preliminarily screened out by Cytoscape software. After identifying the hub genes in the TCGA database, we predicted the potential DEM transcription factors, constructed a miRNA-mRNA regulatory network, and analyzed the relationship between the hub genes and clinical data. RESULTS: A total of 28 DEMs and 764 DEGs were screened out, which were composed of 285 up-regulated genes and 479 down-regulated genes. Enrichment analysis showed that up-regulation of DEGs were mainly enriched in extracellular matrix organization and cancer-related pathway, while down-regulation of DEGs were mainly enriched in muscular system process and adrenaline signal transduction. After preliminary screening by PPI network and identification in TCGA, the up-regulated FN1, COL1A1, COL1A2, AURKA, CCNB1, CCNA2, SPP1, CDC6, and down-regulated ACTN2, TTN, IGF1, CAV3, MYL2, DMD, LDB3, CSRP3, ACTA1, PPARG were identified as hub genes. The miRNA-mRNA regulation network showed that hsa-miR-513b was the DEM with the most regulation, and COL1A1 was the DEG with the most regulation. In addition, CDC6, AURKA, CCNB1 and CCNA2 were related to overall survival and tumor differentiation. CONCLUSIONS: The regulatory relationship of hsa-miR-513b/ CDC6, CCNB1, CCNA2 and the regulatory relationship of hsa-miR-342-5p /AURKA were not only verified in the miRNA-mRNA regulatory network but also related to overall survival and tumor differentiation. These results indicated that they participated in the cellular regulatory process, and provided a molecular mechanism model for the study of pathogenesis.
- Klíčová slova
- bioinformatics, oral squamous cell carcinoma (OSCC), regulatory network,
- MeSH
- adrenalin MeSH
- aurora kinasa A genetika metabolismus MeSH
- dlaždicobuněčné karcinomy hlavy a krku * genetika MeSH
- genové regulační sítě MeSH
- kvalita života MeSH
- lidé MeSH
- messenger RNA genetika metabolismus MeSH
- mikro RNA * genetika MeSH
- nádory úst * genetika MeSH
- PPAR gama genetika metabolismus MeSH
- regulace genové exprese u nádorů MeSH
- stanovení celkové genové exprese MeSH
- transkripční faktory genetika MeSH
- výpočetní biologie metody MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- Názvy látek
- adrenalin MeSH
- aurora kinasa A MeSH
- messenger RNA MeSH
- mikro RNA * MeSH
- PPAR gama MeSH
- transkripční faktory MeSH
BACKGROUND: Integration of multi-omics data can provide a more complex view of the biological system consisting of different interconnected molecular components, the crucial aspect for developing novel personalised therapeutic strategies for complex diseases. Various tools have been developed to integrate multi-omics data. However, an efficient multi-omics framework for regulatory network inference at the genome level that incorporates prior knowledge is still to emerge. RESULTS: We present IntOMICS, an efficient integrative framework based on Bayesian networks. IntOMICS systematically analyses gene expression, DNA methylation, copy number variation and biological prior knowledge to infer regulatory networks. IntOMICS complements the missing biological prior knowledge by so-called empirical biological knowledge, estimated from the available experimental data. Regulatory networks derived from IntOMICS provide deeper insights into the complex flow of genetic information on top of the increasing accuracy trend compared to a published algorithm designed exclusively for gene expression data. The ability to capture relevant crosstalks between multi-omics modalities is verified using known associations in microsatellite stable/instable colon cancer samples. Additionally, IntOMICS performance is compared with two algorithms for multi-omics regulatory network inference that can also incorporate prior knowledge in the inference framework. IntOMICS is also applied to detect potential predictive biomarkers in microsatellite stable stage III colon cancer samples. CONCLUSIONS: We provide IntOMICS, a framework for multi-omics data integration using a novel approach to biological knowledge discovery. IntOMICS is a powerful resource for exploratory systems biology and can provide valuable insights into the complex mechanisms of biological processes that have a vital role in personalised medicine.
- Klíčová slova
- Bayesian networks, Integrative analysis, Knowledge discovery, Multimodal omics, Regulatory networks,
- MeSH
- algoritmy MeSH
- Bayesova věta MeSH
- genové regulační sítě MeSH
- lidé MeSH
- nádory tračníku * MeSH
- systémová biologie metody MeSH
- variabilita počtu kopií segmentů DNA * MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
BACKGROUND: Boolean networks (BNs) provide an effective modelling formalism for various complex biochemical phenomena. Their long term behaviour is represented by attractors-subsets of the state space towards which the BN eventually converges. These are then typically linked to different biological phenotypes. Depending on various logical parameters, the structure and quality of attractors can undergo a significant change, known as a bifurcation. We present a methodology for analysing bifurcations in asynchronous parametrised Boolean networks. RESULTS: In this paper, we propose a computational framework employing advanced symbolic graph algorithms that enable the analysis of large networks with hundreds of Boolean variables. To visualise the results of this analysis, we developed a novel interactive presentation technique based on decision trees, allowing us to quickly uncover parameters crucial to the changes in the attractor landscape. As a whole, the methodology is implemented in our tool AEON. We evaluate the method's applicability on a complex human cell signalling network describing the activity of type-1 interferons and related molecules interacting with SARS-COV-2 virion. In particular, the analysis focuses on explaining the potential suppressive role of the recently proposed drug molecule GRL0617 on replication of the virus. CONCLUSIONS: The proposed method creates a working analogy to the concept of bifurcation analysis widely used in kinetic modelling to reveal the impact of parameters on the system's stability. The important feature of our tool is its unique capability to work fast with large-scale networks with a relatively large extent of unknown information. The results obtained in the case study are in agreement with the recent biological findings.
- Klíčová slova
- Attractor bifurcation, Boolean networks, Software tool, Symbolic computation, type-1 interferons,
- MeSH
- algoritmy MeSH
- aniliny MeSH
- benzamidy MeSH
- COVID-19 * MeSH
- genové regulační sítě * MeSH
- lidé MeSH
- modely genetické MeSH
- naftaleny MeSH
- SARS-CoV-2 MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- Názvy látek
- 5-amino-2-methyl-N-((R)-1-(1-naphthyl)ethyl)benzamide MeSH Prohlížeč
- aniliny MeSH
- benzamidy MeSH
- naftaleny MeSH
