Evaluating statistical methods using plasmode data sets in the age of massive public databases: an illustration using false discovery rates

Plasmode is a term coined several years ago to describe data sets that are derived from real data but for which some truth is known. Omic techniques, most especially microarray and genomewide association studies, have catalyzed a new zeitgeist of data sharing that is making data and data sets publicly available on an unprecedented scale. Coupling such data resources with a science of plasmode use would allow statistical methodologists to vet proposed techniques empirically (as opposed to only theoretically) and with data that are by definition realistic and representative. We illustrate the technique of empirical statistics by consideration of a common task when analyzing high dimensional data: the simultaneous testing of hundreds or thousands of hypotheses to determine which, if any, show statistical significance warranting follow-on research. The now-common practice of multiple testing in high dimensional experiment (HDE) settings has generated new methods for detecting statistically significant results. Although such methods have heretofore been subject to comparative performance analysis using simulated data, simulating data that realistically reflect data from an actual HDE remains a challenge. We describe a simulation procedure using actual data from an HDE where some truth regarding parameters of interest is known. We use the procedure to compare estimates for the proportion of true null hypotheses, the false discovery rate (FDR), and a local version of FDR obtained from 15 different statistical methods.     

TRPV1 channels mediate long-term depression at synapses on hippocampal interneurons

TRPV1 receptors have classically been defined as heat-sensitive, ligand-gated, nonselective cation channels that integrate nociceptive stimuli in sensory neurons. TRPV1 receptors have also been identified in the brain, but their physiological role is poorly understood. Here we report that TRPV1 channel activation is necessary and sufficient to trigger long-term synaptic depression (LTD). Excitatory synapses onto hippocampal interneurons were depressed by either capsaicin, a potent TRPV1 channel activator, or the endogenously released eicosanoid, 12-(S)-HPETE, whereas neighboring excitatory synapses onto CA1 pyramidal cells were unaffected. TRPV1 receptor antagonists also prevented interneuron LTD. In brain slices from TRPV1-/- mice, LTD was absent, and neither capsaicin nor 12-(S)-HPETE elicited synaptic depression. Our results suggest that, in the hippocampus, TRPV1 receptor activation selectively modifies synapses onto interneurons. Like other forms of hippocampal synaptic plasticity, TRPV1-mediated LTD may have a role in long-term changes in physiological and pathological circuit behavior during learning and epileptic activity.     

A radiation hybrid mapping panel for the rhesus macaque

The genomes of nonhuman primates have recently become highly visible candidates for full genome analysis, as they provide powerful models of human disease and a better understanding of the evolution of the human genome. We describe the creation of a 5000 rad radiation hybrid (RH) mapping panel for the rhesus macaque. Duplicate genotypes of 84 microsatellite and coding gene sequence tagged sites from six macaque chromosomes produced an estimated whole genome retention frequency of 0.33. To test the mapping ability of the panel, we constructed RH maps for macaque chromosomes 7 and 9 and compared them to orthologous locus orders in existing human and baboon maps derived from different methodologies. Concordant marker order between all three species maps suggests that the current panel represents a powerful mapping resource for generating high-density comparative maps of the rhesus macaque and other species genomes.     

Peak plasma interleukin-6 and other peripheral markers of inflammation in the first week of ischaemic stroke correlate with brain infarct volume, stroke severity and long-term outcome

Background

Cerebral ischaemia initiates an inflammatory response in the brain and periphery. We assessed the relationship between peak values of plasma interleukin-6 (IL-6) in the first week after ischaemic stroke, with measures of stroke severity and outcome.

Methods

Thirty-seven patients with ischaemic stroke were prospectively recruited. Plasma IL-6, and other markers of peripheral inflammation, were measured at pre-determined timepoints in the first week after stroke onset. Primary analyses were the association between peak plasma IL-6 concentration with both modified Rankin score (mRS) at 3 months and computed tomography (CT) brain infarct volume.

Results

Peak plasma IL-6 concentration correlated significantly (p < 0.001) with CT brain infarct volume (r = 0.75) and mRS at 3 months (r = 0.72). It correlated similarly with clinical outcome at 12 months or stroke severity. Strong associations were also noted between either peak plasma C-reactive protein (CRP) concentration or white blood cell (WBC) count, and all outcome measures.

Conclusions

These data provide evidence that the magnitude of the peripheral inflammatory response is related to the severity of acute ischaemic stroke, and clinical outcome.     

A mode of arthropod brain evolution suggested by Drosophila commissure development

The evolutionary origin of the tritocerebral neuromere, which is a brain segment located at the junction between the supra- and subesophageal ganglia in most mandibulates (arthropods such as crustaceans and insects), is a subject rich in contentious debate. Various models have argued that the tritocerebrum came from a segmental nerve cord ganglia that was recruited into the head during the course of arthropod evolution. However, despite much thought on the subject, the origin of the tritocerebrum remains obscure. Here I describe the development of the tritocerebral commissure in Drosophila and demonstrate that the tritocerebral and mandibular commissures actually form as one commissure and then separate in a manner very similar to how the anterior and posterior commissures of a ventral nerve cord neuromere form. I propose that the tritocerebral neuromere originated from the splitting of an ancestral neuromere located in the anterior subesophageal ganglion into distinct tritocerebral and mandibular neuromeres. Also, I discuss the problem of arthropod brain neuromere homology in reference to this hypothesis.     

Cysteinyl peptide inhibitors of Bacillus cereus zinc beta-lactamase

Several cysteinyl peptides have been synthesised and shown to be reversible competitive inhibitors of the Bacillus cereus metallo-beta-lactamase. The pH dependence of pKi indicates that the thiol anion displaces hydroxide ion from the active site zinc(II). D,D-Peptides bind to the enzyme better than other diastereoisomers, which is compatible with the predicted stereochemistry of the active site.     

A model of primitive streak initiation in the chick embryo

Initiation of the primitive streak in avian embryos provides a well-studied example of a pattern-forming event that displays a striking capacity for regulation. The mechanisms underlying the regulative properties are, however, poorly understood and are not easily accounted for by traditional models of pattern formation, such as reaction-diffusion models. In this paper, we propose a new activator-inhibitor model for streak initiation. We show that the model is consistent with experimental observations, both in its pattern-forming properties and in its ability to form these patterns on the correct time-scales for biologically realistic parameter values. A key component of the model is a travelling wave of inhibition. We present a mathematical analysis of the speed of such waves in both diffusive and juxtacrine relay systems. We use the streak initiation model to make testable predictions. By varying parameters of the model, two very different types of patterning can be obtained, suggesting that our model may be applicable to other processes in addition to streak initiation.     

Mapping of hybrid incompatibility loci in Nasonia

According to theory, F(2) hybrid breakdown (lethality or sterility) is due to incompatibilities between interacting genes of the different species (i.e., the breaking up of coadapted gene complexes). Detection of such incompatibilities is particularly straightforward in haplodiploid species, because virgin F(1) hybrid females will produce haploid recombinant F(2) males. This feature allows for screening of the complete genome for recessive genetic incompatibilities. Crosses were performed between Nasonia vitripennis (v) and its sibling species N. giraulti (g). First, a linkage map was produced using RAPD markers. RAPD markers showed an overall bias toward vitripennis alleles, a pattern not predicted by the basic two-interactor Dobzhansky-Muller model. Recovery patterns of visible markers were consistent with those of linked RAPD markers. If particular genetic interactions between two loci are causing hybrid lethality, then those genotypes should be underrepresented or absent among adult F(2) males. Four sets of significant incompatibilities were detected by performing pairwise comparisons of markers on different chromosomes. Likely explanations for the observed patterns are maternal effect-zygotic gene incompatibilities or clustering of incompatibility loci. Due to the short generation time, advantages of haplodiploidy, and availability of markers, Nasonia promises to be a productive system for investigating the genetics of hybrid inviability.