SimHap GUI: An intuitive graphical user interface for genetic association analysis

<p>Abstract</p> <p>Background</p> <p>Researchers wishing to conduct genetic association analysis involving single nucleotide polymorphisms (SNPs) or haplotypes are often confronted with the lack of user-friendly graphical analysis tools, requiring sophisticated statistical and informatics expertise to perform relatively straightforward tasks. Tools, such as the <it>SimHap </it>package for the R statistics language, provide the necessary statistical operations to conduct sophisticated genetic analysis, but lacks a graphical user interface that allows anyone but a professional statistician to effectively utilise the tool.</p> <p>Results</p> <p>We have developed SimHap GUI, a cross-platform integrated graphical analysis tool for conducting epidemiological, single SNP and haplotype-based association analysis. SimHap GUI features a novel workflow interface that guides the user through each logical step of the analysis process, making it accessible to both novice and advanced users. This tool provides a seamless interface to the <it>SimHap </it>R package, while providing enhanced functionality such as sophisticated data checking, automated data conversion, and real-time estimations of haplotype simulation progress.</p> <p>Conclusion</p> <p>SimHap GUI provides a novel, easy-to-use, cross-platform solution for conducting a range of genetic and non-genetic association analyses. This provides a free alternative to commercial statistics packages that is specifically designed for genetic association analysis.</p

No iron fertilization in the equatorial Pacific Ocean during the last ice age

The equatorial Pacific Ocean is one of the major high-nutrient, low-chlorophyll regions in the global ocean. In such regions, the consumption of the available macro-nutrients such as nitrate and phosphate is thought to be limited in part by the low abundance of the critical micro-nutrient iron1. Greater atmospheric dust deposition2 could have fertilized the equatorial Pacific with iron during the last ice age—the Last Glacial Period (LGP) but the effect of increased ice-age dust fluxes on primary productivity in the equatorial Pacific remains uncertain. Here we present meridional transects of dust (derived from the 232Th proxy), phytoplankton productivity (using opal, 231Pa/230Th and excess Ba), and the degree of nitrate consumption (using foraminifera-bound δ15N) from six cores in the central equatorial Pacific for the Holocene (0–10,000 years ago) and the LGP (17,000–27,000 years ago). We find that, although dust deposition in the central equatorial Pacific was two to three times greater in the LGP than in the Holocene, productivity was the same or lower, and the degree of nitrate consumption was the same. These biogeochemical findings suggest that the relatively greater ice-age dust fluxes were not large enough to provide substantial iron fertilization to the central equatorial Pacific. This may have been because the absolute rate of dust deposition in the LGP (although greater than the Holocene rate) was very low. The lower productivity coupled with unchanged nitrate consumption suggests that the subsurface major nutrient concentrations were lower in the central equatorial Pacific during the LGP. As these nutrients are today dominantly sourced from the Subantarctic Zone of the Southern Ocean, we propose that the central equatorial Pacific data are consistent with more nutrient consumption in the Subantarctic Zone, possibly owing to iron fertilization as a result of higher absolute dust fluxes in this region7,8. Thus, ice-age iron fertilization in the Subantarctic Zone would have ultimately worked to lower, not raise, equatorial Pacific productivity

Is Cortisol Excretion Independent of Menstrual Cycle Day? A Longitudinal Evaluation of First Morning Urinary Specimens

Background Cortisol is frequently used as a marker of physiologic stress levels. Using cortisol for that purpose, however, requires a thorough understanding of its normal longitudinal variability. The current understanding of longitudinal variability of basal cortisol secretion in women is very limited. It is often assumed, for example, that basal cortisol profiles do not vary across the menstrual cycle. This is a critical assumption: if cortisol were to follow a time dependent pattern during the menstrual cycle, then ignoring this cyclic variation could lead to erroneous imputation of physiologic stress. Yet, the assumption that basal cortisol levels are stable across the menstrual cycle rests on partial and contradictory evidence. Here we conduct a thorough test of that assumption using data collected for up to a year from 25 women living in rural Guatemala. Methodology We apply a linear mixed model to describe longitudinal first morning urinary cortisol profiles, accounting for differences in both mean and standard deviation of cortisol among women. To that aim we evaluate the fit of two alternative models. The first model assumes that cortisol does not vary with menstrual cycle day. The second assumes that cortisol mean varies across the menstrual cycle. Menstrual cycles are aligned on ovulation day (day 0). Follicular days are assigned negative numbers and luteal days positive numbers. When we compared Models 1 and 2 restricting our analysis to days between −14 (follicular) and day 14 (luteal) then day of the menstrual cycle did not emerge as a predictor of urinary cortisol levels (p-value &gt;0.05). Yet, when we extended our analyses beyond that central 28-day-period then day of the menstrual cycle become a statistically significant predictor of cortisol levels. Significance The observed trend suggests that studies including cycling women should account for day dependent variation in cortisol in cycles with long follicular and luteal phases

Gradients and Modulation of K+ Channels Optimize Temporal Accuracy in Networks of Auditory Neurons

Accurate timing of action potentials is required for neurons in auditory brainstem nuclei to encode the frequency and phase of incoming sound stimuli. Many such neurons express “high threshold” Kv3-family channels that are required for firing at high rates (>∼200 Hz). Kv3 channels are expressed in gradients along the medial-lateral tonotopic axis of the nuclei. Numerical simulations of auditory brainstem neurons were used to calculate the input-output relations of ensembles of 1–50 neurons, stimulated at rates between 100–1500 Hz. Individual neurons with different levels of potassium currents differ in their ability to follow specific rates of stimulation but all perform poorly when the stimulus rate is greater than the maximal firing rate of the neurons. The temporal accuracy of the combined synaptic output of an ensemble is, however, enhanced by the presence of gradients in Kv3 channel levels over that measured when neurons express uniform levels of channels. Surprisingly, at high rates of stimulation, temporal accuracy is also enhanced by the occurrence of random spontaneous activity, such as is normally observed in the absence of sound stimulation. For any pattern of stimulation, however, greatest accuracy is observed when, in the presence of spontaneous activity, the levels of potassium conductance in all of the neurons is adjusted to that found in the subset of neurons that respond better than their neighbors. This optimization of response by adjusting the K+ conductance occurs for stimulus patterns containing either single and or multiple frequencies in the phase-locking range. The findings suggest that gradients of channel expression are required for normal auditory processing and that changes in levels of potassium currents across the nuclei, by mechanisms such as protein phosphorylation and rapid changes in channel synthesis, adapt the nuclei to the ongoing auditory environment

Dengue Virus Infection Perturbs Lipid Homeostasis in Infected Mosquito Cells

Dengue virus causes ∼50–100 million infections per year and thus is considered one of the most aggressive arthropod-borne human pathogen worldwide. During its replication, dengue virus induces dramatic alterations in the intracellular membranes of infected cells. This phenomenon is observed both in human and vector-derived cells. Using high-resolution mass spectrometry of mosquito cells, we show that this membrane remodeling is directly linked to a unique lipid repertoire induced by dengue virus infection. Specifically, 15% of the metabolites detected were significantly different between DENV infected and uninfected cells while 85% of the metabolites detected were significantly different in isolated replication complex membranes. Furthermore, we demonstrate that intracellular lipid redistribution induced by the inhibition of fatty acid synthase, the rate-limiting enzyme in lipid biosynthesis, is sufficient for cell survival but is inhibitory to dengue virus replication. Lipids that have the capacity to destabilize and change the curvature of membranes as well as lipids that change the permeability of membranes are enriched in dengue virus infected cells. Several sphingolipids and other bioactive signaling molecules that are involved in controlling membrane fusion, fission, and trafficking as well as molecules that influence cytoskeletal reorganization are also up regulated during dengue infection. These observations shed light on the emerging role of lipids in shaping the membrane and protein environments during viral infections and suggest membrane-organizing principles that may influence virus-induced intracellular membrane architecture

Impact of Dendritic Size and Dendritic Topology on Burst Firing in Pyramidal Cells

Neurons display a wide range of intrinsic firing patterns. A particularly relevant pattern for neuronal signaling and synaptic plasticity is burst firing, the generation of clusters of action potentials with short interspike intervals. Besides ion-channel composition, dendritic morphology appears to be an important factor modulating firing pattern. However, the underlying mechanisms are poorly understood, and the impact of morphology on burst firing remains insufficiently known. Dendritic morphology is not fixed but can undergo significant changes in many pathological conditions. Using computational models of neocortical pyramidal cells, we here show that not only the total length of the apical dendrite but also the topological structure of its branching pattern markedly influences inter- and intraburst spike intervals and even determines whether or not a cell exhibits burst firing. We found that there is only a range of dendritic sizes that supports burst firing, and that this range is modulated by dendritic topology. Either reducing or enlarging the dendritic tree, or merely modifying its topological structure without changing total dendritic length, can transform a cell's firing pattern from bursting to tonic firing. Interestingly, the results are largely independent of whether the cells are stimulated by current injection at the soma or by synapses distributed over the dendritic tree. By means of a novel measure called mean electrotonic path length, we show that the influence of dendritic morphology on burst firing is attributable to the effect both dendritic size and dendritic topology have, not on somatic input conductance, but on the average spatial extent of the dendritic tree and the spatiotemporal dynamics of the dendritic membrane potential. Our results suggest that alterations in size or topology of pyramidal cell morphology, such as observed in Alzheimer's disease, mental retardation, epilepsy, and chronic stress, could change neuronal burst firing and thus ultimately affect information processing and cognition

Measurement of the gamma ray background in the Davis cavern at the Sanford Underground Research Facility

Deep underground environments are ideal for low background searches due to the attenuation of cosmic rays by passage through the earth. However, they are affected by backgrounds from γ-rays emitted by 40K and the 238U and 232Th decay chains in the surrounding rock. The LUX-ZEPLIN (LZ) experiment will search for dark matter particle interactions with a liquid xenon TPC located within the Davis campus at the Sanford Underground Research Facility, Lead, South Dakota, at the 4850-foot level. In order to characterise the cavern background, in-situ γ-ray measurements were taken with a sodium iodide detector in various locations and with lead shielding. The integral count rates (0–3300 keV) varied from 596 Hz to 1355 Hz for unshielded measurements, corresponding to a total flux from the cavern walls of 1.9 ± 0.4 γ cm−2s−1. The resulting activity in the walls of the cavern can be characterised as 220 ± 60 Bq/kg of 40K, 29 ± 15 Bq/kg of 238U, and 13 ± 3 Bq/kg of 232Th

Biophysical Basis for Three Distinct Dynamical Mechanisms of Action Potential Initiation

Transduction of graded synaptic input into trains of all-or-none action potentials (spikes) is a crucial step in neural coding. Hodgkin identified three classes of neurons with qualitatively different analog-to-digital transduction properties. Despite widespread use of this classification scheme, a generalizable explanation of its biophysical basis has not been described. We recorded from spinal sensory neurons representing each class and reproduced their transduction properties in a minimal model. With phase plane and bifurcation analysis, each class of excitability was shown to derive from distinct spike initiating dynamics. Excitability could be converted between all three classes by varying single parameters; moreover, several parameters, when varied one at a time, had functionally equivalent effects on excitability. From this, we conclude that the spike-initiating dynamics associated with each of Hodgkin's classes represent different outcomes in a nonlinear competition between oppositely directed, kinetically mismatched currents. Class 1 excitability occurs through a saddle node on invariant circle bifurcation when net current at perithreshold potentials is inward (depolarizing) at steady state. Class 2 excitability occurs through a Hopf bifurcation when, despite net current being outward (hyperpolarizing) at steady state, spike initiation occurs because inward current activates faster than outward current. Class 3 excitability occurs through a quasi-separatrix crossing when fast-activating inward current overpowers slow-activating outward current during a stimulus transient, although slow-activating outward current dominates during constant stimulation. Experiments confirmed that different classes of spinal lamina I neurons express the subthreshold currents predicted by our simulations and, further, that those currents are necessary for the excitability in each cell class. Thus, our results demonstrate that all three classes of excitability arise from a continuum in the direction and magnitude of subthreshold currents. Through detailed analysis of the spike-initiating process, we have explained a fundamental link between biophysical properties and qualitative differences in how neurons encode sensory input