Statistical evaluation of alternative models of human evolution

An appropriate model of recent human evolution is not only important to understand our own history, but it is necessary to disentangle the effects of demography and selection on genome diversity. Although most genetic data support the view that our species originated recently in Africa, it is still unclear if it completely replaced former members of the Homo genus, or if some interbreeding occurred during its range expansion. Several scenarios of modern human evolution have been proposed on the basis of molecular and paleontological data, but their likelihood has never been statistically assessed. Using DNA data from 50 nuclear loci sequenced in African, Asian and Native American samples, we show here by extensive simulations that a simple African replacement model with exponential growth has a higher probability (78%) as compared with alternative multiregional evolution or assimilation scenarios. A Bayesian analysis of the data under this best supported model points to an origin of our species approximately 141 thousand years ago (Kya), an exit out-of-Africa approximately 51 Kya, and a recent colonization of the Americas approximately 10.5 Kya. We also find that the African replacement model explains not only the shallow ancestry of mtDNA or Y-chromosomes but also the occurrence of deep lineages at some autosomal loci, which has been formerly interpreted as a sign of interbreeding with Homo erectus.

Elevated levels of placental growth factor represent an adaptive host response in sepsis

Recently, we demonstrated that circulating levels of vascular endothelial growth factor (VEGF) and placental growth factor (PlGF) are increased in sepsis (Yano, K., P.C. Liaw, J.M. Mullington, S.C. Shih, H. Okada, N. Bodyak, P.M. Kang, L. Toltl, B. Belikoff, J. Buras, et al. 2006. J. Exp. Med. 203:1447-1458). Moreover, enhanced VEGF/Flk-1 signaling was shown to contribute to sepsis morbidity and mortality. We tested the hypothesis that PlGF also contributes to sepsis outcome. In mouse models of endotoxemia and cecal ligation puncture, the genetic absence of PlGF or the systemic administration of neutralizing anti-PlGF antibodies resulted in higher mortality compared with wild-type or immunoglobulin G-injected controls, respectively. The increased mortality associated with genetic deficiency of PlGF was reversed by adenovirus (Ad)-mediated overexpression of PlGF. In the endotoxemia model, PlGF deficiency was associated with elevated circulating levels of VEGF, induction of VEGF expression in the liver, impaired cardiac function, and organ-specific accentuation of barrier dysfunction and inflammation. Mortality of endotoxemic PlGF-deficient mice was increased by Ad-mediated overexpression of VEGF and was blocked by expression of soluble Flt-1. Collectively, these data suggest that up-regulation of PlGF in sepsis is an adaptive host response that exerts its benefit, at least in part, by attenuating VEGF signaling.

Erupted magma volume estimates at Santiaguito and Pacaya Volcanoes, Guatemala using digital elevation models

High resolution digital elevation models (DEMs) of Santiaguito and Pacaya volcanoes, Guatemala, were used to estimate volume changes and eruption rates between 1954 and 2001. The DEMs were generated from contour maps and aerial photography, which were analyzed in ArcGIS 9.0®. Because both volcanoes were growing substantially over the five decade period, they provide a good data set for exploring effective methodology for estimating volume changes. The analysis shows that the Santiaguito dome complex grew by 0.78 ± 0.07 km3 (0.52 ± 0.05 m3 s-1) over the 1954-2001 period with nearly all the growth occurring on the El Brujo (1958-75) and Caliente domes (1971-2001). Adding information from field data prior to 1954, the total volume extruded from Santiaguito since 1922 is estimated at 1.48 ± 0.19 km3. Santiaguito’s growth rate is lower than most other volcanic domes, but it has been sustained over a much longer period and has undergone a change toward more exogenous and progressively slower extrusion with time. At Santiaguito some of the material being added at the dome is subsequently transported downstream by block and ash flows, mudflows and floods, creating channel shifting and areas of aggradation and erosion. At Pacaya volcano a total volume of 0.21 ± 0.05 km3 was erupted between 1961 and 2001 for an average extrusion rate of 0.17 ± 0.04 m3 s-1. Both the Santiaguito and Pacaya eruption rate estimates reported here are minima, because they do not include estimates of materials which are transported downslope after eruption and data on ashfall which may result in significant volumes of material spread over broad areas. Regular analysis of high resolution DEMs using the methods outlined here, would help quantify the effects of fluvial changes to downstream populated areas, as well as assist in tracking hazards related to dome collapse and eruption.

Models with Plants, Microorganisms and Viruses for Basic Research in Homeopathy

Most criticism about homeopathy concerns the lack of a scientific basis and theoretical models. In order to be accepted as a valid part of medical practice, a wellstructured research strategy for homeopathy is needed. This is often hampered by methodological problems as well as by gross underinvestment in the required academic resources. Fundamental research could make important contributions to our understanding of the homeopathic and high dilutions mechanisms of action. Since the pioneering works of Kolisko on wheat germination (Kolisko, 1923) and Junker on growth of microorganisms (paramecium, yeast, fungi) (Junker, 1928), a number of experiments have been performed either with healthy organisms (various physiological aspects of growth) or with artificially diseased organisms, which may react more markedly to homeopathic treatments than healthy ones. In the latter case, the preliminary stress may be either abiotic, e.g. heavy metals, or biotic, e.g. fungal and viral pathogens or nematode infection. Research has also been carried out into the applicability of homeopathic principles to crop growth and disease control (agrohomeopathy): because of the extreme dilutions used, the environmental impact is low and such treatments are well suited to the holistic approach of sustainable agriculture (Betti et al., 2006). Unfortunately, as Scofield reported in an extensive critical review (Scofield, 1984), there is little firm evidence to support the reliability of the reported results, due to poor experimental methodology and inadequate statistical analysis. Moreover, since there is no agricultural homeopathic pharmacopoeia, much work is required to find suitable remedies, potencies and dose levels.

The importance of growth hormone in the regulation of erythropoiesis, red cell mass, and plasma volume in adults with growth hormone deficiency

Total body water (TBW) is reduced in adult GH deficiency (GHD) largely due to a reduction of extracellular water. It is unknown whether total blood volume (TBV) contributes to the reduced extracellular water in GHD. GH and insulin-like growth factor I (IGF-I) have been demonstrated to stimulate erythropoiesis in vitro, in animal models, and in growing children. Whether GH has a regulatory effect on red cell mass (RCM) in adults is not known. We analyzed body composition by bioelectrical impedance and used standard radionuclide dilution methods to measure RCM and plasma volume (PV) along with measuring full blood count, ferritin, vitamin B12, red cell folate, IGF-I, IGF-binding protein-3, and erythropoietin in 13 adult patients with GHD as part of a 3-month, double blind, placebo-controlled trial of GH (0.036 U/kg.day). TBW and lean body mass significantly increased by 2.5 +/- 0.53 kg (mean +/- SEM; P < 0.004) and 3.4 +/- 0.73 kg (P < 0.004), respectively, and fat mass significantly decreased by 2.4 +/- 0.32 kg (P < 0.001) in the GH-treated group. The baseline RCM of all patients with GHD was lower than the predicted normal values (1635 +/- 108 vs. 1850 +/- 104 mL; P < 0.002). GH significantly increased RCM, PV, and TBV by 183 +/- 43 (P < 0.006), 350 +/- 117 (P < 0.03), and 515 +/- 109 (P < 0.004) mL, respectively. The red cell count increased by 0.36 +/- 0.116 x 10(12)/L (P < 0.03) with a decrease in ferritin levels by 39.1 +/- 4.84 micrograms/L (P < 0.001) after GH treatment. Serum IGF-I and IGF-binding protein-3 concentrations increased by 3.0 +/- 0.43 (P < 0.001) and 1.3 +/- 0.15 (P < 0.001) SD, respectively, but the erythropoietin concentration was unchanged after GH treatment. No significant changes in body composition or blood volume were recorded in the placebo group. Significant positive correlations could be established between changes in TBW and TBV, lean body mass and TBV (r = 0.78; P < 0.04 and r = 0.77; P < 0.04, respectively), and a significant negative correlation existed between changes in fat mass and changes in TBV in the GH-treated group (r = -0.95; P < 0.02). We conclude that 1) erythropoiesis is impaired in GHD; 2) GH stimulates erythropoiesis in adult GHD; and 3) GH increases PV and TBV, which may contribute to the increased exercise performance seen in these patients.

Climate-induced interannual variability of marine primary and export production in three global coupled climate carbon cycle models

Fully coupled climate carbon cycle models are sophisticated tools that are used to predict future climate change and its impact on the land and ocean carbon cycles. These models should be able to adequately represent natural variability, requiring model validation by observations. The present study focuses on the ocean carbon cycle component, in particular the spatial and temporal variability in net primary productivity (PP) and export production (EP) of particulate organic carbon (POC). Results from three coupled climate carbon cycle models (IPSL, MPIM, NCAR) are compared with observation-based estimates derived from satellite measurements of ocean colour and results from inverse modelling (data assimilation). Satellite observations of ocean colour have shown that temporal variability of PP on the global scale is largely dominated by the permanently stratified, low-latitude ocean (Behrenfeld et al., 2006) with stronger stratification (higher sea surface temperature; SST) being associated with negative PP anomalies. Results from all three coupled models confirm the role of the low-latitude, permanently stratified ocean for anomalies in globally integrated PP, but only one model (IPSL) also reproduces the inverse relationship between stratification (SST) and PP. An adequate representation of iron and macronutrient co-limitation of phytoplankton growth in the tropical ocean has shown to be the crucial mechanism determining the capability of the models to reproduce observed interactions between climate and PP.

Graft remodeling during growth following anterior cruciate ligament reconstruction in skeletally immature sheep

INTRODUCTION: Ruptures of the anterior cruciate ligament are being diagnosed with increasing frequency in skeletally immature individuals. It was our aim to investigate the graft remodelling process following an autologous, transphyseal reconstruction of the anterior cruciate ligament (ACL) in skeletally immature sheep. We hypothesized that the ligamentisation process in immature sheep is quicker and more complete when compared to adult sheep. MATERIALS AND METHODS: Skeletally immature sheep with an age of 4 months underwent a fully transphyseal ACL reconstruction using an autologous tendon. The animals were subsequently sacrificed at 3, 6, 12 and 24 weeks following surgery. Each group was characterised histomorphometrically, by immunostaining (VEGF, SMA), by transmission electron microscopy (TEM) and biomechanically (UFS Roboter). RESULTS: The histomorphometric analysis and presence of VEGF and SMA positive cells demonstrated a rapid return to a ligament like structure. The biomechanical analysis revealed an anteroposterior translation that was still increased even 6 months following surgery. CONCLUSION: As in adult sheep models, the remodeling of a soft tissue graft used for ACL reconstruction results in a biomechanically inferior substitute. However, the immature tissue seems to remodel faster and more complete when compared to adults.

Planet formation models: the interplay with the planetesimal disc

Context. According to the sequential accretion model (or core-nucleated accretion model), giant planet formation is based first on the formation of a solid core which, when massive enough, can gravitationally bind gas from the nebula to form the envelope. The most critical part of the model is the formation time of the core: to trigger the accretion of gas, the core has to grow up to several Earth masses before the gas component of the protoplanetary disc dissipates. Aims: We calculate planetary formation models including a detailed description of the dynamics of the planetesimal disc, taking into account both gas drag and excitation of forming planets. Methods: We computed the formation of planets, considering the oligarchic regime for the growth of the solid core. Embryos growing in the disc stir their neighbour planetesimals, exciting their relative velocities, which makes accretion more difficult. Here we introduce a more realistic treatment for the evolution of planetesimals' relative velocities, which directly impact on the formation timescale. For this, we computed the excitation state of planetesimals, as a result of stirring by forming planets, and gas-solid interactions. Results: We find that the formation of giant planets is favoured by the accretion of small planetesimals, as their random velocities are more easily damped by the gas drag of the nebula. Moreover, the capture radius of a protoplanet with a (tiny) envelope is also larger for small planetesimals. However, planets migrate as a result of disc-planet angular momentum exchange, with important consequences for their survival: due to the slow growth of a protoplanet in the oligarchic regime, rapid inward type I migration has important implications on intermediate-mass planets that have not yet started their runaway accretion phase of gas. Most of these planets are lost in the central star. Surviving planets have masses either below 10 M⊕ or above several Jupiter masses. Conclusions: To form giant planets before the dissipation of the disc, small planetesimals (~0.1 km) have to be the major contributors of the solid accretion process. However, the combination of oligarchic growth and fast inward migration leads to the absence of intermediate-mass planets. Other processes must therefore be at work to explain the population of extrasolar planets that are presently known.

Antitumor effect of SIRT1 inhibition in human HCC tumor models in vitro and in vivo

Sirtuins (SIRT1-7) are a highly conserved family of NAD(+)-dependent enzymes that control the activity of histone and nonhistone regulatory proteins. SIRT1 is purposed to promote longevity and to suppress the initiation of some cancers. Nevertheless, SIRT1 is reported to function as a tumor suppressor as well as an oncogenic protein. Our data show that compared with normal liver or surrounding tumor tissue, SIRT1 is strongly overexpressed in human hepatocellular carcinoma (HCC). In addition, human HCC cell lines (Hep3B, HepG2, HuH7, HLE, HLF, HepKK1, skHep1) were screened for the expression of the sirtuin family members and only SIRT1 was consistently overexpressed compared with normal hepatocytes. To determine its effect on HCC growth, SIRT1 activity was inhibited either with lentiviruses expressing short hairpin RNAs or with the small molecule inhibitor, cambinol. Knockdown or inhibition of SIRT1 activity had a cytostatic effect, characterized by an altered morphology, impaired proliferation, an increased expression of differentiation markers, and cellular senescence. In an orthotopic xenograft model, knockdown of SIRT1 resulted in 50% fewer animals developing tumors and cambinol treatment resulted in an overall lower tumor burden. Taken together, our data show that inhibition of SIRT1 in HCC cells impairs their proliferation in vitro and tumor formation in vivo. These data suggest that SIRT1 expression positively influences the growth of HCC and support further studies aimed to block its activity alone or in combination as a novel treatment strategy.

Cardiomyocyte growth and sarcomerogenesis at the intercalated disc

Cardiomyocytes grow during heart maturation or disease-related cardiac remodeling. We present evidence that the intercalated disc (ID) is integral to both longitudinal and lateral growth: increases in width are accommodated by lateral extension of the plicate tread regions and increases in length by sarcomere insertion within the ID. At the margin between myofibril and the folded membrane of the ID lies a transitional junction through which the thin filaments from the last sarcomere run to the ID membrane and it has been suggested that this junction acts as a proto Z-disc for sarcomere addition. In support of this hypothesis, we have investigated the ultrastructure of the ID in mouse hearts from control and dilated cardiomyopathy (DCM) models, the MLP-null and a cardiac-specific β-catenin mutant, cΔex3, as well as in human left ventricle from normal and DCM samples. We find that the ID amplitude can vary tenfold from 0.2 μm up to a maximum of ~2 μm allowing gradual expansion during heart growth. At the greatest amplitude, equivalent to a sarcomere length, A-bands and thick filaments are found within the ID membrane loops together with a Z-disc, which develops at the transitional junction position. Here, also, the tops of the membrane folds, which are rich in αII spectrin, become enlarged and associated with junctional sarcoplasmic reticulum. Systematically larger ID amplitudes are found in DCM samples. Other morphological differences between mouse DCM and normal hearts suggest that sarcomere inclusion is compromised in the diseased hearts.

Computer simulation of glioma growth and morphology.

Despite major advances in the study of glioma, the quantitative links between intra-tumor molecular/cellular properties, clinically observable properties such as morphology, and critical tumor behaviors such as growth and invasiveness remain unclear, hampering more effective coupling of tumor physical characteristics with implications for prognosis and therapy. Although molecular biology, histopathology, and radiological imaging are employed in this endeavor, studies are severely challenged by the multitude of different physical scales involved in tumor growth, i.e., from molecular nanoscale to cell microscale and finally to tissue centimeter scale. Consequently, it is often difficult to determine the underlying dynamics across dimensions. New techniques are needed to tackle these issues. Here, we address this multi-scalar problem by employing a novel predictive three-dimensional mathematical and computational model based on first-principle equations (conservation laws of physics) that describe mathematically the diffusion of cell substrates and other processes determining tumor mass growth and invasion. The model uses conserved variables to represent known determinants of glioma behavior, e.g., cell density and oxygen concentration, as well as biological functional relationships and parameters linking phenomena at different scales whose specific forms and values are hypothesized and calculated based on in vitro and in vivo experiments and from histopathology of tissue specimens from human gliomas. This model enables correlation of glioma morphology to tumor growth by quantifying interdependence of tumor mass on the microenvironment (e.g., hypoxia, tissue disruption) and on the cellular phenotypes (e.g., mitosis and apoptosis rates, cell adhesion strength). Once functional relationships between variables and associated parameter values have been informed, e.g., from histopathology or intra-operative analysis, this model can be used for disease diagnosis/prognosis, hypothesis testing, and to guide surgery and therapy. In particular, this tool identifies and quantifies the effects of vascularization and other cell-scale glioma morphological characteristics as predictors of tumor-scale growth and invasion.

Delta like ligand 4 is a critical regulator of bone marrow cell differentiation into pericytes/vascular smooth muscle cells and is essential for the vasculogenesis that supports the growth of Ewing’s sarcoma

We have previously shown that vasculogenesis, the process by which bone marrow-derived cells are recruited to the tumor and organized to form a blood vessel network de novo, is essential for the growth of Ewing’s sarcoma. We further demonstrated that these bone marrow cells differentiate into pericytes/vascular smooth muscle cells(vSMC) and contribute to the formation of the functional vascular network. The molecular mechanisms that control bone marrow cell differentiation into pericytes/vSMC in Ewing’s sarcoma are poorly understood. Here, we demonstrate that the Notch ligand Delta like ligand 4 (DLL4) plays a critical role in this process. DLL4 is essential for the formation of mature blood vessels during development and in several tumor models. Inhibition of DLL4 causes increased vascular sprouting, decreased pericyte coverage, and decreased vessel functionality. We demonstrate for the first time that DLL4 is expressed by bone marrow-derived pericytes/vascular smooth muscle cells in two Ewing’s sarcoma xenograft models and by perivascular cells in 12 out of 14 patient samples. Using dominant negative mastermind to inhibit Notch, we demonstrate that Notch signaling is essential for bone marrow cell participation in vasculogenesis. Further, inhibition of DLL4 using either shRNA or the monoclonal DLL4 neutralizing antibody YW152F led to dramatic changes in blood vessel morphology and function. Vessels in tumors where DLL4 was inhibited were smaller, lacked lumens, had significantly reduced numbers of bone marrow-derived pericyte/vascular smooth muscle cells, and were less functional. Importantly, growth of TC71 and A4573 tumors was significantly inhibited by treatment with YW152F. Additionally, we provide in vitro evidence that DLL4-Notch signaling is involved in bone marrow-derived pericyte/vascular smooth muscle cell formation outside of the Ewing’s sarcoma environment. Pericyte/vascular smooth muscle cell marker expression by whole bone marrow cells cultured with mouse embryonic stromal cells was reduced when DLL4 was inhibited by YW152F. For the first time, our findings demonstrate a role for DLL4 in bone marrow-derived pericyte/vascular smooth muscle differentiation as well as a critical role for DLL4 in Ewing’s sarcoma tumor growth.

Computer simulation of glioma growth and morphology.

Despite major advances in the study of glioma, the quantitative links between intra-tumor molecular/cellular properties, clinically observable properties such as morphology, and critical tumor behaviors such as growth and invasiveness remain unclear, hampering more effective coupling of tumor physical characteristics with implications for prognosis and therapy. Although molecular biology, histopathology, and radiological imaging are employed in this endeavor, studies are severely challenged by the multitude of different physical scales involved in tumor growth, i.e., from molecular nanoscale to cell microscale and finally to tissue centimeter scale. Consequently, it is often difficult to determine the underlying dynamics across dimensions. New techniques are needed to tackle these issues. Here, we address this multi-scalar problem by employing a novel predictive three-dimensional mathematical and computational model based on first-principle equations (conservation laws of physics) that describe mathematically the diffusion of cell substrates and other processes determining tumor mass growth and invasion. The model uses conserved variables to represent known determinants of glioma behavior, e.g., cell density and oxygen concentration, as well as biological functional relationships and parameters linking phenomena at different scales whose specific forms and values are hypothesized and calculated based on in vitro and in vivo experiments and from histopathology of tissue specimens from human gliomas. This model enables correlation of glioma morphology to tumor growth by quantifying interdependence of tumor mass on the microenvironment (e.g., hypoxia, tissue disruption) and on the cellular phenotypes (e.g., mitosis and apoptosis rates, cell adhesion strength). Once functional relationships between variables and associated parameter values have been informed, e.g., from histopathology or intra-operative analysis, this model can be used for disease diagnosis/prognosis, hypothesis testing, and to guide surgery and therapy. In particular, this tool identifies and quantifies the effects of vascularization and other cell-scale glioma morphological characteristics as predictors of tumor-scale growth and invasion.

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.

Antibodies targeting hepatoma-derived growth factor as a novel strategy in treating lung cancer.

Hepatoma-derived growth factor (HDGF) is overexpressed in lung cancer and the overexpression correlates with aggressive biological behaviors and poor clinical outcomes. We developed anti-HDGF monoclonal antibodies and tested their antitumor activity in lung cancer xenograft models. We also determined biological effects in tumors treated with the antibody alone or in combination with bevacizumab/avastin (an anti-vascular endothelial growth factor antibody) and/or gemcitabine (a chemotherapeutic agent). We found the anti-HDGF was effective to inhibit tumor growth in non-small cell lung cancer xenograft models. In the A549 model, compared with control IgG, tumor growth was substantially inhibited in animals treated with anti-HDGF antibodies, particularly HDGF-C1 (P = 0.002) and HDGF-H3 (P = 0.005). When HDGF-H3 was combined with either bevacizumab or gemcitabine, we observed enhanced tumor growth inhibition, particularly when the three agents were used together. HDGF-H3-treated tumors exhibited significant reduction of microvessel density with a pattern distinctive from the microvessel reduction pattern observed in bevacizumab-treated tumors. HDGF-H3-treated but not bevacizumab-treated tumors also showed a significant increase of apoptosis. Interestingly, many of the apoptotic cells in HDGF-H3-treated tumors are stroma cells, suggesting that the mechanism of the antitumor activity is, at least in part, through disrupting formation of tumor-stroma structures. Our results show that HDGF is a novel therapeutic target for lung cancer and can be effectively targeted by an antibody-based approach.