Comparison of burst detectors for spike trains

Accurate identification of bursting activity is an essential element in the characterization of neuronal network activity. Despite this, no one technique for identifying bursts in spike trains has been widely adopted. Instead, many methods have been developed for the analysis of bursting activity, often on an ad hoc basis. Here we provide an unbiased assessment of the effectiveness of eight of these methods at detecting bursts in a range of spike trains. We suggest a list of features that an ideal burst detection technique should possess and use synthetic data to assess each method in regard to these properties. We further employ each of the methods to reanalyze microelectrode array (MEA) recordings from mouse retinal ganglion cells and examine their coherence with bursts detected by a human observer. We show that several common burst detection techniques perform poorly at analyzing spike trains with a variety of properties. We identify four promising burst detection techniques, which are then applied to MEA recordings of networks of human induced pluripotent stem cell-derived neurons and used to describe the ontogeny of bursting activity in these networks over several months of development. We conclude that no current method can provide "perfect" burst detection results across a range of spike trains; however, two burst detection techniques, the MaxInterval and logISI methods, outperform compared with others. We provide recommendations for the robust analysis of bursting activity in experimental recordings using current techniques.Experimental data collection was supported by the BBSRC (PC, OP, grant number BB/H008608/1). EC was supported by a Wellcome Trust PhD Studentship and NIHR Cambridge Biomedical Research Centre Studentship. CWT was supported by a bursary from the Bridgwater Summer Undergraduate Research programme.This is the final version of the article. It first appeared from the American Physiological Society via https://doi.org/10.1152/jn.00093.201

The purification and characterization of ATP synthase complexes from the mitochondria of four fungal species.

The ATP synthases have been isolated by affinity chromatography from the mitochondria of the fungal species Yarrowia lipolytica, Pichia pastoris, Pichia angusta and Saccharomyces cerevisiae. The subunit compositions of the purified enzyme complexes depended on the detergent used to solubilize and purify the complex, and the presence or absence of exogenous phospholipids. All four enzymes purified in the presence of n-dodecyl-β-D-maltoside had a complete complement of core subunits involved directly in the synthesis of ATP, but they were deficient to different extents in their supernumerary membrane subunits. In contrast, the enzymes from P. angusta and S. cerevisiae purified in the presence of n-decyl-β-maltose neopentyl glycol and the phospholipids 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine, cardiolipin (diphosphatidylglycerol) and 1-palmitoyl-2-oleoyl-sn-glycero-3-[phospho-rac-(1-glycerol)] had a complete complement of core subunits and also contained all of the known supernumerary membrane subunits, e, f, g, j, k and ATP8 (or Aap1), plus an additional new membrane component named subunit l, related in sequence to subunit k. The catalytic domain of the enzyme from P. angusta was more resistant to thermal denaturation than the enzyme from S. cerevisiae, but less stable than the catalytic domain of the bovine enzyme, but the stator and the integrity of the transmembrane proton pathway were most stable in the enzyme from P. angusta. The P. angusta enzyme provides a suitable source of enzyme for studying the structure of the membrane domain and properties associated with that sector of the enzyme complex.This work was funded by the intramural programme of the Medical Research Council (MRC). T.J.C was the recipient of an MRC graduate studentship

Wound healing and hyper-hydration - a counter intuitive model

Winters seminal work in the 1960s relating to providing an optimal level of moisture to aid wound healing (granulation and re-epithelialisation) has been the single most effective advance in wound care over many decades. As such the development of advanced wound dressings that manage the fluidic wound environment have provided significant benefits in terms of healing to both patient and clinician. Although moist wound healing provides the guiding management principle confusion may arise between what is deemed to be an adequate level of tissue hydration and the risk of developing maceration. In addition, the counter-intuitive model ‘hyper-hydration’ of tissue appears to frustrate the moist wound healing approach and advocate a course of intervention whereby tissue is hydrated beyond what is a normally acceptable therapeutic level. This paper discusses tissue hydration, the cause and effect of maceration and distinguishes these from hyper-hydration of tissue. The rationale is to provide the clinician with a knowledge base that allows optimisation of treatment and outcomes and explains the reasoning behind wound healing using hyper-hydration

Genetic signal maximization using environmental regression

Joint analyses of correlated phenotypes in genetic epidemiology studies are common. However, these analyses primarily focus on genetic correlation between traits and do not take into account environmental correlation. We describe a method that optimizes the genetic signal by accounting for stochastic environmental noise through joint analysis of a discrete trait and a correlated quantitative marker. We conducted bivariate analyses where heritability and the environmental correlation between the discrete and quantitative traits were calculated using Genetic Analysis Workshop 17 (GAW17) family data. The resulting inverse value of the environmental correlation between these traits was then used to determine a new β coefficient for each quantitative trait and was constrained in a univariate model. We conducted genetic association tests on 7,087 nonsynonymous SNPs in three GAW17 family replicates for Affected status with the β coefficient fixed for three quantitative phenotypes and compared these to an association model where the β coefficient was allowed to vary. Bivariate environmental correlations were 0.64 (± 0.09) for Q1, 0.798 (± 0.076) for Q2, and −0.169 (± 0.18) for Q4. Heritability of Affected status improved in each univariate model where a constrained β coefficient was used to account for stochastic environmental effects. No genome-wide significant associations were identified for either method but we demonstrated that constraining β for covariates slightly improved the genetic signal for Affected status. This environmental regression approach allows for increased heritability when the β coefficient for a highly correlated quantitative covariate is constrained and increases the genetic signal for the discrete trait

Strong signature of natural selection within an FHIT intron implicated in prostate cancer risk

Previously, a candidate gene linkage approach on brother pairs affected with prostate cancer identified a locus of prostate cancer susceptibility at D3S1234 within the fragile histidine triad gene (FHIT), a tumor suppressor that induces apoptosis. Subsequent association tests on 16 SNPs spanning approximately 381 kb surrounding D3S1234 in Americans of European descent revealed significant evidence of association for a single SNP within intron 5 of FHIT. In the current study, resequencing and genotyping within a 28.5 kb region surrounding this SNP further delineated the association with prostate cancer risk to a 15 kb region. Multiple SNPs in sequences under evolutionary constraint within intron 5 of FHIT defined several related haplotypes with an increased risk of prostate cancer in European-Americans. Strong associations were detected for a risk haplotype defined by SNPs 138543, 142413, and 152494 in all cases (Pearson's χ2 = 12.34, df 1, P = 0.00045) and for the homozygous risk haplotype defined by SNPs 144716, 142413, and 148444 in cases that shared 2 alleles identical by descent with their affected brothers (Pearson's χ2 = 11.50, df 1, P = 0.00070). In addition to highly conserved sequences encompassing SNPs 148444 and 152413, population studies revealed strong signatures of natural selection for a 1 kb window covering the SNP 144716 in two human populations, the European American (π = 0.0072, Tajima's D= 3.31, 14 SNPs) and the Japanese (π = 0.0049, Fay & Wu's H = 8.05, 14 SNPs), as well as in chimpanzees (Fay & Wu's H = 8.62, 12 SNPs). These results strongly support the involvement of the FHIT intronic region in an increased risk of prostate cancer. © 2008 Ding et al

Bacterial microevolution and the Pangenome

The comparison of multiple genome sequences sampled from a bacterial population reveals considerable diversity in both the core and the accessory parts of the pangenome. This diversity can be analysed in terms of microevolutionary events that took place since the genomes shared a common ancestor, especially deletion, duplication, and recombination. We review the basic modelling ingredients used implicitly or explicitly when performing such a pangenome analysis. In particular, we describe a basic neutral phylogenetic framework of bacterial pangenome microevolution, which is not incompatible with evaluating the role of natural selection. We survey the different ways in which pangenome data is summarised in order to be included in microevolutionary models, as well as the main methodological approaches that have been proposed to reconstruct pangenome microevolutionary history

Cytogerontology since 1881: A reappraisal of August Weismann and a review of modern progress

Cytogerontology, the science of cellular ageing, originated in 1881 with the prediction by August Weismann that the somatic cells of higher animals have limited division potential. Weismann's prediction was derived by considering the role of natural selection in regulating the duration of an organism's life. For various reasons, Weismann's ideas on ageing fell into neglect following his death in 1914, and cytogerontology has only reappeared as a major research area following the demonstration by Hayflick and Moorhead in the early 1960s that diploid human fibroblasts are restricted to a finite number of divisions in vitro. In this review we give a detailed account of Weismann's theory, and we reveal that his ideas were both more extensive in their scope and more pertinent to current research than is generally recognised. We also appraise the progress which has been made over the past hundred years in investigating the causes of ageing, with particular emphasis being given to (i) the evolution of ageing, and (ii) ageing at the cellular level. We critically assess the current state of knowledge in these areas and recommend a series of points as primary targets for future research

Improved Heterosis Prediction by Combining Information on DNA- and Metabolic Markers

Background: Hybrids represent a cornerstone in the success story of breeding programs. The fundamental principle underlying this success is the phenomenon of hybrid vigour, or heterosis. It describes an advantage of the offspring as compared to the two parental lines with respect to parameters such as growth and resistance against abiotic or biotic stress. Dominance, overdominance or epistasis based models are commonly used explanations. Conclusion/Significance: The heterosis level is clearly a function of the combination of the parents used for offspring production. This results in a major challenge for plant breeders, as usually several thousand combinations of parents have to be tested for identifying the best combinations. Thus, any approach to reliably predict heterosis levels based on properties of the parental lines would be highly beneficial for plant breeding. Methodology/Principal Findings: Recently, genetic data have been used to predict heterosis. Here we show that a combination of parental genetic and metabolic markers, identified via feature selection and minimum-description-length based regression methods, significantly improves the prediction of biomass heterosis in resulting offspring. These findings will help furthering our understanding of the molecular basis of heterosis, revealing, for instance, the presence of nonlinear genotype-phenotype relationships. In addition, we describe a possible approach for accelerated selection in plant breeding

The Pattern of R2 Retrotransposon Activity in Natural Populations of Drosophila simulans Reflects the Dynamic Nature of the rDNA Locus

The pattern and frequency of insertions that enable transposable elements to remain active in a population are poorly understood. The retrotransposable element R2 exclusively inserts into the 28S rRNA genes where it establishes long-term, stable relationships with its animal hosts. Previous studies with laboratory stocks of Drosophila simulans have suggested that control over R2 retrotransposition resides within the rDNA loci. In this report, we sampled 180 rDNA loci of animals collected from two natural populations of D. simulans. The two populations were found to have similar patterns of R2 activity. About half of the rDNA loci supported no or very low levels of R2 transcripts with no evidence of R2 retrotransposition. The remaining half of the rDNA loci had levels of R2 transcripts that varied in a continuous manner over almost a 100-fold range and did support new retrotransposition events. Structural analysis of the rDNA loci in 18 lines that spanned the range of R2 transcript levels in these populations revealed that R2 number and rDNA locus size varied 2-fold; however, R2 activity was not readily correlated with either of these parameters. Instead R2 activity was best correlated with the distribution of elements within the rDNA locus. Loci with no activity had larger contiguous blocks of rDNA units free of R2-insertions. These data suggest a model in which frequent recombination within the rDNA locus continually redistributes R2-inserted units resulting in changing levels of R2 activity within individual loci and persistent R2 activity within the population

Genetic architecture and evolution of the S locus supergene in Primula vulgaris

Darwin’s studies on heterostyly in Primula described two floral morphs, pin and thrum, with reciprocal anther and stigma heights that promote insect-mediated cross-pollination. This key innovation evolved independently in several angiosperm families. Subsequent studies on heterostyly in Primula contributed to the foundation of modern genetic theory and the neo-Darwinian synthesis. The established genetic model for Primula heterostyly involves a diallelic S locus comprising several genes, with rare recombination events that result in self-fertile homostyle flowers with anthers and stigma at the same height. Here we reveal the S locus supergene as a tightly-linked cluster of thrum-specific genes that are absent in pins. We show that thrums are hemizygous not heterozygous for the S locus, which suggests that homostyles do not arise by recombination between S locus haplotypes as previously proposed. Duplication of a floral homeotic gene 51.7 MYA, followed by its neofunctionalisation, created the current S locus assemblage which led to floral heteromorphy in Primula. Our findings provide new insights into the structure, function and evolution of this archetypal supergene