Emerging infectious diseases and animal social systems

Emerging infectious diseases threaten a wide diversity of animals, and important questions remain concerning disease emergence in socially structured populations. We developed a spatially explicit simulation model to investigate whether—and under what conditions—disease-related mortality can impact rates of pathogen spread in populations of polygynous groups. Specifically, we investigated whether pathogen-mediated dispersal (PMD) can occur when females disperse after the resident male dies from disease, thus carrying infections to new groups. We also examined the effects of incubation period and virulence, host mortality and rates of background dispersal, and we used the model to investigate the spread of the virus responsible for Ebola hemorrhagic fever, which currently is devastating African ape populations. Output was analyzed using regression trees, which enable exploration of hierarchical and non-linear relationships. Analyses revealed that the incidence of disease in single-male (polygynous) groups was significantly greater for those groups containing an average of more than six females, while the total number of infected hosts in the population was most sensitive to the number of females per group. Thus, as expected, PMD occurs in polygynous groups and its effects increase as harem size (the number of females) increases. Simulation output further indicated that population-level effects of Ebola are likely to differ among multi-male–multi-female chimpanzees and polygynous gorillas, with larger overall numbers of chimpanzees infected, but more gorilla groups becoming infected due to increased dispersal when the resident male dies. Collectively, our results highlight the importance of social system on the spread of disease in wild mammals.     

Kinetochore life histories reveal an Aurora-B-dependent error correction mechanism in anaphase

1. Download : Download high-res image (280KB)
2. Download : Download full-size image

     

Simulation of Y-chromosomal haplotype data

The non-recombining nature of the Y-chromosome determines the non-independence of alleles between loci. The evolution of short tandem repeat (STR) loci in the Y-chromosome is the result of different factors such as differential mutation rates, mutation modes, gene conversion, selection and demographic processes. The degree of correlation between loci is dependent on the magnitude of these processes. The simulation of data is a routine tool used for testing hypotheses in population and evolutionary studies. The most basic parameters hitherto used in lineage haplotype simulations are the allele frequency distributions and mutation rates, assuming either full independence or linkage between loci. In this study we introduce use of the Spearman correlation coefficient to estimate the degree of dependence between non-recombining loci. Then, both the interdependence between loci and the allele frequency distributions at multi-allelic loci are incorporated in an algorithm for simulating haplotypes. We illustrate the method using published and unpublished Y-chromosome STR data.     

Conserved water-mediated H-bonding dynamics of catalytic Asn 175 in plant thiol protease

The role of invariant water molecules in the activity of plant cysteine protease is ubiquitous in nature. On analysing the 11 different Protein DataBank (PDB) structures of plant thiol proteases, the two invariant water molecules W1 and W2 (W220 and W222 in the template 1PPN structure) were observed to form H-bonds with the O_b atom of Asn 175. Extensive energy minimization and molecular dynamics simulation studies up to 2 ns on all the PDB and solvated structures clearly revealed the involvement of the H-bonding association of the two water molecules in fixing the orientation of the asparagine residue of the catalytic triad. From this study, it is suggested that H-bonding of the water molecule at the W1 invariant site better stabilizes the Asn residue at the active site of the catalytic triad.     

A Molecular Dynamics Investigation of Chain Conformational Changes In Compressed Bilayers

Abstract

The conformations of the chains constituting the hydrophilic component of alkyl monolayers and bilayers are investigated by performing molecular dynamics atomistic simulations on these systems at different temperatures. Several monitoring techniques are used to reveal the chain conformations, including atom pair radial distribution functions, evolutions of the torsional angles over thousands of timesteps, frequency distributions of the torsionl angles and ‘snapshot’ plots of the atomic configurations. These methods consistently testify to the stability of the trans (fully extended) character of the strain-free alkyl chains up to room temperature. The chains retain much of this conformation even when the layers are compressed by the application of pressure, to which the chains respond by ‘folding’ at the ends attaching them to the substrate planes while maintaining directions which are mainly normal to these planes. A non-zero gap between the layers is also maintained. A pressure of about 50 kbar abruptly causes all motion in the chains to cease, resulting in a highly ordered lattice structure.     

A genetic and spatial Bayesian analysis of mastitis resistance

A nationwide health card recording system for dairy cattle was introduced in Norway in 1975 (the Norwegian Cattle Health Services). The data base holds information on mastitis occurrences on an individual cow basis. A reduction in mastitis frequency across the population is desired, and for this purpose risk factors are investigated. In this paper a Bayesian proportional hazards model is used for modelling the time to first veterinary treatment of clinical mastitis, including both genetic and environmental covariates. Sire effects were modelled as shared random components, and veterinary district was included as an environmental effect with prior spatial smoothing. A non-informative smoothing prior was assumed for the baseline hazard, and Markov chain Monte Carlo methods (MCMC) were used for inference. We propose a new measure of quality for sires, in terms of their posterior probability of being among the, say 10% best sires. The probability is an easily interpretable measure that can be directly used to rank sires. Estimating these complex probabilities is straightforward in an MCMC setting. The results indicate considerable differences between sires with regards to their daughters disease resistance. A regional effect was also discovered with the lowest risk of disease in the south-eastern parts of Norway.     

Novel Monte Carlo dose calculation algorithm for robotic radiosurgery with multi leaf collimator: Dosimetric evaluation

PurposeAt introduction in 2014, dose calculation for the first MLC on a robotic SRS/SBRT platform was limited to a correction-based Finite-Size Pencil Beam (FSPB) algorithm. We report on the dosimetric accuracy of a novel Monte Carlo (MC) dose calculation algorithm for this MLC, included in the Precision™ treatment planning system.MethodsA phantom was built of one slab (5.0 cm) of lung-equivalent material (0.09…0.29 g/cc) enclosed by 3.5 cm (above) and 5 cm (below) slabs of solid water (1.045 g/cc). This was irradiated using rectangular (15.4 × 15.4 mm² to 53.8 × 53.7 mm²) and two irregular MLC-fields. Radiochromic film (EBT3) was positioned perpendicular to the slabs and parallel to the beam. Calculated dose distributions were compared to film measurements using line scans and 2D gamma analysis.ResultsMeasured and MC calculated percent depth dose curves showed a characteristic dose drop within the low-density region, which was not correctly reproduced by FSPB. Superior average gamma pass rates (2%/1 mm) were found for MC (91.2 ± 1.5%) compared to FSPB (55.4 ± 2.7%). However, MC calculations exhibited localized anomalies at mass density transitions around 0.15 g/cc, which were traced to a simplification in electron transport. Absence of these anomalies was confirmed in a modified build of the MC engine, which increased gamma pass rates to 96.6 ± 1.2%.ConclusionsThe novel MC algorithm greatly improves dosimetric accuracy in heterogeneous tissue, potentially expanding the clinical use of robotic radiosurgery with MLC. In-depth, independent validation is paramount to identify and reduce the residual uncertainties in any software solution.     

Computational Oral Absorption Simulation for Low‐Solubility Compounds

Bile micelles play an important role in oral absorption of low‐solubility compounds. Bile micelles can affect solubility, dissolution rate, and permeability. For the pH–solubility profile in bile micelles, the Henderson–Hasselbalch equation should be modified to take bile‐micelle partition into account. For the dissolution rate, in the Nernst–Brunner equation, the effective diffusion coefficient in bile‐micelle media should be used instead of the monomer diffusion coefficient. The diffusion coefficient of bile micelles is 8‐ to 18‐fold smaller than that of monomer molecules. For permeability, the effective diffusion coefficient in the unstirred water layer adjacent to the epithelial membrane, and the free fraction at the epithelial membrane surface should be taken into account. The importance of these aspects is demonstrated here using several in vivo and clinical oral‐absorption data of low‐solubility model compounds. Using the theoretical equations, the food effect on oral absorption is further discussed.     

Application of a photosystem II model for analysis of fluorescence induction curves in the 100 ns to 10 s time domain after excitation with a saturating light pulse

A mathematical model of photosystem II (PSII) events was used to analyze chlorophyll fluorescence transients in the time domain from 100 ns to 10 s after excitation with a saturating 10-ns flash, applied as a part of specialized illumination protocol, using preparations of a thermophilic strain of the unicellular green alga, Chlorella pyrenoidosa Chick (using both intact and diuron-treated cells). Analysis of simulation results has proven that particular attention should be given to flash-induced recombination processes, including nonradiative recombination in PSII, while subsequent charge transfer along the electron transport chain of thylakoid membrane can be adequately described by a single reaction of quinone reoxidation. The PSII model was extended by taking inhibition by diuron of the electron transport in the acceptor side of PSII into account, which allowed simulation of fluorescence induction curves observed in the presence of this inhibitor. The model parameters were determined (stromal pH, rate constants of nonradiative recombination, and the initial reduction state of the quinone pool) which provided adequate simulation of experimentally observed ratios of the maximal and initial fluorescence levels (F _m/F ₀).