Improved prediction of hiking speeds using a data driven approach.

Hikers and hillwalkers typically use the gradient in the direction of travel (walking slope) as the main variable in established methods for predicting walking time (via the walking speed) along a route. Research into fell-running has suggested further variables which may improve speed algorithms in this context; the gradient of the terrain (hill slope) and the level of terrain obstruction. Recent improvements in data availability, as well as widespread use of GPS tracking now make it possible to explore these variables in a walking speed model at a sufficient scale to test statistical significance. We tested various established models used to predict walking speed against public GPS data from almost 88,000 km of UK walking / hiking tracks. Tracks were filtered to remove breaks and non-walking sections. A new generalised linear model (GLM) was then used to predict walking speeds. Key differences between the GLM and established rules were that the GLM considered the gradient of the terrain (hill slope) irrespective of walking slope, as well as the terrain type and level of terrain obstruction in off-road travel. All of these factors were shown to be highly significant, and this is supported by a lower root-mean-square-error compared to existing functions. We also observed an increase in RMSE between the GLM and established methods as hill slope increases, further supporting the importance of this variable.

A phospho-dawn of protein modification anticipates light onset in the picoeukaryote Ostreococcus tauri.

Diel regulation of protein levels and protein modification had been less studied than transcript rhythms. Here, we compare transcriptome data under light-dark cycles with partial proteome and phosphoproteome data, assayed using shotgun MS, from the alga Ostreococcus tauri, the smallest free-living eukaryote. A total of 10% of quantified proteins but two-thirds of phosphoproteins were rhythmic. Mathematical modelling showed that light-stimulated protein synthesis can account for the observed clustering of protein peaks in the daytime. Prompted by night-peaking and apparently dark-stable proteins, we also tested cultures under prolonged darkness, where the proteome changed less than under the diel cycle. Among the dark-stable proteins were prasinophyte-specific sequences that were also reported to accumulate when O. tauri formed lipid droplets. In the phosphoproteome, 39% of rhythmic phospho-sites reached peak levels just before dawn. This anticipatory phosphorylation suggests that a clock-regulated phospho-dawn prepares green cells for daytime functions. Acid-directed and proline-directed protein phosphorylation sites were regulated in antiphase, implicating the clock-related casein kinases 1 and 2 in phase-specific regulation, alternating with the CMGC protein kinase family. Understanding the dynamic phosphoprotein network should be facilitated by the minimal kinome and proteome of O. tauri. The data are available from ProteomeXchange, with identifiers PXD001734, PXD001735, and PXD002909.

Comparison of rule- and ordinary differential equation-based dynamic model of DARPP-32 signalling network.

Dynamic modelling has considerably improved our understanding of complex molecular mechanisms. Ordinary differential equations (ODEs) are the most detailed and popular approach to modelling the dynamics of molecular systems. However, their application in signalling networks, characterised by multi-state molecular complexes, can be prohibitive. Contemporary modelling methods, such as rule- based (RB) modelling, have addressed these issues. The advantages of RB modelling over ODEs have been presented and discussed in numerous reviews. In this study, we conduct a direct comparison of the time courses of a molecular system founded on the same reaction network but encoded in the two frameworks. To make such a comparison, a set of reactions that underlie an ODE model was manually encoded in the Kappa language, one of the RB implementations. A comparison of the models was performed at the level of model specification and dynamics, acquired through model simulations. In line with previous reports, we confirm that the Kappa model recapitulates the general dynamics of its ODE counterpart with minor differences. These occur when molecules have multiple sites binding the same interactor. Furthermore, activation of these molecules in the RB model is slower than in the ODE one. As reported for other molecular systems, we find that, also for the DARPP-32 reaction network, the RB representation offers a more expressive and flexible syntax that facilitates access to fine details of the model, easing model reuse. In parallel with these analyses, we report a refactored model of the DARPP-32 interaction network that can serve as a canvas for the development of more complex dynamic models to study this important molecular system.

SFARI genes and where to find them; modelling Autism Spectrum Disorder specific gene expression dysregulation with RNA-seq data.

Autism Spectrum Disorders (ASD) have a strong, yet heterogeneous, genetic component. Among the various methods that are being developed to help reveal the underlying molecular aetiology of the disease one approach that is gaining popularity is the combination of gene expression and clinical genetic data, often using the SFARI-gene database, which comprises lists of curated genes considered to have causative roles in ASD when mutated in patients. We build a gene co-expression network to study the relationship between ASD-specific transcriptomic data and SFARI genes and then analyse it at different levels of granularity. No significant evidence is found of association between SFARI genes and differential gene expression patterns when comparing ASD samples to a control group, nor statistical enrichment of SFARI genes in gene co-expression network modules that have a strong correlation with ASD diagnosis. However, classification models that incorporate topological information from the whole ASD-specific gene co-expression network can predict novel SFARI candidate genes that share features of existing SFARI genes and have support for roles in ASD in the literature. A statistically significant association is also found between the absolute level of gene expression and SFARI's genes and Scores, which can confound the analysis if uncorrected. We propose a novel approach to correct for this that is general enough to be applied to other problems affected by continuous sources of bias. It was found that only co-expression network analyses that integrate information from the whole network are able to reveal signatures linked to ASD diagnosis and novel candidate genes for the study of ASD, which individual gene or module analyses fail to do. It was also found that the influence of SFARI genes permeates not only other ASD scoring systems, but also lists of genes believed to be involved in other neurodevelopmental disorders.

Presence of DNA from Chlamydia-like organisms in the nasal cavities of grey seal pups (Halichoerus grypus) and three different substrates present in a breeding colony.

BACKGROUND: Chlamydia-like organisms (CLO) have been found to be present in many environmental niches, including human sewage and agricultural run-off, as well as in a number of aquatic species worldwide. Therefore, monitoring their presence in sentinel wildlife species may be useful in assessing the wider health of marine food webs in response to habitat loss, pollution and disease. We used nasal swabs from live (n = 42) and dead (n = 50) pre-weaned grey seal pups and samples of differing natal substrates (n = 8) from an off-shore island devoid of livestock and permanent human habitation to determine if CLO DNA is present in these mammals and to identify possible sources. RESULTS: We recovered CLO DNA from 32/92 (34.7%) nasal swabs from both live (n = 17) and dead (n = 15) seal pups that clustered most closely with currently recognised species belonging to three chlamydial families: Parachlamydiaceae (n = 22), Rhabdochlamydiaceae (n = 6), and Simkaniaceae (n = 3). All DNA positive sediment samples (n = 7) clustered with the Rhabdochlamydiaceae. No difference was found in rates of recovery of CLO DNA in live versus dead pups suggesting the organisms are commensal but their potential as opportunistic secondary pathogens could not be determined. CONCLUSION: This is the first report of CLO DNA being found in marine mammals. This identification warrants further investigation in other seal populations around the coast of the UK and in other areas of the world to determine if this finding is unique or more common than shown by this data. Further investigation would also be warranted to determine if they are present as purely commensal organisms or whether they could also be opportunistic pathogens in seals, as well as to investigate possible sources of origin, including whether they originated as a result of anthropogenic impacts, including human waste and agricultural run-off.

A unified resource and configurable model of the synapse proteome and its role in disease.

Genes encoding synaptic proteins are highly associated with neuronal disorders many of which show clinical co-morbidity. We integrated 58 published synaptic proteomic datasets that describe over 8000 proteins and combined them with direct protein-protein interactions and functional metadata to build a network resource that reveals the shared and unique protein components that underpin multiple disorders. All the data are provided in a flexible and accessible format to encourage custom use.

MicroRNA-34a Acutely Regulates Synaptic Efficacy in the Adult Dentate Gyrus In Vivo.

Activity-dependent synaptic plasticity involves rapid regulation of neuronal protein synthesis on a time-scale of minutes. miRNA function in synaptic plasticity and memory formation has been elucidated by stable experimental manipulation of miRNA expression and activity using transgenic approaches and viral vectors. However, the impact of rapid miRNA modulation on synaptic efficacy is unknown. Here, we examined the effect of acute (12 min), intrahippocampal infusion of a miR-34a antagonist (antimiR) on medial perforant path-evoked synaptic transmission in the dentate gyrus of adult anesthetised rats. AntimiR-34a infusion acutely depressed medial perforant path-evoked field excitatory post-synaptic potentials (fEPSPs). The fEPSP decrease was detected within 9 min of infusion, lasted for hours, and was associated with knockdown of antimiR-34a levels. AntimiR-34a-induced synaptic depression was sequence-specific; no changes were elicited by infusion of scrambled or mismatch control. The rapid modulation suggests that a target, or set of targets, is regulated by miR-34a. Western blot analysis of dentate gyrus lysates revealed enhanced expression of Arc, a known miR-34a target, and four novel predicted targets (Ctip2, PKI-1α, TCF4 and Ube2g1). Remarkably, antimiR-34a had no effect when infused during the maintenance phase of long-term potentiation. We conclude that miR-34a regulates basal synaptic efficacy in the adult dentate gyrus in vivo. To our knowledge, these in vivo findings are the first to demonstrate acute (< 9 min) regulation of synaptic efficacy in the adult brain by a miRNA.

Mixed-species RNA-seq for elucidation of non-cell-autonomous control of gene transcription.

Transcriptomic changes induced in one cell type by another mediate many biological processes in the brain and elsewhere; however, achieving artifact-free physical separation of cell types to study them is challenging and generally allows for analysis of only a single cell type. We describe an approach using a co-culture of distinct cell types from different species that enables physical cell sorting to be replaced by in silico RNA sequencing (RNA-seq) read sorting, which is possible because of evolutionary divergence of messenger RNA (mRNA) sequences. As an exemplary experiment, we describe the co-culture of purified neurons, astrocytes, and microglia from different species (12-14 d). We describe how to use our Python tool, Sargasso, to separate the reads from conventional RNA-seq according to species and to eliminate any artifacts borne of imperfect genome annotation (10 h). We show how this procedure, which requires no special skills beyond those that might normally be expected of wet lab and bioinformatics researchers, enables the simultaneous transcriptomic profiling of different cell types, revealing the distinct influence of microglia on astrocytic and neuronal transcriptomes under inflammatory conditions.

Author Correction: Neurons and neuronal activity control gene expression in astrocytes to regulate their development and metabolism.

This corrects the article DOI: 10.1038/ncomms15132.

Neurons and neuronal activity control gene expression in astrocytes to regulate their development and metabolism.

The influence that neurons exert on astrocytic function is poorly understood. To investigate this, we first developed a system combining cortical neurons and astrocytes from closely related species, followed by RNA-seq and in silico species separation. This approach uncovers a wide programme of neuron-induced astrocytic gene expression, involving Notch signalling, which drives and maintains astrocytic maturity and neurotransmitter uptake function, is conserved in human development, and is disrupted by neurodegeneration. Separately, hundreds of astrocytic genes are acutely regulated by synaptic activity via mechanisms involving cAMP/PKA-dependent CREB activation. This includes the coordinated activity-dependent upregulation of major astrocytic components of the astrocyte-neuron lactate shuttle, leading to a CREB-dependent increase in astrocytic glucose metabolism and elevated lactate export. Moreover, the groups of astrocytic genes induced by neurons or neuronal activity both show age-dependent decline in humans. Thus, neurons and neuronal activity regulate the astrocytic transcriptome with the potential to shape astrocyte-neuron metabolic cooperation.

Evidence for evolutionary divergence of activity-dependent gene expression in developing neurons.

Evolutionary differences in gene regulation between humans and lower mammalian experimental systems are incompletely understood, a potential translational obstacle that is challenging to surmount in neurons, where primary tissue availability is poor. Rodent-based studies show that activity-dependent transcriptional programs mediate myriad functions in neuronal development, but the extent of their conservation in human neurons is unknown. We compared activity-dependent transcriptional responses in developing human stem cell-derived cortical neurons with those induced in developing primary- or stem cell-derived mouse cortical neurons. While activity-dependent gene-responsiveness showed little dependence on developmental stage or origin (primary tissue vs. stem cell), notable species-dependent differences were observed. Moreover, differential species-specific gene ortholog regulation was recapitulated in aneuploid mouse neurons carrying human chromosome-21, implicating promoter/enhancer sequence divergence as a factor, including human-specific activity-responsive AP-1 sites. These findings support the use of human neuronal systems for probing transcriptional responses to physiological stimuli or indeed pharmaceutical agents.

Analysis of compound heterozygotes reveals that the mouse floxed Pax6 (tm1Ued) allele produces abnormal eye phenotypes.

Analysis of abnormal phenotypes produced by different types of mutations has been crucial for our understanding of gene function. Some floxed alleles that retain a neomycin-resistance selection cassette (neo cassette) are not equivalent to wild-type alleles and provide useful experimental resources. Pax6 is an important developmental gene and the aim of this study was to determine whether the floxed Pax6 (tm1Ued) (Pax6 (fl) ) allele, which has a retained neo cassette, produced any abnormal eye phenotypes that would imply that it differs from the wild-type allele. Homozygous Pax6 (fl/fl) and heterozygous Pax6 (fl/+) mice had no overt qualitative eye abnormalities but morphometric analysis showed that Pax6 (fl/fl) corneas tended be thicker and smaller in diameter. To aid identification of weak effects, we produced compound heterozygotes with the Pax6 (Sey-Neu) (Pax6 (-)) null allele. Pax6 (fl/-) compound heterozygotes had more severe eye abnormalities than Pax6 (+/-) heterozygotes, implying that Pax6 (fl) differs from the wild-type Pax6 (+) allele. Immunohistochemistry showed that the Pax6 (fl/-) corneal epithelium was positive for keratin 19 and negative for keratin 12, indicating that it was abnormally differentiated. This Pax6 (fl) allele provides a useful addition to the existing Pax6 allelic series and this study demonstrates the utility of using compound heterozygotes with null alleles to unmask cryptic effects of floxed alleles.

Fluid consumption and taste novelty determines transcription temporal dynamics in the gustatory cortex.

BACKGROUND: Novel taste memories, critical for animal survival, are consolidated to form long term memories which are dependent on translation regulation in the gustatory cortex (GC) hours following acquisition. However, the role of transcription regulation in the process is unknown. RESULTS: Here, we report that transcription in the GC is necessary for taste learning in rats, and that drinking and its consequences, as well as the novel taste experience, affect transcription in the GC during taste memory consolidation. We show differential effects of learning on temporal dynamics in set of genes in the GC, including Arc/Arg3.1, known to regulate the homeostasis of excitatory synapses. CONCLUSIONS: We demonstrate that in taste learning, transcription programs were activated following the physiological responses (i.e., fluid consumption following a water restriction regime, reward, arousal of the animal, etc.) and the specific information about a given taste (i.e., taste novelty). Moreover, the cortical differential prolonged kinetics of mRNA following novel versus familiar taste learning may represent additional novelty related molecular response, where not only the total amount, but also the temporal dynamics of transcription is modulated by sensory experience of novel information.

Expression of mRNA Encoding Mcu and Other Mitochondrial Calcium Regulatory Genes Depends on Cell Type, Neuronal Subtype, and Ca2+ Signaling.

Uptake of Ca2+ into the mitochondrial matrix controls cellular metabolism and survival-death pathways. Several genes are implicated in controlling mitochondrial Ca2+ uptake (mitochondrial calcium regulatory genes, MCRGs), however, less is known about the factors which influence their expression level. Here we have compared MCRG mRNA expression, in neural cells of differing type (cortical neurons vs. astrocytes), differing neuronal subtype (CA3 vs. CA1 hippocampus) and in response to Ca2+ influx, using a combination of qPCR and RNA-seq analysis. Of note, we find that the Mcu-regulating Micu gene family profile differs substantially between neurons and astrocytes, while expression of Mcu itself is markedly different between CA3 and CA1 regions in the adult hippocampus. Moreover, dynamic control of MCRG mRNA expression in response to membrane depolarization-induced Ca2+ influx is also apparent, resulting in repression of Letm1, as well as Mcu. Thus, the mRNA expression profile of MCRGs is not fixed, which may cause differences in the coupling between cytoplasmic and mitochondrial Ca2+, as well as diversity of mitochondrial Ca2+ uptake mechanisms.

Network analyses reveal novel aspects of ALS pathogenesis.

Amyotrophic Lateral Sclerosis (ALS) is a fatal neurodegenerative disease characterized by selective loss of motor neurons, muscle atrophy and paralysis. Mutations in the human VAMP-associated protein B (hVAPB) cause a heterogeneous group of motor neuron diseases including ALS8. Despite extensive research, the molecular mechanisms underlying ALS pathogenesis remain largely unknown. Genetic screens for key interactors of hVAPB activity in the intact nervous system, however, represent a fundamental approach towards understanding the in vivo function of hVAPB and its role in ALS pathogenesis. Targeted expression of the disease-causing allele leads to neurodegeneration and progressive decline in motor performance when expressed in the adult Drosophila, eye or in its entire nervous system, respectively. By using these two phenotypic readouts, we carried out a systematic survey of the Drosophila genome to identify modifiers of hVAPB-induced neurotoxicity. Modifiers cluster in a diverse array of biological functions including processes and genes that have been previously linked to hVAPB function, such as proteolysis and vesicular trafficking. In addition to established mechanisms, the screen identified endocytic trafficking and genes controlling proliferation and apoptosis as potent modifiers of ALS8-mediated defects. Surprisingly, the list of modifiers was mostly enriched for proteins linked to lipid droplet biogenesis and dynamics. Computational analysis reveals that most modifiers can be linked into a complex network of interacting genes, and that the human genes homologous to the Drosophila modifiers can be assembled into an interacting network largely overlapping with that in flies. Identity markers of the endocytic process were also found to abnormally accumulate in ALS patients, further supporting the relevance of the fly data for human biology. Collectively, these results not only lead to a better understanding of hVAPB function but also point to potentially relevant targets for therapeutic intervention.

The reduced kinome of Ostreococcus tauri: core eukaryotic signalling components in a tractable model species.

BACKGROUND: The current knowledge of eukaryote signalling originates from phenotypically diverse organisms. There is a pressing need to identify conserved signalling components among eukaryotes, which will lead to the transfer of knowledge across kingdoms. Two useful properties of a eukaryote model for signalling are (1) reduced signalling complexity, and (2) conservation of signalling components. The alga Ostreococcus tauri is described as the smallest free-living eukaryote. With less than 8,000 genes, it represents a highly constrained genomic palette. RESULTS: Our survey revealed 133 protein kinases and 34 protein phosphatases (1.7% and 0.4% of the proteome). We conducted phosphoproteomic experiments and constructed domain structures and phylogenies for the catalytic protein-kinases. For each of the major kinases families we review the completeness and divergence of O. tauri representatives in comparison to the well-studied kinomes of the laboratory models Arabidopsis thaliana and Saccharomyces cerevisiae, and of Homo sapiens. Many kinase clades in O. tauri were reduced to a single member, in preference to the loss of family diversity, whereas TKL and ABC1 clades were expanded. We also identified kinases that have been lost in A. thaliana but retained in O. tauri. For three, contrasting eukaryotic pathways - TOR, MAPK, and the circadian clock - we established the subset of conserved components and demonstrate conserved sites of substrate phosphorylation and kinase motifs. CONCLUSIONS: We conclude that O. tauri satisfies our two central requirements. Several of its kinases are more closely related to H. sapiens orthologs than S. cerevisiae is to H. sapiens. The greatly reduced kinome of O. tauri is therefore a suitable model for signalling in free-living eukaryotes.

NMDA receptor-dependent regulation of miRNA expression and association with Argonaute during LTP in vivo.

microRNAs (miRNAs) are major regulators of protein synthesis in the brain. A major goal is to identify changes in miRNA expression underlying protein synthesis-dependent forms of synaptic plasticity such as long-term potentiation (LTP). Previous analyses focused on changes in miRNA levels in total lysate samples. Here, we asked whether changes in total miRNA accurately reflect changes in the amount of miRNA bound to Argonaute protein within the miRNA-induced silencing complex (miRISC). Ago2 immunoprecipitation was used to isolate RISC-associated miRNAs following high-frequency stimulation (HFS)-induced LTP in the dentate gyrus of anesthetized rats. Using locked-nucleic acid-based PCR cards for high-throughput screening and independent validation by quantitative TaqMan RT-PCR, we identified differential regulation of Ago2-associated and total miRNA expression. The ratio of Ago2/total miRNA expression was regulated bidirectionally in a miRNA-specific manner and was largely dependent on N-methyl-D-aspartate receptor (NMDA) activation during LTP induction. The present results identify miRNA association with Ago2 as a potential control point in activity-dependent synaptic plasticity in the adult brain. Finally, novel computational analysis for targets of the Ago2-associated miRNAs identifies 21 pathways that are enriched and differentially targeted by the miRNAs including axon guidance, mTOR, MAPK, Ras, and LTP.

Pax6 exerts regional control of cortical progenitor proliferation via direct repression of Cdk6 and hypophosphorylation of pRb.

The mechanisms by which early spatiotemporal expression patterns of transcription factors such as Pax6 regulate cortical progenitors in a region-specific manner are poorly understood. Pax6 is expressed in a gradient across the developing cortex and is essential for normal corticogenesis. We found that constitutive or conditional loss of Pax6 increases cortical progenitor proliferation by amounts that vary regionally with normal Pax6 levels. We compared the gene expression profiles of equivalent Pax6-expressing progenitors isolated from Pax6⁺/⁺ and Pax6⁻/⁻ cortices and identified many negatively regulated cell-cycle genes, including Cyclins and Cdks. Biochemical assays indicated that Pax6 directly represses Cdk6 expression. Cyclin/Cdk repression inhibits retinoblastoma protein (pRb) phosphorylation, thereby limiting the transcription of genes that directly promote the mechanics of the cell cycle, and we found that Pax6 inhibits pRb phosphorylation and represses genes involved in DNA replication. Our results indicate that Pax6's modulation of cortical progenitor cell cycles is regional and direct.

Consolidation and translation regulation.

mRNA translation, or protein synthesis, is a major component of the transformation of the genetic code into any cellular activity. This complicated, multistep process is divided into three phases: initiation, elongation, and termination. Initiation is the step at which the ribosome is recruited to the mRNA, and is regarded as the major rate-limiting step in translation, while elongation consists of the elongation of the polypeptide chain; both steps are frequent targets for regulation, which is defined as a change in the rate of translation of an mRNA per unit time. In the normal brain, control of translation is a key mechanism for regulation of memory and synaptic plasticity consolidation, i.e., the off-line processing of acquired information. These regulation processes may differ between different brain structures or neuronal populations. Moreover, dysregulation of translation leads to pathological brain function such as memory impairment. Both normal and abnormal function of the translation machinery is believed to lead to translational up-regulation or down-regulation of a subset of mRNAs. However, the identification of these newly synthesized proteins and determination of the rates of protein synthesis or degradation taking place in different neuronal types and compartments at different time points in the brain demand new proteomic methods and system biology approaches. Here, we discuss in detail the relationship between translation regulation and memory or synaptic plasticity consolidation while focusing on a model of cortical-dependent taste learning task and hippocampal-dependent plasticity. In addition, we describe a novel systems biology perspective to better describe consolidation.