Predictive oncology: a review of multidisciplinary, multiscale in silico modeling linking phenotype, morphology and growth.

Empirical evidence and theoretical studies suggest that the phenotype, i.e., cellular- and molecular-scale dynamics, including proliferation rate and adhesiveness due to microenvironmental factors and gene expression that govern tumor growth and invasiveness, also determine gross tumor-scale morphology. It has been difficult to quantify the relative effect of these links on disease progression and prognosis using conventional clinical and experimental methods and observables. As a result, successful individualized treatment of highly malignant and invasive cancers, such as glioblastoma, via surgical resection and chemotherapy cannot be offered and outcomes are generally poor. What is needed is a deterministic, quantifiable method to enable understanding of the connections between phenotype and tumor morphology. Here, we critically assess advantages and disadvantages of recent computational modeling efforts (e.g., continuum, discrete, and cellular automata models) that have pursued this understanding. Based on this assessment, we review a multiscale, i.e., from the molecular to the gross tumor scale, mathematical and computational "first-principle" approach based on mass conservation and other physical laws, such as employed in reaction-diffusion systems. Model variables describe known characteristics of tumor behavior, and parameters and functional relationships across scales are informed from in vitro, in vivo and ex vivo biology. We review the feasibility of this methodology that, once coupled to tumor imaging and tumor biopsy or cell culture data, should enable prediction of tumor growth and therapy outcome through quantification of the relation between the underlying dynamics and morphological characteristics. In particular, morphologic stability analysis of this mathematical model reveals that tumor cell patterning at the tumor-host interface is regulated by cell proliferation, adhesion and other phenotypic characteristics: histopathology information of tumor boundary can be inputted to the mathematical model and used as a phenotype-diagnostic tool to predict collective and individual tumor cell invasion of surrounding tissue. This approach further provides a means to deterministically test effects of novel and hypothetical therapy strategies on tumor behavior.

Interplay between temperature, fish partial migration and trophic dynamics

Whereas many studies have addressed the mechanisms driving partial migration, few have focused on the consequences of partial migration on trophic dynamics, and integrated studies combining the two approaches are virtually nonexistent. Here we show that temperature affects seasonal partial migration of cyprinid fish from lakes to predation refuges in streams during winter and that this migration in combination with temperature affects the characteristics and phenology of lower trophic levels in the lake ecosystem. Specifically, our six-year study showed that the proportion of fish migrating was positively related to lake temperature during the pre-migration growth period, i.e. during summer. Migration from the lake occurred later when autumn water temperatures were high, and timing of return migration to the lake occurred earlier at higher spring water temperatures. Moreover, the winter mean size of zooplankton in the lake increased with the proportion of fish being away from the lake, likely as a consequence of decreased predation pressure. Peak biomass of phytoplankton in spring occurred earlier at higher spring water temperatures and with less fish being away from the lake. Accordingly, peak zooplankton biomass occurred earlier at higher spring water temperature, but relatively later if less fish were away from the lake. Hence, the time between phyto- and zooplankton peaks depended only on the amount of fish being away from the lake, and not on temperature. The intensity of fish migration thereby had a major effect on plankton spring dynamics. These results significantly contribute to our understanding of the interplay between partial migration and trophic dynamics, and suggest that ongoing climate change may significantly affect such dynamics.

A meta-analysis of cambium phenology and growth: linear and non-linear patterns in conifers of the northern hemisphere

Background and Aims Ongoing global warming has been implicated in shifting phenological patterns such as the timing and duration of the growing season across a wide variety of ecosystems. Linear models are routinely used to extrapolate these observed shifts in phenology into the future and to estimate changes in associated ecosystem properties such as net primary productivity. Yet, in nature, linear relationships may be special cases. Biological processes frequently follow more complex, non-linear patterns according to limiting factors that generate shifts and discontinuities, or contain thresholds beyond which responses change abruptly. This study investigates to what extent cambium phenology is associated with xylem growth and differentiation across conifer species of the northern hemisphere. Methods Xylem cell production is compared with the periods of cambial activity and cell differentiation assessed on a weekly time scale on histological sections of cambium and wood tissue collected from the stems of nine species in Canada and Europe over 1–9 years per site from 1998 to 2011. Key Results The dynamics of xylogenesis were surprisingly homogeneous among conifer species, although dispersions from the average were obviously observed. Within the range analysed, the relationships between the phenological timings were linear, with several slopes showing values close to or not statistically different from 1. The relationships between the phenological timings and cell production were distinctly non-linear, and involved an exponential pattern. Conclusions The trees adjust their phenological timings according to linear patterns. Thus, shifts of one phenological phase are associated with synchronous and comparable shifts of the successive phases. However, small increases in the duration of xylogenesis could correspond to a substantial increase in cell production. The findings suggest that the length of the growing season and the resulting amount of growth could respond differently to changes in environmental conditions.

Recent trends in Inner Asian forest dynamics to temperature and precipitation indicate high sensitivity to climate change

Semi-arid ecosystems play an important role in regulating global climate with the fate of these ecosystems in the Anthropocene depending upon interactions among temperature, precipitation, and CO2. However, in cool-arid environments, precipitation is not the only limitation to forest productivity. Interactions between changes in precipitation and air temperature may enhance soil moisture stress while simultaneously extending growing season length, with unclear consequences for net carbon uptake. This study evaluates recent trends in productivity and phenology of Inner Asian forests (in Mongolia and Northern China) using satellite remote sensing, dendrochronology, and dynamic global vegetation model (DGVM) simulations to quantify the sensitivity of forest dynamics to decadal climate variability and trends. Trends in photosynthetically active radiation fraction (FPAR) between 1982 and 2010 show a greening of about 7% of the region in spring (March, April, May), and 3% of the area ‘browning’ during summertime (June, July, August). These satellite observations of FPAR are corroborated by trends in NPP simulated by the LPJ DGVM. Spring greening trends in FPAR are mainly explained by long-term trends in precipitation whereas summer browning trends are correlated with decreasing precipitation. Tree ring data from 25 sites confirm annual growth increments are mainly limited by summer precipitation (June, July, August) in Mongolia, and spring precipitation in northern China (March, April, May), with relatively weak prior-year lag effects. An ensemble of climate projections from the IPCC CMIP3 models indicates that warming temperatures (spring, summer) are expected to be associated with higher summer precipitation, which combined with CO2 causes large increases in NPP and possibly even greater forest cover in the Mongolian steppe. In the absence of a strong direct CO2 fertilization effect on plant growth (e.g., due to nutrient limitation), water stress or decreased carbon gain from higher autotrophic respiration results in decreased productivity and loss of forest cover. The fate of these semi-arid ecosystems thus appears to hinge upon the magnitude and subtleties of CO2 fertilization effects, for which experimental observations in arid systems are needed to test and refine vegetation models.

Assessing dynamics of forced livestock movements, livelihoods and future development options for pastoralists/agro-pastoralists in Ruvuma and Lindi Regions, in the Southern Tanzania

The study that aimed at understanding the dynamics of forced livestock movements and pastoral livelihood and development options was conducted in Lindi and Ruvuma regions, using both formal and informal approaches. Data were collected from 60 randomly selected Agro-pastoralists/Pastoralists and native farmers using a structured questionnaire. Four villages were involved; two in Lindi region (Matandu and Mkwajuni) and the other two in Ruvuma region (Gumbiro and Muhuwesi). Data were analyzed using descriptive statistics of SPSS to generate means and frequencies. The results indicate that a large number of animals moved into the study area following the eviction order of the government in Ihefu wetlands in 2006/2007. Lindi region was earmarked by the government to receive all the evicted pastoralists. However, by 2008 only 30% of the total cattle that were expected to move into the region had been received. Deaths of many animals on transit, selling of the animals to pay for transportation and other costs while on transit and many pastoralists settling in Coastal and Ruvuma regions before reaching their destinations were reported to be the reasons for the discrepancy observed. To mitigate anticipated conflicts between farmers and pastoralists, Participatory Land Use Management (PLUM) plans were developed in all the study villages in order to demarcate village land area into different uses, including grazing, cropping, settlement and forests. Land units for grazing were supposed to be provided with all necessary livestock infrastructures (dips, charcoal dams, livestock markets and stock routes). However, the land use plans were not able to prevent the anticipated conflicts because most of the livestock infrastructures were lacking, the land use boundaries were not clearly demarcated and there was limited enforcement of village by-laws, since most had not been enacted by the respective district councils. Similarly, the areas allocated for grazing were inadequate for the number of livestock available and thus the carrying capacity exceeded. Thus, land resource-based conflicts between farmers and pastoralists were emerging in the study areas for the reason that most of the important components in the PLUM plans were not in place. Nevertheless, the arrival of pastoralists in the study areas had positive effects on food security and growth of social interactions between pastoralists and farmers including marriages between them. Environmental degradations due to the arrival of livestock were also not evident. Thus, there is a need for the government to purposely set aside enough grazing land with all necessary infrastructures in place for the agro-pastoral/pastoral communities in the country.

Herbivores differentially limit the seedling growth and sapling recruitment of two dominant rain forest trees

Resource heterogeneity may influence how plants are attacked and respond to consumers in multiple ways. Perhaps a better understanding of how this interaction might limit sapling recruitment in tree populations may be achieved by examining species’ functional responses to herbivores on a continuum of resource availability. Here, we experimentally reduced herbivore pressure on newly established seedlings of two dominant masting trees in 40 canopy gaps, across c. 80 ha of tropical rain forest in central Africa (Korup, Cameroon). Mesh cages were built to protect individual seedlings, and their leaf production and changes in height were followed for 22 months. With more light, herbivores increasingly prevented the less shade-tolerant Microberlinia bisulcata from growing as tall as it could and producing more leaves, indicating an undercompensation. The more shade-tolerant Tetraberlinia bifoliolata was much less affected by herbivores, showing instead near to full compensation for leaf numbers, and a negligible to weak impact of herbivores on its height growth. A stage-matrix model that compared control and caged populations lent evidence for a stronger impact of herbivores on the long-term population dynamics of M. bisulcata than T. bifoliolata. Our results suggest that insect herbivores can contribute to the local coexistence of two abundant tree species at Korup by disproportionately suppressing sapling recruitment of the faster-growing dominant via undercompensation across the light gradient created by canopy disturbances. The functional patterns we have documented here are consistent with current theory, and, because gap formations are integral to forest regeneration, they may be more widely applicable in other tropical forest communities. If so, the interaction between life-history and herbivore impact across light gradients may play a substantial role in tree species coexistence.

Climate sensitivity of Mediterranean pine growth reveals distinct east-west dipole

The European Mediterranean region is governed by a characteristic climate of summer drought that is likely to increase in duration and intensity under predicted climate change. However, large-scale network analyses investigating spatial aspects of pre-instrumental drought variability for this biogeographic zone are still scarce. In this study we introduce 54 mid- to high-elevation tree-ring width (TRW) chronologies comprising 2186 individual series from pine trees (Pinus spp.). This compilation spans a 4000-km east–west transect from Spain to Turkey, and was subjected to quality control and standardization prior to the development of site chronologies. A principal component analysis (PCA) was applied to identify spatial growth patterns during the network's common period 1862–1976, and new composite TRW chronologies were developed and investigated. The PCA reveals a common variance of 19.7% over the 54 Mediterranean pine chronologies. More interestingly, a dipole pattern in growth variability is found between the western (15% explained variance) and eastern (9.6%) sites, persisting back to 1330 AD. Pine growth on the Iberian Peninsula and Italy favours warm early growing seasons, but summer drought is most critical for ring width formation in the eastern Mediterranean region. Synoptic climate dynamics that have been in operation for the last seven centuries have been identified as the driving mechanism of a distinct east–west dipole in the growth variability of Mediterranean pines.

Maize Domestication and Anti-Herbivore Defences: Leaf-Specific Dynamics during Early Ontogeny of Maize and Its Wild Ancestors

As a consequence of artificial selection for specific traits, crop plants underwent considerable genotypic and phenotypic changes during the process of domestication. These changes may have led to reduced resistance in the cultivated plant due to shifts in resource allocation from defensive traits to increased growth rates and yield. Modern maize (Zea mays ssp. mays) was domesticated from its ancestor Balsas teosinte (Z. mays ssp. parviglumis) approximately 9000 years ago. Although maize displays a high genetic overlap with its direct ancestor and other annual teosintes, several studies show that maize and its ancestors differ in their resistance phenotypes with teosintes being less susceptible to herbivore damage. However, the underlying mechanisms are poorly understood. Here we addressed the question to what extent maize domestication has affected two crucial chemical and one physical defence traits and whether differences in their expression may explain the differences in herbivore resistance levels. The ontogenetic trajectories of 1,4-benzoxazin-3-ones, maysin and leaf toughness were monitored for different leaf types across several maize cultivars and teosinte accessions during early vegetative growth stages. We found significant quantitative and qualitative differences in 1,4-benzoxazin-3-one accumulation in an initial pairwise comparison, but we did not find consistent differences between wild and cultivated genotypes during a more thorough examination employing several cultivars/accessions. Yet, 1,4-benzoxazin-3-one levels tended to decline more rapidly with plant age in the modern maize cultivars. Foliar maysin levels and leaf toughness increased with plant age in a leaf-specific manner, but were also unaffected by domestication. Based on our findings we suggest that defence traits other than the ones that were investigated are responsible for the observed differences in herbivore resistance between teosinte and maize. Furthermore, our results indicate that single pairwise comparisons may lead to false conclusions regarding the effects of domestication on defensive and possibly other traits.

A novel approach for assessing shifting cultivation dynamics in a regional conservation hotspot - insights from nort-eastern Madagascar

The north-eastern escarpment of Madagascar has been deemed a global hotspot of biodiversity due to its high levels of endemic speciesbeing heavily threatened by accelerated deforestation rates and landscape changes. The main concern for conservation of the remaining humid primary forests is the shifting cultivation practices of local smallholder farmers for rice production. According to the mainstream narrative, human population growth leads to a shortening of crop-fallow cycles and thus to the accelerated conversion of forests to agricultural land. However, little is currently known about the dynamic changes between forest and shifting cultivation systems at the regional level. Existing land cover change analyses in this area have so far only focused on binary forest to non-forest changes and have therefore failed to account for the dynamic nature of the change processes between forest and different agriculture land use systems. This can be partly explained by the significant challenge to delineate shifting cultivation systems on land cover maps using traditional remote sensing classification approaches. To address this gap we therefore applied a novel GIS approach, that was originally developed for the assessment of shifting cultivation dynamics in Laos and has so far never been applied elsewhere, to map shifting cultivation of different crop-fallow lengths as well as permanent agriculture land use at the regional level. Change analyses of land use maps between 1995 and 2011 allowed us to comprehend the general trends of land use trajectories and their spatial variation. This more detailed understanding of land use change dynamics is key to plan for successful interventions to slow forest loss while at the same time improving local livelihoods. We further believe that this approach holds great potential for conservation monitoring in this resource-rich but povertyprone conservation hotspot.

Density-dependent dynamics of a dominant rain forest tree change with juvenile stage and time of masting

Although negative density dependence (NDD) can facilitate tree species coexistence in forests, the underlying mechanisms can differ, and rarely are the dynamics of seedlings and saplings studied together. Herein we present and discuss a novel mechanism based on our investigation of NDD predictions for the large, grove-forming ectomycorrhizal mast fruiting tree, Microberlinia bisulcata (Caesalpiniaceae), in an 82.5-ha plot at Korup, Cameroon. We tested whether juvenile density, size, growth and survival decreases with increasing conspecific adult basal area for 3245 ‘new’ seedlings and 540 ‘old’ seedlings (< 75-cm tall) during an approximately 4-year study period (2008–2012) and for 234 ‘saplings’ (≥ 75-cm tall) during an approximately 6-year study period (2008–2014). We found that the respective densities of new seedlings, old seedlings and saplings were positively, not and negatively related to increasing BA. Maximum leaf numbers and heights of old seedlings were negatively correlated with increasing basal areas, as were sapling heights and stem diameters. Whereas survivorship of new seedlings decreased by more than one-half with increasing basal area over its range in 2010–2012, that of old seedlings decreased by almost two-thirds, but only in 2008–2010, and was generally unrelated to conspecific seedling density. In 2010–2012 relative growth rates in new seedlings’ heights decreased with increasing basal area, as well as with increasing seedling density, together with increasing leaf numbers, whereas old seedlings’ growth was unrelated to either conspecific density or basal area. Saplings of below-average height had reduced survivorship with increasing basal area (probability decreasing from approx. 0.4 to 0.05 over the basal area range tested), but only sapling growth in terms of leaf numbers decreased with increasing basal area. These static and dynamic results indicate that NDD is operating within this system, possibly stabilizing the M. bisulcata population. However, these NDD patterns are unlikely to be caused by symmetric competition or by consumers. Instead, an alternative mechanism for conspecific adult–juvenile negative feedback is proposed, one which involves the interaction between tree phenology and ectomycorrhizal linkages.

Comparative carbon cycle dynamics of the present and last interglacial

Changes in temperature and carbon dioxide during glacial cycles recorded in Antarctic ice cores are tightly coupled. However, this relationship does not hold for interglacials. While climate cooled towards the end of both the last (Eemian) and present (Holocene) interglacials, CO₂ remained stable during the Eemian while rising in the Holocene. We identify and review twelve biogeochemical mechanisms of terrestrial (vegetation dynamics and CO₂ fertilization, land use, wild fire, accumulation of peat, changes in permafrost carbon, subaerial volcanic outgassing) and marine origin (changes in sea surface temperature, carbonate compensation to deglaciation and terrestrial biosphere regrowth, shallow-water carbonate sedimentation, changes in the soft tissue pump, and methane hydrates), which potentially may have contributed to the CO₂ dynamics during interglacials but which remain not well quantified. We use three Earth System Models (ESMs) of intermediate complexity to compare effects of selected mechanisms on the interglacial CO₂ and δ¹³ CO₂ changes, focusing on those with substantial potential impacts: namely carbonate sedimentation in shallow waters, peat growth, and (in the case of the Holocene) human land use. A set of specified carbon cycle forcings could qualitatively explain atmospheric CO₂ dynamics from 8ka BP to the pre-industrial. However, when applied to Eemian boundary conditions from 126 to 115 ka BP, the same set of forcings led to disagreement with the observed direction of CO₂ changes after 122 ka BP. This failure to simulate late-Eemian CO₂ dynamics could be a result of the imposed forcings such as prescribed CaCO₃ accumulation and/or an incorrect response of simulated terrestrial carbon to the surface cooling at the end of the interglacial. These experiments also reveal that key natural processes of interglacial CO₂ dynamics eshallow water CaCO₃ accumulation, peat and permafrost carbon dynamics are not well represented in the current ESMs. Global-scale modeling of these long-term carbon cycle components started only in the last decade, and uncertainty in parameterization of these mechanisms is a main limitation in the successful modeling of interglacial CO₂ dynamics.

Spatio-temporal co-ordination of RhoA, Rac1 and Cdc42 activation during prototypical edge protrusion and retraction dynamics

The three canonical Rho GTPases RhoA, Rac1 and Cdc42 co-ordinate cytoskeletal dynamics. Recent studies indicate that all three Rho GTPases are activated at the leading edge of motile fibroblasts, where their activity fluctuates at subminute time and micrometer length scales. Here, we use a microfluidic chip to acutely manipulate fibroblast edge dynamics by applying pulses of platelet-derived growth factor (PDGF) or the Rho kinase inhibitor Y-27632 (which lowers contractility). This induces acute and robust membrane protrusion and retraction events, that exhibit stereotyped cytoskeletal dynamics, allowing us to fairly compare specific morphodynamic states across experiments. Using a novel Cdc42, as well as previously described, second generation RhoA and Rac1 biosensors, we observe distinct spatio-temporal signaling programs that involve all three Rho GTPases, during protrusion/retraction edge dynamics. Our results suggest that Rac1, Cdc42 and RhoA regulate different cytoskeletal and adhesion processes to fine tune the highly plastic edge protrusion/retraction dynamics that power cell motility.

Frequency modulation of ERK activation dynamics rewires cell fate.

Transient versus sustained ERK MAP kinase (MAPK) activation dynamics induce proliferation versus differentiation in response to epidermal (EGF) or nerve (NGF) growth factors in PC-12 cells. Duration of ERK activation has therefore been proposed to specify cell fate decisions. Using a biosensor to measure ERK activation dynamics in single living cells reveals that sustained EGF/NGF application leads to a heterogeneous mix of transient and sustained ERK activation dynamics in distinct cells of the population, different than the population average. EGF biases toward transient, while NGF biases toward sustained ERK activation responses. In contrast, pulsed growth factor application can repeatedly and homogeneously trigger ERK activity transients across the cell population. These datasets enable mathematical modeling to reveal salient features inherent to the MAPK network. Ultimately, this predicts pulsed growth factor stimulation regimes that can bypass the typical feedback activation to rewire the system toward cell differentiation irrespective of growth factor identity.

Primary forest dynamics in lowland dipterocarp forest at Danum Valley, Sabah, Malaysia, and the role of the understorey

Changes in species composition in two 4–ha plots of lowland dipterocarp rainforest at Danum, Sabah, were measured over ten years (1986 to 1996) for trees greater than or equal to 10 cm girth at breast height (gbh). Each included a lower–slope to ridge gradient. The period lay between two drought events of moderate intensity but the forest showed no large lasting responses, suggesting that its species were well adapted to this regime. Mortality and recruitment rates were not unusual in global or regional comparisons. The forest continued to aggrade from its relatively (for Sabah) low basal area in 1986 and, together with the very open upper canopy structure and an abundance of lianas, this suggests a forest in a late stage of recovery from a major disturbance, yet one continually affected by smaller recent setbacks. Mortality and recruitment rates were not related to population size in 1986, but across subplots recruitment was positively correlated with the density and basal area of small trees (10 to <50 cm gbh) forming the dense understorey. Neither rate was related to topography. While species with larger mean gbh had greater relative growth rates (rgr) than smaller ones, subplot mean recruitment rates were correlated with rgr among small trees. Separating understorey species (typically the Euphorbiaceae) from the overstorey (Dipterocarpaceae) showed marked differences in change in mortality with increasing gbh: in the former it increased, in the latter it decreased. Forest processes are centred on this understorey quasi–stratum. The two replicate plots showed a high correspondence in the mortality, recruitment, population changes and growth rates of small trees for the 49 most abundant species in common to both. Overstorey species had higher rgrs than understorey ones, but both showed considerable ranges in mortality and recruitment rates. The supposed trade–off in traits, viz slower rgr, shade tolerance and lower population turnover in the understorey group versus faster potential growth rate, high light responsiveness and high turnover in the overstorey group, was only partly met, as some understorey species were also very dynamic. The forest at Danum, under such a disturbance–recovery regime, can be viewed as having a dynamic equilibrium in functional and structural terms. A second trade–off in shade–tolerance versus drought–tolerance is suggested for among the understorey species. A two–storey (or vertical component) model is proposed where the understorey–overstorey species’ ratio of small stems (currently 2:1) is maintained by a major feedback process. The understorey appears to be an important part of this forest, giving resilience against drought and protecting the overstorey saplings in the long term. This view could be valuable for understanding forest responses to climate change where drought frequency in Borneo is predicted to intensify in the coming decades.

Job Matching and Wage Growth in the U.S. and Germany

This paper examines the contribution of job matching to wage growth in the U.S. and Germany using data drawn from the Panel Study of Income Dynamics and the German Socio-Economic Panel from 1984 through 1992. Using a symmetrical set of variables and data handling procedures, real wage growth is found to be higher in the U.S. than in Germany during this period. Also, using two different estimators, job matches are found to enhance wage growth in the U.S. and retard it in Germany. The relationship of general skills to employment in each country appears responsible for this result.