Finding flies in the mushroom soup : Host specificity of fungus-associated communities revisited with a novel molecular method

Fruiting bodies of fungi constitute an important resource for thousands of other taxa. The structure of these diverse assemblages has traditionally been studied with labour-intensive methods involving cultivation and morphology-based species identification, to which molecular information might offer convenient complements. To overcome challenges in DNA extraction and PCR associated with the complex chemical properties of fruiting bodies, we developed a pipeline applicable for extracting amplifiable total DNA from soft fungal samples of any size. Our protocol purifies DNA in two sequential steps: (a) initial salt-isopropanol extraction of all nucleic acids in the sample is followed by (b) an extra clean-up step using solid-phase reversible immobilization (SPRI) magnetic beads. The protocol proved highly efficient, with practically all of our samples-regardless of biomass or other properties-being successfully PCR-amplified using metabarcoding primers and subsequently sequenced. As a proof of concept, we apply our methods to address a topical ecological question: is host specificity a major characteristic of fungus-associated communities, that is, do different fungus species harbour different communities of associated organisms? Based on an analysis of 312 fungal fruiting bodies representing 10 species in five genera from three orders, we show that molecular methods are suitable for studying this rich natural microcosm. Comparing to previous knowledge based on rearing and morphology-based identifications, we find a species-rich assemblage characterized by a low degree of host specialization. Our method opens up new horizons for molecular analyses of fungus-associated interaction webs and communities.Peer reviewe

Mathematical analysis of the Escherichia coli chemotaxis signalling pathway

We undertake a detailed mathematical analysis of a recent nonlinear ordinary differential equation (ODE) model describing the chemotactic signalling cascade within an {\it Escherichia coli} cell. The model includes a detailed description of the cell signalling cascade and an average approximation of the receptor activity. A steady-state stability analysis reveals the system exhibits one positive real steady-state which is shown to be asymptotically stable. Given the occurrence of a negative feedback between phosphorylated CheB (CheB-P) and the receptor state, we ask under what conditions, the system may exhibit oscillatory type behaviour. A detailed analysis of parameter space reveals that whilst variation in kinetic rate parameters within known biological limits is unlikely to lead to such behaviour, changes in the total concentration of the signalling proteins does. We postulate that experimentally observed overshoot behaviour can actually be described by damped oscillatory dynamics and consider the relationship between overshoot amplitude, total cell protein concentration and the magnitude of the external ligand stimulus. Model reductions of the full ODE model allow us to understand the link between phosphorylation events and the negative feedback between CheB-P and receptor methylation, as well as elucidate why some mathematical models exhibit overshoot and others do not. Our manuscript closes by discussing intercell variability of total protein concentration as means of ensuring the overall survival of a population as cells are subjected to different environments

Dependence of Bacterial Chemotaxis on Gradient Shape and Adaptation Rate

Simulation of cellular behavior on multiple scales requires models that are sufficiently detailed to capture central intracellular processes but at the same time enable the simulation of entire cell populations in a computationally cheap way. In this paper we present RapidCell, a hybrid model of chemotactic Escherichia coli that combines the Monod-Wyman-Changeux signal processing by mixed chemoreceptor clusters, the adaptation dynamics described by ordinary differential equations, and a detailed model of cell tumbling. Our model dramatically reduces computational costs and allows the highly efficient simulation of E. coli chemotaxis. We use the model to investigate chemotaxis in different gradients, and suggest a new, constant-activity type of gradient to systematically study chemotactic behavior of virtual bacteria. Using the unique properties of this gradient, we show that optimal chemotaxis is observed in a narrow range of CheA kinase activity, where concentration of the response regulator CheY-P falls into the operating range of flagellar motors. Our simulations also confirm that the CheB phosphorylation feedback improves chemotactic efficiency by shifting the average CheY-P concentration to fit the motor operating range. Our results suggest that in liquid media the variability in adaptation times among cells may be evolutionary favorable to ensure coexistence of subpopulations that will be optimally tactic in different gradients. However, in a porous medium (agar) such variability appears to be less important, because agar structure poses mainly negative selection against subpopulations with low levels of adaptation enzymes. RapidCell is available from the authors upon request

Cell-signalling dynamics in time and space

The specificity of cellular responses to receptor stimulation is encoded by the spatial and temporal dynamics of downstream signalling networks. Computational models provide insights into the intricate relationships between stimuli and responses and reveal mechanisms that enable networks to amplify signals, reduce noise and generate discontinuous bistable dynamics or oscillations. These temporal dynamics are coupled to precipitous spatial gradients of signalling activities, which guide pivotal intracellular processes, but also necessitate mechanisms to facilitate signal propagation across a cell