Phage-induced efflux down-regulation boosts antibiotic efficacy.

The interactions between a virus and its host vary in space and time and are affected by the presence of molecules that alter the physiology of either the host or the virus. Determining the molecular mechanisms at the basis of these interactions is paramount for predicting the fate of bacterial and phage populations and for designing rational phage-antibiotic therapies. We study the interactions between stationary phase Burkholderia thailandensis and the phage ΦBp-AMP1. Although heterogeneous genetic resistance to phage rapidly emerges in B. thailandensis, the presence of phage enhances the efficacy of three major antibiotic classes, the quinolones, the beta-lactams and the tetracyclines, but antagonizes tetrahydrofolate synthesis inhibitors. We discovered that enhanced antibiotic efficacy is facilitated by reduced antibiotic efflux in the presence of phage. This new phage-antibiotic therapy allows for eradication of stationary phase bacteria, whilst requiring reduced antibiotic concentrations, which is crucial for treating infections in sites where it is difficult to achieve high antibiotic concentrations.

Combined multiplex panel test results are a poor estimate of disease prevalence without adjustment for test error.

Multiplex panel tests identify many individual pathogens at once, using a set of component tests. In some panels the number of components can be large. If the panel is detecting causative pathogens for a single syndrome or disease then we might estimate the burden of that disease by combining the results of the panel, for example determining the prevalence of pneumococcal pneumonia as caused by many individual pneumococcal serotypes. When we are dealing with multiplex test panels with many components, test error in the individual components of a panel, even when present at very low levels, can cause significant overall error. Uncertainty in the sensitivity and specificity of the individual tests, and statistical fluctuations in the numbers of false positives and false negatives, will cause large uncertainty in the combined estimates of disease prevalence. In many cases this can be a source of significant bias. In this paper we develop a mathematical framework to characterise this issue, we determine expressions for the sensitivity and specificity of panel tests. In this we identify a counter-intuitive relationship between panel test sensitivity and disease prevalence that means panel tests become more sensitive as prevalence increases. We present novel statistical methods that adjust for bias and quantify uncertainty in prevalence estimates from panel tests, and use simulations to test these methods. As multiplex testing becomes more commonly used for screening in routine clinical practice, accumulation of test error due to the combination of large numbers of test results needs to be identified and corrected for.

HormoneBayes: A novel Bayesian framework for the analysis of pulsatile hormone dynamics.

The hypothalamus is the central regulator of reproductive hormone secretion. Pulsatile secretion of gonadotropin releasing hormone (GnRH) is fundamental to physiological stimulation of the pituitary gland to release luteinizing hormone (LH) and follicle stimulating hormone (FSH). Furthermore, GnRH pulsatility is altered in common reproductive disorders such as polycystic ovary syndrome (PCOS) and hypothalamic amenorrhea (HA). LH is measured routinely in clinical practice using an automated chemiluminescent immunoassay method and is the gold standard surrogate marker of GnRH. LH can be measured at frequent intervals (e.g., 10 minutely) to assess GnRH/LH pulsatility. However, this is rarely done in clinical practice because it is resource intensive, and there is no open-access, graphical interface software for computational analysis of the LH data available to clinicians. Here we present hormoneBayes, a novel open-access Bayesian framework that can be easily applied to reliably analyze serial LH measurements to assess LH pulsatility. The framework utilizes parsimonious models to simulate hypothalamic signals that drive LH dynamics, together with state-of-the-art (sequential) Monte-Carlo methods to infer key parameters and latent hypothalamic dynamics. We show that this method provides estimates for key pulse parameters including inter-pulse interval, secretion and clearance rates and identifies LH pulses in line with the widely used deconvolution method. We show that these parameters can distinguish LH pulsatility in different clinical contexts including in reproductive health and disease in men and women (e.g., healthy men, healthy women before and after menopause, women with HA or PCOS). A further advantage of hormoneBayes is that our mathematical approach provides a quantified estimation of uncertainty. Our framework will complement methods enabling real-time in-vivo hormone monitoring and therefore has the potential to assist translation of personalized, data-driven, clinical care of patients presenting with conditions of reproductive hormone dysfunction.

Smoking during pregnancy and its effect on placental weight: a Mendelian randomization study.

BACKGROUND: The causal relationship between maternal smoking in pregnancy and reduced offspring birth weight is well established and is likely due to impaired placental function. However, observational studies have given conflicting results on the association between smoking and placental weight. We aimed to estimate the causal effect of newly pregnant mothers quitting smoking on their placental weight at the time of delivery. METHODS: We used one-sample Mendelian randomization, drawing data from the Avon Longitudinal Study of Parents and Children (ALSPAC) (N = 690 to 804) and the Norwegian Mother, Father and Child Cohort Study (MoBa) (N = 4267 to 4606). The sample size depends on the smoking definition used for different analyses. The analysis was performed in pre-pregnancy smokers only, due to the specific role of the single-nucleotide polymorphism (SNP) rs1051730 (CHRNA5 - CHRNA3 - CHRNB4) in affecting smoking cessation but not initiation. RESULTS: Fixed effect meta-analysis showed a 182 g [95%CI: 29,335] higher placental weight for pre-pregnancy smoking mothers who continued smoking at the beginning of pregnancy, compared with those who stopped smoking. Using the number of cigarettes smoked per day in the first trimester as the exposure, the causal effect on placental weight was 11 g [95%CI: 1,21] per cigarette per day. Similarly, smoking at the end of pregnancy was causally associated with higher placental weight. Using the residuals of birth weight regressed on placental weight as the outcome, we showed evidence of lower offspring birth weight relative to the placental weight, both for continuing smoking at the start of pregnancy as well as continuing smoking throughout pregnancy (change in z-score birth weight adjusted for z-score placental weight: -0.8 [95%CI: -1.6,-0.1]). CONCLUSION: Our results suggest that continued smoking during pregnancy causes higher placental weights.

Microbial Dysregulation of the Gut-Lung Axis in Bronchiectasis.

Rationale: Emerging data support the existence of a microbial "gut-lung" axis that remains unexplored in bronchiectasis. Methods: Prospective and concurrent sampling of gut (stool) and lung (sputum) was performed in a cohort of n = 57 individuals with bronchiectasis and subjected to bacteriome (16S rRNA) and mycobiome (18S Internal Transcribed Spacer) sequencing (total, 228 microbiomes). Shotgun metagenomics was performed in a subset (n = 15; 30 microbiomes). Data from gut and lung compartments were integrated by weighted similarity network fusion, clustered, and subjected to co-occurrence analysis to evaluate gut-lung networks. Murine experiments were undertaken to validate specific Pseudomonas-driven gut-lung interactions. Results: Microbial communities in stable bronchiectasis demonstrate a significant gut-lung interaction. Multibiome integration followed by unsupervised clustering reveals two patient clusters, differing by gut-lung interactions and with contrasting clinical phenotypes. A high gut-lung interaction cluster, characterized by lung Pseudomonas, gut Bacteroides, and gut Saccharomyces, is associated with increased exacerbations and greater radiological and overall bronchiectasis severity, whereas the low gut-lung interaction cluster demonstrates an overrepresentation of lung commensals, including Prevotella, Fusobacterium, and Porphyromonas with gut Candida. The lung Pseudomonas-gut Bacteroides relationship, observed in the high gut-lung interaction bronchiectasis cluster, was validated in a murine model of lung Pseudomonas aeruginosa infection. This interaction was abrogated after antibiotic (imipenem) pretreatment in mice confirming the relevance and therapeutic potential of targeting the gut microbiome to influence the gut-lung axis. Metagenomics in a subset of individuals with bronchiectasis corroborated our findings from targeted analyses. Conclusions: A dysregulated gut-lung axis, driven by lung Pseudomonas, associates with poorer clinical outcomes in bronchiectasis.

Sensitisation to recombinant Aspergillus fumigatus allergens and clinical outcomes in COPD.

BACKGROUND: Variable clinical outcomes are reported with fungal sensitisation in chronic obstructive pulmonary disease (COPD), and it remains unclear which fungi and what allergens associate with the poorest outcomes. The use of recombinant as opposed to crude allergens for such assessment is unknown. METHODS: A prospective multicentre assessment of stable COPD (n=614) was undertaken in five hospitals across three countries: Singapore, Malaysia and Hong Kong. Clinical and serological assessment was performed against a panel of 35 fungal allergens including crude and recombinant Aspergillus and non-Aspergillus allergens. Unsupervised clustering and topological data analysis (TDA) approaches were employed using the measured sensitisation responses to elucidate if sensitisation subgroups exist and their related clinical outcomes. RESULTS: Aspergillus fumigatus sensitisation was associated with increased exacerbations in COPD. Unsupervised cluster analyses revealed two "fungal sensitisation" groups. The first was characterised by Aspergillus sensitisation and increased exacerbations, poorer lung function and worse prognosis. Polysensitisation in this group conferred even poorer outcome. The second group, characterised by Cladosporium sensitisation, was more symptomatic. Significant numbers of individuals demonstrated sensitisation responses to only recombinant (as opposed to crude) A. fumigatus allergens f 1, 3, 5 and 6, and exhibited increased exacerbations, poorer lung function and an overall worse prognosis. TDA validated these findings and additionally identified a subgroup within Aspergillus-sensitised COPD of patients with frequent exacerbations. CONCLUSION: Aspergillus sensitisation is a treatable trait in COPD. Measuring sensitisation responses to recombinant Aspergillus allergens identifies an important patient subgroup with poor COPD outcomes that remains overlooked by assessment of only crude Aspergillus allergens.

Using immersive virtual reality to remotely examine performance differences between dominant and non-dominant hands.

Circle drawing may be a useful task to study upper-limb function in patient populations. However, previous studies rely on expensive and bulky robotics to measure performance. For clinics or hospitals with limited budgets and space, this may be unfeasible. Virtual reality (VR) provides a portable and low-cost tool with integrated motion capture. It offers potentially a more feasible medium by which to assess upper-limb motor function. Prior to use with patient populations, it is important to validate and test the capabilities of VR with healthy users. This study examined whether a VR-based circle drawing task, completed remotely using participant's own devices, could capture differences between movement kinematics of the dominant and non-dominant hands in healthy individuals. Participants (n = 47) traced the outline of a circle presented on their VR head-mounted displays with each hand, while the positions of the hand-held controllers were continuously recorded. Although there were no differences observed in the size or roundness of circles drawn with each hand, consistent with prior literature our results did show that the circles drawn with the dominant hand were completed faster than those with the non-dominant hand. This provides preliminary evidence that a VR-based circle drawing task may be a feasible method for detecting subtle differences in function in clinical populations. Supplementary Information: The online version contains supplementary material available at 10.1007/s10055-023-00794-z.

Separable actions of acetylcholine and noradrenaline on neuronal ensemble formation in hippocampal CA3 circuits.

In the hippocampus, episodic memories are thought to be encoded by the formation of ensembles of synaptically coupled CA3 pyramidal cells driven by sparse but powerful mossy fiber inputs from dentate gyrus granule cells. The neuromodulators acetylcholine and noradrenaline are separately proposed as saliency signals that dictate memory encoding but it is not known if they represent distinct signals with separate mechanisms. Here, we show experimentally that acetylcholine, and to a lesser extent noradrenaline, suppress feed-forward inhibition and enhance Excitatory-Inhibitory ratio in the mossy fiber pathway but CA3 recurrent network properties are only altered by acetylcholine. We explore the implications of these findings on CA3 ensemble formation using a hierarchy of models. In reconstructions of CA3 pyramidal cells, mossy fiber pathway disinhibition facilitates postsynaptic dendritic depolarization known to be required for synaptic plasticity at CA3-CA3 recurrent synapses. We further show in a spiking neural network model of CA3 how acetylcholine-specific network alterations can drive rapid overlapping ensemble formation. Thus, through these distinct sets of mechanisms, acetylcholine and noradrenaline facilitate the formation of neuronal ensembles in CA3 that encode salient episodic memories in the hippocampus but acetylcholine selectively enhances the density of memory storage.

Algorithmic hospital catchment area estimation using label propagation.

BACKGROUND: Hospital catchment areas define the primary population of a hospital and are central to assessing the potential demand on that hospital, for example, due to infectious disease outbreaks. METHODS: We present a novel algorithm, based on label propagation, for estimating hospital catchment areas, from the capacity of the hospital and demographics of the nearby population, and without requiring any data on hospital activity. RESULTS: The algorithm is demonstrated to produce a mapping from fine grained geographic regions to larger scale catchment areas, providing contiguous and realistic subdivisions of geographies relating to a single hospital or to a group of hospitals. In validation against an alternative approach predicated on activity data gathered during the COVID-19 outbreak in the UK, the label propagation algorithm is found to have a high level of agreement and perform at a similar level of accuracy. RESULTS: The algorithm can be used to make estimates of hospital catchment areas in new situations where activity data is not yet available, such as in the early stages of a infections disease outbreak.

scReQTL: an approach to correlate SNVs to gene expression from individual scRNA-seq datasets.

BACKGROUND: Recently, pioneering expression quantitative trait loci (eQTL) studies on single cell RNA sequencing (scRNA-seq) data have revealed new and cell-specific regulatory single nucleotide variants (SNVs). Here, we present an alternative QTL-related approach applicable to transcribed SNV loci from scRNA-seq data: scReQTL. ScReQTL uses Variant Allele Fraction (VAFRNA) at expressed biallelic loci, and corelates it to gene expression from the corresponding cell. RESULTS: Our approach employs the advantage that, when estimated from multiple cells, VAFRNA can be used to assess effects of SNVs in a single sample or individual. In this setting scReQTL operates in the context of identical genotypes, where it is likely to capture RNA-mediated genetic interactions with cell-specific and transient effects. Applying scReQTL on scRNA-seq data generated on the 10 × Genomics Chromium platform using 26,640 mesenchymal cells derived from adipose tissue obtained from three healthy female donors, we identified 1272 unique scReQTLs. ScReQTLs common between individuals or cell types were consistent in terms of the directionality of the relationship and the effect size. Comparative assessment with eQTLs from bulk sequencing data showed that scReQTL analysis identifies a distinct set of SNV-gene correlations, that are substantially enriched in known gene-gene interactions and significant genome-wide association studies (GWAS) loci. CONCLUSION: ScReQTL is relevant to the rapidly growing source of scRNA-seq data and can be applied to outline SNVs potentially contributing to cell type-specific and/or dynamic genetic interactions from an individual scRNA-seq dataset. AVAILABILITY: https://github.com/HorvathLab/NGS/tree/master/scReQTL.

ReQTL: identifying correlations between expressed SNVs and gene expression using RNA-sequencing data.

MOTIVATION: By testing for associations between DNA genotypes and gene expression levels, expression quantitative trait locus (eQTL) analyses have been instrumental in understanding how thousands of single nucleotide variants (SNVs) may affect gene expression. As compared to DNA genotypes, RNA genetic variation represents a phenotypic trait that reflects the actual allele content of the studied system. RNA genetic variation at expressed SNV loci can be estimated using the proportion of alleles bearing the variant nucleotide (variant allele fraction, VAFRNA). VAFRNA is a continuous measure which allows for precise allele quantitation in loci where the RNA alleles do not scale with the genotype count. We describe a method to correlate VAFRNA with gene expression and assess its ability to identify genetically regulated expression solely from RNA-sequencing (RNA-seq) datasets. RESULTS: We introduce ReQTL, an eQTL modification which substitutes the DNA allele count for the variant allele fraction at expressed SNV loci in the transcriptome (VAFRNA). We exemplify the method on sets of RNA-seq data from human tissues obtained though the Genotype-Tissue Expression (GTEx) project and demonstrate that ReQTL analyses are computationally feasible and can identify a subset of expressed eQTL loci. AVAILABILITY AND IMPLEMENTATION: A toolkit to perform ReQTL analyses is available at https://github.com/HorvathLab/ReQTL. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Domino-like transient dynamics at seizure onset in epilepsy.

The International League Against Epilepsy (ILAE) groups seizures into "focal", "generalized" and "unknown" based on whether the seizure onset is confined to a brain region in one hemisphere, arises in several brain region simultaneously, or is not known, respectively. This separation fails to account for the rich diversity of clinically and experimentally observed spatiotemporal patterns of seizure onset and even less so for the properties of the brain networks generating them. We consider three different patterns of domino-like seizure onset in Idiopathic Generalized Epilepsy (IGE) and present a novel approach to classification of seizures. To understand how these patterns are generated on networks requires understanding of the relationship between intrinsic node dynamics and coupling between nodes in the presence of noise, which currently is unknown. We investigate this interplay here in the framework of domino-like recruitment across a network. In particular, we use a phenomenological model of seizure onset with heterogeneous coupling and node properties, and show that in combination they generate a range of domino-like onset patterns observed in the IGE seizures. We further explore the individual contribution of heterogeneous node dynamics and coupling by interpreting in-vitro experimental data in which the speed of onset can be chemically modulated. This work contributes to a better understanding of possible drivers for the spatiotemporal patterns observed at seizure onset and may ultimately contribute to a more personalized approach to classification of seizure types in clinical practice.

Reduced Models of Cardiomyocytes Excitability: Comparing Karma and FitzHugh-Nagumo.

Since Noble adapted in 1962 the model of Hodgkin and Huxley to fit Purkinje fibres, the refinement of models for cardiomyocytes has continued. Most of these models are high-dimensional systems of coupled equations so that the possible mathematical analysis is quite limited, even numerically. This has inspired the development of reduced, phenomenological models that preserve qualitatively the main feature of cardiomyocyte's dynamics. In this paper, we present a systematic comparison of the dynamics between two notable low-dimensional models, the FitzHugh-Nagumo model (FitzHugh in Bull Math Biophys 17:257-269, 1955, J Gen Physiol 43:867-896, 1960, Biophys J 1:445-466, 1961) as a prototype of excitable behaviour and a polynomial version of the Karma model (Karma in Phys Rev Lett 71(7):16, 1993, Chaos 4:461, 1994) which is specifically developed to fit cardiomyocyte's behaviour well. We start by introducing the models and considering their pure ODE versions. We analyse the ODEs employing the main ideas and steps used in the setting of geometric singular perturbation theory. Next, we turn to the spatially extended models, where we focus on travelling wave solutions in 1D. Finally, we perform numerical simulations of the 1D PDE Karma model varying model parameters in order to systematically investigate the impact on wave propagation velocity and shape. In summary, our study provides a reference regarding key similarities as well as key differences of the two models.

Metagenomics Reveals a Core Macrolide Resistome Related to Microbiota in Chronic Respiratory Disease.

Rationale: Long-term antibiotic use for managing chronic respiratory disease is increasing; however, the role of the airway resistome and its relationship to host microbiomes remains unknown.Objectives: To evaluate airway resistomes and relate them to host and environmental microbiomes using ultradeep metagenomic shotgun sequencing.Methods: Airway specimens from 85 individuals with and without chronic respiratory disease (severe asthma, chronic obstructive pulmonary disease, and bronchiectasis) were subjected to metagenomic sequencing to an average depth exceeding 20 million reads. Respiratory and device-associated microbiomes were evaluated on the basis of taxonomical classification and functional annotation including the Comprehensive Antibiotic Resistance Database to determine airway resistomes. Co-occurrence networks of gene-microbe association were constructed to determine potential microbial sources of the airway resistome. Paired patient-inhaler metagenomes were compared (n = 31) to assess for the presence of airway-environment overlap in microbiomes and/or resistomes.Measurements and Main Results: Airway metagenomes exhibit taxonomic and metabolic diversity and distinct antimicrobial resistance patterns. A "core" airway resistome dominated by macrolide but with high prevalence of ß-lactam, fluoroquinolone, and tetracycline resistance genes exists and is independent of disease status or antibiotic exposure. Streptococcus and Actinomyces are key potential microbial reservoirs of macrolide resistance including the ermX, ermF, and msrD genes. Significant patient-inhaler overlap in airway microbiomes and their resistomes is identified where the latter may be a proxy for airway microbiome assessment in chronic respiratory disease.Conclusions: Metagenomic analysis of the airway reveals a core macrolide resistome harbored by the host microbiome.

Decoding identity from motion: how motor similarities colour our perception of self and others.

For more than 4 decades, it has been shown that humans are particularly sensitive to biological motion and extract socially relevant information from it such as gender, intentions, emotions or a person's identity. A growing number of findings, however, indicate that identity perception is not always highly accurate, especially due to large inter-individual differences and a fuzzy self-recognition advantage compared to the recognition of others. Here, we investigated the self-other identification performance and sought to relate this performance to the metric properties of perceptual/physical representations of individual motor signatures. We show that identity perception ability varies substantially across individuals and is associated to the perceptual/physical motor similarities between self and other stimuli. Specifically, we found that the perceptual representations of postural signatures are veridical in the sense that closely reflects the physical postural trajectories and those similarities between people' actions elicit numerous misattributions. While, on average, people can well recognize their self-generated actions, they more frequently attribute to themselves the actions of those acting in a similar way. These findings are consistent with the common coding theory and support that perception and action are tightly linked and may modulate each other by virtue of similarity.

The Origin of GnRH Pulse Generation: An Integrative Mathematical-Experimental Approach.

Fertility critically depends on the gonadotropin-releasing hormone (GnRH) pulse generator, a neural construct comprised of hypothalamic neurons coexpressing kisspeptin, neurokoinin-B and dynorphin. Here, using mathematical modeling and in vivo optogenetics we reveal for the first time how this neural construct initiates and sustains the appropriate ultradian frequency essential for reproduction. Prompted by mathematical modeling, we show experimentally using female estrous mice that robust pulsatile release of luteinizing hormone, a proxy for GnRH, emerges abruptly as we increase the basal activity of the neuronal network using continuous low-frequency optogenetic stimulation. Further increase in basal activity markedly increases pulse frequency and eventually leads to pulse termination. Additional model predictions that pulsatile dynamics emerge from nonlinear positive and negative feedback interactions mediated through neurokinin-B and dynorphin signaling respectively are confirmed neuropharmacologically. Our results shed light on the long-elusive GnRH pulse generator offering new horizons for reproductive health and wellbeing.SIGNIFICANCE STATEMENT The gonadotropin-releasing hormone (GnRH) pulse generator controls the pulsatile secretion of the gonadotropic hormones LH and FSH and is critical for fertility. The hypothalamic arcuate kisspeptin neurons are thought to represent the GnRH pulse generator, since their oscillatory activity is coincident with LH pulses in the blood; a proxy for GnRH pulses. However, the mechanisms underlying GnRH pulse generation remain elusive. We developed a mathematical model of the kisspeptin neuronal network and confirmed its predictions experimentally, showing how LH secretion is frequency-modulated as we increase the basal activity of the arcuate kisspeptin neurons in vivo using continuous optogenetic stimulation. Our model provides a quantitative framework for understanding the reproductive neuroendocrine system and opens new horizons for fertility regulation.

Oscillatory stimuli differentiate adapting circuit topologies.

Biology emerges from interactions between molecules, which are challenging to elucidate with current techniques. An orthogonal approach is to probe for 'response signatures' that identify specific circuit motifs. For example, bistability, hysteresis, or irreversibility are used to detect positive feedback loops. For adapting systems, such signatures are not known. Only two circuit motifs generate adaptation: negative feedback loops (NFLs) and incoherent feed-forward loops (IFFLs). On the basis of computational testing and mathematical proofs, we propose differential signatures: in response to oscillatory stimulation, NFLs but not IFFLs show refractory-period stabilization (robustness to changes in stimulus duration) or period skipping. Applying this approach to yeast, we identified the circuit dominating cell cycle timing. In Caenorhabditis elegans AWA neurons, which are crucial for chemotaxis, we uncovered a Ca2+ NFL leading to adaptation that would be difficult to find by other means. These response signatures allow direct access to the outlines of the wiring diagrams of adapting systems.

Neurologically Motivated Coupling Functions in Models of Motor Coordination.

We present an analysis of two Haken-Kelso-Bunz (HKB) oscillators coupled by a neurologically motivated function. We study the effect of time delay and weighted self-feedback and mutual feedback on the synchronization behavior of the model. We focus on identifying parameter regimes supporting experimentally observed decrease in oscillation amplitude and loss of anti-phase stability that has inspired the development of the HKB model. We show that a combination of cross-talk and nonlinearity in the coupling, along with physiologically relevant time delay, is able to quantitatively account for both drop in oscillation amplitude and loss of anti-phase stability in a frequency dependent manner. Furthermore, we demonstrate that the transition between discrete and rhythmic movements could be captured by this model. To this end, we carry out theoretical and numerical analysis of the emergence of in-phase and anti-phase oscillations.

Shaw and Shal voltage-gated potassium channels mediate circadian changes in Drosophila clock neuron excitability.

As in mammals, Drosophila circadian clock neurons display rhythms of activity with higher action potential firing rates and more positive resting membrane potentials during the day. This rhythmic excitability has been widely observed but, critically, its regulation remains unresolved. We have characterized and modelled the changes underlying these electrical activity rhythms in the lateral ventral clock neurons (LNvs). We show that currents mediated by the voltage-gated potassium channels Shaw (Kv3) and Shal (Kv4) oscillate in a circadian manner. Disruption of these channels, by expression of dominant negative (DN) subunits, leads to changes in circadian locomotor activity and shortens lifespan. LNv whole-cell recordings then show that changes in Shaw and Shal currents drive changes in action potential firing rate and that these rhythms are abolished when the circadian molecular clock is stopped. A whole-cell biophysical model using Hodgkin-Huxley equations can recapitulate these changes in electrical activity. Based on this model and by using dynamic clamp to manipulate clock neurons directly, we can rescue the pharmacological block of Shaw and Shal, restore the firing rhythm, and thus demonstrate the critical importance of Shaw and Shal. Together, these findings point to a key role for Shaw and Shal in controlling circadian firing of clock neurons and show that changes in clock neuron currents can account for this. Moreover, with dynamic clamp we can switch the LNvs between morning-like and evening-like states of electrical activity. We conclude that changes in Shaw and Shal underlie the daily oscillation in LNv firing rate.