Comparing methods for modelling spreading cell fronts

Spreading cell fronts play an essential role in many physiological processes. Classically, models of this process are based on the Fisher-Kolmogorov equation; however, such continuum representations are not always suitable as they do not explicitly represent behaviour at the level of individual cells. Additionally, many models examine only the large time asymptotic behaviour, where a travelling wave front with a constant speed has been established. Many experiments, such as a scratch assay, never display this asymptotic behaviour, and in these cases the transient behaviour must be taken into account. We examine the transient and asymptotic behaviour of moving cell fronts using techniques that go beyond the continuum approximation via a volume-excluding birth-migration process on a regular one-dimensional lattice. We approximate the averaged discrete results using three methods: (i) mean-field, (ii) pair-wise, and (iii) one-hole approximations. We discuss the performace of these methods, in comparison to the averaged discrete results, for a range of parameter space, examining both the transient and asymptotic behaviours. The one-hole approximation, based on techniques from statistical physics, is not capable of predicting transient behaviour but provides excellent agreement with the asymptotic behaviour of the averaged discrete results, provided that cells are proliferating fast enough relative to their rate of migration. The mean-field and pair-wise approximations give indistinguishable asymptotic results, which agree with the averaged discrete results when cells are migrating much more rapidly than they are proliferating. The pair-wise approximation performs better in the transient region than does the mean-field, despite having the same asymptotic behaviour. Our results show that each approximation only works in specific situations, thus we must be careful to use a suitable approximation for a given system, otherwise inaccurate predictions could be made.

Distinguishing between mechanisms of cell aggregation using pair-correlation functions

Many cell types form clumps or aggregates when cultured in vitro through a variety of mechanisms including rapid cell proliferation, chemotaxis, or direct cell-to-cell contact. In this paper we develop an agent-based model to explore the formation of aggregates in cultures where cells are initially distributed uniformly, at random, on a two-dimensional substrate. Our model includes unbiased random cell motion, together with two mechanisms which can produce cell aggregates: (i) rapid cell proliferation, and (ii) a biased cell motility mechanism where cells can sense other cells within a finite range, and will tend to move towards areas with higher numbers of cells. We then introduce a pair-correlation function which allows us to quantify aspects of the spatial patterns produced by our agent-based model. In particular, these pair-correlation functions are able to detect differences between domains populated uniformly at random (i.e. at the exclusion complete spatial randomness (ECSR) state) and those where the proliferation and biased motion rules have been employed - even when such differences are not obvious to the naked eye. The pair-correlation function can also detect the emergence of a characteristic inter-aggregate distance which occurs when the biased motion mechanism is dominant, and is not observed when cell proliferation is the main mechanism of aggregate formation. This suggests that applying the pair-correlation function to experimental images of cell aggregates may provide information about the mechanism associated with observed aggregates. As a proof of concept, we perform such analysis for images of cancer cell aggregates, which are known to be associated with rapid proliferation. The results of our analysis are consistent with the predictions of the proliferation-based simulations, which supports the potential usefulness of pair correlation functions for providing insight into the mechanisms of aggregate formation.

Interpreting scratch assays using pair density dynamics and approximate Bayesian computation

Quantifying the impact of biochemical compounds on collective cell spreading is an essential element of drug design, with various applications including developing treatments for chronic wounds and cancer. Scratch assays are a technically simple and inexpensive method used to study collective cell spreading; however, most previous interpretations of scratch assays are qualitative and do not provide estimates of the cell diffusivity, D, or the cell proliferation rate,l. Estimating D and l is important for investigating the efficacy of a potential treatment and provides insight into the mechanism through which the potential treatment acts. While a few methods for estimating D and l have been proposed, these previous methods lead to point estimates of D and l, and provide no insight into the uncertainty in these estimates. Here, we compare various types of information that can be extracted from images of a scratch assay, and quantify D and l using discrete computational simulations and approximate Bayesian computation. We show that it is possible to robustly recover estimates of D and l from synthetic data, as well as a new set of experimental data. For the first time, our approach also provides a method to estimate the uncertainty in our estimates of D and l. We anticipate that our approach can be generalized to deal with more realistic experimental scenarios in which we are interested in estimating D and l, as well as additional relevant parameters such as the strength of cell-to-cell adhesion or the strength of cell-to-substrate adhesion.

Expression of mammalian glutathione S-transferase 5-5 in Salmonella typhimurium TA1535 leads to base-pair mutations upon exposure to dihalomethanes

Dihalomethanes can produce liver tumors in mice but not in rats, and concern exists about the risk of these compounds to humans. Glutathione (GSH) conjugation of dihalomethanes has been considered to be a critical event in the bioactivation process, and risk assessment is based upon this premise; however, there is little experimental support for this view or information about the basis of genotoxicity. A plasmid vector containing rat GSH S-transferase 5-5 was transfected into the Salmonella typhimurium tester strain TA1535, which then produced active enzyme. The transfected bacteria produced base-pair revertants in the presence of ethylene dihalides or dihalomethanes, in the order CH2Br2 > CH2BrCl > CH2Cl2. However, revertants were not seen when cells were exposed to GSH, CH2Br2, and an amount of purified GSH S-transferase 5-5 (20-fold excess in amount of that expressed within the cells). HCHO, which is an end product of the reaction of GSH with dihalomethanes, also did not produce mutations. S-(1-Acetoxymethyl)GSH was prepared as an analog of the putative S-(1-halomethyl)GSH reactive intermediates. This analog did not produce revertants, consistent with the view that activation of dihalomethanes must occur within the bacteria to cause genetic damage, presenting a model to be considered in studies with mammalian cells. S-(1-Acetoxymethyl)GSH reacted with 2′-deoxyguanosine to yield a major adduct, identified as S-[1-(N2-deoxyguanosinyl)methyl]GSH. Demonstration of the activation of dihalomethanes by this mammalian GSH S-transferase theta class enzyme should be of use in evaluating the risk of these chemicals, particularly in light of reports of the polymorphic expression of a similar activity in humans.

Disability-adjusted life years (DALYs) for 291 diseases and injuries in 21 regions, 1990-2010: A systematic analysis for the Global Burden of Disease Study 2010

BACKGROUND Measuring disease and injury burden in populations requires a composite metric that captures both premature mortality and the prevalence and severity of ill-health. The 1990 Global Burden of Disease study proposed disability-adjusted life years (DALYs) to measure disease burden. No comprehensive update of disease burden worldwide incorporating a systematic reassessment of disease and injury-specific epidemiology has been done since the 1990 study. We aimed to calculate disease burden worldwide and for 21 regions for 1990, 2005, and 2010 with methods to enable meaningful comparisons over time. METHODS We calculated DALYs as the sum of years of life lost (YLLs) and years lived with disability (YLDs). DALYs were calculated for 291 causes, 20 age groups, both sexes, and for 187 countries, and aggregated to regional and global estimates of disease burden for three points in time with strictly comparable definitions and methods. YLLs were calculated from age-sex-country-time-specific estimates of mortality by cause, with death by standardised lost life expectancy at each age. YLDs were calculated as prevalence of 1160 disabling sequelae, by age, sex, and cause, and weighted by new disability weights for each health state. Neither YLLs nor YLDs were age-weighted or discounted. Uncertainty around cause-specific DALYs was calculated incorporating uncertainty in levels of all-cause mortality, cause-specific mortality, prevalence, and disability weights. FINDINGS Global DALYs remained stable from 1990 (2·503 billion) to 2010 (2·490 billion). Crude DALYs per 1000 decreased by 23% (472 per 1000 to 361 per 1000). An important shift has occurred in DALY composition with the contribution of deaths and disability among children (younger than 5 years of age) declining from 41% of global DALYs in 1990 to 25% in 2010. YLLs typically account for about half of disease burden in more developed regions (high-income Asia Pacific, western Europe, high-income North America, and Australasia), rising to over 80% of DALYs in sub-Saharan Africa. In 1990, 47% of DALYs worldwide were from communicable, maternal, neonatal, and nutritional disorders, 43% from non-communicable diseases, and 10% from injuries. By 2010, this had shifted to 35%, 54%, and 11%, respectively. Ischaemic heart disease was the leading cause of DALYs worldwide in 2010 (up from fourth rank in 1990, increasing by 29%), followed by lower respiratory infections (top rank in 1990; 44% decline in DALYs), stroke (fifth in 1990; 19% increase), diarrhoeal diseases (second in 1990; 51% decrease), and HIV/AIDS (33rd in 1990; 351% increase). Major depressive disorder increased from 15th to 11th rank (37% increase) and road injury from 12th to 10th rank (34% increase). Substantial heterogeneity exists in rankings of leading causes of disease burden among regions. INTERPRETATION Global disease burden has continued to shift away from communicable to non-communicable diseases and from premature death to years lived with disability. In sub-Saharan Africa, however, many communicable, maternal, neonatal, and nutritional disorders remain the dominant causes of disease burden. The rising burden from mental and behavioural disorders, musculoskeletal disorders, and diabetes will impose new challenges on health systems. Regional heterogeneity highlights the importance of understanding local burden of disease and setting goals and targets for the post-2015 agenda taking such patterns into account. Because of improved definitions, methods, and data, these results for 1990 and 2010 supersede all previously published Global Burden of Disease results.

Psychological health and the risk of diabetes mellitus in Australian women: A 21-year prospective study

Background Symptoms of depression can be recurrent or limited to one episode. This study discusses the prospective association between psychological health, measured as change in depression symptoms, and the risk of diabetes mellitus in Australian women. Methods Data obtained from the Mater-University of Queensland Study of Pregnancy. Depression was measured using the Delusions-Symptoms: States Inventory. To examine possible transitions over time, depression was grouped into four categories and assessed at different phases over the 21-year period. Multiple logistic regression models and sensitivity analysis to assess the robustness of our analytical strategy were performed. Results Three hundred and one women reported diabetes 21 years after the index pregnancy. Almost one-third of the women who reported depression symptoms continued to report these at a subsequent follow-up (FU) phase. About 1 in 20 women who had not reported depression symptoms at the 5-year FU did so at the subsequent 14-year FU. In prospective analyses, we did not find a significant association between diabetes and negative change (not depressed to depressed, at subsequent phase); however, for women with positive history of symptoms of depression and women with persistent symptoms, there was a 1.97-fold (95% confidence interval [CI]: 1.14–3.40) to 2.23-fold (95% CI: 1.09–4.57) greater risk of diabetes. Conclusions Our study suggests that an increased risk of diabetes is significantly associated with persistent depression symptoms. It highlights the importance of recognizing depression symptoms in terms of women's psychological wellbeing and thus provides a basis for targeting those most at risk.

A comparative risk assessment of burden of disease and injury attributable to 67 risk factors and risk factor clusters in 21 regions, 1990–2010: A systematic analysis for the Global Burden of Disease Study 2010

BACKGROUND Quantification of the disease burden caused by different risks informs prevention by providing an account of health loss different to that provided by a disease-by-disease analysis. No complete revision of global disease burden caused by risk factors has been done since a comparative risk assessment in 2000, and no previous analysis has assessed changes in burden attributable to risk factors over time. METHODS We estimated deaths and disability-adjusted life years (DALYs; sum of years lived with disability [YLD] and years of life lost [YLL]) attributable to the independent effects of 67 risk factors and clusters of risk factors for 21 regions in 1990 and 2010. We estimated exposure distributions for each year, region, sex, and age group, and relative risks per unit of exposure by systematically reviewing and synthesising published and unpublished data. We used these estimates, together with estimates of cause-specific deaths and DALYs from the Global Burden of Disease Study 2010, to calculate the burden attributable to each risk factor exposure compared with the theoretical-minimum-risk exposure. We incorporated uncertainty in disease burden, relative risks, and exposures into our estimates of attributable burden. FINDINGS In 2010, the three leading risk factors for global disease burden were high blood pressure (7·0% [95% uncertainty interval 6·2-7·7] of global DALYs), tobacco smoking including second-hand smoke (6·3% [5·5-7·0]), and alcohol use (5·5% [5·0-5·9]). In 1990, the leading risks were childhood underweight (7·9% [6·8-9·4]), household air pollution from solid fuels (HAP; 7·0% [5·6-8·3]), and tobacco smoking including second-hand smoke (6·1% [5·4-6·8]). Dietary risk factors and physical inactivity collectively accounted for 10·0% (95% UI 9·2-10·8) of global DALYs in 2010, with the most prominent dietary risks being diets low in fruits and those high in sodium. Several risks that primarily affect childhood communicable diseases, including unimproved water and sanitation and childhood micronutrient deficiencies, fell in rank between 1990 and 2010, with unimproved water and sanitation accounting for 0·9% (0·4-1·6) of global DALYs in 2010. However, in most of sub-Saharan Africa childhood underweight, HAP, and non-exclusive and discontinued breastfeeding were the leading risks in 2010, while HAP was the leading risk in south Asia. The leading risk factor in Eastern Europe, most of Latin America, and southern sub-Saharan Africa in 2010 was alcohol use; in most of Asia, North Africa and Middle East, and central Europe it was high blood pressure. Despite declines, tobacco smoking including second-hand smoke remained the leading risk in high-income north America and western Europe. High body-mass index has increased globally and it is the leading risk in Australasia and southern Latin America, and also ranks high in other high-income regions, North Africa and Middle East, and Oceania. INTERPRETATION Worldwide, the contribution of different risk factors to disease burden has changed substantially, with a shift away from risks for communicable diseases in children towards those for non-communicable diseases in adults. These changes are related to the ageing population, decreased mortality among children younger than 5 years, changes in cause-of-death composition, and changes in risk factor exposures. New evidence has led to changes in the magnitude of key risks including unimproved water and sanitation, vitamin A and zinc deficiencies, and ambient particulate matter pollution. The extent to which the epidemiological shift has occurred and what the leading risks currently are varies greatly across regions. In much of sub-Saharan Africa, the leading risks are still those associated with poverty and those that affect children.

Genomic architecture of ecologically divergent body shape in a pair of sympatric crater lake cichlid fishes

Determining the genetic bases of adaptations and their roles in speciation is a prominent issue in evolutionary biology. Cichlid fish species flocks are a prime example of recent rapid radiations, often associated with adaptive phenotypic divergence from a common ancestor within a short period of time. In several radiations of freshwater fishes, divergence in ecomorphological traits - including body shape, colour, lips and jaws - is thought to underlie their ecological differentiation, specialization and, ultimately, speciation. The Midas cichlid species complex (Amphilophus spp.) of Nicaragua provides one of the few known examples of sympatric speciation where species have rapidly evolved different but parallel morphologies in young crater lakes. This study identified significant QTL for body shape using SNPs generated via ddRAD sequencing and geometric morphometric analyses of a cross between two ecologically and morphologically divergent, sympatric cichlid species endemic to crater Lake Apoyo: an elongated limnetic species (Amphilophus zaliosus) and a high-bodied benthic species (Amphilophus astorquii). A total of 453 genome-wide informative SNPs were identified in 240 F-2 hybrids. These markers were used to construct a genetic map in which 25 linkage groups were resolved. Seventy-two segregating SNPs were linked to 11 QTL. By annotating the two most highly supported QTL-linked genomic regions, genes that might contribute to divergence in body shape along the benthic-limnetic axis in Midas cichlid sympatric adaptive radiations were identified. These results suggest that few genomic regions of large effect contribute to early stage divergence in Midas cichlids.

Genome-wide association study for radiographic vertebral fractures: A potential role for the 16q24 BMD locus

Vertebral fracture risk is a heritable complex trait. The aim of this study was to identify genetic susceptibility factors for osteoporotic vertebral fractures applying a genome-wide association study (GWAS) approach. The GWAS discovery was based on the Rotterdam Study, a population-based study of elderly Dutch individuals aged >55years; and comprising 329 cases and 2666 controls with radiographic scoring (McCloskey-Kanis) and genetic data. Replication of one top-associated SNP was pursued by de-novo genotyping of 15 independent studies across Europe, the United States, and Australia and one Asian study. Radiographic vertebral fracture assessment was performed using McCloskey-Kanis or Genant semi-quantitative definitions. SNPs were analyzed in relation to vertebral fracture using logistic regression models corrected for age and sex. Fixed effects inverse variance and Han-Eskin alternative random effects meta-analyses were applied. Genome-wide significance was set at p<5×10-8. In the discovery, a SNP (rs11645938) on chromosome 16q24 was associated with the risk for vertebral fractures at p=4.6×10-8. However, the association was not significant across 5720 cases and 21,791 controls from 14 studies. Fixed-effects meta-analysis summary estimate was 1.06 (95% CI: 0.98-1.14; p=0.17), displaying high degree of heterogeneity (I2=57%; Qhet p=0.0006). Under Han-Eskin alternative random effects model the summary effect was significant (p=0.0005). The SNP maps to a region previously found associated with lumbar spine bone mineral density (LS-BMD) in two large meta-analyses from the GEFOS consortium. A false positive association in the GWAS discovery cannot be excluded, yet, the low-powered setting of the discovery and replication settings (appropriate to identify risk effect size >1.25) may still be consistent with an effect size <1.10, more of the type expected in complex traits. Larger effort in studies with standardized phenotype definitions is needed to confirm or reject the involvement of this locus on the risk for vertebral fractures.

Mimicry of 21. 5 KDa myelin basic protein peptide by a Maedi Visna virus polymerase peptide does not contribute to the pathogenesis of encephalitis in sheep

Epitope mimicry is the theory that an infectious agent such as a virus causes pathological effects via mimicry of host proteins and thus elicits a cross-reactive immune response to host tissues. Weise and Carnegie (1988) found a region of sequence similarity between the pol gene of the Maedi Visna virus (MVV), which induces demyelinating encephalitis in sheep, and myelin basic protein (MBP), which is known to induce experimental allergic encephalitis (EAE) in laboratory animals. In this study, cross-reactions between sera raised in sheep against synthetic peptides of MVV (TGKIPWILLPGR) and 21.5 kDa MBP (SGKVPWLKRPGR) were demonstrated using enzyme-linked immunosorbant assay (ELISA) and thin layer chromatography (TLC) immunoprobing. The antibody responses of MVV-infected sheep were investigated using ELISA against the peptides, and MBP protein, immunoprobing of the peptides on TPC plates and Western blotting against MBP. Slight significant reactions to the 21.5 kDa MBP peptide (P < 0.001) and to a lesser extent sheep MBP (P < 0.004) were detected in ELISA. The MBP peptide evoked stronger responses from more sera than the MVV peptide on immunoprobed TLC plates. On the Western blots, eight of the 23 sheep with Visna had serum reactivity to MBP. This slight reaction to MBP in MVV-infected sheep is of interest because of the immune responses to MBP evident in multiple sclerosis and EAE, but its relevance in Visna is limited since no correlation with disease severity was observed. The cell-mediated immune responses of MVV-infected sheep against similar peptides was assessed. The peptides did not stimulate proliferation of peripheral blood lymphocytes of MVV-infected sheep. Since the MVV peptide was not recognised by antibodies or T lymphocytes from MVV-infected and encephalic sheep, it was concluded that epitope mimicry of this 21.5 kDa MBP peptide by the similar MVV pol peptide was not contributing to the immunopathogenesis of Visna. The slight antibody response to MBP and the MBP peptide can be attributed to by-stander effects of the immunopathology of MVV-induced encephalitis.

Common variants at the MHC locus and at chromosome 16q24.1 predispose to Barrett's Esophagus

Barrett's esophagus is an increasingly common disease that is strongly associated with reflux of stomach acid and usually a hiatus hernia, and it strongly predisposes to esophageal adenocarcinoma (EAC), a tumor with a very poor prognosis. We report the first genome-wide association study on Barrett's esophagus, comprising 1,852 UK cases and 5,172 UK controls in the discovery stage and 5,986 cases and 12,825 controls in the replication stage. Variants at two loci were associated with disease risk: chromosome 6p21, rs9257809 (P combined = 4.09 × 10-9; odds ratio (OR) = 1.21, 95% confidence interval (CI) =1.13-1.28), within the major histocompatibility complex locus, and chromosome 16q24, rs9936833 (P combined = 2.74 × 10-10; OR = 1.14, 95% CI = 1.10-1.19), for which the closest protein-coding gene is FOXF1, which is implicated in esophageal development and structure. We found evidence that many common variants of small effect contribute to genetic susceptibility to Barrett's esophagus and that SNP alleles predisposing to obesity also increase risk for Barrett's esophagus. © 2012 Nature America, Inc. All rights reserved.

Single 2×1 domain orientation on Si(001) surfaces using aperiodic Bi line arrays

Studies of Bi heteroepitaxy on Si(001) have shown that lines grow to lengths of up to 500nm if the substrate is heated to above the Bi desorption temperature (500°C) during or after Bi deposition. Unlike many other nanoline systems, the lines formed by this nonequilibrium growth process have no detectable width dispersion. Although much attention has been given to the atomic geometery of the line, in this paper, we focus on how the lines can be used to create a majority 2×1 domain orientation. It is demonstrated that the Bi lines can be used to produce a single-domain orientation on Si(001) if the lines are grown on Si(001) surfaces with a regular distribution of single height steps. This is a compelling example of how a nanoscale motif can be used to modify mesoscopic surface structure on Si(001).

pi-conjugation and conformation in a semiconducting polymer: small angle x-ray scattering study

Small angle x-ray scattering (SAXS) in a poly[2-methoxy-5-(2'-ethyl-hexyloxy)-1,4-phenylene vinylene] (MEH-PPV) solution has shown the important role of pi-electron conjugation in controlling the chain conformation and assembly. By increasing the extent of conjugation from 30 to 100%, the persistence length (l(p)) increases from 20 to 66 angstrom. Moreover, a pronounced second peak in the pair distribution function has been observed in a fully conjugated chain, at larger length scales. This feature indicates that the chain segments tend to self-assemble as the conjugation along the chain increases. Xylene enhances the rigidity of the PPV backbone to yield extended structures, while tetrahydrofuran solvates the side groups to form compact coils in which the lp is much shorter.

Practical Algorithms For Mean Velocity Estimation In Pulse Doppler Weather Radars Using A Small Number Of Samples

Doppler weather radars with fast scanning rates must estimate spectral moments based on a small number of echo samples. This paper concerns the estimation of mean Doppler velocity in a coherent radar using a short complex time series. Specific results are presented based on 16 samples. A wide range of signal-to-noise ratios are considered, and attention is given to ease of implementation. It is shown that FFT estimators fare poorly in low SNR and/or high spectrum-width situations. Several variants of a vector pulse-pair processor are postulated and an algorithm is developed for the resolution of phase angle ambiguity. This processor is found to be better than conventional processors at very low SNR values. A feasible approximation to the maximum entropy estimator is derived as well as a technique utilizing the maximization of the periodogram. It is found that a vector pulse-pair processor operating with four lags for clear air observation and a single lag (pulse-pair mode) for storm observation may be a good way to estimate Doppler velocities over the entire gamut of weather phenomena.

Reliable stabilization using a multi-controller configuration

Given a plant P, we consider the problem of designing a pair of controllers C1 and C2 such that their sum stabilizes P, and in addition, each of them also stabilizes P should the other one fail. This is referred to as the reliable stabilization problem. It is shown that every strongly stabilizable plant can be reliably stabilized; moreover, one of the two controllers can be specified arbitrarily, subject only to the constraint that it should be stable. The stabilization technique is extended to reliable regulation.