Search

The availability of COVID-19 vaccines promised a reduction in the severity of disease and relief from the strict public health and social measures (PHSMs) imposed in many countries to limit spread and burden of COVID-19. We were asked to define vaccine coverage thresholds for Australia’s transition to easing restrictions and reopening international borders. Using evidence of vaccine effectiveness against the then-circulating Delta variant, we used a mathematical model to determine coverage targets. The absence of any COVID-19 infections in many sub-national jurisdictions in Australia posed particular methodological challenges. We used a novel metric called Transmission Potential (TP) as a proxy measure of the population-level effective reproduction number. We estimated TP of the Delta variant under a range of PHSMs, test-trace-isolate-quarantine (TTIQ) efficiencies, vaccination coverage thresholds, and age-based vaccine allocation strategies. We found that high coverage across all ages (>_70%) combined with ongoing TTIQ and minimal PHSMs was sufficient to avoid lockdowns. At lesser coverage (<_60%) rapid case escalation risked overwhelming of the health sector or the need to reimpose stricter restrictions. Maintaining low case numbers was most beneficial for health and the economy, and at higher coverage levels (>_80%) further easing of restrictions was deemed possible. These results directly informed easing of COVID-19 restrictions in Australia.

Disease surveillance data was critical in supporting public health decisions throughout the coronavirus disease 2019 (COVID-19) pandemic. At the same time, the unprecedented circumstances of the pandemic revealed many shortcomings of surveillance systems for viral respiratory pathogens. Strengthening of surveillance systems was identified as a priority for the recently established Australian Centre for Disease Control, which represents a critical opportunity to review pre-pandemic and pandemic surveillance practices, and to decide on future priorities, during both pandemic and inter-pandemic periods. On 20 October 2022, we ran a workshop with experts from the academic and government sectors who had contributed to the COVID-19 response in Australia on 'The role of surveillance in epidemic response', at the University of New South Wales, Sydney, Australia. Following the workshop, we developed five recommendations to strengthen respiratory virus surveillance systems in Australia, which we present here. Our recommendations are not intended to be exhaustive. We instead chose to focus on data types that are highly valuable yet typically overlooked by surveillance planners. Three of the recommendations focus on data collection activities that support the monitoring and prediction of disease impact and the effectiveness of interventions (what to measure) and two focus on surveillance methods and capabilities (how to measure). Implementation of our recommendations would enable more robust, timely, and impactful epidemic analysis.

Current malaria elimination targets must withstand a colossal challenge-resistance to the current gold standard antimalarial drug, namely artemisinin derivatives. If artemisinin resistance significantly expands to Africa or India, cases and malaria-related deaths are set to increase substantially. Spatial information on the changing levels of artemisinin resistance in Southeast Asia is therefore critical for health organisations to prioritise malaria control measures, but available data on artemisinin resistance are sparse. We use a comprehensive database from the WorldWide Antimalarial Resistance Network on the prevalence of non-synonymous mutations in the Kelch 13 (K13) gene, which are known to be associated with artemisinin resistance, and a Bayesian geostatistical model to produce spatio-temporal predictions of artemisinin resistance. Our maps of estimated prevalence show an expansion of the K13 mutation across the Greater Mekong Subregion from 2000 to 2022. Moreover, the period between 2010 and 2015 demonstrated the most spatial change across the region. Our model and maps provide important insights into the spatial and temporal trends of artemisinin resistance in a way that is not possible using data alone, thereby enabling improved spatial decision support systems on an unprecedented fine-scale spatial resolution. By predicting for the first time spatio-temporal patterns and extents of artemisinin resistance at the subcontinent level, this study provides critical information for supporting malaria elimination goals in Southeast Asia.

BACKGROUND: Australian states and territories used test-trace-isolate-quarantine (TTIQ) systems extensively in their response to the COVID-19 pandemic in 2020-2021. We report on an analysis of Australian case data to estimate the impact of test-trace-isolate-quarantine systems on SARS-CoV-2 transmission. METHODS: Our analysis uses a novel mathematical modelling framework and detailed surveillance data on COVID-19 cases including dates of infection and dates of isolation. First, we directly translate an empirical distribution of times from infection to isolation into reductions in potential for onward transmission during periods of relatively low caseloads (tens to hundreds of reported cases per day). We then apply a simulation approach, validated against case data, to assess the impact of case-initiated contact tracing on transmission during a period of relatively higher caseloads and system stress (up to thousands of cases per day). RESULTS: We estimate that under relatively low caseloads in the state of New South Wales (tens of cases per day), TTIQ contributed to a 54% reduction in transmission. Under higher caseloads in the state of Victoria (hundreds of cases per day), TTIQ contributed to a 42% reduction in transmission. Our results also suggest that case-initiated contact tracing can support timely quarantine in times of system stress (thousands of cases per day). CONCLUSION: Contact tracing systems for COVID-19 in Australia were highly effective and adaptable in supporting the national suppression strategy from 2020-21, prior to the emergence of the Omicron variant in November 2021. TTIQ systems were critical to the maintenance of the strong suppression strategy and were more effective when caseloads were (relatively) low.

Since the emergence of SARS-CoV-2 in 2019 through to mid-2021, much of the Australian population lived in a COVID-19-free environment. This followed the broadly successful implementation of a strong suppression strategy, including international border closures. With the availability of COVID-19 vaccines in early 2021, the national government sought to transition from a state of minimal incidence and strong suppression activities to one of high vaccine coverage and reduced restrictions but with still-manageable transmission. This transition is articulated in the national 're-opening' plan released in July 2021. Here, we report on the dynamic modelling study that directly informed policies within the national re-opening plan including the identification of priority age groups for vaccination, target vaccine coverage thresholds and the anticipated requirements for continued public health measures-assuming circulation of the Delta SARS-CoV-2 variant. Our findings demonstrated that adult vaccine coverage needed to be at least 60% to minimize public health and clinical impacts following the establishment of community transmission. They also supported the need for continued application of test-trace-isolate-quarantine and social measures during the vaccine roll-out phase and beyond.

BACKGROUND: Geostatistical analysis of health data is increasingly used to model spatial variation in malaria prevalence, burden, and other metrics. Traditional inference methods for geostatistical modelling are notoriously computationally intensive, motivating the development of newer, approximate methods for geostatistical analysis or, more broadly, computational modelling of spatial processes. The appeal of faster methods is particularly great as the size of the region and number of spatial locations being modelled increases. METHODS: This work presents an applied comparison of four proposed 'fast' computational methods for spatial modelling and the software provided to implement them-Integrated Nested Laplace Approximation (INLA), tree boosting with Gaussian processes and mixed effect models (GPBoost), Fixed Rank Kriging (FRK) and Spatial Random Forests (SpRF). The four methods are illustrated by estimating malaria prevalence on two different spatial scales-country and continent. The performance of the four methods is compared on these data in terms of accuracy, computation time, and ease of implementation. RESULTS: Two of these methods-SpRF and GPBoost-do not scale well as the data size increases, and so are likely to be infeasible for larger-scale analysis problems. The two remaining methods-INLA and FRK-do scale well computationally, however the resulting model fits are very sensitive to the user's modelling assumptions and parameter choices. The binomial observation distribution commonly used for disease prevalence mapping with INLA fails to account for small-scale overdispersion present in the malaria prevalence data, which can lead to poor predictions. Selection of an appropriate alternative such as the Beta-binomial distribution is required to produce a reliable model fit. The small-scale random effect term in FRK overcomes this pitfall, but FRK model estimates are very reliant on providing a sufficient number and appropriate configuration of basis functi

As of January 2021, Australia had effectively controlled local transmission of COVID-19 despite a steady influx of imported cases and several local, but contained, outbreaks in 2020. Throughout 2020, state and territory public health responses were informed by weekly situational reports that included an ensemble forecast of daily COVID-19 cases for each jurisdiction. We present here an analysis of one forecasting model included in this ensemble across the variety of scenarios experienced by each jurisdiction from May to October 2020. We examine how successfully the forecasts characterised future case incidence, subject to variations in data timeliness and completeness, showcase how we adapted these forecasts to support decisions of public health priority in rapidly-evolving situations, evaluate the impact of key model features on forecast skill, and demonstrate how to assess forecast skill in real-time before the ground truth is known. Conditioning the model on the most recent, but incomplete, data improved the forecast skill, emphasising the importance of developing strong quantitative models of surveillance system characteristics, such as ascertainment delay distributions. Forecast skill was highest when there were at least 10 reported cases per day, the circumstances in which authorities were most in need of forecasts to aid in planning and response.

BACKGROUND: Buruli ulcer (BU) is a neglected tropical disease caused by infection of subcutaneous tissue with Mycobacterium ulcerans. BU is commonly reported across rural regions of Central and West Africa but has been increasing dramatically in temperate southeast Australia around the major metropolitan city of Melbourne, with most disease transmission occurring in the summer months. Previous research has shown that Australian native possums are reservoirs of M. ulcerans and that they shed the bacteria in their fecal material (excreta). Field surveys show that locales where possums harbor M. ulcerans overlap with human cases of BU, raising the possibility of using possum excreta surveys to predict the risk of disease occurrence in humans. METHODS: We thus established a highly structured 12 month possum excreta surveillance program across an area of 350 km2 in the Mornington Peninsula area 70 km south of Melbourne, Australia. The primary objective of our study was to assess using statistical modeling if M. ulcerans surveillance of possum excreta provided useful information for predicting future human BU case locations. RESULTS: Over two sampling campaigns in summer and winter, we collected 2,282 possum excreta specimens of which 11% were PCR positive for M. ulcerans-specific DNA. Using the spatial scanning statistical tool SaTScan, we observed non-random, co-correlated clustering of both M. ulcerans positive possum excreta and human BU cases. We next trained a statistical model with the Mornington Peninsula excreta survey data to predict the future likelihood of human BU cases occurring in the region. By observing where human BU cases subsequently occurred, we show that the excreta model performance was superior to a null model trained using the previous year's human BU case incidence data (AUC 0.66 vs 0.55). We then used data unseen by the excreta-informed model from a new survey of 661 possum excreta specimens in Geelong, a geographically separate BU endemic area

Against a backdrop of widespread global transmission, a number of countries have successfully brought large outbreaks of COVID-19 under control and maintained near-elimination status. A key element of epidemic response is the tracking of disease transmissibility in near real-time. During major outbreaks, the effective reproduction number can be estimated from a time-series of case, hospitalisation or death counts. In low or zero incidence settings, knowing the potential for the virus to spread is a response priority. Absence of case data means that this potential cannot be estimated directly. We present a semi-mechanistic modelling framework that draws on time-series of both behavioural data and case data (when disease activity is present) to estimate the transmissibility of SARS-CoV-2 from periods of high to low - or zero - case incidence, with a coherent transition in interpretation across the changing epidemiological situations. Of note, during periods of epidemic activity, our analysis recovers the effective reproduction number, while during periods of low - or zero - case incidence, it provides an estimate of transmission risk. This enables tracking and planning of progress towards the control of large outbreaks, maintenance of virus suppression, and monitoring the risk posed by re-introduction of the virus. We demonstrate the value of our methods by reporting on their use throughout 2020 in Australia, where they have become a central component of the national COVID-19 response.

By mapping land use under projections of socio-economic change, ecological changes can be predicted to inform conservation decision-making. We present a land use model that enables the fine-scale mapping of land use change under future scenarios. Its predictions can be used as input to virtually all existing spatially-explicit ecological models. Our model maps the fractional cover of land use within each grid cell, providing higher information content than discrete classes at the same spatial resolution. The method accurately reproduced land use patterns observed in the Amazon, both in terms of the allocated fractional amounts and also the direction of predicted land use changes. A small case study showcases the application of our model to reproduce patterns of agricultural expansion and natural habitat declines. The model source code is provided as an open-source R package, making this new method available to bridge the gap between socio-economic, land use and biodiversity modelling.

Due to global climate change-induced shifts in species distributions, estimating changes in community composition through the use of Species Distribution Models has become a key management tool. Being able to determine how species associations change along environmental gradients is likely to be pivotal in exploring the magnitude of future changes in species' distributions. This is particularly important in connectivity-limited ecosystems, such as freshwater ecosystems, where increased human translocation is creating species associations over previously unseen environmental gradients. Here, we use a large-scale presence-absence dataset of freshwater fish from lakes across the Fennoscandian region in a Joint Species Distribution Model, to measure the effect of temperature on species associations. We identified a trend of negative associations between species tolerant of cold waters and those tolerant of warmer waters, as well as positive associations between several more warm-tolerant species, with these associations often shifting depending on local temperatures. Our results confirm that freshwater ecosystems can expect to see a large-scale shift towards communities dominated by more warm-tolerant species. While there remains much work to be done to predict exactly where and when local extinctions may take place, the model implemented provides a starting-point for the exploration of climate-driven community trends. This approach is especially informative in regards to determining which species associations are most central in shaping future community composition, and which areas are most vulnerable to local extinctions.

In an effort to effectively model observed patterns in the spatial configuration of individuals of multiple species in nature, we introduce the saturated pairwise interaction Gibbs point process. Its main strength lies in its ability to model both attraction and repulsion within and between species, over different scales. As such, it is particularly well-suited to the study of associations in complex ecosystems. Based on the existing literature, we provide an easy to implement fitting procedure as well as a technique to make inference for the model parameters. We also prove that under certain hypotheses the point process is locally stable, which allows us to use the well-known ‘coupling from the past’ algorithm to draw samples from the model. Different numerical experiments show the robustness of the model. We study three different ecological data sets, demonstrating in each one that our model helps disentangle competing ecological effects on species' distribution.

Joint species distribution models (JSDMs) simultaneously model the distributions of multiple species, while accounting for residual co‐occurrence patterns. Despite increasing adoption of JSDMs in the literature, the question of how to define and evaluate JSDM predictions has only begun to be explored. We define four different JSDM prediction types that correspond to different aspects of species distribution and community assemblage processes. Marginal predictions are environment‐only predictions akin to predictions from single‐species models; joint predictions simultaneously predict entire community assemblages; and conditional marginal and conditional joint predictions are made at the species or assemblage level, conditional on the known occurrence state of one or more species at a site. We define five different classes of metrics that can be used to evaluate these types of predictions: threshold‐dependent, threshold‐independent, community dissimilarity, species richness and likelihood metrics. We illustrate different prediction types and evaluation metrics using a case study in which we fit a JSDM to a frog occurrence dataset collected in Melbourne, Australia. Joint species distribution models present opportunities to investigate the facets of species distribution and community assemblage processes that are not possible to explore with single‐species models. We show that there are a variety of different metrics available to evaluate JSDM predictions, and that choice of prediction type and evaluation metric should closely match the questions being investigated.

Light‐level geolocators are popular bio‐logging tools, with advantageous sizes, longevity and affordability. Biologists tracking seabirds often presume geolocator spatial accuracies between 186 and 202 km from previously innovative, yet taxonomically, spatially and computationally limited, studies. Using recently developed methods, we investigated whether assumed uncertainty norms held across a larger‐scale, multispecies study. We field‐tested geolocator spatial accuracy by synchronously deploying these with GPS loggers on scores of seabirds across five species and 11 Mediterranean Sea, east Atlantic and south Pacific breeding colonies. We first interpolated geolocations using the geolocation package FLightR without prior knowledge of GPS tracked routes. We likewise applied another package, probGLS, additionally testing whether sea‐surface temperatures could improve route accuracy. Geolocator spatial accuracy was lower than the ~200 km often assumed. probGLS produced the best accuracy (mean ± SD = 304 ± 413 km, n = 185 deployments) with 84.5% of GPS‐derived latitudes and 88.8% of longitudes falling within resulting uncertainty estimates. FLightR produced lower spatial accuracy (408 ± 473 km, n = 171 deployments) with 38.6% of GPS‐derived latitudes and 23.7% of longitudes within package‐specific uncertainty estimates. Expected inter‐twilight period (from GPS position and date) was the strongest predictor of accuracy, with increasingly equatorial solar profiles (i.e. closer temporally to equinoxes and/or spatially to the Equator) inducing more error. Individuals, species and geolocator model also significantly affected accuracy, while the impact of distance travelled between successive twilights depended on the geolocation package. Geolocation accuracy is not uniform among seabird species and can be considerably lower than assumed. Individual idiosyncrasies and spatiotemporal dynamics (i.e. shallower inter‐twilight shifts by date and latitude) mean that pract

Climate change threatens biodiversity directly by influencing biophysical variables that drive species' geographic distributions and indirectly through socio-economic changes that influence land use patterns, driven by global consumption, production and climate. To date, no detailed analyses have been produced that assess the relative importance of, or interaction between, these direct and indirect climate change impacts on biodiversity at large scales. Here, we apply a new integrated modelling framework to quantify the relative influence of biophysical and socio-economically mediated impacts on avian species in Vietnam and Australia and we find that socio-economically mediated impacts on suitable ranges are largely outweighed by biophysical impacts. However, by translating economic futures and shocks into spatially explicit predictions of biodiversity change, we now have the power to analyse in a consistent way outcomes for nature and people of any change to policy, regulation, trading conditions or consumption trend at any scale from sub-national to global.

Understanding human movement patterns at local, national and international scales is critical in a range of fields, including transportation, logistics and epidemiology. Data on human movement is increasingly available, and when combined with statistical models, enables predictions of movement patterns across broad regions. Movement characteristics, however, strongly depend on the scale and type of movement captured for a given study. The models that have so far been proposed for human movement are best suited to specific spatial scales and types of movement. Selecting both the scale of data collection, and the appropriate model for the data remains a key challenge in predicting human movements. We used two different data sources on human movement in Australia, at different spatial scales, to train a range of statistical movement models and evaluate their ability to predict movement patterns for each data type and scale. Whilst the five commonly-used movement models we evaluated varied markedly between datasets in their predictive ability, we show that an ensemble modelling approach that combines the predictions of these models consistently outperformed all individual models against hold-out data.

Approximately 150 triatomine species are suspected to be infected with the Chagas parasite, Trypanosoma cruzi, but they differ in the risk they pose to human populations. The largest risk comes from species that have a domestic life cycle and these species have been targeted by indoor residual spraying campaigns, which have been successful in many locations. It is now important to consider residual transmission that may be linked to persistent populations of dominant vectors, or to secondary or minor vectors. The aim of this project was to define the geographical distributions of the community of triatomine species across the Chagas endemic region. Presence-only data with over 12, 000 observations of triatomine vectors were extracted from a public database and target-group background data were generated to account for sampling bias in the presence data. Geostatistical regression was then applied to estimate species distributions and fine-scale distribution maps were generated for thirty triatomine vector species including those found within one or two countries and species that are more widely distributed from northern Argentina to Guatemala, Bolivia to southern Mexico, and Mexico to the southern United States of America. The results for Rhodnius pictipes, Panstrongylus geniculatus, Triatoma dimidiata, Triatoma gerstaeckeri, and Triatoma infestans are presented in detail, including model predictions and uncertainty in these predictions, and the model validation results for each of the 30 species are presented in full. The predictive maps for all species are made publicly available so that they can be used to assess the communities of vectors present within different regions of the endemic zone. The maps are presented alongside key indicators for the capacity of each species to transmit T. cruzi to humans. These indicators include infection prevalence, evidence for human blood meals, and colonisation or invasion of homes. A summary of the published evidence for these

With the expansion in the quantity and types of biodiversity data being collected, there is a need to find ways to combine these different sources to provide cohesive summaries of species' potential and realized distributions in space and time. Recently, model-based data integration has emerged as a means to achieve this by combining datasets in ways that retain the strengths of each. We describe a flexible approach to data integration using point process models, which provide a convenient way to translate across ecological currencies. We highlight recent examples of large-scale ecological models based on data integration and outline the conceptual and technical challenges and opportunities that arise.

Globally, protected areas are being established to protect biodiversity and to promote ecosystem resilience. The typical spatial conservation planning process leading to the creation of these protected areas focuses on representation and replication of ecological features, often using decision support tools such as Marxan. Yet, despite the important role ecological connectivity has in metapopulation persistence and resilience, Marxan currently requires manual input or specialized scripts to explicitly consider connectivity. ‘Marxan Connect’ is a new open source, open access Graphical User Interface (GUI) tool designed to assist conservation planners with the appropriate use of data on ecological connectivity in protected area network planning. Marxan Connect can facilitate the use of estimates of demographic connectivity (e.g. derived from animal tracking data, dispersal models, or genetic tools) or structural landscape connectivity (e.g. isolation by resistance). This is accomplished by calculating metapopulation‐relevant connectivity metrics (e.g. eigenvector centrality) and treating those as conservation features or by including the connectivity data as a spatial dependency amongst sites in the prioritization process. Marxan Connect allows a wide group of users to incorporate directional ecological connectivity into conservation planning with Marxan. The solutions provided by Marxan Connect, combined with ecologically relevant post‐hoc testing, are more likely to support persistent and resilient metapopulations (e.g. fish stocks) and provide better protection for biodiversity.

Species population dynamics are driven by spatial and temporal changes in the environment, anthropogenic activities and conservation management actions. Understanding how populations will change in response to these drivers is fundamental to a wide range of ecological applications, but there are few open‐source software options accessible to researchers and managers that allow them to predict these changes in a flexible and transparent way. We introduce an open‐source, multi‐platformrpackage,steps, that models spatial changes in species populations as a function of drivers of distribution and abundance, such as climate, disturbance, landscape dynamics and species ecological and physiological requirements. To illustrate the functionality ofsteps, we model the population dynamics of the greater gliderPetauroides volans, an arboreal Australian mammal. We demonstrate howstepscan be used to simulate population responses of the glider to forest dynamics and management with the types of data commonly used in ecological analyses. stepsexpands on the features found in existing software packages, can easily incorporate a range of spatial layers (e.g. habitat suitability, vegetation dynamics and disturbances), facilitates integrated and transparent analyses within a single platform and produces interpretable outputs of changes in species' populations through space and time. Further,stepsoffers both ready‐to‐use, built‐in functionality, as well as the ability for advanced users to define their own modules for custom analyses. Thus, we anticipate thatstepswill be of significant value to environment and wildlife managers and researchers from a broad range of disciplines.

BACKGROUND: Asymptomatic or subclinical SARS-CoV-2 infections are often unreported, which means that confirmed case counts may not accurately reflect underlying epidemic dynamics. Understanding the level of ascertainment (the ratio of confirmed symptomatic cases to the true number of symptomatic individuals) and undetected epidemic progression is crucial to informing COVID-19 response planning, including the introduction and relaxation of control measures. Estimating case ascertainment over time allows for accurate estimates of specific outcomes such as seroprevalence, which is essential for planning control measures. METHODS: Using reported data on COVID-19 cases and fatalities globally, we estimated the proportion of symptomatic cases (i.e. any person with any of fever ≥ 37.5 °C, cough, shortness of breath, sudden onset of anosmia, ageusia or dysgeusia illness) that were reported in 210 countries and territories, given those countries had experienced more than ten deaths. We used published estimates of the baseline case fatality ratio (CFR), which was adjusted for delays and under-ascertainment, then calculated the ratio of this baseline CFR to an estimated local delay-adjusted CFR to estimate the level of under-ascertainment in a particular location. We then fit a Bayesian Gaussian process model to estimate the temporal pattern of under-ascertainment. RESULTS: Based on reported cases and deaths, we estimated that, during March 2020, the median percentage of symptomatic cases detected across the 84 countries which experienced more than ten deaths ranged from 2.4% (Bangladesh) to 100% (Chile). Across the ten countries with the highest number of total confirmed cases as of 6 July 2020, we estimated that the peak number of symptomatic cases ranged from 1.4 times (Chile) to 18 times (France) larger than reported. Comparing our model with national and regional seroprevalence data where available, we find that our estimates are consistent with observed values. Finally

As of 1 May 2020, there had been 6808 confirmed cases of COVID-19 in Australia. Of these, 98 had died from the disease. The epidemic had been in decline since mid-March, with 308 cases confirmed nationally since 14 April. This suggests that the collective actions of the Australian public and government authorities in response to COVID-19 were sufficiently early and assiduous to avert a public health crisis - for now. Analysing factors that contribute to individual country experiences of COVID-19, such as the intensity and timing of public health interventions, will assist in the next stage of response planning globally. We describe how the epidemic and public health response unfolded in Australia up to 13 April. We estimate that the effective reproduction number was likely below one in each Australian state since mid-March and forecast that clinical demand would remain below capacity thresholds over the forecast period (from mid-to-late April).

Background

Methods

Results

Understanding the conditions under which small native Australian mammals can persist in the presence of introduced predators remains a key challenge to conservation ecologists. Bettong‐specific one‐way gates were used at a predator‐free reserve in South Australia to allow the burrowing bettong (Bettongia lesueur) – a small potoroid, listed as ‘vulnerable’ nationally – to disperse out of the reserve. We conducted a field experiment to explore the conditions affecting residence time of bettongs that left the reserve. We monitored bettong and mammalian predator activity outside the fence using track surveys across 18 sites over two seasons. We examined the effect of supplementary feeding as a strategy for increasing residence time, as well as the influence of predator presence and habitat quality, using linear mixed models. Bettong activity was positively associated with supplementary feeding, midstorey vegetation cover and shelter availability. After gates were closed, bettong activity near gates declined to almost zero the following weeks, likely either due to death from predation or due to movement away from the sites. To a small extent, mammalian predators were more likely to be present at sites with high bettong activity. Further research on conditions to support persistence of burrowing bettongs and other small mammals, including understanding minimum necessary predator control effort, is required before successful establishment of populations outside of fences can occur.

Dengue has enormous health impacts globally. A novel approach to decrease dengue incidence involves the introduction of Wolbachia endosymbionts that block dengue virus transmission into populations of the primary vector mosquito, Aedes aegypti. The wMel Wolbachia strain has previously been trialed in open releases of Ae. aegypti; however, the wAlbB strain has been shown to maintain higher density than wMel at high larval rearing temperatures. Releases of Ae. aegypti mosquitoes carrying wAlbB were carried out in 6 diverse sites in greater Kuala Lumpur, Malaysia, with high endemic dengue transmission. The strain was successfully established and maintained at very high population frequency at some sites or persisted with additional releases following fluctuations at other sites. Based on passive case monitoring, reduced human dengue incidence was observed in the release sites when compared to control sites. The wAlbB strain of Wolbachia provides a promising option as a tool for dengue control, particularly in very hot climates.

This Article was mistakenly not made Open Access when originally published; this has now been amended, and information about the Creative Commons Attribution 4.0 International License has been added into the 'Additional information' section.

The global population at risk from mosquito-borne diseases-including dengue, yellow fever, chikungunya and Zika-is expanding in concert with changes in the distribution of two key vectors: Aedes aegypti and Aedes albopictus. The distribution of these species is largely driven by both human movement and the presence of suitable climate. Using statistical mapping techniques, we show that human movement patterns explain the spread of both species in Europe and the United States following their introduction. We find that the spread of Ae. aegypti is characterized by long distance importations, while Ae. albopictus has expanded more along the fringes of its distribution. We describe these processes and predict the future distributions of both species in response to accelerating urbanization, connectivity and climate change. Global surveillance and control efforts that aim to mitigate the spread of chikungunya, dengue, yellow fever and Zika viruses must consider the so far unabated spread of these mosquitos. Our maps and predictions offer an opportunity to strategically target surveillance and control programmes and thereby augment efforts to reduce arbovirus burden in human populations globally.

Dengue is a mosquito-borne viral infection that has spread throughout the tropical world over the past 60 years and now affects over half the world's population. The geographical range of dengue is expected to further expand due to ongoing global phenomena including climate change and urbanization. We applied statistical mapping techniques to the most extensive database of case locations to date to predict global environmental suitability for the virus as of 2015. We then made use of climate, population and socioeconomic projections for the years 2020, 2050 and 2080 to project future changes in virus suitability and human population at risk. This study is the first to consider the spread of Aedes mosquito vectors to project dengue suitability. Our projections provide a key missing piece of evidence for the changing global threat of vector-borne disease and will help decision-makers worldwide to better prepare for and respond to future changes in dengue risk.

1. Joint species distribution models (JSDMs) account for biotic interactions and missing environmental predictors in correlative species distribution models. Several different JSDMs have been proposed in the literature, but the use of different or conflicting nomenclature and statistical notation potentially obscures similarities and differences among them. Furthermore, new JSDM implementations have been illustrated with different case studies, preventing direct comparisons of computational and statistical performance. 2. We aim to resolve these outstanding issues by (a) highlighting similarities among seven presence–absence JSDMs using a clearly defined, singular notation; and (b) evaluating the computational and statistical performance of each JSDM using six datasets that vary widely in numbers of sites, species, and environmental covariates considered. 3. Our singular notation shows that many of the JSDMs are very similar, and in turn parameter estimates of different JSDMs are moderate to strongly, positively correlated. In contrast, the different JSDMs clearly differ in computational efficiency and memory limitations. 4. Our framework will allow ecologists to make educated decisions about the JSDM that best suits their objective, and enable wider uptake of JSDM methods among the ecological community.

Human mobility is an important driver of geographic spread of infectious pathogens. Detailed information about human movements during outbreaks are, however, difficult to obtain and may not be available during future epidemics. The Ebola virus disease (EVD) outbreak in West Africa between 2014-16 demonstrated how quickly pathogens can spread to large urban centers following one cross-species transmission event. Here we describe a flexible transmission model to test the utility of generalised human movement models in estimating EVD cases and spatial spread over the course of the outbreak. A transmission model that includes a general model of human mobility significantly improves prediction of EVD's incidence compared to models without this component. Human movement plays an important role not only to ignite the epidemic in locations previously disease free, but over the course of the entire epidemic. We also demonstrate important differences between countries in population mixing and the improved prediction attributable to movement metrics. Given their relative rareness, locally derived mobility data are unlikely to exist in advance of future epidemics or pandemics. Our findings show that transmission patterns derived from general human movement models can improve forecasts of spatio-temporal transmission patterns in places where local mobility data is unavailable.

Since 2000, many countries have achieved considerable success in improving child survival, but localized progress remains unclear. To inform efforts towards United Nations Sustainable Development Goal 3.2—to end preventable child deaths by 2030—we need consistently estimated data at the subnational level regarding child mortality rates and trends. Here we quantified, for the period 2000–2017, the subnational variation in mortality rates and number of deaths of neonates, infants and children under 5 years of age within 99 low- and middle-income countries using a geostatistical survival model. We estimated that 32% of children under 5 in these countries lived in districts that had attained rates of 25 or fewer child deaths per 1,000 live births by 2017, and that 58% of child deaths between 2000 and 2017 in these countries could have been averted in the absence of geographical inequality. This study enables the identification of high-mortality clusters, patterns of progress and geographical inequalities to inform appropriate investments and implementations that will help to improve the health of all populations. © 2019, The Author(s).

BACKGROUND: Yellow fever cases are under-reported and the exact distribution of the disease is unknown. An effective vaccine is available but more information is needed about which populations within risk zones should be targeted to implement interventions. Substantial outbreaks of yellow fever in Angola, Democratic Republic of the Congo, and Brazil, coupled with the global expansion of the range of its main urban vector, Aedes aegypti, suggest that yellow fever has the propensity to spread further internationally. The aim of this study was to estimate the disease's contemporary distribution and potential for spread into new areas to help inform optimal control and prevention strategies. METHODS: We assembled 1155 geographical records of yellow fever virus infection in people from 1970 to 2016. We used a Poisson point process boosted regression tree model that explicitly incorporated environmental and biological explanatory covariates, vaccination coverage, and spatial variability in disease reporting rates to predict the relative risk of apparent yellow fever virus infection at a 5 × 5 km resolution across all risk zones (47 countries across the Americas and Africa). We also used the fitted model to predict the receptivity of areas outside at-risk zones to the introduction or reintroduction of yellow fever transmission. By use of previously published estimates of annual national case numbers, we used the model to map subnational variation in incidence of yellow fever across at-risk countries and to estimate the number of cases averted by vaccination worldwide. FINDINGS: Substantial international and subnational spatial variation exists in relative risk and incidence of yellow fever as well as varied success of vaccination in reducing incidence in several high-risk regions, including Brazil, Cameroon, and Togo. Areas with the highest predicted average annual case numbers include large parts of Nigeria, the Democratic Republic of the Congo, and South Sudan, where vac

BACKGROUND: Podoconiosis is a non-filarial elephantiasis, which causes massive swelling of the lower legs. It was identified as a neglected tropical disease by WHO in 2011. Understanding of the geographical distribution of the disease is incomplete. As part of a global mapping of podoconiosis, this study was conducted in Cameroon to map the distribution of the disease. This mapping work will help to generate data on the geographical distribution of podoconiosis in Cameroon and contribute to the global atlas of podoconiosis. METHODS: We used a multi-stage sampling design with stratification of the country by environmental risk of podoconiosis. We sampled 76 villages from 40 health districts from the ten Regions of Cameroon. All individuals of 15-years old or older in the village were surveyed house-to-house and screened for lymphedema. A clinical algorithm was used to reliably diagnose podoconiosis, excluding filarial-associated lymphedema. Individuals with lymphoedema were tested for circulating Wuchereria bancrofti antigen and specific IgG4 using the Alere Filariasis Test Strips (FTS) test and the Standard Diagnostics (SD) BIOLINE lymphatic filariasis IgG4 test (Wb123) respectively, in addition to thick blood films. Presence of DNA specific to W. bancrofti was checked on night blood using a qPCR technique. PRINCIPAL FINDINGS: Overall, 10,178 individuals from 4,603 households participated in the study. In total, 83 individuals with lymphedema were identified. Of the 83 individuals with lymphedema, two were found to be FTS positive and all were negative using the Wb123 test. No microfilaria of W. bancrofti were found in the night blood of any individual with clinical lymphedema. None were found to be positive for W. bancrofti using qPCR. Of the two FTS positive cases, one was positive for Mansonella perstans DNA, while the other harbored Loa loa microfilaria. Overall, 52 people with podoconiosis were identified after applying the clinical algorithm. The overall preva