A divide-and-conquer method for analyzing high-dimensional noisy gene expression networks

Abstract: Intracellular gene expression systems are inevitably random due to low molecular counts. Consequently, mechanistic models for gene expression should be stochastic, and central to the analysis and inference of such models is solving the Chemical Master Equation (CME), which characterizes the probability evolution of the randomly evolving copy-numbers of the reacting species. While conventional methods such as Monte-Carlo simulations and finite state projections exist for estimating CME solutions, they suffer from the curse of dimensionality, significantly decreasing their efficacy for high-dimensional systems. Here, we propose a new computational method that resolves this issue through a novel divide-and-conquer approach. Our method divides the system into a leader system and several conditionally independent follower subsystems. The solution of the CME is then constructed by combining Monte Carlo estimation for the leader system with stochastic filtering procedures for the follower subsystems. We develop an optimized system decomposition, which ensures the low-dimensionality of the sub-problems, thereby allowing for improved scalability with increasing system dimension. The efficiency and accuracy of the method are demonstrated through several biologically relevant examples in high-dimensional estimation and inference problems. We demonstrate that our method can successfully identify a yeast transcription system at the single-cell resolution, leveraging mRNA time-course microscopy data, allowing us to rigorously examine the heterogeneity in rate parameters among isogenic cells cultured under identical conditions. Furthermore, we validate this finding using a novel noise decomposition technique introduced in this study. This technique exploits experimental time-course data to quantify intrinsic and extrinsic noise components, without requiring supplementary components, such as dual-reporter systems.

algebraic geometrydynamical systemsprobability

Audience: researchers in the topic

( video )

---

Seminar on the Mathematics of Reaction Networks

Series comments: Subscription link: list.ku.dk/postorius/lists/morn.list.ku.dk/

This seminar series focuses on progress in mathematical theory for the study of reaction networks, mainly in biology and chemistry. The scope is broad and accommodates works arising from dynamical systems, stochastics, algebra, topology and beyond.

We aim at providing a common forum for sharing knowledge and encouraging discussion across subfields. In particular we aim at facilitating interactions between junior and established researchers. These considerations will be represented in the choice of invited speakers and we will strive to create an excellent, exciting and diverse schedule.

The seminar runs twice a month, typically on the 2nd and 4th Thursday of the month, at 17:00 Brussels time (observe that this webpage shows the schedule in your current time zone). Each session consists of two 25-minute talks followed by 5-minute questions. After the two talks, longer discussions will take place for those interested. To this end, we will use breakout rooms. For this to work well, you need to have the latest version of Zoom installed (version 5.3.0 or higher), and use the desktop client or mobile app (not supported on ChromeOS).

We look forward hearing about new work and meeting many of you over zoom! Many of the talks are recorded; to see the recording, from Past Talks, open details of the listed talk for a video link.

The organizers.

---