Recasting the theory of mosquito-borne pathogen transmission dynamics and control

Mosquito-borne diseases pose some of the greatest challenges in public health, especially in tropical and sub-tropical regions of theworld. Efforts to control these diseases have been underpinned by a theoretical framework developed for malaria by Ross and Macdonald, including models, metrics for measuring transmission, and theory of control that identifies key vulnerabilities in the transmission cycle. That framework, especially Macdonald’s formula for R0 and its entomological derivative, vectorial capacity, are nowused to study dynamics and design interventions for many mosquito-borne diseases. A systematic review of 388 models published between 1970 and 2010 found that the vast majority adopted the Ross–Macdonald assumption of homogeneous transmission in a well-mixed population. Studies comparing models and data question these assumptions and point to the capacity to model heterogeneous, focal transmission as the most important but relatively unexplored component in current theory. Fine-scale heterogeneity causes transmission dynamics to be nonlinear, and poses problems for modeling, epidemiology and measurement. Novel mathematical approaches show how heterogeneity arises from the biology and the landscape on which the processes of mosquito biting and pathogen transmission unfold. Emerging theory focuses attention on the ecological and social context formosquito blood feeding, themovement of both hosts and mosquitoes, and the relevant spatial scales for measuring transmission and for modeling dynamics and control.

Characteristics of spatial variation, size distribution, formation and growth of particles in urban environments

This thesis is the first comprehensive study of important parameters relating to aerosols' impact on climate and human health, namely spatial variation, particle size distribution and new particle formation. We determined the importance of spatial variation of particle number concentration in microscale environments, developed a method for particle size parameterisation and provided knowledge about the chemistry of new particle formation. This is a significant contribution to our understanding of processes behind the transformation and dynamics of urban aerosols. This PhD project included extensive measurements of air quality parameters using state of the art instrumentation at each of the 25 sites within the Brisbane metropolitan area and advanced statistical analysis.

Modelling the dynamics of Plasmodium falciparum histidine-rich protein 2 in human malaria to better understand malaria rapid diagnostic test performance

BACKGROUND: Effective diagnosis of malaria is a major component of case management. Rapid diagnostic tests (RDTs) based on Plasmodium falciparumhistidine-rich protein 2 (PfHRP2) are popular for diagnosis of this most virulent malaria infection. However, concerns have been raised about the longevity of the PfHRP2 antigenaemia following curative treatment in endemic regions. METHODS: A model of PfHRP2 production and decay was developed to mimic the kinetics of PfHRP2 antigenaemia during infections. Data from two human infection studies was used to fit the model, and to investigate PfHRP2 kinetics. Four malaria RDTs were assessed in the laboratory to determine the minimum detectable concentration of PfHRP2. RESULTS: Fitting of the PfHRP2 dynamics model indicated that in malaria naive hosts, P. falciparum parasites of the 3D7 strain produce 1.4 x 10(-)(1)(3) g of PfHRP2 per parasite per replication cycle. The four RDTs had minimum detection thresholds between 6.9 and 27.8 ng/mL. Combining these detection thresholds with the kinetics of PfHRP2, it is predicted that as few as 8 parasites/muL may be required to maintain a positive RDT in a chronic infection. CONCLUSIONS: The results of the model indicate that good quality PfHRP2-based RDTs should be able to detect parasites on the first day of symptoms, and that the persistence of the antigen will cause the tests to remain positive for at least seven days after treatment. The duration of a positive test result following curative treatment is dependent on the duration and density of parasitaemia prior to treatment and the presence and affinity of anti-PfHRP2 antibodies.

The dynamics of cellular automata on 2-manifolds is affected by topology

In this paper, we demonstrate that the distribution of Wolfram classes within a cellular automata rule space in the triangular tessellation is not consistent across different topological general. Using a statistical mechanics approach, cellular automata dynamical classes were approximated for cellular automata defined on genus-0, genus-1 and genus-2 2-manifolds. A distribution-free equality test for empirical distributions was applied to identify cases in which Wolfram classes were distributed differently across topologies. This result implies that global structure and local dynamics contribute to the long term evolution of cellular automata.

Particle size distribution effects on preferential deposition areas in metal foam wrapped tube bundle

This paper presents a numerical model for understanding particle transport and deposition in metal foam heat exchangers. Two-dimensional steady and unsteady numerical simulations of a standard single row metal foam-wrapped tube bundle are performed for different particle size distributions, i.e. uniform and normal distributions. Effects of different particle sizes and fluid inlet velocities on the overall particle transport inside and outside the foam layer are also investigated. It was noted that the simplification made in the previously-published numerical works in the literature, e.g. uniform particle deposition in the foam, is not necessarily accurate at least for the cases considered here. The results highlight the preferential particle deposition areas both along the tube walls and inside the foam using a developed particle deposition likelihood matrix. This likelihood matrix is developed based on three criteria being particle local velocity, time spent in the foam, and volume fraction. It was noted that the particles tend to deposit near both front and rear stagnation points. The former is explained by the higher momentum and direct exposure of the particles to the foam while the latter only accommodate small particles which can be entrained in the recirculation region formed behind the foam-wrapped tubes.

Formation of carbon nanoscrolls from graphene nanoribbons : a molecular dynamics study

Carbon nanoscrolls (CNSs) are one of the carbon-based nanomaterials similar to carbon nanotubes (CNTs) but are not widely studied in spite of their great potential applications. Their practical applications are hindered by the challenging fabrication of the CNSs. A physical approach has been proposed recently to fabricate the CNS by rolling up a monolayer graphene nanoribbon (GNR) around a CNT driven by the interaction energy between them. In this study, we perform extensive molecular dynamics (MD) simulations to investigate the various factors that impact the formation of the CNS from GNR. Our simulation results show that the formation of the CNS is sensitive to the length of the CNT and temperature. When the GNR is functionalized with hydrogen, the formation of the CNS is determined by the density and distribution of the hydrogen atoms. Graphyne, the allotrope of graphene, is inferior to graphene in the formation of the CNS due to the weaker bonds and the associated smaller atom density. The mechanism behind the rolling of GNR into CNS lies in the balance between the GNR–CNT van der Waals (vdW) interactions and the strain energy of GNR. The present work reveals new important insights and provides useful guidelines for the fabrication of the CNS.

Stability analysis of STATCOM in distribution networks

This chapter presents the stability analysis based on bifurcation theory of the distribution static compensator (DSTATCOM) operating both in current control mode as in voltage control mode. The bifurcation analysis allows delimiting the operating zones of nonlinear power systems and hence the computation of these boundaries is of interest for practical design and planning purposes. Suitable mathematical representations of the DSTATCOM are proposed to carry out the bifurcation analyses efficiently. The stability regions in the Thevenin equivalent plane are computed for different power factors at the Point of Common Coupling (PCC). In addition, the stability regions in the control gain space are computed, and the DC capacitor and AC capacitor impact on the stability are analyzed in detail. It is shown through bifurcation analysis that the loss of stability in the DSTATCOM is in general due to the emergence of oscillatory dynamics. The observations are verified through detailed simulation studies.

A simulation model of the within-host dynamics of Plasmodium vivax infection

Background The benign reputation of Plasmodium vivax is at odds with the burden and severity of the disease. This reputation, combined with restricted in vitro techniques, has slowed efforts to gain an understanding of the parasite biology and interaction with its human host. Methods A simulation model of the within-host dynamics of P. vivax infection is described, incorporating distinctive characteristics of the parasite such as the preferential invasion of reticulocytes and hypnozoite production. The developed model is fitted using digitized time-series’ from historic neurosyphilis studies, and subsequently validated against summary statistics from a larger study of the same population. The Chesson relapse pattern was used to demonstrate the impact of released hypnozoites. Results The typical pattern for dynamics of the parasite population is a rapid exponential increase in the first 10 days, followed by a gradual decline. Gametocyte counts follow a similar trend, but are approximately two orders of magnitude lower. The model predicts that, on average, an infected naïve host in the absence of treatment becomes infectious 7.9 days post patency and is infectious for a mean of 34.4 days. In the absence of treatment, the effect of hypnozoite release was not apparent as newly released parasites were obscured by the existing infection. Conclusions The results from the model provides useful insights into the dynamics of P. vivax infection in human hosts, in particular the timing of host infectiousness and the role of the hypnozoite in perpetuating infection.

Thermal behavior of kaolinite–urea intercalation complex and molecular dynamics simulation for urea molecule orientation

The thermal behavior of kaolinite–urea intercalation complex was investigated by thermogravimetry–differential scanning calorimetry (TG–DSC), X-ray diffraction (XRD), and fourier transform infrared spectroscopy (FTIR). In addition, the interaction mode of urea molecules intercalated into the kaolinite gallery was studied by means of molecular dynamics simulation. Three main mass losses were observed at 136 °C, in the range of 210–270 °C, and at 500 °C in the TG–DSC curves, which were, respectively, attributed to (1) melting of the surface-adsorbed urea, (2) removal of the intercalated urea, and (3) dehydroxylation of the deintercalated kaolinite. The three DSC endothermic peaks at 218, 250, and 261 °C were related to the successive removals of intercalated urea with three different distribution structures. Based on the angle between the dipole moment vector of urea and the basal surface of kaolinite, the three urea models could be described as follows: (1) Type A, the dipole moment vector is nearly parallel to the basal surface of kaolinite; (2) Type B, the dipole moment vector points to the silica tetrahedron with the angle between it and the basal surface of kaolinite ranging from 20°to 40°; and (3) Type C, the dipole moment vector is nearly perpendicular to the basal surface of kaolinite. The three distribution structures of urea molecules were validated by the results of the molecular dynamics simulation. Furthermore, the thermal behavior of the kaolinite–urea intercalation complex investigated by TG–DSC was also supported by FTIR and XRD analyses.

Dynamics of gray matter loss in Alzheimer's disease

We detected and mapped a dynamically spreading wave of gray matter loss in the brains of patients with Alzheimer's disease (AD). The loss pattern was visualized in four dimensions as it spread over time from temporal and limbic cortices into frontal and occipital brain regions, sparing sensorimotor cortices. The shifting deficits were asymmetric (left hemisphere > right hemisphere) and correlated with progressively declining cognitive status (p < 0.0006). Novel brain mapping methods allowed us to visualize dynamic patterns of atrophy in 52 high-resolution magnetic resonance image scans of 12 patients with AD (age 68.4 ± 1.9 years) and 14 elderly matched controls (age 71.4 ± 0.9 years) scanned longitudinally (two scans; interscan interval 2.1 ± 0.4 years). A cortical pattern matching technique encoded changes in brain shape and tissue distribution across subjects and time. Cortical atrophy occurred in a well defined sequence as the disease progressed, mirroring the sequence of neurofibrillary tangle accumulation observed in cross sections at autopsy. Advancing deficits were visualized as dynamic maps that change over time. Frontal regions, spared early in the disease, showed pervasive deficits later (< 15% loss). The maps distinguished different phases of AD and differentiated AD from normal aging. Local gray matter loss rates (5.3 ± 2.3% per year in AD v 0.9 ± 0.9% per year in controls) were faster in the left hemisphere (p < 0.029) than the right. Transient barriers to disease progression appeared at limbic/frontal boundaries. This degenerative sequence, observed in vivo as it developed, provides the first quantitative, dynamic visualization of cortical atrophic rates in normal elderly populations and in those with dementia.

Spatiotemporal transmission dynamics of hemorrhagic fever with renal syndrome in China, 2005-2012

Background Hemorrhagic fever with renal syndrome (HFRS) is a rodent-borne disease caused by many serotypes of hantaviruses. In China, HFRS has been recognized as a severe public health problem with 90% of the total reported cases in the world. This study describes the spatiotemporal dynamics of HFRS cases in China and identifies the regions, time, and populations at highest risk, which could help the planning and implementation of key preventative measures. Methods Data on all reported HFRS cases at the county level from January 2005 to December 2012 were collected from Chinese Center for Disease Control and Prevention. Geographic Information System-based spatiotemporal analyses including Local Indicators of Spatial Association and Kulldorff's space-time scan statistic were performed to detect local high-risk space-time clusters of HFRS in China. In addition, cases from high-risk and low-risk counties were compared to identify significant demographic differences. Results A total of 100,868 cases were reported during 2005–2012 in mainland China. There were significant variations in the spatiotemporal dynamics of HFRS. HFRS cases occurred most frequently in June, November, and December. There was a significant positive spatial autocorrelation of HFRS incidence during the study periods, with Moran's I values ranging from 0.46 to 0.56 (P<0.05). Several distinct HFRS cluster areas were identified, mainly concentrated in northeastern, central, and eastern of China. Compared with cases from low-risk areas, a higher proportion of cases were younger, non-farmer, and floating residents in high-risk counties. Conclusions This study identified significant space-time clusters of HFRS in China during 2005–2012 indicating that preventative strategies for HFRS should be particularly focused on the northeastern, central, and eastern of China to achieve the most cost-effective outcomes.

Stochastic dynamics of interacting haematopoietic stem cell niche lineages

Since we still know very little about stem cells in their natural environment, it is useful to explore their dynamics through modelling and simulation, as well as experimentally. Most models of stem cell systems are based on deterministic differential equations that ignore the natural heterogeneity of stem cell populations. This is not appropriate at the level of individual cells and niches, when randomness is more likely to affect dynamics. In this paper, we introduce a fast stochastic method for simulating a metapopulation of stem cell niche lineages, that is, many sub-populations that together form a heterogeneous metapopulation, over time. By selecting the common limiting timestep, our method ensures that the entire metapopulation is simulated synchronously. This is important, as it allows us to introduce interactions between separate niche lineages, which would otherwise be impossible. We expand our method to enable the coupling of many lineages into niche groups, where differentiated cells are pooled within each niche group. Using this method, we explore the dynamics of the haematopoietic system from a demand control system perspective. We find that coupling together niche lineages allows the organism to regulate blood cell numbers as closely as possible to the homeostatic optimum. Furthermore, coupled lineages respond better than uncoupled ones to random perturbations, here the loss of some myeloid cells. This could imply that it is advantageous for an organism to connect together its niche lineages into groups. Our results suggest that a potential fruitful empirical direction will be to understand how stem cell descendants communicate with the niche and how cancer may arise as a result of a failure of such communication.

Fractal and complex network analyses of protein molecular dynamics

Based on protein molecular dynamics, we investigate the fractal properties of energy, pressure and volume time series using the multifractal detrended fluctuation analysis (MF-DFA) and the topological and fractal properties of their converted horizontal visibility graphs (HVGs). The energy parameters of protein dynamics we considered are bonded potential, angle potential, dihedral potential, improper potential, kinetic energy, Van der Waals potential, electrostatic potential, total energy and potential energy. The shape of the h(q)h(q) curves from MF-DFA indicates that these time series are multifractal. The numerical values of the exponent h(2)h(2) of MF-DFA show that the series of total energy and potential energy are non-stationary and anti-persistent; the other time series are stationary and persistent apart from series of pressure (with H≈0.5H≈0.5 indicating the absence of long-range correlation). The degree distributions of their converted HVGs show that these networks are exponential. The results of fractal analysis show that fractality exists in these converted HVGs. For each energy, pressure or volume parameter, it is found that the values of h(2)h(2) of MF-DFA on the time series, exponent λλ of the exponential degree distribution and fractal dimension dBdB of their converted HVGs do not change much for different proteins (indicating some universality). We also found that after taking average over all proteins, there is a linear relationship between 〈h(2)〉〈h(2)〉 (from MF-DFA on time series) and 〈dB〉〈dB〉 of the converted HVGs for different energy, pressure and volume.

Hot spot detection and spatio-temporal dynamics of dengue in Queensland, Australia

Dengue has been a major public health concern in Australia since it re-emerged in Queensland in 1992-1993. This study explored spatio-temporal distribution and clustering of locally-acquired dengue cases in Queensland State, Australia and identified target areas for effective interventions. A computerised locally-acquired dengue case dataset was collected from Queensland Health for Queensland from 1993 to 2012. Descriptive spatial and temporal analyses were conducted using geographic information system tools and geostatistical techniques. Dengue hot spots were detected using SatScan method. Descriptive spatial analysis showed that a total of 2,398 locally-acquired dengue cases were recorded in central and northern regions of tropical Queensland. A seasonal pattern was observed with most of the cases occurring in autumn. Spatial and temporal variation of dengue cases was observed in the geographic areas affected by dengue over time. Tropical areas are potential high-risk areas for mosquito-borne diseases such as dengue. This study demonstrated that the locally-acquired dengue cases have exhibited a spatial and temporal variation over the past twenty years in tropical Queensland, Australia. There is a clear evidence for the existence of statistically significant clusters of dengue and these clusters varied over time. These findings enabled us to detect and target dengue clusters suggesting that the use of geospatial information can assist the health authority in planning dengue control activities and it would allow for better design and implementation of dengue management programs.

Dynamics of a dilute sheared inelastic fluid. II. The effect of correlations

The effect of correlations on the viscosity of a dilute sheared inelastic fluid is analyzed using the ring-kinetic equation for the two-particle correlation function. The leading-order contribution to the stress in an expansion in epsilon=(1-e)(1/2) is calculated, and it is shown that the leading-order viscosity is identical to that obtained from the Green-Kubo formula, provided the stress autocorrelation function in a sheared steady state is used in the Green-Kubo formula. A systemmatic extension of this to higher orders is also formulated, and the higher-order contributions to the stress from the ring-kinetic equation are determined in terms of the terms in the Chapman-Enskog solution for the Boltzmann equation. The series is resummed analytically to obtain a renormalized stress equation. The most dominant contributions to the two-particle correlation function are products of the eigenvectors of the conserved hydrodynamic modes of the two correlated particles. In Part I, it was shown that the long-time tails of the velocity autocorrelation function are not present in a sheared fluid. Using those results, we show that correlations do not cause a divergence in the transport coefficients; the viscosity is not divergent in two dimensions, and the Burnett coefficients are not divergent in three dimensions. The equations for three-particle and higher correlations are analyzed diagrammatically. It is found that the contributions due to the three-particle and higher correlation functions to the renormalized viscosity are smaller than those due to the two-particle distribution function in the limit epsilon -> 0. This implies that the most dominant correlation effects are due to the two-particle correlations.