Survival analysis of time-to-event data in respiratory health research studies

This article provides a review of techniques for the analysis of survival data arising from respiratory health studies. Popular techniques such as the Kaplan–Meier survival plot and the Cox proportional hazards model are presented and illustrated using data from a lung cancer study. Advanced issues are also discussed, including parametric proportional hazards models, accelerated failure time models, time-varying explanatory variables, simultaneous analysis of multiple types of outcome events and the restricted mean survival time, a novel measure of the effect of treatment.

A boat hitchhiker's guide to survival: Cabomba caroliniana desiccation resistance and survival ability

Cabomba caroliniana is a submersed macrophyte that has become a serious invader. Cabomba predominantly spreads by stem fragments, in particular through unintentional transport on boat trailers ('hitch hiking'). Desiccation resistance affects the potential dispersal radius. Therefore, knowledge of maximum survival times allows predicting future dispersal. Experiments were conducted to assess desiccation resistance and survival ability of cabomba fragments under various environmental scenarios. Cabomba fragments were highly tolerant of desiccation. However, even relatively low wind speeds resulted in rapid mass loss, indicating a low survival rate of fragments exposed to air currents, such as fragments transported on a boat trailer. The experiments indicated that cabomba could survive at least 3 h of overland transport if exposed to wind. However, even small clumps of cabomba could potentially survive up to 42 h. Thus, targeting the transport of clumps of macrophytes should receive high priority in management. The high resilience of cabomba to desiccation demonstrates the risk of continuing spread. Because of the high probability of fragment viability on arrival, preventing fragment uptake on boat trailers is paramount to reduce the risk of further spread. These findings will assist improving models that predict the spread of aquatic invasive macrophytes.

A boat hitchhiker’s guide to survival: Cabomba caroliniana desiccation resistance and survival ability

Cabomba caroliniana is a submersed macrophyte that has become a serious invader. Cabomba predominantly spreads by stem fragments, in particular through unintentional transport on boat trailers (‘hitch hiking’). Desiccation resistance affects the potential dispersal radius. Therefore, knowledge of maximum survival times allows predicting future dispersal. Experiments were conducted to assess desiccation resistance and survival ability of cabomba fragments under various environmental scenarios. Cabomba fragments were highly tolerant of desiccation. However, even relatively low wind speeds resulted in rapid mass loss, indicating a low survival rate of fragments exposed to air currents, such as fragments transported on a boat trailer. The experiments indicated that cabomba could survive at least 3 h of overland transport if exposed to wind. However, even small clumps of cabomba could potentially survive up to 42 h. Thus, targeting the transport of clumps of macrophytes should receive high priority in management. The high resilience of cabomba to desiccation demonstrates the risk of continuing spread. Because of the high probability of fragment viability on arrival, preventing fragment uptake on boat trailers is paramount to reduce the risk of further spread. These findings will assist improving models that predict the spread of aquatic invasive macrophytes.

Hazard Function Models to Estimate Mortality Rates Affecting Fish Populations with Application to the Sea Mullet (Mugil cephalus) Fishery on the Queensland Coast (Australia)

Fisheries management agencies around the world collect age data for the purpose of assessing the status of natural resources in their jurisdiction. Estimates of mortality rates represent a key information to assess the sustainability of fish stocks exploitation. Contrary to medical research or manufacturing where survival analysis is routinely applied to estimate failure rates, survival analysis has seldom been applied in fisheries stock assessment despite similar purposes between these fields of applied statistics. In this paper, we developed hazard functions to model the dynamic of an exploited fish population. These functions were used to estimate all parameters necessary for stock assessment (including natural and fishing mortality rates as well as gear selectivity) by maximum likelihood using age data from a sample of catch. This novel application of survival analysis to fisheries stock assessment was tested by Monte Carlo simulations to assert that it provided unbiased estimations of relevant quantities. The method was applied to the data from the Queensland (Australia) sea mullet (Mugil cephalus) commercial fishery collected between 2007 and 2014. It provided, for the first time, an estimate of natural mortality affecting this stock: 0.22±0.08 year −1 .

The mesenchymal-to-epithelial reverting transition is enriched in metastatic castration resistant prostate cancer and correlates with poor survival

Despite recent recognition that the epithelial-mesenchymal transition (EMT) program acts in a dynamic manner (termed Epithelial to Mesenchymal Plasticity or EMP) during carcinoma metastasis, it has largely been ignored in the discovery and development of EMT-targeted therapies. In part, this has stemmed from a lack of preclinical models that can mimic the full dynamic nature of EMP and the perception that the EMT-reverting transition [or mesenchymal-epithelial reverting transition; (MErT)] is a mere antithesis of EMT. The objective of this study was to develop the first PCa model capable of recapitulating the dynamic nature of EMP.

Modelling initial survival and growth of radiata pine in New Zealand.

A sensitive framework has been developed for modelling young radiata pine survival, its growth and its size class distribution, from time of planting to age 5 or 6 years. The data and analysis refer to the Central North Island region of New Zealand. The survival function is derived from a Weibull probability density function, to reflect diminishing mortality with the passage of time in young stands. An anamorphic family of trends was used, as very little between-tree competition can be expected in young stands. An exponential height function was found to fit best the lower portion of its sigmoid form. The most appropriate basal area/ha exponential function included an allometric adjustment which resulted in compatible mean height and basal area/ha models. Each of these equations successfully represented the effects of several establishment practices by making coefficients linear functions of site factors, management activities and their interactions. Height and diameter distribution modelling techniques that ensured compatibility with stand values were employed to represent the effects of management practices on crop variation. Model parameters for this research were estimated using data from site preparation experiments in the region and were tested with some independent data sets.

Putative neuroprotective compounds in in vitro models of dopaminergic neurodegeneration

Parkinson´s disease (PD) is a debilitating age-related neurological disorder that affects various motor skills and can lead to a loss of cognitive functions. The motor symptoms are the result of the progressive degeneration of dopaminergic neurons within the substantia nigra. The factors that influence the pathogenesis and the progression of the neurodegeneration remain mostly unclear. This study investigated the role of various programmed cell death (PCD) pathways, oxidative stress, and glial cells both in dopaminergic neurodegeneration and in the protective action of various drugs. To this end, we exposed dopaminergic neuroblastoma cells (SH-SY5Y cells) to 6-OHDA, which produces oxidative stress and activates various PCD modalities that result in neuronal degeneration. Additionally, to explore the role of glia, we prepared rat midbrain primary mixed-cell cultures containing both neurons and glial cell types such as microglia and astroglia and then exposed the cultures to either MPP plus or lipopolysaccharide. Our results revealed that 6-OHDA activated several PCD pathways including apoptosis, autophagic stress, lysosomal membrane permeabilization, and perhaps paraptosis in SH-SY5Y cells. Furthermore, we found that minocycline protected SH-SY5Y cells from 6-OHDA by inhibiting both apoptotic and non-apoptotic PCD modalities. We also observed an inconsistent neuroprotective effect of various dietary anti-oxidant compounds against 6-OHDA toxicity in vitro in SH-SY5Y cells. Specifically, quercetin and curcumin exerted neuroprotection only within a narrow concentration range and a limited time frame, whereas resveratrol and epigallocatechin 3-gallate provided no protection whatsoever. Lastly, we found that molecules such as amantadine may delay or even halt the neurodegeneration in primary cell cultures by inhibiting the release of neurotoxic factors from overactivated microglia and by enhancing the pro-survival actions of astroglia. Together these data suggest that the strategy of dampening oxidative species with anti-oxidants is less effective than preventing the production of toxic factors such as oxidative and pro-inflammatory molecules by pathologically activated microglia. This would subsequently prevent the activation of various PCD modalities that cause neuronal degeneration.

Non-exponentiality in electron transfer kinetics: Static versus dynamic disorder models

Non-exponential electron transfer kinetics in complex systems are often analyzed in terms of a quenched, static disorder model. In this work we present an alternative analysis in terms of a simple dynamic disorder model where the solvent is characterized by highly non-exponential dynamics. We consider both low and high barrier reactions. For the former, the main result is a simple analytical expression for the survival probability of the reactant. In this case, electron transfer, in the long time, is controlled by the solvent polarization relaxation-in agreement with the analyses of Rips and Jortner and of Nadler and Marcus. The short time dynamics is also non-exponential, but for different reasons. The high barrier reactions, on the other hand, show an interesting dynamic dependence on the electronic coupling element, V-el.

Large rapidity gaps survival probabilities at LHC

We calculate the probability of large rapidity gaps in high energy hadronic collisions using a model based on QCD mini-jets and soft gluon emission down into the infrared region. Comparing with other models we find a remarkable agreement among most predictions.

STAT-1 expression is regulated by IGFBP-3 in malignant glioma cells and is a strong predictor of poor survival in patients with glioblastoma Laboratory investigation

Object. Insulin-like growth factor binding proteins (IGEBPs) have been implicated in the pathogenesis of glioma. In a previous study the authors demonstrated that IGFBP-3 is a novel glioblastoma biomarker associated with poor survival. Since signal transducer and activator of transcription 1 (STAT-1) has been shown to be regulated by IGFBP-3 during chondrogenesis and is a prosurvival and radioresistant molecule in different tumors, the aim in the present study was to explore the functional significance of IGFBP-3 in malignant glioma cells, to determine if STAT-1 is indeed regulated by IGFBP-3, and to study the potential of STAT-1 as a biomarker in glioblastoma. Methods. The functional significance of IGFBP-3 was investigated using the short hairpin (sh)RNA gene knockdown approach on U251MG cells. STAT-1 regulation by IGFBP-3 was tested on U251MG and U87MG cells by shRNA gene knockdown and exogenous treatment with recombinant IGFBP-3 protein. Subsequently, the expression of STAT-1 was analyzed with real-time reverse transcription polymerase chain reaction (RT-PCR) and immunohistochemistry (IHC) in glioblastoma and control brain tissues. Survival analyses were done on a uniformly treated prospective cohort of adults with newly diagnosed glioblastoma (136 patients) using Kaplan-Meier and Cox regression models. Results. IGFBP-3 knockdown significantly impaired proliferation, motility, migration, and invasive capacity of U251MG cells in vitro (p < 0.005). Exogenous overexpression of IGFBP-3 in U251MG and U87MG cells demonstrated STAT-1 regulation. The mean transcript levels (by real-time RT-PCR) and the mean labeling index of STAT-1 (by IHC) were significantly higher in glioblastoma than in control brain tissues (p = 0.0239 and p < 0.001, respectively). Multivariate survival analysis revealed that STAT-1 protein expression (HR 1.015, p = 0.033, 95% CI 1.001-1.029) along with patient age (HR 1.025, p = 0.005, 95% CI 1.008-1.042) were significant predictors of shorter survival in patients with glioblastoma. Conclusions. IGFBP-3 influences tumor cell proliferation, migration, and invasion and regulates STAT-1 expression in malignant glioma cells. STAT-1 is overexpressed in human glioblastoma tissues and emerges as a novel prognostic biomarker.

A Mycobacterial Phosphoribosyltransferase Promotes Bacillary Survival by Inhibiting Oxidative Stress and Autophagy Pathways in Macrophages and Zebrafish

Mycobacterium tuberculosis employs various strategies to modulate host immune responses to facilitate its persistence in macrophages. The M. tuberculosis cell wall contains numerous glycoproteins with unknown roles in pathogenesis. Here, by using Concanavalin A and LC-MS analysis, we identified a novel mannosylated glycoprotein phosphoribosyltransferase, encoded by Rv3242c from M. tuberculosis cell walls. Homology modeling, bioinformatic analyses, and an assay of phosphoribosyltransferase activity in Mycobacterium smegmatis expressing recombinant Rv3242c (MsmRv3242c) confirmed the mass spectrometry data. Using Mycobacterium marinum-zebrafish and the surrogate MsmRv3242c infection models, we proved that phosphoribosyltransferase is involved in mycobacterial virulence. Histological and infection assays showed that the M. marinum mimG mutant, an Rv3242c orthologue in a pathogenic M. marinum strain, was strongly attenuated in adult zebrafish and also survived less in macrophages. In contrast, infection with wild type and the complemented Delta mimG: Rv3242c M. marinum strains showed prominent pathological features, such as severe emaciation, skin lesions, hemorrhaging, and more zebrafish death. Similarly, recombinant Msm Rv3242c bacteria showed increased invasion in non-phagocytic epithelial cells and longer intracellular survival in macrophages as compared with wild type and vector control M. smegmatis strains. Further mechanistic studies revealed that the Rv3242c- and mimG-mediated enhancement of intramacrophagic survival was due to inhibition of autophagy, reactive oxygen species, and reduced activities of superoxide dismutase and catalase enzymes. Infection with MsmRv3242c also activated the MAPK pathway, NF-kappa B, and inflammatory cytokines. In summary, we show that a novel mycobacterial mannosylated phosphoribosyltransferase acts as a virulence and immunomodulatory factor, suggesting that it may constitute a novel target for antimycobacterial drugs.

Structural dynamic eutrophication models

This article discusses problems of modelling the seasonal succession of algal species in lakes and reservoirs, and the adaptive selection of certain groups of algae in response to changes in the inputs and relative concentrations of nutrients and other environmental variables. A new generation of quantitative models is being developed which attempts to translate some important biological properties of species (survival, variation, inheritance, reproductive rates and population growth) into predictions about the survival of the fittest, where ”fitness” is measured or estimated in thermodynamic terms. The concept of ”exergy” and its calculation is explored to examine maximal exergy as a measure of fitness in ecosystems, and its use for calculating changes in species composition by means of structural dynamic models. These models accomodate short-term changes in parameters that affect the adaptive responses (species selection) of algae.

Is there a subgroup of long-term evolution among patients with advanced lung cancer?: Hints from the analysis of survival curves from cancer registry data

Background: Recently, with the access of low toxicity biological and targeted therapies, evidence of the existence of a long-term survival subpopulation of cancer patients is appearing. We have studied an unselected population with advanced lung cancer to look for evidence of multimodality in survival distribution, and estimate the proportion of long-term survivors. Methods: We used survival data of 4944 patients with non-small-cell lung cancer (NSCLC) stages IIIb-IV at diagnostic, registered in the National Cancer Registry of Cuba (NCRC) between January 1998 and December 2006. We fitted one-component survival model and two-component mixture models to identify short-and long-term survivors. Bayesian information criterion was used for model selection. Results: For all of the selected parametric distributions the two components model presented the best fit. The population with short-term survival (almost 4 months median survival) represented 64% of patients. The population of long-term survival included 35% of patients, and showed a median survival around 12 months. None of the patients of short-term survival was still alive at month 24, while 10% of the patients of long-term survival died afterwards. Conclusions: There is a subgroup showing long-term evolution among patients with advanced lung cancer. As survival rates continue to improve with the new generation of therapies, prognostic models considering short-and long-term survival subpopulations should be considered in clinical research.

Life-stage recruitment models for Atlantic cod (Gadus morhua) and haddock (Melanogrammus aeglefinus) on Georges Bank

Ichthyoplankton surveys have been used to provide an independent estimate of adult spawning biomass of commercially exploited species and to further our understanding of the recruitment processes in the early life stages. However, predicting recruitment has been difficult because of the complex interaction of physical and biological processes operating at different spatial and temporal scales that can occur at the different life stages. A model of first-year life-stage recruitment was applied to Georges Bank Atlantic cod (Gadus morhua) and haddock (Melanogrammus aeglefinus) stocks over the years 1977–2004 by using environmental and densitydependent relationships. The best lifestage mortality relationships for eggs, larvae, pelagic juveniles, and demersal juveniles were first determined by hindcasting recruitment estimates based on egg and larval abundance and mortality rates derived from two intensive sampling periods, 1977–87 and 1995–99. A wind-driven egg mortality relationship was used to estimate losses due to transport off the bank, and a wind-stress larval mortality relationship was derived from feeding and survival studies. A simple metric for the density-dependent effects of Atlantic cod was used for both Atlantic cod and haddock. These life stage proxies were then applied to the virtual population analysis (VPA) derived annual egg abundances to predict age-1 recruitment. Best models were determined from the correlation of predicted and VPA-derived age-1 abundance. The larval stage was the most quantifiable of any stage from surveys, whereas abundance estimates of the demersal juvenile stage were not available because of undersampling. Attempts to forecast recruitment from spawning stock biomass or egg abundance, however, will always be poor because of variable egg survival.

The effects of temperature on hatching and survival of northern rock sole larvae (Lepidopsetta polyxystra)

Northern rock sole (Lepidopsetta polyxystra) is a commercially important flatfish in Alaska and was recently classified as a distinct species from southern rock sole (L. bilineata). Taxonomic and vital rate data for northern rock sole are still not fully described, notably at early egg and larval stages. In this study, we provide new taxonomic descriptions of late-stage eggs and newly hatched larvae, as well as temperature-response models of hatching (timing, duration, success), and larval size-at-hatch and posthatch survival at four temperatures (2°, 5°, 9°, and 12°C). Time-to-first-hatch, hatch cycle duration, and overall hatching success showed a negative relationship with temperature. Early hatching larvae within each temperature treatment were smaller and had larger yolk sacs, but larvae incubated at higher temperatures (9° and 12°C) had the largest yolk reserves overall. Despite having smaller yolks, size-at-hatch and the maximum size achieved during the hatching cycle was highest for larvae reared at cold temperatures (2° and 5°C), indicating that endogenous reserves are more efficiently used for growth at these temperatures. In addition, larvae reared at high temperatures died more rapidly in the absence of food despite having more yolk reserves than cold-incubated larvae. Overall, northern rock sole eggs and larvae display early life history traits consistent with coldwater adaptation for winter spawning in the North Pacific.