The SARS-CoV-2 Delta (Pango lineage B.1.617.2) variant of concern spread globally, causing resurgences of COVID-19 worldwide^1,2. The emergence of the Delta variant in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions. Here we analyse 52,992 SARS-CoV-2 genomes from England together with 93,649 genomes from the rest of the world to reconstruct the emergence of Delta and quantify its introduction to and regional dissemination across England in the context of changing travel and social restrictions. Using analysis of human movement, contact tracing and virus genomic data, we find that the geographic focus of the expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced more than 1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers reduced onward transmission from importations; however, the transmission chains that later dominated the Delta wave in England were seeded before travel restrictions were introduced. Increasing inter-regional travel within England drove the nationwide dissemination of Delta, with some cities receiving more than 2,000 observable lineage introductions from elsewhere. Subsequently, increased levels of local population mixing—and not the number of importations—were associated with the faster relative spread of Delta. The invasion dynamics of Delta depended on spatial heterogeneity in contact patterns, and our findings will inform optimal spatial interventions to reduce the transmission of current and future variants of concern, such as Omicron (Pango lineage B.1.1.529).

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) now arise in the context of heterogeneous human connectivity and population immunity. Through a large-scale phylodynamic analysis of 115,622 Omicron BA.1 genomes, we identified >6,000 introductions of the antigenically distinct VOC into England and analyzed their local transmission and dispersal history. We find that six of the eight largest English Omicron lineages were already transmitting when Omicron was first reported in southern Africa (22 November 2021). Multiple datasets show that importation of Omicron continued despite subsequent restrictions on travel from southern Africa as a result of export from well-connected secondary locations. Initiation and dispersal of Omicron transmission lineages in England was a two-stage process that can be explained by models of the country’s human geography and hierarchical travel network. Our results enable a comparison of the processes that drive the invasion of Omicron and other VOCs across multiple spatial scales.

Mosquito control with bednets, residual sprays or fumigation remains the most effective tool for preventing vector-borne diseases such as malaria, dengue and Zika, though there are no widely used entomological methods for directly assessing its efficacy. Mosquito age is the most informative metric for evaluating interventions that kill adult mosquitoes but there is no simple or reliable way of measuring it in the field. Near-Infrared Spectroscopy (NIRS) has been shown to be a promising, high-throughput method that can estimate the age of mosquitoes. Currently the ability of NIRS to measure mosquito age is biased, and has relatively high individual mosquito measurement error, though its capacity to rigorously monitor mosquito populations in the field has never been assessed. In this study, we use machine learning methods from the chemometric literature to generate more accurate, unbiased estimates of individual mosquito age. These unbiased estimates produce precise population-level measurements, which are relatively insensitive to further increases in NIRS accuracy when feasible numbers of mosquitoes are sampled. The utility of NIRS to directly measure the impact of pyrethroid resistance on mosquito control is illustrated, showing how the technology has potential as a highly valuable tool for directly assessing the efficacy of mosquito control interventions.

The adult mammalian kidney has a complex, highly-branched collecting duct epithelium that arises as a ureteric bud sidebranch from an epithelial tube known as the nephric duct. Subsequent branching of the ureteric bud to form the collecting duct tree is regulated by subcellular interactions between the epithelium and a population of mesenchymal cells that surround the tips of outgrowing branches. The mesenchymal cells produce glial cell-line derived neurotrophic factor (GDNF), that binds with RET receptors on the surface of the epithelial cells to stimulate several subcellular pathways in the epithelium. Such interactions are known to be a prerequisite for normal branching development, although competing theories exist for their role in morphogenesis. Here we introduce the first agent-based model of ex vivo kidney uretic branching. Through comparison with experimental data, we show that growth factor-regulated growth mechanisms can explain early epithelial cell branching, but only if epithelial cell division depends in a switch-like way on the local growth factor concentration; cell division occurring only if the driving growth factor level exceeds a threshold. We also show how a recently-developed method, “Approximate Approximate Bayesian Computation”, can be used to infer key model parameters, and reveal the dependency between the parameters controlling a growth factor-dependent growth switch. These results are consistent with a requirement for signals controlling proliferation and chemotaxis, both of which are previously identified roles for GDNF.

During sporogony, malaria-causing parasites infect a mosquito, reproduce and migrate to the mosquito salivary glands where they can be transmitted the next time blood feeding occurs. The time required for sporogony, known as the extrinsic incubation period (EIP), is an important determinant of malaria transmission intensity. The EIP is typically estimated as the time for a given percentile, x, of infected mosquitoes to develop salivary gland sporozoites (the infectious parasite life stage), which is denoted by EIP_x. Many mechanisms, however, affect the observed sporozoite prevalence including the human-to-mosquito transmission probability and possibly differences in mosquito mortality according to infection status. To account for these various mechanisms, we present a mechanistic mathematical model, which explicitly models key processes at the parasite, mosquito and observational scales. Fitting this model to experimental data, we find greater variation in the EIP than previously thought: we estimated the range between EIP₁₀ and EIP₉₀ (at 27°C) as 4.5 days compared to 0.9 days using existing statistical methods. This pattern holds over the range of study temperatures included in the dataset. Increasing temperature from 21°C to 34°C decreased the EIP₅₀ from 16.1 to 8.8 days. Our work highlights the importance of mechanistic modelling of sporogony to (1) improve estimates of malaria transmission under different environmental conditions or disease control programs and (2) evaluate novel interventions that target the mosquito life stages of the parasite.

Flapping flight is the most energetically demanding form of sustained forwards locomotion that vertebrates perform. Flock dynamics therefore have significant implications for energy expenditure. Despite this, no studies have quantified the biomechanical consequences of flying in a cluster flock or pair relative to flying solo. Here, we compared the flight characteristics of homing pigeons (Columba livia) flying solo and in pairs released from a site 7 km from home, using high-precision 5 Hz global positioning system (GPS) and 200 Hz tri-axial accelerometer bio-loggers. As expected, paired individuals benefitted from improved homing route accuracy, which reduced flight distance by 7% and time by 9%. However, realising these navigational gains involved substantial changes in flight kinematics and energetics. Both individuals in a pair increased their wingbeat frequency by 18% by decreasing the duration of their upstroke. This sharp increase in wingbeat frequency caused just a 3% increase in airspeed but reduced the oscillatory displacement of the body by 22%, which we hypothesise relates to an increased requirement for visual stability and manoeuvrability when flying in a flock or pair. The combination of the increase in airspeed and a higher wingbeat frequency would result in a minimum 2.2% increase in the total aerodynamic power requirements if the wingbeats were fully optimised. Overall, the enhanced navigational performance will offset any additional energetic costs as long as the metabolic power requirements are not increased above 9%. Our results demonstrate that the increases in wingbeat frequency when flying together have previously been underestimated by an order of magnitude and force reinterpretation of their mechanistic origin. We show that, for pigeons flying in pairs, two heads are better than one but keeping a steady head necessitates energetically costly kinematics.