No Arabic abstract
Cell division timing is critical for cell fate specification and morphogenesis during embryogenesis. How division timings are regulated among cells during development is poorly understood. Here we focus on the comparison of asynchrony of division between sister cells (ADS) between wild-type and mutant individuals of Caenorhabditis elegans. Since the replicate number of mutant individuals of each mutated gene, usually one, is far smaller than that of wild-type, direct comparison of two distributions of ADS between wild-type and mutant type, such as Kolmogorov- Smirnov test, is not feasible. On the other hand, we find that sometimes ADS is correlated with the life span of corresponding mother cell in wild-type. Hence, we apply a semiparametric Bayesian quantile regression method to estimate the 95% confidence interval curve of ADS with respect to life span of mother cell of wild-type individuals. Then, mutant-type ADSs outside the corresponding confidence interval are selected out as abnormal one with a significance level of 0.05. Simulation study demonstrates the accuracy of our method and Gene Enrichment Analysis validates the results of real data sets.
Aging in Caenorhabditis elegans is controlled, in part, by the insulin-like signaling and heat shock response pathways. Following thermal stress, expression levels of small heat shock protein 16.2 show a spatial patterning across the 20 intestinal cells that reside along the length of the worm. Here, we present a hypothesized mechanism that could lead to this patterned response and develop a mathematical model of this system to test our hypothesis. We propose that the patterned expression of heat shock protein is caused by a diffusion-driven instability within the pseudocoelom, or fluid-filled cavity, that borders the intestinal cells in C. elegans. This instability is due to the interactions between two classes of insulin like peptides that serve antagonistic roles. We examine output from the developed model and compare it to experimental data on heat shock protein expression. Furthermore, we use the model to gain insight on possible biological parameters in the system. The model presented is capable of producing patterns similar to what is observed experimentally and provides a first step in mathematically modeling aging-related mechanisms in C. elegans.
To study how a zygote develops into an embryo with different tissues, large-scale 4D confocal movies of C. elegans embryos have been produced recently by experimental biologists. However, the lack of principled statistical methods for the highly noisy data has hindered the comprehensive analysis of these data sets. We introduced a probabilistic change point model on the cell lineage tree to estimate the embryonic gene expression onset time. A Bayesian approach is used to fit the 4D confocal movies data to the model. Subsequent classification methods are used to decide a model selection threshold and further refine the expression onset time from the branch level to the specific cell time level. Extensive simulations have shown the high accuracy of our method. Its application on real data yields both previously known results and new findings.
Nematodes have evolved to swim in highly viscous environments. Artificial mechanisms that mimic the locomotory functions of nematodes can be efficient viscous pumps. We experimentally simulate the motion of the head segment of Caenorhabditis elegans by introducing a reciprocating and rocking blade. We show that the bio-inspired blades motion not only induces a flow structure similar to that of the worm, but also mixes the surrounding fluid by generating a circulatory flow. When confined between two parallel walls, the blade causes a steady Poiseuille flow through closed circuits. The pumping efficiency is comparable with the swimming efficiency of the worm. If implanted in a sealed chamber and actuated remotely, the blade can provide pumping and mixing functions for microprocessors cooled by polymeric flows and microfluidic devices.
The connectome, or the entire connectivity of a neural system represented by network, ranges various scales from synaptic connections between individual neurons to fibre tract connections between brain regions. Although the modularity they commonly show has been extensively studied, it is unclear whether connection specificity of such networks can already be fully explained by the modularity alone. To answer this question, we study two networks, the neuronal network of C. elegans and the fibre tract network of human brains yielded through diffusion spectrum imaging (DSI). We compare them to their respective benchmark networks with varying modularities, which are generated by link swapping to have desired modularity values but otherwise maximally random. We find several network properties that are specific to the neural networks and cannot be fully explained by the modularity alone. First, the clustering coefficient and the characteristic path length of C. elegans and human connectomes are both higher than those of the benchmark networks with similar modularity. High clustering coefficient indicates efficient local information distribution and high characteristic path length suggests reduced global integration. Second, the total wiring length is smaller than for the alternative configurations with similar modularity. This is due to lower dispersion of connections, which means each neuron in C. elegans connectome or each region of interest (ROI) in human connectome reaches fewer ganglia or cortical areas, respectively. Third, both neural networks show lower algorithmic entropy compared to the alternative arrangements. This implies that fewer rules are needed to encode for the organisation of neural systems.
Graph theoretical analyses of nervous systems usually omit the aspect of connection polarity, due to data insufficiency. The chemical synapse network of Caenorhabditis elegans is a well-reconstructed directed network, but the signs of its connections are yet to be elucidated. Here, we present the gene expression-based sign prediction of the ionotropic chemical synapse connectome of C. elegans (3,638 connections and 20,589 synapses total), incorporating available presynaptic neurotransmitter and postsynaptic receptor gene expression data for three major neurotransmitter systems. We made predictions for more than two-thirds of these chemical synapses and observed an excitatory-inhibitory (E:I) ratio close to 4:1 which was found similar to that observed in many real-world networks. Our open source tool (http://EleganSign.linkgroup.hu) is simple but efficient in predicting polarities by integrating neuronal connectome and gene expression data.