@article{7463,
  abstract     = {Resting-state brain activity is characterized by the presence of neuronal avalanches showing absence of characteristic size. Such evidence has been interpreted in the context of criticality and associated with the normal functioning of the brain. A distinctive attribute of systems at criticality is the presence of long-range correlations. Thus, to verify the hypothesis that the brain operates close to a critical point and consequently assess deviations from criticality for diagnostic purposes, it is of primary importance to robustly and reliably characterize correlations in resting-state brain activity. Recent works focused on the analysis of narrow-band electroencephalography (EEG) and magnetoencephalography (MEG) signal amplitude envelope, showing evidence of long-range temporal correlations (LRTC) in neural oscillations. However, brain activity is a broadband phenomenon, and a significant piece of information useful to precisely discriminate between normal (critical) and pathological behavior (non-critical), may be encoded in the broadband spatio-temporal cortical dynamics. Here we propose to characterize the temporal correlations in the broadband brain activity through the lens of neuronal avalanches. To this end, we consider resting-state EEG and long-term MEG recordings, extract the corresponding neuronal avalanche sequences, and study their temporal correlations. We demonstrate that the broadband resting-state brain activity consistently exhibits long-range power-law correlations in both EEG and MEG recordings, with similar values of the scaling exponents. Importantly, although we observe that the avalanche size distribution depends on scale parameters, scaling exponents characterizing long-range correlations are quite robust. In particular, they are independent of the temporal binning (scale of analysis), indicating that our analysis captures intrinsic characteristics of the underlying dynamics. Because neuronal avalanches constitute a fundamental feature of neural systems with universal characteristics, the proposed approach may serve as a general, systems- and experiment-independent procedure to infer the existence of underlying long-range correlations in extended neural systems, and identify pathological behaviors in the complex spatio-temporal interplay of cortical rhythms.},
  author       = {Lombardi, Fabrizio and Shriki, Oren and Herrmann, Hans J and de Arcangelis, Lucilla},
  issn         = {1872-8286},
  journal      = {Neurocomputing},
  pages        = {657--666},
  publisher    = {Elsevier},
  title        = {{Long-range temporal correlations in the broadband resting state activity of the human brain revealed by neuronal avalanches}},
  doi          = {10.1016/j.neucom.2020.05.126},
  volume       = {461},
  year         = {2021},
}

@article{7553,
  abstract     = {Normative theories and statistical inference provide complementary approaches for the study of biological systems. A normative theory postulates that organisms have adapted to efficiently solve essential tasks, and proceeds to mathematically work out testable consequences of such optimality; parameters that maximize the hypothesized organismal function can be derived ab initio, without reference to experimental data. In contrast, statistical inference focuses on efficient utilization of data to learn model parameters, without reference to any a priori notion of biological function, utility, or fitness. Traditionally, these two approaches were developed independently and applied separately. Here we unify them in a coherent Bayesian framework that embeds a normative theory into a family of maximum-entropy “optimization priors.” This family defines a smooth interpolation between a data-rich inference regime (characteristic of “bottom-up” statistical models), and a data-limited ab inito prediction regime (characteristic of “top-down” normative theory). We demonstrate the applicability of our framework using data from the visual cortex, and argue that the flexibility it affords is essential to address a number of fundamental challenges relating to inference and prediction in complex, high-dimensional biological problems.},
  author       = {Mlynarski, Wiktor F and Hledik, Michal and Sokolowski, Thomas R and Tkačik, Gašper},
  journal      = {Neuron},
  number       = {7},
  pages        = {1227--1241.e5},
  publisher    = {Cell Press},
  title        = {{Statistical analysis and optimality of neural systems}},
  doi          = {10.1016/j.neuron.2021.01.020},
  volume       = {109},
  year         = {2021},
}

@unpublished{10077,
  abstract     = {Although much is known about how single neurons in the hippocampus represent an animal’s position, how cell-cell interactions contribute to spatial coding remains poorly understood. Using a novel statistical estimator and theoretical modeling, both developed in the framework of maximum entropy models, we reveal highly structured cell-to-cell interactions whose statistics depend on familiar vs. novel environment. In both conditions the circuit interactions optimize the encoding of spatial information, but for regimes that differ in the signal-to-noise ratio of their spatial inputs. Moreover, the topology of the interactions facilitates linear decodability, making the information easy to read out by downstream circuits. These findings suggest that the efficient coding hypothesis is not applicable only to individual neuron properties in the sensory periphery, but also to neural interactions in the central brain.},
  author       = {Nardin, Michele and Csicsvari, Jozsef L and Tkačik, Gašper and Savin, Cristina},
  booktitle    = {bioRxiv},
  publisher    = {Cold Spring Harbor Laboratory},
  title        = {{The structure of hippocampal CA1 interactions optimizes spatial coding across experience}},
  doi          = {10.1101/2021.09.28.460602},
  year         = {2021},
}

@article{10535,
  abstract     = {Realistic models of biological processes typically involve interacting components on multiple scales, driven by changing environment and inherent stochasticity. Such models are often analytically and numerically intractable. We revisit a dynamic maximum entropy method that combines a static maximum entropy with a quasi-stationary approximation. This allows us to reduce stochastic non-equilibrium dynamics expressed by the Fokker-Planck equation to a simpler low-dimensional deterministic dynamics, without the need to track microscopic details. Although the method has been previously applied to a few (rather complicated) applications in population genetics, our main goal here is to explain and to better understand how the method works. We demonstrate the usefulness of the method for two widely studied stochastic problems, highlighting its accuracy in capturing important macroscopic quantities even in rapidly changing non-stationary conditions. For the Ornstein-Uhlenbeck process, the method recovers the exact dynamics whilst for a stochastic island model with migration from other habitats, the approximation retains high macroscopic accuracy under a wide range of scenarios in a dynamic environment.},
  author       = {Bod'ová, Katarína and Szep, Eniko and Barton, Nicholas H},
  issn         = {1553-7358},
  journal      = {PLoS Computational Biology},
  number       = {12},
  publisher    = {Public Library of Science},
  title        = {{Dynamic maximum entropy provides accurate approximation of structured population dynamics}},
  doi          = {10.1371/journal.pcbi.1009661},
  volume       = {17},
  year         = {2021},
}

@unpublished{10579,
  abstract     = {We consider a totally asymmetric simple exclusion process (TASEP) consisting of particles on a lattice that require binding by a "token" to move. Using a combination of theory and simulations, we address the following questions: (i) How token binding kinetics affects the current-density relation; (ii) How the current-density relation depends on the scarcity of tokens; (iii) How tokens propagate the effects of the locally-imposed disorder (such a slow site) over the entire lattice; (iv) How a shared pool of tokens couples concurrent TASEPs running on multiple lattices; (v) How our results translate to TASEPs with open boundaries that exchange particles with the reservoir. Since real particle motion (including in systems that inspired the standard TASEP model, e.g., protein synthesis or movement of molecular motors) is often catalyzed, regulated, actuated, or otherwise mediated, the token-driven TASEP dynamics analyzed in this paper should allow for a better understanding of real systems and enable a closer match between TASEP theory and experimental observations.},
  author       = {Kavcic, Bor and Tkačik, Gašper},
  booktitle    = {arXiv},
  title        = {{Token-driven totally asymmetric simple exclusion process}},
  doi          = {10.48550/arXiv.2112.13558},
  year         = {2021},
}

@article{9439,
  abstract     = {The ability to adapt to changes in stimulus statistics is a hallmark of sensory systems. Here, we developed a theoretical framework that can account for the dynamics of adaptation from an information processing perspective. We use this framework to optimize and analyze adaptive sensory codes, and we show that codes optimized for stationary environments can suffer from prolonged periods of poor performance when the environment changes. To mitigate the adversarial effects of these environmental changes, sensory systems must navigate tradeoffs between the ability to accurately encode incoming stimuli and the ability to rapidly detect and adapt to changes in the distribution of these stimuli. We derive families of codes that balance these objectives, and we demonstrate their close match to experimentally observed neural dynamics during mean and variance adaptation. Our results provide a unifying perspective on adaptation across a range of sensory systems, environments, and sensory tasks.},
  author       = {Mlynarski, Wiktor F and Hermundstad, Ann M.},
  issn         = {1546-1726},
  journal      = {Nature Neuroscience},
  pages        = {998--1009},
  publisher    = {Springer Nature},
  title        = {{Efficient and adaptive sensory codes}},
  doi          = {10.1038/s41593-021-00846-0},
  volume       = {24},
  year         = {2021},
}

@article{9822,
  abstract     = {Attachment of adhesive molecules on cell culture surfaces to restrict cell adhesion to defined areas and shapes has been vital for the progress of in vitro research. In currently existing patterning methods, a combination of pattern properties such as stability, precision, specificity, high-throughput outcome, and spatiotemporal control is highly desirable but challenging to achieve. Here, we introduce a versatile and high-throughput covalent photoimmobilization technique, comprising a light-dose-dependent patterning step and a subsequent functionalization of the pattern via click chemistry. This two-step process is feasible on arbitrary surfaces and allows for generation of sustainable patterns and gradients. The method is validated in different biological systems by patterning adhesive ligands on cell-repellent surfaces, thereby constraining the growth and migration of cells to the designated areas. We then implement a sequential photopatterning approach by adding a second switchable patterning step, allowing for spatiotemporal control over two distinct surface patterns. As a proof of concept, we reconstruct the dynamics of the tip/stalk cell switch during angiogenesis. Our results show that the spatiotemporal control provided by our “sequential photopatterning” system is essential for mimicking dynamic biological processes and that our innovative approach has great potential for further applications in cell science.},
  author       = {Zisis, Themistoklis and Schwarz, Jan and Balles, Miriam and Kretschmer, Maibritt and Nemethova, Maria and Chait, Remy P and Hauschild, Robert and Lange, Janina and Guet, Calin C and Sixt, Michael K and Zahler, Stefan},
  issn         = {19448252},
  journal      = {ACS Applied Materials and Interfaces},
  number       = {30},
  pages        = {35545–35560},
  publisher    = {American Chemical Society},
  title        = {{Sequential and switchable patterning for studying cellular processes under spatiotemporal control}},
  doi          = {10.1021/acsami.1c09850},
  volume       = {13},
  year         = {2021},
}

@article{9828,
  abstract     = {Amplitude demodulation is a classical operation used in signal processing. For a long time, its effective applications in practice have been limited to narrowband signals. In this work, we generalize amplitude demodulation to wideband signals. We pose demodulation as a recovery problem of an oversampled corrupted signal and introduce special iterative schemes belonging to the family of alternating projection algorithms to solve it. Sensibly chosen structural assumptions on the demodulation outputs allow us to reveal the high inferential accuracy of the method over a rich set of relevant signals. This new approach surpasses current state-of-the-art demodulation techniques apt to wideband signals in computational efficiency by up to many orders of magnitude with no sacrifice in quality. Such performance opens the door for applications of the amplitude demodulation procedure in new contexts. In particular, the new method makes online and large-scale offline data processing feasible, including the calculation of modulator-carrier pairs in higher dimensions and poor sampling conditions, independent of the signal bandwidth. We illustrate the utility and specifics of applications of the new method in practice by using natural speech and synthetic signals.},
  author       = {Gabrielaitis, Mantas},
  issn         = {1941-0476},
  journal      = {IEEE Transactions on Signal Processing},
  pages        = {4039 -- 4054},
  publisher    = {Institute of Electrical and Electronics Engineers},
  title        = {{Fast and accurate amplitude demodulation of wideband signals}},
  doi          = {10.1109/TSP.2021.3087899},
  volume       = {69},
  year         = {2021},
}

@unpublished{10912,
  abstract     = {Brain dynamics display collective phenomena as diverse as neuronal oscillations and avalanches. Oscillations are rhythmic, with fluctuations occurring at a characteristic scale, whereas avalanches are scale-free cascades of neural activity. Here we show that such antithetic features can coexist in a very generic class of adaptive neural networks. In the most simple yet fully microscopic model from this class we make direct contact with human brain resting-state activity recordings via tractable inference of the model's two essential parameters. The inferred model quantitatively captures the dynamics over a broad range of scales, from single sensor fluctuations, collective behaviors of nearly-synchronous extreme events on multiple sensors, to neuronal avalanches unfolding over multiple sensors across multiple time-bins. Importantly, the inferred parameters correlate with model-independent signatures of "closeness to criticality", suggesting that the coexistence of scale-specific (neural oscillations) and scale-free (neuronal avalanches) dynamics in brain activity occurs close to a non-equilibrium critical point at the onset of self-sustained oscillations.},
  author       = {Lombardi, Fabrizio and Pepic, Selver and Shriki, Oren and Tkačik, Gašper and De Martino, Daniele},
  pages        = {37},
  publisher    = {arXiv},
  title        = {{Quantifying the coexistence of neuronal oscillations and avalanches}},
  doi          = {10.48550/ARXIV.2108.06686},
  year         = {2021},
}

@article{8250,
  abstract     = {Antibiotics that interfere with translation, when combined, interact in diverse and difficult-to-predict ways. Here, we explain these interactions by “translation bottlenecks”: points in the translation cycle where antibiotics block ribosomal progression. To elucidate the underlying mechanisms of drug interactions between translation inhibitors, we generate translation bottlenecks genetically using inducible control of translation factors that regulate well-defined translation cycle steps. These perturbations accurately mimic antibiotic action and drug interactions, supporting that the interplay of different translation bottlenecks causes these interactions. We further show that growth laws, combined with drug uptake and binding kinetics, enable the direct prediction of a large fraction of observed interactions, yet fail to predict suppression. However, varying two translation bottlenecks simultaneously supports that dense traffic of ribosomes and competition for translation factors account for the previously unexplained suppression. These results highlight the importance of “continuous epistasis” in bacterial physiology.},
  author       = {Kavcic, Bor and Tkačik, Gašper and Bollenbach, Tobias},
  issn         = {2041-1723},
  journal      = {Nature Communications},
  publisher    = {Springer Nature},
  title        = {{Mechanisms of drug interactions between translation-inhibiting antibiotics}},
  doi          = {10.1038/s41467-020-17734-z},
  volume       = {11},
  year         = {2020},
}

@phdthesis{8657,
  abstract     = {Synthesis of proteins – translation – is a fundamental process of life. Quantitative studies anchor translation into the context of bacterial physiology and reveal several mathematical relationships, called “growth laws,” which capture physiological feedbacks between protein synthesis and cell growth. Growth laws describe the dependency of the ribosome abundance as a function of growth rate, which can change depending on the growth conditions. Perturbations of translation reveal that bacteria employ a compensatory strategy in which the reduced translation capability results in increased expression of the translation machinery.
Perturbations of translation are achieved in various ways; clinically interesting is the application of translation-targeting antibiotics – translation inhibitors. The antibiotic effects on bacterial physiology are often poorly understood. Bacterial responses to two or more simultaneously applied antibiotics are even more puzzling. The combined antibiotic effect determines the type of drug interaction, which ranges from synergy (the effect is stronger than expected) to antagonism (the effect is weaker) and suppression (one of the drugs loses its potency).
In the first part of this work, we systematically measure the pairwise interaction network for translation inhibitors that interfere with different steps in translation. We find that the interactions are surprisingly diverse and tend to be more antagonistic. To explore the underlying mechanisms, we begin with a minimal biophysical model of combined antibiotic action. We base this model on the kinetics of antibiotic uptake and binding together with the physiological response described by the growth laws. The biophysical model explains some drug interactions, but not all; it specifically fails to predict suppression.
In the second part of this work, we hypothesize that elusive suppressive drug interactions result from the interplay between ribosomes halted in different stages of translation. To elucidate this putative mechanism of drug interactions between translation inhibitors, we generate translation bottlenecks genetically using in- ducible control of translation factors that regulate well-defined translation cycle steps. These perturbations accurately mimic antibiotic action and drug interactions, supporting that the interplay of different translation bottlenecks partially causes these interactions.
We extend this approach by varying two translation bottlenecks simultaneously. This approach reveals the suppression of translocation inhibition by inhibited translation. We rationalize this effect by modeling dense traffic of ribosomes that move on transcripts in a translation factor-mediated manner. This model predicts a dissolution of traffic jams caused by inhibited translocation when the density of ribosome traffic is reduced by lowered initiation. We base this model on the growth laws and quantitative relationships between different translation and growth parameters.
In the final part of this work, we describe a set of tools aimed at quantification of physiological and translation parameters. We further develop a simple model that directly connects the abundance of a translation factor with the growth rate, which allows us to extract physiological parameters describing initiation. We demonstrate the development of tools for measuring translation rate.
This thesis showcases how a combination of high-throughput growth rate mea- surements, genetics, and modeling can reveal mechanisms of drug interactions. Furthermore, by a gradual transition from combinations of antibiotics to precise genetic interventions, we demonstrated the equivalency between genetic and chemi- cal perturbations of translation. These findings tile the path for quantitative studies of antibiotic combinations and illustrate future approaches towards the quantitative description of translation.},
  author       = {Kavcic, Bor},
  isbn         = {978-3-99078-011-4},
  issn         = {2663-337X},
  pages        = {271},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Perturbations of protein synthesis: from antibiotics to genetics and physiology}},
  doi          = {10.15479/AT:ISTA:8657},
  year         = {2020},
}

@article{8698,
  abstract     = {The brain represents and reasons probabilistically about complex stimuli and motor actions using a noisy, spike-based neural code. A key building block for such neural computations, as well as the basis for supervised and unsupervised learning, is the ability to estimate the surprise or likelihood of incoming high-dimensional neural activity patterns. Despite progress in statistical modeling of neural responses and deep learning, current approaches either do not scale to large neural populations or cannot be implemented using biologically realistic mechanisms. Inspired by the sparse and random connectivity of real neuronal circuits, we present a model for neural codes that accurately estimates the likelihood of individual spiking patterns and has a straightforward, scalable, efficient, learnable, and realistic neural implementation. This model’s performance on simultaneously recorded spiking activity of >100 neurons in the monkey visual and prefrontal cortices is comparable with or better than that of state-of-the-art models. Importantly, the model can be learned using a small number of samples and using a local learning rule that utilizes noise intrinsic to neural circuits. Slower, structural changes in random connectivity, consistent with rewiring and pruning processes, further improve the efficiency and sparseness of the resulting neural representations. Our results merge insights from neuroanatomy, machine learning, and theoretical neuroscience to suggest random sparse connectivity as a key design principle for neuronal computation.},
  author       = {Maoz, Ori and Tkačik, Gašper and Esteki, Mohamad Saleh and Kiani, Roozbeh and Schneidman, Elad},
  issn         = {10916490},
  journal      = {Proceedings of the National Academy of Sciences of the United States of America},
  number       = {40},
  pages        = {25066--25073},
  publisher    = {National Academy of Sciences},
  title        = {{Learning probabilistic neural representations with randomly connected circuits}},
  doi          = {10.1073/pnas.1912804117},
  volume       = {117},
  year         = {2020},
}

@misc{8930,
  abstract     = {Phenomenological relations such as Ohm’s or Fourier’s law have a venerable history in physics but are still scarce in biology. This situation restrains predictive theory. Here, we build on bacterial “growth laws,” which capture physiological feedback between translation and cell growth, to construct a minimal biophysical model for the combined action of ribosome-targeting antibiotics. Our model predicts drug interactions like antagonism or synergy solely from responses to individual drugs. We provide analytical results for limiting cases, which agree well with numerical results. We systematically refine the model by including direct physical interactions of different antibiotics on the ribosome. In a limiting case, our model provides a mechanistic underpinning for recent predictions of higher-order interactions that were derived using entropy maximization. We further refine the model to include the effects of antibiotics that mimic starvation and the presence of resistance genes. We describe the impact of a starvation-mimicking antibiotic on drug interactions analytically and verify it experimentally. Our extended model suggests a change in the type of drug interaction that depends on the strength of resistance, which challenges established rescaling paradigms. We experimentally show that the presence of unregulated resistance genes can lead to altered drug interaction, which agrees with the prediction of the model. While minimal, the model is readily adaptable and opens the door to predicting interactions of second and higher-order in a broad range of biological systems.},
  author       = {Kavcic, Bor},
  keywords     = {Escherichia coli, antibiotic combinations, translation, growth laws, drug interactions, bacterial physiology, translation inhibitors},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Analysis scripts and research data for the paper "Minimal biophysical model of combined antibiotic action"}},
  doi          = {10.15479/AT:ISTA:8930},
  year         = {2020},
}

@article{8955,
  abstract     = {Skeletal muscle activity is continuously modulated across physiologic states to provide coordination, flexibility and responsiveness to body tasks and external inputs. Despite the central role the muscular system plays in facilitating vital body functions, the network of brain-muscle interactions required to control hundreds of muscles and synchronize their activation in relation to distinct physiologic states has not been investigated. Recent approaches have focused on general associations between individual brain rhythms and muscle activation during movement tasks. However, the specific forms of coupling, the functional network of cortico-muscular coordination, and how network structure and dynamics are modulated by autonomic regulation across physiologic states remains unknown. To identify and quantify the cortico-muscular interaction network and uncover basic features of neuro-autonomic control of muscle function, we investigate the coupling between synchronous bursts in cortical rhythms and peripheral muscle activation during sleep and wake. Utilizing the concept of time delay stability and a novel network physiology approach, we find that the brain-muscle network exhibits complex dynamic patterns of communication involving multiple brain rhythms across cortical locations and different electromyographic frequency bands. Moreover, our results show that during each physiologic state the cortico-muscular network is characterized by a specific profile of network links strength, where particular brain rhythms play role of main mediators of interaction and control. Further, we discover a hierarchical reorganization in network structure across physiologic states, with high connectivity and network link strength during wake, intermediate during REM and light sleep, and low during deep sleep, a sleep-stage stratification that demonstrates a unique association between physiologic states and cortico-muscular network structure. The reported empirical observations are consistent across individual subjects, indicating universal behavior in network structure and dynamics, and high sensitivity of cortico-muscular control to changes in autonomic regulation, even at low levels of physical activity and muscle tone during sleep. Our findings demonstrate previously unrecognized basic principles of brain-muscle network communication and control, and provide new perspectives on the regulatory mechanisms of brain dynamics and locomotor activation, with potential clinical implications for neurodegenerative, movement and sleep disorders, and for developing efficient treatment strategies.},
  author       = {Rizzo, Rossella and Zhang, Xiyun and Wang, Jilin W.J.L. and Lombardi, Fabrizio and Ivanov, Plamen Ch},
  issn         = {1664042X},
  journal      = {Frontiers in Physiology},
  publisher    = {Frontiers},
  title        = {{Network physiology of cortico–muscular interactions}},
  doi          = {10.3389/fphys.2020.558070},
  volume       = {11},
  year         = {2020},
}

@article{9000,
  abstract     = {In prokaryotes, thermodynamic models of gene regulation provide a highly quantitative mapping from promoter sequences to gene-expression levels that is compatible with in vivo and in vitro biophysical measurements. Such concordance has not been achieved for models of enhancer function in eukaryotes. In equilibrium models, it is difficult to reconcile the reported short transcription factor (TF) residence times on the DNA with the high specificity of regulation. In nonequilibrium models, progress is difficult due to an explosion in the number of parameters. Here, we navigate this complexity by looking for minimal nonequilibrium enhancer models that yield desired regulatory phenotypes: low TF residence time, high specificity, and tunable cooperativity. We find that a single extra parameter, interpretable as the “linking rate,” by which bound TFs interact with Mediator components, enables our models to escape equilibrium bounds and access optimal regulatory phenotypes, while remaining consistent with the reported phenomenology and simple enough to be inferred from upcoming experiments. We further find that high specificity in nonequilibrium models is in a trade-off with gene-expression noise, predicting bursty dynamics—an experimentally observed hallmark of eukaryotic transcription. By drastically reducing the vast parameter space of nonequilibrium enhancer models to a much smaller subspace that optimally realizes biological function, we deliver a rich class of models that could be tractably inferred from data in the near future.},
  author       = {Grah, Rok and Zoller, Benjamin and Tkačik, Gašper},
  issn         = {10916490},
  journal      = {PNAS},
  number       = {50},
  pages        = {31614--31622},
  publisher    = {National Academy of Sciences},
  title        = {{Nonequilibrium models of optimal enhancer function}},
  doi          = {10.1073/pnas.2006731117},
  volume       = {117},
  year         = {2020},
}

@misc{7383,
  abstract     = {Organisms cope with change by employing transcriptional regulators. However, when faced with rare environments, the evolution of transcriptional regulators and their promoters may be too slow. We ask whether the intrinsic instability of gene duplication and amplification provides a generic alternative to canonical gene regulation. By real-time monitoring of gene copy number mutations in E. coli, we show that gene duplications and amplifications enable adaptation to fluctuating environments by rapidly generating copy number, and hence expression level, polymorphism. This ‘amplification-mediated gene expression tuning’ occurs on timescales similar to canonical gene regulation and can deal with rapid environmental changes. Mathematical modeling shows that amplifications also tune gene expression in stochastic environments where transcription factor-based schemes are hard to evolve or maintain. The fleeting nature of gene amplifications gives rise to a generic population-level mechanism that relies on genetic heterogeneity to rapidly tune expression of any gene, without leaving any genomic signature.},
  author       = {Grah, Rok},
  keywords     = {Matlab scripts, analysis of microfluidics, mathematical model},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Matlab scripts for the Paper: Gene Amplification as a Form of Population-Level Gene Expression regulation}},
  doi          = {10.15479/AT:ISTA:7383},
  year         = {2020},
}

@article{7490,
  abstract     = {In plants, clathrin mediated endocytosis (CME) represents the major route for cargo internalisation from the cell surface. It has been assumed to operate in an evolutionary conserved manner as in yeast and animals. Here we report characterisation of ultrastructure, dynamics and mechanisms of plant CME as allowed by our advancement in electron microscopy and quantitative live imaging techniques. Arabidopsis CME appears to follow the constant curvature model and the bona fide CME population generates vesicles of a predominantly hexagonal-basket type; larger and with faster kinetics than in other models. Contrary to the existing paradigm, actin is dispensable for CME events at the plasma membrane but plays a unique role in collecting endocytic vesicles, sorting of internalised cargos and directional endosome movement that itself actively promote CME events. Internalized vesicles display a strongly delayed and sequential uncoating. These unique features highlight the independent evolution of the plant CME mechanism during the autonomous rise of multicellularity in eukaryotes.},
  author       = {Narasimhan, Madhumitha and Johnson, Alexander J and Prizak, Roshan and Kaufmann, Walter and Tan, Shutang and Casillas Perez, Barbara E and Friml, Jiří},
  issn         = {2050-084X},
  journal      = {eLife},
  publisher    = {eLife Sciences Publications},
  title        = {{Evolutionarily unique mechanistic framework of clathrin-mediated endocytosis in plants}},
  doi          = {10.7554/eLife.52067},
  volume       = {9},
  year         = {2020},
}

@article{7569,
  abstract     = {Genes differ in the frequency at which they are expressed and in the form of regulation used to control their activity. In particular, positive or negative regulation can lead to activation of a gene in response to an external signal. Previous works proposed that the form of regulation of a gene correlates with its frequency of usage: positive regulation when the gene is frequently expressed and negative regulation when infrequently expressed. Such network design means that, in the absence of their regulators, the genes are found in their least required activity state, hence regulatory intervention is often necessary. Due to the multitude of genes and regulators, spurious binding and unbinding events, called “crosstalk”, could occur. To determine how the form of regulation affects the global crosstalk in the network, we used a mathematical model that includes multiple regulators and multiple target genes. We found that crosstalk depends non-monotonically on the availability of regulators. Our analysis showed that excess use of regulation entailed by the formerly suggested network design caused high crosstalk levels in a large part of the parameter space. We therefore considered the opposite ‘idle’ design, where the default unregulated state of genes is their frequently required activity state. We found, that ‘idle’ design minimized the use of regulation and thus minimized crosstalk. In addition, we estimated global crosstalk of S. cerevisiae using transcription factors binding data. We demonstrated that even partial network data could suffice to estimate its global crosstalk, suggesting its applicability to additional organisms. We found that S. cerevisiae estimated crosstalk is lower than that of a random network, suggesting that natural selection reduces crosstalk. In summary, our study highlights a new type of protein production cost which is typically overlooked: that of regulatory interference caused by the presence of excess regulators in the cell. It demonstrates the importance of whole-network descriptions, which could show effects missed by single-gene models.},
  author       = {Grah, Rok and Friedlander, Tamar},
  issn         = {1553-7358},
  journal      = {PLOS Computational Biology},
  number       = {2},
  publisher    = {Public Library of Science},
  title        = {{The relation between crosstalk and gene regulation form revisited}},
  doi          = {10.1371/journal.pcbi.1007642},
  volume       = {16},
  year         = {2020},
}

@article{7652,
  abstract     = {Organisms cope with change by taking advantage of transcriptional regulators. However, when faced with rare environments, the evolution of transcriptional regulators and their promoters may be too slow. Here, we investigate whether the intrinsic instability of gene duplication and amplification provides a generic alternative to canonical gene regulation. Using real-time monitoring of gene-copy-number mutations in Escherichia coli, we show that gene duplications and amplifications enable adaptation to fluctuating environments by rapidly generating copy-number and, therefore, expression-level polymorphisms. This amplification-mediated gene expression tuning (AMGET) occurs on timescales that are similar to canonical gene regulation and can respond to rapid environmental changes. Mathematical modelling shows that amplifications also tune gene expression in stochastic environments in which transcription-factor-based schemes are hard to evolve or maintain. The fleeting nature of gene amplifications gives rise to a generic population-level mechanism that relies on genetic heterogeneity to rapidly tune the expression of any gene, without leaving any genomic signature.},
  author       = {Tomanek, Isabella and Grah, Rok and Lagator, M. and Andersson, A. M. C. and Bollback, Jonathan P and Tkačik, Gašper and Guet, Calin C},
  issn         = {2397-334X},
  journal      = {Nature Ecology & Evolution},
  number       = {4},
  pages        = {612--625},
  publisher    = {Springer Nature},
  title        = {{Gene amplification as a form of population-level gene expression regulation}},
  doi          = {10.1038/s41559-020-1132-7},
  volume       = {4},
  year         = {2020},
}

@article{7656,
  abstract     = {We propose that correlations among neurons are generically strong enough to organize neural activity patterns into a discrete set of clusters, which can each be viewed as a population codeword. Our reasoning starts with the analysis of retinal ganglion cell data using maximum entropy models, showing that the population is robustly in a frustrated, marginally sub-critical, or glassy, state. This leads to an argument that neural populations in many other brain areas might share this structure. Next, we use latent variable models to show that this glassy state possesses well-defined clusters of neural activity. Clusters have three appealing properties: (i) clusters exhibit error correction, i.e., they are reproducibly elicited by the same stimulus despite variability at the level of constituent neurons; (ii) clusters encode qualitatively different visual features than their constituent neurons; and (iii) clusters can be learned by downstream neural circuits in an unsupervised fashion. We hypothesize that these properties give rise to a “learnable” neural code which the cortical hierarchy uses to extract increasingly complex features without supervision or reinforcement.},
  author       = {Berry, Michael J. and Tkačik, Gašper},
  issn         = {16625188},
  journal      = {Frontiers in Computational Neuroscience},
  publisher    = {Frontiers},
  title        = {{Clustering of neural activity: A design principle for population codes}},
  doi          = {10.3389/fncom.2020.00020},
  volume       = {14},
  year         = {2020},
}

