Shifting patterns of natural variation in the nuclear genome of caenorhabditis elegans

<p>Abstract</p> <p>Background</p> <p>Genome wide analysis of variation within a species can reveal the evolution of fundamental biological processes such as mutation, recombination, and natural selection. We compare genome wide sequence differences between two independent isolates of the nematode <it>Caenorhabditis elegans </it>(CB4856 and CB4858) and the reference genome (N2).</p> <p>Results</p> <p>The base substitution pattern when comparing N2 against CB4858 reveals a transition over transversion bias (1.32:1) that is not present in CB4856. In CB4856, there is a significant bias in the direction of base substitution. The frequency of A or T bases in N2 that are G or C bases in CB4856 outnumber the opposite frequencies for transitions as well as transversions. These differences were not observed in the N2/CB4858 comparison. Similarly, we observed a strong bias for deletions over insertions in CB4856 (1.44: 1) that is not present in CB4858. In both CB4856 and CB4858, there is a significant correlation between SNP rate and recombination rate on the autosomes but not on the X chromosome. Furthermore, we identified numerous significant hotspots of variation in the CB4856-N2 comparison.</p> <p>In both CB4856 and CB4858, based on a measure of the strength of selection (k_a/k_s), all the chromosomes are under negative selection and in CB4856, there is no difference in the strength of natural selection in either the autosomes versus X or between any of the chromosomes. By contrast, in CB4858, k_a/k_svalues are smaller in the autosomes than in the X chromosome. In addition, in CB4858, k_a/k_svalues differ between chromosomes.</p> <p>Conclusions</p> <p>The clear bias of deletions over insertions in CB4856 suggests that either the CB4856 genome is becoming smaller or the N2 genome is getting larger. We hypothesize the hotspots found represent alleles that are shared between CB4856 and CB4858 but not N2. Because the k_a/k_sratio in the X chromosome is higher than the autosomes on average in CB4858, purifying selection is reduced on the X chromosome.</p

Simple, Fast and Accurate Implementation of the Diffusion Approximation Algorithm for Stochastic Ion Channels with Multiple States

The phenomena that emerge from the interaction of the stochastic opening and closing of ion channels (channel noise) with the non-linear neural dynamics are essential to our understanding of the operation of the nervous system. The effects that channel noise can have on neural dynamics are generally studied using numerical simulations of stochastic models. Algorithms based on discrete Markov Chains (MC) seem to be the most reliable and trustworthy, but even optimized algorithms come with a non-negligible computational cost. Diffusion Approximation (DA) methods use Stochastic Differential Equations (SDE) to approximate the behavior of a number of MCs, considerably speeding up simulation times. However, model comparisons have suggested that DA methods did not lead to the same results as in MC modeling in terms of channel noise statistics and effects on excitability. Recently, it was shown that the difference arose because MCs were modeled with coupled activation subunits, while the DA was modeled using uncoupled activation subunits. Implementations of DA with coupled subunits, in the context of a specific kinetic scheme, yielded similar results to MC. However, it remained unclear how to generalize these implementations to different kinetic schemes, or whether they were faster than MC algorithms. Additionally, a steady state approximation was used for the stochastic terms, which, as we show here, can introduce significant inaccuracies. We derived the SDE explicitly for any given ion channel kinetic scheme. The resulting generic equations were surprisingly simple and interpretable - allowing an easy and efficient DA implementation. The algorithm was tested in a voltage clamp simulation and in two different current clamp simulations, yielding the same results as MC modeling. Also, the simulation efficiency of this DA method demonstrated considerable superiority over MC methods.Comment: 32 text pages, 10 figures, 1 supplementary text + figur

On Conduction in a Bacterial Sodium Channel

Voltage-gated Na+-channels are transmembrane proteins that are responsible for the fast depolarizing phase of the action potential in nerve and muscular cells. Selective permeability of Na+ over Ca2+ or K+ ions is essential for the biological function of Na+-channels. After the emergence of the first high-resolution structure of a Na+-channel, an anionic coordination site was proposed to confer Na+ selectivity through partial dehydration of Na+ via its direct interaction with conserved glutamate side chains. By combining molecular dynamics simulations and free-energy calculations, a low-energy permeation pathway for Na+ ion translocation through the selectivity filter of the recently determined crystal structure of a prokaryotic sodium channel from Arcobacter butzleri is characterised. The picture that emerges is that of a pore preferentially occupied by two ions, which can switch between different configurations by crossing low free-energy barriers. In contrast to K+-channels, the movements of the ions appear to be weakly coupled in Na+-channels. When the free-energy maps for Na+ and K+ ions are compared, a selective site is characterised in the narrowest region of the filter, where a hydrated Na+ ion, and not a hydrated K+ ion, is energetically stable

Heritable determinants of male fertilization success in the nematode Caenorhabditis elegans

<p>Abstract</p> <p>Background</p> <p>Sperm competition is a driving force in the evolution of male sperm characteristics in many species. In the nematode <it>Caenorhabditis elegans</it>, larger male sperm evolve under experimentally increased sperm competition and larger male sperm outcompete smaller hermaphrodite sperm for fertilization within the hermaphrodite reproductive tract. To further elucidate the relative importance of sperm-related traits that contribute to differential reproductive success among males, we quantified within- and among-strain variation in sperm traits (size, rate of production, number transferred, competitive ability) for seven male genetic backgrounds known previously to differ with respect to some sperm traits. We also quantified male mating ability in assays for rates of courtship and successful copulation, and then assessed the roles of these pre- and post-mating traits in first- and second-male fertilization success.</p> <p>Results</p> <p>We document significant variation in courtship ability, mating ability, sperm size and sperm production rate. Sperm size and production rate were strong indicators of early fertilization success for males that mated second, but male genetic backgrounds conferring faster sperm production make smaller sperm, despite virgin males of all genetic backgrounds transferring indistinguishable numbers of sperm to mating partners.</p> <p>Conclusions</p> <p>We have demonstrated that sperm size and the rate of sperm production represent dominant factors in determining male fertilization success and that <it>C. elegans </it>harbors substantial heritable variation for traits contributing to male reproductive success. <it>C. elegans </it>provides a powerful, tractable system for studying sexual selection and for dissecting the genetic basis and evolution of reproduction-related traits.</p

Caenorhabditis briggsae Recombinant Inbred Line Genotypes Reveal Inter-Strain Incompatibility and the Evolution of Recombination

The nematode Caenorhabditis briggsae is an emerging model organism that allows evolutionary comparisons with C. elegans and exploration of its own unique biological attributes. To produce a high-resolution C. briggsae recombination map, recombinant inbred lines were generated from reciprocal crosses between two strains and genotyped at over 1,000 loci. A second set of recombinant inbred lines involving a third strain was also genotyped at lower resolution. The resulting recombination maps exhibit discrete domains of high and low recombination, as in C. elegans, indicating these are a general feature of Caenorhabditis species. The proportion of a chromosome's physical size occupied by the central, low-recombination domain is highly correlated between species. However, the C. briggsae intra-species comparison reveals striking variation in the distribution of recombination between domains. Hybrid lines made with the more divergent pair of strains also exhibit pervasive marker transmission ratio distortion, evidence of selection acting on hybrid genotypes. The strongest effect, on chromosome III, is explained by a developmental delay phenotype exhibited by some hybrid F2 animals. In addition, on chromosomes IV and V, cross direction-specific biases towards one parental genotype suggest the existence of cytonuclear epistatic interactions. These interactions are discussed in relation to surprising mitochondrial genome polymorphism in C. briggsae, evidence that the two strains diverged in allopatry, the potential for local adaptation, and the evolution of Dobzhansky-Muller incompatibilities. The genetic and genomic resources resulting from this work will support future efforts to understand inter-strain divergence as well as facilitate studies of gene function, natural variation, and the evolution of recombination in Caenorhabditis nematodes

Stochastically Gating Ion Channels Enable Patterned Spike Firing through Activity-Dependent Modulation of Spike Probability

The transformation of synaptic input into patterns of spike output is a fundamental operation that is determined by the particular complement of ion channels that a neuron expresses. Although it is well established that individual ion channel proteins make stochastic transitions between conducting and non-conducting states, most models of synaptic integration are deterministic, and relatively little is known about the functional consequences of interactions between stochastically gating ion channels. Here, we show that a model of stellate neurons from layer II of the medial entorhinal cortex implemented with either stochastic or deterministically gating ion channels can reproduce the resting membrane properties of stellate neurons, but only the stochastic version of the model can fully account for perithreshold membrane potential fluctuations and clustered patterns of spike output that are recorded from stellate neurons during depolarized states. We demonstrate that the stochastic model implements an example of a general mechanism for patterning of neuronal output through activity-dependent changes in the probability of spike firing. Unlike deterministic mechanisms that generate spike patterns through slow changes in the state of model parameters, this general stochastic mechanism does not require retention of information beyond the duration of a single spike and its associated afterhyperpolarization. Instead, clustered patterns of spikes emerge in the stochastic model of stellate neurons as a result of a transient increase in firing probability driven by activation of HCN channels during recovery from the spike afterhyperpolarization. Using this model, we infer conditions in which stochastic ion channel gating may influence firing patterns in vivo and predict consequences of modifications of HCN channel function for in vivo firing patterns

Non-Lytic, Actin-Based Exit of Intracellular Parasites from C. elegans Intestinal Cells

The intestine is a common site for invasion by intracellular pathogens, but little is known about how pathogens restructure and exit intestinal cells in vivo. The natural microsporidian parasite N. parisii invades intestinal cells of the nematode C. elegans, progresses through its life cycle, and then exits cells in a transmissible spore form. Here we show that N. parisii causes rearrangements of host actin inside intestinal cells as part of a novel parasite exit strategy. First, we show that N. parisii infection causes ectopic localization of the normally apical-restricted actin to the basolateral side of intestinal cells, where it often forms network-like structures. Soon after this actin relocalization, we find that gaps appear in the terminal web, a conserved cytoskeletal structure that could present a barrier to exit. Reducing actin expression creates terminal web gaps in the absence of infection, suggesting that infection-induced actin relocalization triggers gap formation. We show that terminal web gaps form at a distinct stage of infection, precisely timed to precede spore exit, and that all contagious animals exhibit gaps. Interestingly, we find that while perturbations in actin can create these gaps, actin is not required for infection progression or spore formation, but actin is required for spore exit. Finally, we show that despite large numbers of spores exiting intestinal cells, this exit does not cause cell lysis. These results provide insight into parasite manipulation of the host cytoskeleton and non-lytic escape from intestinal cells in vivo

Serotonergic chemosensory neurons modify the <i>C. elegans</i> immune response by regulating G-protein signaling in epithelial cells

The nervous and immune systems influence each other, allowing animals to rapidly protect themselves from changes in their internal and external environment. However, the complex nature of these systems in mammals makes it difficult to determine how neuronal signaling influences the immune response. Here we show that serotonin, synthesized in Caenorhabditis elegans chemosensory neurons, modulates the immune response. Serotonin released from these cells acts, directly or indirectly, to regulate G-protein signaling in epithelial cells. Signaling in these cells is required for the immune response to infection by the natural pathogen Microbacterium nematophilum. Here we show that serotonin signaling suppresses the innate immune response and limits the rate of pathogen clearance. We show that C. elegans uses classical neurotransmitters to alter the immune response. Serotonin released from sensory neurons may function to modify the immune system in response to changes in the animal's external environment such as the availability, or quality, of food

Macro-level Modeling of the Response of C. elegans Reproduction to Chronic Heat Stress

A major goal of systems biology is to understand how organism-level behavior arises from a myriad of molecular interactions. Often this involves complex sets of rules describing interactions among a large number of components. As an alternative, we have developed a simple, macro-level model to describe how chronic temperature stress affects reproduction in C. elegans. Our approach uses fundamental engineering principles, together with a limited set of experimentally derived facts, and provides quantitatively accurate predictions of performance under a range of physiologically relevant conditions. We generated detailed time-resolved experimental data to evaluate the ability of our model to describe the dynamics of C. elegans reproduction. We find considerable heterogeneity in responses of individual animals to heat stress, which can be understood as modulation of a few processes and may represent a strategy for coping with the ever-changing environment. Our experimental results and model provide quantitative insight into the breakdown of a robust biological system under stress and suggest, surprisingly, that the behavior of complex biological systems may be determined by a small number of key components