Bayesian calibration of a stochastic,multiscale agent-based model for predicting in vitro tumor growth

Hybrid multiscale agent-based models (ABMs) are unique in their ability to simulate individual cell interactions and microenvironmental dynamics. Unfortunately, the high computational cost of modeling individual cells, the inherent stochasticity of cell dynamics, and numerous model parameters are fundamental limitations of applying such models to predict tumor dynamics. To overcome these challenges, we have developed a coarse-grained two-scale ABM (cgABM) with a reduced parameter space that allows for an accurate and efficient calibration using a set of time-resolved microscopy measurements of cancer cells grown with different initial conditions. The multiscale model consists of a reaction-diffusion type model capturing the spatio-temporal evolution of glucose and growth factors in the tumor microenvironment (at tissue scale), coupled with a lattice-free ABM to simulate individual cell dynamics (at cellular scale). The experimental data consists of BT474 human breast carcinoma cells initialized with different glucose concentrations and tumor cell confluences. The confluence of live and dead cells was measured every three hours over four days. Given this model, we perform a time-dependent global sensitivity analysis to identify the relative importance of the model parameters. The subsequent cgABM is calibrated within a Bayesian framework to the experimental data to estimate model parameters, which are then used to predict the temporal evolution of the living and dead cell populations. To this end, a moment-based Bayesian inference is proposed to account for the stochasticity of the cgABM while quantifying uncertainties due to limited temporal observational data. The cgABM reduces the computational time of ABM simulations by 93% to 97% while staying within a 3% difference in prediction compared to ABM. Additionally, the cgABM can reliably predict the temporal evolution of breast cancer cells observed by the microscopy data with an average error and standard deviation for live and dead cells being 7.61±2.01 and 5.78±1.13, respectively.     

Revisiting the particular role of host shifts in initiating insect speciation

The notion that shifts to new hosts can initiate insect speciation is more than 150 years old, yet widespread conflation with paradigms of sympatric speciation has led to confusion about how much support exists for this hypothesis. Here, we review 85 insect systems and evaluate the relationship between host shifting, reproductive isolation, and speciation. We sort insects into five categories: (1) systems in which a host shift has initiated speciation; (2) systems in which a host shift has made a contribution to speciation; (3) systems in which a host shift has caused the evolution of new reproductive isolating barriers; (4) systems with host‐associated genetic differences; and (5) systems with no evidence of host‐associated genetic differences. We find host‐associated genetic structure in 65 systems, 43 of which show that host shifts have resulted in the evolution of new reproductive barriers. Twenty‐six of the latter also support a role for host shifts in speciation, including eight studies that definitively support the hypothesis that a host shift has initiated speciation. While this review is agnostic as to the fraction of all insect speciation events to which host shifts have contributed, it clarifies that host shifts absolutely can and do initiate speciation.     

No difference in 1RM strength and muscle activation during the barbell chest press on a stable and unstable surface

Exercise or Swiss balls are increasingly being used with conventional resistance exercises. There is little evidence supporting the efficacy of this approach compared to traditional resistance training on a stable surface. Previous studies have shown that force output may be reduced with no change in muscle electromyography (EMG) activity while others have shown increased muscle EMG activity when performing resistance exercises on an unstable surface. This study compared 1RM strength, and upper body and trunk muscle EMG activity during the barbell chest press exercise on a stable (flat bench) and unstable surface (exercise ball). After familiarization, 13 subjects underwent testing for 1RM strength for the barbell chest press on both a stable bench and an exercise ball, each separated by at least 7 days. Surface EMG was recorded for 5 upper body muscles and one trunk muscle from which average root mean square of the muscle activity was calculated for the whole 1RM lift and the concentric and eccentric phases. Elbow angle during each lift was recorded to examine any range-of-motion differences between the two surfaces. The results show that there was no difference in 1RM strength or muscle EMG activity for the stable and unstable surfaces. In addition, there was no difference in elbow range-of-motion between the two surfaces. Taken together, these results indicate that there is no reduction in 1RM strength or any differences in muscle EMG activity for the barbell chest press exercise on an unstable exercise ball when compared to a stable flat surface. Moreover, these results do not support the notion that resistance exercises performed on an exercise ball are more efficacious than traditional stable exercises.     

Forging Links between Human Mental Retardation–Associated CNVs and Mouse Gene Knockout Models

Rare copy number variants (CNVs) are frequently associated with common neurological disorders such as mental retardation (MR; learning disability), autism, and schizophrenia. CNV screening in clinical practice is limited because pathological CNVs cannot be distinguished routinely from benign CNVs, and because genes underlying patients'' phenotypes remain largely unknown. Here, we present a novel, statistically robust approach that forges links between 148 MR–associated CNVs and phenotypes from ∼5,000 mouse gene knockout experiments. These CNVs were found to be significantly enriched in two classes of genes, those whose mouse orthologues, when disrupted, result in either abnormal axon or dopaminergic neuron morphologies. Additional enrichments highlighted correspondences between relevant mouse phenotypes and secondary presentations such as brain abnormality, cleft palate, and seizures. The strength of these phenotype enrichments (>100% increases) greatly exceeded molecular annotations (<30% increases) and allowed the identification of 78 genes that may contribute to MR and associated phenotypes. This study is the first to demonstrate how the power of mouse knockout data can be systematically exploited to better understand genetically heterogeneous neurological disorders.     

Identification of microRNAs associated with hyperthermia-induced cellular stress response

MicroRNAs (miRNAs) are a class of small RNAs that play a critical role in the coordination of fundamental cellular processes. Recent studies suggest that miRNAs participate in the cellular stress response (CSR), but their specific involvement remains unclear. In this study, we identify a group of thermally regulated miRNAs (TRMs) that are associated with the CSR. Using miRNA microarrays, we show that dermal fibroblasts differentially express 123 miRNAs when exposed to hyperthermia. Interestingly, only 27 of these miRNAs are annotated in the current Sanger registry. We validated the expression of the annotated miRNAs using qPCR techniques, and we found that the qPCR and microarray data was in well agreement. Computational target-prediction studies revealed that putative targets for the TRMs are heat shock proteins and Argonaute-2—the core functional unit of RNA silencing. These results indicate that cells express a specific group of miRNAs when exposed to hyperthermia, and these miRNAs may function in the regulation of the CSR. Future studies will be conducted to determine if other cells lines differentially express these miRNAs when exposed to hyperthermia.     

Infection of CD127+ (interleukin-7 receptor+) CD4+ cells and overexpression of CTLA-4 are linked to loss of antigen-specific CD4 T cells during primary human immunodeficiency virus type 1 infection

We recently found that human immunodeficiency virus (HIV)-specific CD4+ T cells express coreceptor CCR5 and activation antigen CD38 during early primary HIV-1 infection (PHI) but then rapidly disappear from the circulation. This cell loss may be due to susceptibility to infection with HIV-1 but could also be due to inappropriate apoptosis, an expansion of T regulatory cells, trafficking out of the circulation, or dysfunction. We purified CD38+++CD4+ T cells from peripheral blood mononuclear cells, measured their level of HIV-1 DNA by PCR, and found that about 10% of this population was infected. However, a small subset of HIV-specific CD4+) T cells also expressed CD127, a marker of long-term memory cells. Purified CD127+CD4+ lymphocytes contained fivefold more copies of HIV-1 DNA per cell than did CD127-negative CD4+ cells, suggesting preferential infection of long-term memory cells. We observed no apoptosis of antigen-specific CD4+ T cells in vitro and only a small increase in CD45RO+CD25+CD127dimCD4+ T regulatory cells during PHI. However, 40% of CCR5+CD38+++ CD4+ T cells expressed gut-homing integrins, suggesting trafficking through gut-associated lymphoid tissue (GALT). Furthermore, 80% of HIV-specific CD4+ T cells expressed high levels of the negative regulator CTLA-4 in response to antigen stimulation in vitro, which was probably contributing to their inability to produce interleukin-2 and proliferate. Taken together, the loss of HIV-specific CD4+ T cells is associated with a combination of an infection of CCR5+ CD127+ memory CD4+ T cells, possibly in GALT, and a high expression of the inhibitory receptor CTLA-4.     

The good, the bad and the recovery in an assisted migration

Background

Assisted migration or translocation of species to ameliorate effects of habitat loss or changing environment is currently under scrutiny as a conservation tool. A large scale experiment of assisted migration over hundreds of kilometres was tested on a morph from a commercial fishery of southern rock lobster Jasus edwardsii, to enhance depleted populations, improve the yield and sustainability of the fishery, and test resilience to a changing climate.

Methodology and Principal Findings

Approximately 10,000 lower-valued, pale-coloured lobsters were moved from deep water to inshore sites (2 in Tasmania [TAS] and 2 in South Australia [SA]) where the high-value, red morph occurs. In TAS this was a northwards movement of 1° latitude. Growth was measured only in TAS lobsters, and reproductive status was recorded in lobsters from all locations. Pale females (TAS) grew 4 times faster than resident pale lobsters from the original site and twice as fast as red lobsters at their new location. Approximately 30% of translocated pale lobsters deferred reproduction for one year after release (SA and TAS), and grew around 1 mm yr⁻¹ less compared to translocated pale lobsters that did not defer reproduction. In spite of this stress response to translocation, females that deferred reproduction still grew 2–6 mm yr⁻¹ more than lobsters at the source site. Lobsters have isometric growth whereby volume increases as a cube of length. Consequently despite the one-year hiatus in reproduction, increased growth increases fecundity of translocated lobsters, as the increase in size provided a larger volume for producing and incubating eggs in future years.

Conclusions and Significance

Assisted migration improved egg production and growth, despite a temporary stress response, and offers a tool to improve the production, sustainability and resilience of the fishery.     

Cell pairing and methylation in Tetrahymena thermophila are altered by exogenous homocysteine

Homocysteine is causally associated with birth defects such as spina bifida, and with premature vascular disease. We have investigated the effects of homocysteine on a cell-cell interaction in a fundamental eukaryotic system, the free-living ciliate Tetrahymena. Exogenously added homocysteine inhibits cell pairing in a dose-dependent manner. These effects are exacerbated by adenosine, which by itself has little demonstrable influence on pairing. S-adenosylhomocysteine (SAH) is a product of the reaction between adenosine and homocysteine, and is an inhibitor of methyl transferases. We therefore predicted that protein methylation would be significantly inhibited by homocysteine. A direct test of that hypothesis involved a demonstration that incorporation of an isotopically labeled methyl group from methionine into proteins was significantly reduced by homocysteine. The undermethylated proteins are of low molecular weight, and might correspond to known methylatable signaling proteins. We show that vanadate, an inhibitor of protein phosphatase, also inhibits cell pairing, and that the effects of vanadate and homocysteine are additive. This is the first demonstration that methylation and possibly phosphorylation play a regulatory role in cell-cell interactions in ciliates.