Tracking pathogen transmissibility during infectious disease outbreaks is essential for assessing the effectiveness of public health measures and planning future control strategies. A key measure of transmissibility is the time-dependent reproduction number, which has been estimated in real-time during outbreaks of a range of pathogens from disease incidence time series data. While commonly used approaches for estimating the time-dependent reproduction number can be reliable when disease incidence is recorded frequently, such incidence data are often aggregated temporally (for example, numbers of cases may be reported weekly rather than daily). As we show, commonly used methods for estimating transmissibility can be unreliable when the timescale of transmission is shorter than the timescale of data recording. To address this, here we develop a simulation-based approach involving Approximate Bayesian Computation for estimating the time-dependent reproduction number from temporally aggregated disease incidence time series data. We first use a simulated dataset representative of a situation in which daily disease incidence data are unavailable and only weekly summary values are reported, demonstrating that our method provides accurate estimates of the time-dependent reproduction number under such circumstances. We then apply our method to two outbreak datasets consisting of weekly influenza case numbers in 2019-20 and 2022-23 in Wales (in the United Kingdom). Our simple-to-use approach will allow accurate estimates of time-dependent reproduction numbers to be obtained from temporally aggregated data during future infectious disease outbreaks.

Surface water nutrient pollution, the primary cause of eutrophication, remains a major environmental concern in Western Lake Erie despite intergovernmental efforts to regulate nutrient sources. The Maumee River Basin has been the largest nutrient contributor. The two primary nutrient sources are inorganic fertilizer and livestock manure applied to croplands, which are later carried to the streams via runoff and soil erosion. Prior studies of nutrient source attribution have focused on large watersheds or counties at annual time scales. Source attribution at finer spatiotemporal scales, which enables more effective nutrient management, remains a substantial challenge. This study aims to address this challenge by developing a generalizable Bayesian network model for phosphorus source attribution at the subwatershed scale (12-digit Hydrologic Unit Code). Since phosphorus release is uncertain, we combine excess phosphorus derived from manure and fertilizer application and crop uptake data, flow information simulated by the SWAT model, and in-stream water quality measurements using Approximate Bayesian Computation to derive a posterior that attributes phosphorus contributions to subwatersheds. Our results show significant variability in subwatershed-scale phosphorus release that is lost in coarse-scale attribution. Phosphorus contributions attributed to the subwatersheds are on average lower than the excess phosphorus estimated by the nutrient balance approach currently adopted by environmental agencies. Fertilizer contributes more soluble reactive phosphorus than manure, while manure contributes most of the unreactive phosphorus. While developed for the specific context of Maumee River Basin, our lightweight and generalizable model framework could be adapted to other regions and pollutants and could help inform targeted environmental regulation and enforcement.

African Swine Fever (ASF) disease transmission parameters are crucial for making response and control decisions when faced with an outbreak, yet they are poorly quantified for smallholder and village contexts within Southeast Asia. Whilst disease-specific factors - such as latent and infectious periods - should remain reasonably consistent, host, environmental and management factors are likely to affect the rate of disease spread. These differences are investigated using Approximate Bayesian Computation with Sequential Monte-Carlo methods to provide disease parameter estimates in four naïve pig populations in villages of Lao People\'s Democratic Republic. The villages represent smallholder pig farmers of the Northern province of Oudomxay and the Southern province of Savannakhet, and the model utilised field mortality data to validate the transmission parameter estimates over the course of multiple model generations. The basic reproductive number between-pigs was estimated to range from 3.08 to 7.80, whilst the latent and infectious periods were consistent with those published in the literature for similar genotypes in the region (4.72 to 6.19 days and 2.63 to 5.50 days, respectively). These findings demonstrate that smallholder village pigs interact similarly to commercial pigs, however the spread of disease may occur slightly slower than in commercial study groups. Furthermore, the findings demonstrated that despite diversity across the study groups, the disease behaved in a consistent manner. This data can be used in disease control programs or for future modelling of ASF in smallholder contexts.

The abnormal aggregation of extracellular amyloid-β(Aβ) in senile plaques resulting in calcium Ca+2 dyshomeostasis is one of the primary symptoms of Alzheimer\'s disease (AD). Significant research efforts have been devoted in the past to better understand the underlying molecular mechanisms driving Aβ deposition and Ca+2 dysregulation. Importantly, synaptic impairments, neuronal loss, and cognitive failure in AD patients are all related to the buildup of intraneuronal Aβ accumulation. Moreover, increasing evidence show a feed-forward loop between Aβ and Ca+2 levels, i.e. Aβ disrupts neuronal Ca+2 levels, which in turn affects the formation of Aβ. To better understand this interaction, we report a novel stochastic model where we analyze the positive feedback loop between Aβ and Ca+2 using ADNI data. A good therapeutic treatment plan for AD requires precise predictions. Stochastic models offer an appropriate framework for modelling AD since AD studies are observational in nature and involve regular patient visits. The etiology of AD may be described as a multi-state disease process using the approximate Bayesian computation method. So, utilizing ADNI data from 2-year visits for AD patients, we employ this method to investigate the interplay between Aβ and Ca+2 levels at various disease development phases. Incorporating the ADNI data in our physics-based Bayesian model, we discovered that a sufficiently large disruption in either Aβ metabolism or intracellular Ca+2 homeostasis causes the relative growth rate in both Ca+2 and Aβ, which corresponds to the development of AD. The imbalance of Ca+2 ions causes Aβ disorders by directly or indirectly affecting a variety of cellular and subcellular processes, and the altered homeostasis may worsen the abnormalities of Ca+2 ion transportation and deposition. This suggests that altering the Ca+2 balance or the balance between Aβ and Ca+2 by chelating them may be able to reduce disorders associated with AD and open up new research possibilities for AD therapy.

The development of mathematical models for studying newly emerging and re-emerging infectious diseases has gained momentum due to global events. The gyrodactylid-fish system, like many host-parasite systems, serves as a valuable resource for ecological, evolutionary, and epidemiological investigations owing to its ease of experimental manipulation and long-term monitoring. Although this system has an existing individual-based model, it falls short in capturing information about species-specific microhabitat preferences and other biological details for different Gyrodactylus strains across diverse fish populations. This current study introduces a new individual-based stochastic simulation model that uses a hybrid τ -leaping algorithm to incorporate this essential data, enhancing our understanding of the complexity of the gyrodactylid-fish system. We compare the infection dynamics of three gyrodactylid strains across three host populations. A modified sequential-type approximate Bayesian computation (ABC) method, based on sequential Monte Carlo and sequential importance sampling, is developed. Additionally, we establish two penalised local-linear regression methods (based on L1 and L2 regularisations) for ABC post-processing analysis to fit our model using existing empirical data. With the support of experimental data and the fitted mathematical model, we address open biological questions for the first time and propose directions for future studies on the gyrodactylid-fish system. The adaptability of the mathematical model extends beyond the gyrodactylid-fish system to other host-parasite systems. Furthermore, the modified ABC methodologies provide efficient calibration for other multi-parameter models characterised by a large set of correlated or independent summary statistics.

Characterizing the processes underlying reproductive isolation between diverging lineages is central to understanding speciation. Here, we present RIDGE-Reproductive Isolation Detection using Genomic polymorphisms-a tool tailored for quantifying gene flow barrier proportion and identifying the relevant genomic regions. RIDGE relies on an Approximate Bayesian Computation with a model-averaging approach to accommodate diverse scenarios of lineage divergence. It captures heterogeneity in effective migration rate along the genome while accounting for variation in linked selection and recombination. The barrier detection test relies on numerous summary statistics to compute a Bayes factor, offering a robust statistical framework that facilitates cross-species comparisons. Simulations revealed RIDGE\'s efficiency in capturing signals of ongoing migration. Model averaging proved particularly valuable in scenarios of high model uncertainty where no migration or migration homogeneity can be wrongly assumed, typically for recent divergence times <0.1 2Ne generations. Applying RIDGE to four published crow data sets, we first validated our tool by identifying a well-known large genomic region associated with mate choice patterns. Second, while we identified a significant overlap of outlier loci using RIDGE and traditional genomic scans, our results suggest that a substantial portion of previously identified outliers are likely false positives. Outlier detection relies on allele differentiation, relative measures of divergence and the count of shared polymorphisms and fixed differences. Our analyses also highlight the value of incorporating multiple summary statistics including our newly developed outlier ones that can be useful in challenging detection conditions.

Fibroblasts in a confluent monolayer are known to adopt elongated morphologies in which cells are oriented parallel to their neighbors. We collected and analyzed new microscopy movies to show that confluent fibroblasts are motile and that neighboring cells often move in anti-parallel directions in a collective motion phenomenon we refer to as \"fluidization\" of the cell population. We used machine learning to perform cell tracking for each movie and then leveraged topological data analysis (TDA) to show that time-varying point-clouds generated by the tracks contain significant topological information content that is driven by fluidization, i.e., the anti-parallel movement of individual neighboring cells and neighboring groups of cells over long distances. We then utilized the TDA summaries extracted from each movie to perform Bayesian parameter estimation for the D\'Orsgona model, an agent-based model (ABM) known to produce a wide array of different patterns, including patterns that are qualitatively similar to fluidization. Although the D\'Orsgona ABM is a phenomenological model that only describes inter-cellular attraction and repulsion, the estimated region of D\'Orsogna model parameter space was consistent across all movies, suggesting that a specific level of inter-cellular repulsion force at close range may be a mechanism that helps drive fluidization patterns in confluent mesenchymal cell populations.

Landlocking of diadromous fish in freshwater systems can have significant genomic consequences. For instance, the loss of the migratory life stage can dramatically reduce gene flow across populations, leading to increased genetic structuring and stronger effects of local adaptation. These genomic consequences have been well-studied in some mainland systems, but the evolutionary impacts of landlocking in island ecosystems are largely unknown. In this study, we used a genotyping-by-sequencing (GBS) approach to examine the evolutionary history of landlocking in common smelt (Retropinna retropinna) on Chatham Island, a small isolated oceanic island 800 kilometres east of mainland New Zealand. We examined the relationship between Chatham Island and mainland smelt and used coalescent analyses to test the number and timing of landlocking events on Chatham Island. Our genomic analysis, based on 21,135 SNPs across 169 individuals, revealed that the Chatham Island smelt was genomically distinct from the mainland New Zealand fish, consistent with a single ancestral colonisation event of Chatham Island in the Pleistocene. Significant genetic structure was also evident within the Chatham Island smelt, with a diadromous Chatham Island smelt group, along with three geographically structured landlocked groups. Coalescent demographic analysis supported three independent landlocking events, with this loss of diadromy significantly pre-dating human colonisation. Our results illustrate how landlocking of diadromous fish can occur repeatedly across a narrow spatial scale, and highlight a unique system to study the genomic basis of repeated adaptation.

Introduction: Atrial fibrillation (AF) is the most common arrhythmia, associated with significant burdens to patients and the healthcare system. The atrioventricular (AV) node plays a vital role in regulating heart rate during AF by filtering electrical impulses from the atria. However, it is often insufficient in regards to maintaining a healthy heart rate, thus the AV node properties are modified using rate-control drugs. Moreover, treatment selection during permanent AF is currently done empirically. Quantifying individual differences in diurnal and short-term variability of AV-nodal function could aid in personalized treatment selection. Methods: This study presents a novel methodology for estimating the refractory period (RP) and conduction delay (CD) trends, and their uncertainty in the two pathways of the AV node during 24 h using non-invasive data. This was achieved by utilizing a network model together with a problem-specific genetic algorithm and an approximate Bayesian computation algorithm. Diurnal variability in the estimated RP and CD was quantified by the difference between the daytime and nighttime estimates, and short-term variability was quantified by the Kolmogorov-Smirnov distance between adjacent 10-min segments in the 24-h trends. Additionally, the predictive value of the derived parameter trends regarding drug outcome was investigated using several machine learning tools. Results: Holter electrocardiograms from 51 patients with permanent AF during baseline were analyzed, and the predictive power of variations in RP and CD on the resulting heart rate reduction after treatment with four rate control drugs was investigated. Diurnal variability yielded no correlation to treatment outcome, and no prediction of drug outcome was possible using the machine learning tools. However, a correlation between the short-term variability for the RP and CD in the fast pathway and resulting heart rate reduction during treatment with metoprolol (ρ = 0.48, p < 0.005 in RP, ρ = 0.35, p < 0.05 in CD) were found. Discussion: The proposed methodology enables non-invasive estimation of the AV node properties during 24 h, which-indicated by the correlation between the short-term variability and heart rate reduction-may have the potential to assist in treatment selection.

The Australian lungfish is a primitive and endangered representative of the subclass Dipnoi. The distribution of this species is limited to south-east Queensland, with some populations considered endemic and others possibly descending from translocations in the late nineteenth century shortly after European discovery. Attempts to resolve the historical distribution of this species have met with conflicting results based on descriptive genetic studies. Understanding if all populations are endemic or some are the result of, or influenced by, translocation events, has implications for conservation management. In this work, we analysed the genetic variation at three types of markers (mtDNA genomes, 11 STRs and 5196 nuclear SNPs) using the approximate Bayesian computation (ABC) algorithm to compare several demographic models. We postulated different contributions of Mary River and Burnett River gene pools into the Brisbane River and North Pine River populations, related to documented translocation events. We ran the analysis for each marker type separately, and we also estimated the posterior probabilities of the models combining the markers. Nuclear SNPs have the highest power to correctly identify the true model among the simulated datasets (where the model was known), but different marker types typically provided similar answers. The most supported demographic model able to explain the real dataset implies that an endemic gene pool is still present in the Brisbane and North Pine Rivers and coexists with the gene pools derived from past documented translocation events. These results support the view that ABC modelling can be useful to reconstruct complex historical translocation events with contemporary implications, and will inform ongoing conservation efforts for the endangered and iconic Australian lungfish.