Evaluation of Viral Load Thresholds for Predicting New World Health Organization Stage 3 and 4 Events in HIV-Infected Children Receiving Highly Active Antiretroviral Therapy

Background: This study evaluated a wide range of viral load (VL) thresholds to identify a cut-point that best predicts new clinical events in children on stable highly active antiretroviral therapy (HAART). Methods: Cox proportional hazards modeling was used to assess the adjusted risk for World Health Organization stage 3 or 4 clinical events (WHO events) as a function of time-varying CD4, VL, and hemoglobin values in a cohort study of Latin American children on HAART >= 6 months. Models were fit using different VL cut-points between 400 and 50,000 copies per milliliter, with model fit evaluated on the basis of the minimum Akaike information criterion value, a standard model fit statistic. Results: Models were based on 67 subjects with WHO events out of 550 subjects on study. The VL cut-points of >2600 and >32,000 copies per milliliter corresponded to the lowest Akaike information criterion values and were associated with the highest hazard ratios (2.0, P = 0.015; and 2.1, P = 0.0058, respectively) for WHO events. Conclusions: In HIV-infected Latin American children on stable HAART, 2 distinct VL thresholds (>2600 and >32,000 copies/mL) were identified for predicting children at significantly increased risk for HIV-related clinical illness, after accounting for CD4 level, hemoglobin level, and other significant factors.

Association between UCP2 A55V polymorphism and risk of cardiovascular events in patients with multi-vessel coronary arterial disease

Background UCP2 (uncoupling protein 2) plays an important role in cardiovascular diseases and recent studies have suggested that the A55V polymorphism can cause UCP2 dysfunction. The main aim was to investigate the association of A55V polymorphism with cardiovascular events in a group of 611 patients enrolled in the Medical, Angioplasty or Surgery Study II (MASS II), a randomized trial comparing treatments for patients with coronary artery disease and preserved left ventricular function. Methods The participants of the MASS II were genotyped for the A55V polymorphism using allele-specific PCR assay. Survival curves were calculated with the Kaplan–Meier method and evaluated with the log-rank statistic. The relationship between baseline variables and the composite end-point of cardiac death, acute myocardial infarction (AMI), refractory angina requiring revascularization and cerebrovascular accident were assessed using a Cox proportional hazards survival model. Results There were no significant differences for baseline variables according genotypes. After 2 years of follow-up, dysglycemic patients harboring the VV genotype had higher occurrence of AMI (p=0.026), Death+AMI (p=0.033), new revascularization intervention (p=0.009) and combined events (p=0.037) as compared with patients carrying other genotypes. This association was not evident in normoglycemic patients. Conclusions These findings support the hypothesis that A55V polymorphism is associated with UCP2 functional alterations that increase the risk of cardiovascular events in patients with previous coronary artery disease and dysglycemia.

Risks of non-accidental mortality by baseline CD4+ T-cell strata in hepatitis-C-positive and -negative individuals initiating highly active antiretroviral therapy

BACKGROUND: Patients coinfected with hepatitis C virus (HCV) and HIV experience higher mortality rates than patients infected with HIV alone. We designed a study to determine whether risks for later mortality are similar for HCV-positive and HCV-negative individuals when subjects are stratified on the basis of baseline CD4+ T-cell counts. METHODS: Antiretroviral-naive individuals, who initiated highly active antiretroviral therapy (HAART) between 1996 and 2002 were included in the study. HCV-positive and HCV-negative individuals were stratified separately by baseline CD4+ T-cell counts of 50 cell/microl increments. Cox-proportional hazards regression was used to model the effect of these strata with other variables on survival. RESULTS: CD4+ T-cell strata below 200 cells/microl, but not above, imparted an increased relative hazard (RH) of mortality for both HCV-positive and HCV-negative individuals. Among HCV-positive individuals, after adjustment for baseline age, HIV RNA levels, history of injection drug use and adherence to therapy, only CD4+ T-cell strata of <50 cells/microl (RH=4.60; 95% confidence interval [CI] 2.72-7.76) and 50-199 cells/microl (RH=2.49; 95% CI 1.63-3.81) were significantly associated with increased mortality when compared with those initiating therapy at cell counts >500 cells/microl. The same baseline CD4+ T-cell strata were found for HCV-negative individuals. CONCLUSION: In a within-groups analysis, the baseline CD4+ T-cell strata that are associated with increased RHs for mortality are the same for HCV-positive and HCV-negative individuals initiating HAART. However, a between-groups analysis reveals a higher absolute mortality risk for HCV-positive individuals.

Inference on Survival Data with Covariate Measurement Error - An Imputation-based Approach

We propose a new method for fitting proportional hazards models with error-prone covariates. Regression coefficients are estimated by solving an estimating equation that is the average of the partial likelihood scores based on imputed true covariates. For the purpose of imputation, a linear spline model is assumed on the baseline hazard. We discuss consistency and asymptotic normality of the resulting estimators, and propose a stochastic approximation scheme to obtain the estimates. The algorithm is easy to implement, and reduces to the ordinary Cox partial likelihood approach when the measurement error has a degenerative distribution. Simulations indicate high efficiency and robustness. We consider the special case where error-prone replicates are available on the unobserved true covariates. As expected, increasing the number of replicate for the unobserved covariates increases efficiency and reduces bias. We illustrate the practical utility of the proposed method with an Eastern Cooperative Oncology Group clinical trial where a genetic marker, c-myc expression level, is subject to measurement error.

Designed Extension of Survival Studies: Application to Clinical Trials with Unrecognized Heterogeneity

It is well known that unrecognized heterogeneity among patients, such as is conferred by genetic subtype, can undermine the power of randomized trial, designed under the assumption of homogeneity, to detect a truly beneficial treatment. We consider the conditional power approach to allow for recovery of power under unexplained heterogeneity. While Proschan and Hunsberger (1995) confined the application of conditional power design to normally distributed observations, we consider more general and difficult settings in which the data are in the framework of continuous time and are subject to censoring. In particular, we derive a procedure appropriate for the analysis of the weighted log rank test under the assumption of a proportional hazards frailty model. The proposed method is illustrated through application to a brain tumor trial.

Additive Hazards Models with Latent Treatment Effectiveness Lag Time

In many clinical trials to evaluate treatment efficacy, it is believed that there may exist latent treatment effectiveness lag times after which medical procedure or chemical compound would be in full effect. In this article, semiparametric regression models are proposed and studied to estimate the treatment effect accounting for such latent lag times. The new models take advantage of the invariance property of the additive hazards model in marginalizing over random effects, so parameters in the models are easy to be estimated and interpreted, while the flexibility without specifying baseline hazard function is kept. Monte Carlo simulation studies demonstrate the appropriateness of the proposed semiparametric estimation procedure. Data collected in the actual randomized clinical trial, which evaluates the effectiveness of biodegradable carmustine polymers for treatment of recurrent brain tumors, are analyzed.

Modeling multilevel sleep transitional data via Poisson log-linear multilevel models

This paper proposes Poisson log-linear multilevel models to investigate population variability in sleep state transition rates. We specifically propose a Bayesian Poisson regression model that is more flexible, scalable to larger studies, and easily fit than other attempts in the literature. We further use hierarchical random effects to account for pairings of individuals and repeated measures within those individuals, as comparing diseased to non-diseased subjects while minimizing bias is of epidemiologic importance. We estimate essentially non-parametric piecewise constant hazards and smooth them, and allow for time varying covariates and segment of the night comparisons. The Bayesian Poisson regression is justified through a re-derivation of a classical algebraic likelihood equivalence of Poisson regression with a log(time) offset and survival regression assuming piecewise constant hazards. This relationship allows us to synthesize two methods currently used to analyze sleep transition phenomena: stratified multi-state proportional hazards models and log-linear models with GEE for transition counts. An example data set from the Sleep Heart Health Study is analyzed.

Time to renal disease and end-stage renal disease in PROFILE: a multiethnic lupus cohort.

BACKGROUND: Renal involvement is a serious manifestation of systemic lupus erythematosus (SLE); it may portend a poor prognosis as it may lead to end-stage renal disease (ESRD). The purpose of this study was to determine the factors predicting the development of renal involvement and its progression to ESRD in a multi-ethnic SLE cohort (PROFILE). METHODS AND FINDINGS: PROFILE includes SLE patients from five different United States institutions. We examined at baseline the socioeconomic-demographic, clinical, and genetic variables associated with the development of renal involvement and its progression to ESRD by univariable and multivariable Cox proportional hazards regression analyses. Analyses of onset of renal involvement included only patients with renal involvement after SLE diagnosis (n = 229). Analyses of ESRD included all patients, regardless of whether renal involvement occurred before, at, or after SLE diagnosis (34 of 438 patients). In addition, we performed a multivariable logistic regression analysis of the variables associated with the development of renal involvement at any time during the course of SLE.In the time-dependent multivariable analysis, patients developing renal involvement were more likely to have more American College of Rheumatology criteria for SLE, and to be younger, hypertensive, and of African-American or Hispanic (from Texas) ethnicity. Alternative regression models were consistent with these results. In addition to greater accrued disease damage (renal damage excluded), younger age, and Hispanic ethnicity (from Texas), homozygosity for the valine allele of FcgammaRIIIa (FCGR3A*GG) was a significant predictor of ESRD. Results from the multivariable logistic regression model that included all cases of renal involvement were consistent with those from the Cox model. CONCLUSIONS: Fcgamma receptor genotype is a risk factor for progression of renal disease to ESRD. Since the frequency distribution of FCGR3A alleles does not vary significantly among the ethnic groups studied, the additional factors underlying the ethnic disparities in renal disease progression remain to be elucidated.

Health-related quality of life and survival in liver transplant candidates.

Health-related quality of life (HRQOL) is an important measure of the effects of chronic liver disease in affected patients that helps guide interventions to improve well-being. However, the relationship between HRQOL and survival in liver transplant candidates remains unclear. We examined whether the Physical Component Summary (PCS) and Mental Component Summary (MCS) scores from the Short Form 36 (SF-36) Health Survey were associated with survival in liver transplant candidates. We administered the SF-36 questionnaire (version 2.0) to patients in the Pulmonary Vascular Complications of Liver Disease study, a multicenter prospective cohort of patients evaluated for liver transplantation in 7 academic centers in the United States between 2003 and 2006. Cox proportional hazards models were used with death as the primary outcome and adjustment for liver transplantation as a time-varying covariate. The mean age of the 252 participants was 54 +/- 10 years, 64% were male, and 94% were white. During the 422 person years of follow-up, 147 patients (58%) were listed, 75 patients (30%) underwent transplantation, 49 patients (19%) died, and 3 patients were lost to follow-up. Lower baseline PCS scores were associated with an increased mortality rate despite adjustments for age, gender, Model for End-Stage Liver Disease score, and liver transplantation (P for the trend = 0.0001). The MCS score was not associated with mortality (P for the trend = 0.53). In conclusion, PCS significantly predicts survival in liver transplant candidates, and interventions directed toward improving the physical status may be helpful in improving outcomes in liver transplant candidates.

Relationships between the alterations of tumor suppressor genes, risk factors and clinical outcomes of gliomas

The relationship between MMAC/PTEN, DMBT1 and the progression and prognosis of glioma, and the association between the alterations of MMAC/PTEN, p53, p16, and Rb and some cancer risk factors, such as smoking, exposure to radiation, family cancer history, and previous cancer history, were assessed in 4 studies. ^ By allelic deletion analysis, MMAC/PTEN locus was shown to be frequently lost in glioblastomas multiforme (GM) but maintained in most lower-grade astrocytic tumors. DMBT1 locus, however, was frequently lost in all grades of gliomas examined. The potential biological significance of these two regions was frontier assessed by examining microcell-hybrids that contained various fragments of 10q. Somatic cell hybrid clones that retained the MMAC/PTEN locus have less transformed phenotypes, exhibiting an inability to grow in soft agarose. On the other hand, the presence or absence of DAMT1 did not correlate with any in vitro phenotype assessed in our model system. Further, Cox proportional hazards regression analysis, adjusted for age at surgery and histologic grades (GM, and non-GM), showed that without LOH at the MMAC/PTEN locus had a significantly better prognosis than did patients with LOH at MMAC/ PTEN (hazard ratio = 0.5; 95% Cl = 0.28–0.89; P = 0.018). Furthermore, status of LOH at MMAC/PTEN was found to be significantly associated with age, while that for DMBT1 was not. These results suggest that the DMBT1 may be involved early in the oncogenesis of gliomas, while alterations in the MMAC /PTEN may be a late event in the oncogenesis related with progression of gliomas and provide a significant prognostic marker for patient survival. ^ The associations between 4 cancer risk factors and 4 tumor suppressor genes were assessed. The expression of p16 was observed to be associated with current smoking (adjusted OR = 1.9, 95% CI = 1.02–3.6) but not the former smoking (adjusted OR = 1.1, 95% Cl = 0.5–3.5). The expression of p53 was found to be associated with the family cancer history (OR = 3.5, 95% Cl = 1.07–11 for patients with first-degree family history of cancer). MMAC/ PTEN was associated with the histologic grade (OR = 2.8, 95% CI = 1.2–6.6) and age (P = 0.035). Also, the OR for LOH around MMAC/PTEN in patients with a family history of cancer was elevated (OR = 1.9, 95% CI = 0.8–4.6 for patients with first-degree family history of cancer). The associations between exposure and the alterations of tumor suppressor genes, between smoking and p16, between family history of cancer and p53 and MMAC/PTEN, provide suggestive evidences that those exposures are related to the development of gliomas. ^

Patients with community acquired pneumonia admitted to European intensive care units: an epidemiological survey of the GenOSept cohort.

INTRODUCTION Community acquired pneumonia (CAP) is the most common infectious reason for admission to the Intensive Care Unit (ICU). The GenOSept study was designed to determine genetic influences on sepsis outcome. Phenotypic data was recorded using a robust clinical database allowing a contemporary analysis of the clinical characteristics, microbiology, outcomes and independent risk factors in patients with severe CAP admitted to ICUs across Europe. METHODS Kaplan-Meier analysis was used to determine mortality rates. A Cox Proportional Hazards (PH) model was used to identify variables independently associated with 28-day and six-month mortality. RESULTS Data from 1166 patients admitted to 102 centres across 17 countries was extracted. Median age was 64 years, 62% were male. Mortality rate at 28 days was 17%, rising to 27% at six months. Streptococcus pneumoniae was the commonest organism isolated (28% of cases) with no organism identified in 36%. Independent risk factors associated with an increased risk of death at six months included APACHE II score (hazard ratio, HR, 1.03; confidence interval, CI, 1.01-1.05), bilateral pulmonary infiltrates (HR1.44; CI 1.11-1.87) and ventilator support (HR 3.04; CI 1.64-5.62). Haematocrit, pH and urine volume on day one were all associated with a worse outcome. CONCLUSIONS The mortality rate in patients with severe CAP admitted to European ICUs was 27% at six months. Streptococcus pneumoniae was the commonest organism isolated. In many cases the infecting organism was not identified. Ventilator support, the presence of diffuse pulmonary infiltrates, lower haematocrit, urine volume and pH on admission were independent predictors of a worse outcome.

Assessing efficacy of different nucleos(t)ide backbones in NNRTI-containing regimens in the Swiss HIV Cohort Study.

BACKGROUND The most recommended NRTI combinations as first-line antiretroviral treatment for HIV-1 infection in resource-rich settings are tenofovir/emtricitabine, abacavir/lamivudine, tenofovir/lamivudine and zidovudine/lamivudine. Efficacy studies of these combinations also considering pill numbers, dosing frequencies and ethnicities are rare. METHODS We included patients starting first-line combination ART (cART) with or switching from first-line cART without treatment failure to tenofovir/emtricitabine, abacavir/lamivudine, tenofovir/lamivudine and zidovudine/lamivudine plus efavirenz or nevirapine. Cox proportional hazards regression was used to investigate the effect of the different NRTI combinations on two primary outcomes: virological failure (VF) and emergence of NRTI resistance. Additionally, we performed a pill burden analysis and adjusted the model for pill number and dosing frequency. RESULTS Failure events per treated patient for the four NRTI combinations were as follows: 19/1858 (tenofovir/emtricitabine), 9/387 (abacavir/lamivudine), 11/344 (tenofovir/lamivudine) and 45/1244 (zidovudine/lamivudine). Compared with tenofovir/emtricitabine, abacavir/lamivudine had an adjusted HR for having VF of 2.01 (95% CI 0.86-4.55), tenofovir/lamivudine 2.89 (1.22-6.88) and zidovudine/lamivudine 2.28 (1.01-5.14), whereas for the emergence of NRTI resistance abacavir/lamivudine had an HR of 1.17 (0.11-12.2), tenofovir/lamivudine 11.3 (2.34-55.3) and zidovudine/lamivudine 4.02 (0.78-20.7). Differences among regimens disappeared when models were additionally adjusted for pill burden. However, non-white patients compared with white patients and higher pill number per day were associated with increased risks of VF and emergence of NRTI resistance: HR of non-white ethnicity for VF was 2.85 (1.64-4.96) and for NRTI resistance 3.54 (1.20-10.4); HR of pill burden for VF was 1.41 (1.01-1.96) and for NRTI resistance 1.72 (0.97-3.02). CONCLUSIONS Although VF and emergence of resistance was very low in the population studied, tenofovir/emtricitabine appears to be superior to abacavir/lamivudine, tenofovir/lamivudine and zidovudine/lamivudine. However, it is unclear whether these differences are due to the substances as such or to an association of tenofovir/emtricitabine regimens with lower pill burden.

Normal Thyroid Function and the Risk of Atrial Fibrillation: the Rotterdam Study

CONTEXT Hyperthyroidism is an established risk factor for atrial fibrillation (AF), but information concerning the association with variations within the normal range of thyroid function and subgroups at risk is lacking. OBJECTIVE This study aimed to investigate the association between normal thyroid function and AF prospectively and explore potential differential risk patterns. DESIGN, SETTING, AND PARTICIPANTS From the Rotterdam Study we included 9166 participants ≥ 45 y with TSH and/or free T4 (FT4) measurements and AF assessment (1997-2012 median followup, 6.8 y), with 399 prevalent and 403 incident AF cases. MAIN OUTCOME MEASURES Outcome measures were 3-fold: 1) hazard ratios (HRs) for the risk of incident AF by Cox proportional-hazards models, 2) 10-year absolute risks taking competing risk of death into account, and 3) discrimination ability of adding FT4 to the CHARGE-AF simple model, an established prediction model for AF. RESULTS Higher FT4 levels were associated with higher risks of AF (HR 1.63, 95% confidence interval, 1.19-2.22), when comparing those in the highest quartile to those in lowest quartile. Absolute 10-year risks increased with higher FT4 in participants ≤ 65 y from 1-9% and from 6-12% in subjects ≥ 65 y. Discrimination of the prediction model improved when adding FT4 to the simple model (c-statistic, 0.722 vs 0.729; P = .039). TSH levels were not associated with AF. CONCLUSIONS There is an increased risk of AF with higher FT4 levels within the normal range, especially in younger subjects. Adding FT4 to the simple model slightly improved discrimination of risk prediction.

Modeling survival and the development of second primary tumors among Hodgkin's disease patients

Hodgkin's disease (HD) is a cancer of the lymphatic system. Survivors of HD face varieties of consequent adverse effects, in which secondary primary tumors (SPT) is one of the most serious consequences. This dissertation is aimed to model time-to-SPT in the presence of death and HD relapses during follow-up.^ The model is designed to handle a mixture phenomenon of SPT and the influence of death. Relapses of HD are adjusted as a covariate. Proportional hazards framework is used to define SPT intensity function, which includes an exponential term to estimate explanatory variables. Death as a competing risk is considered according to different scenarios, depending on which terminal event comes first. Newton-Raphson method is used to estimate the parameter estimates in the end.^ The proposed method is applied to a real data set containing a group of HD patients. Several risk factors for the development of SPT are identified and the findings are noteworthy in the development of healthcare guidelines that may lead to the early detection or prevention of SPT.^

Diabetes genes and risk of prostate cancer in the Atherosclerosis Risk in Communities study

Prostate cancer (PrCa) is a leading cause of morbidity and mortality, yet the etiology remains uncertain. Meta-analyses show that PrCa risk is reduced by 16% in men with type 2 diabetes (T2D), but the mechanism is unknown. Recent genome-wide association studies and meta-analyses have found single nucleotide polymorphisms (SNPs) that consistently predict T2D risk. We evaluated associations of incident PrCa with 14 T2D SNPs in the Atherosclerosis Risk in Communities (ARIC) study. From 1987-2000, there were 397 incident PrCa cases ascertained from state or local cancer registries among 6,642 men (1,560 blacks and 5,082 whites) aged 45-64 years at baseline. Genotypes were determined by TaqMan assay. Cox proportional hazards models were used to assess the association between PrCa and increasing number of T2D risk-raising alleles for individual SNPs and for genetic risk scores (GRS) comprised of the number of T2D risk-raising alleles across SNPs. Two-way gene-gene interactions were evaluated with likelihood ratio tests. Using additive genetic models, the T2D risk-raising allele was associated with significantly reduced risk of PrCa for IGF2BP2 rs4402960 (hazard ratio [HR]=0.79; P=0.07 among blacks only), SLC2A2 rs5400 (race-adjusted HR=0.85; P=0.05) and UCP2 rs660339 (race-adjusted HR=0.84; P=0.02), but significantly increased risk of PrCa for CAPN10 rs3792267 (race-adjusted HR=1.20; P=0.05). No other SNPs were associated with PrCa using an additive genetic model. However, at least one copy of the T2D risk-raising allele for TCF7L2 rs7903146 was associated with reduced PrCa risk using a dominant genetic model (race-adjusted HR=0.79; P=0.03). These results imply that the T2D-PrCa association may be partly due to shared genetic variation, but these results should be verified since multiple tests were performed. When the combined, additive effects of these SNPs were tested using a GRS, there was nearly a 10% reduction in risk of PrCa per T2D risk-raising allele (race-adjusted HR=0.92; P=0.02). SNPs in IGF2BP2, KCNJ11 and SLC2A2 were also involved in multiple synergistic gene-gene interactions on a multiplicative scale. In conclusion, it appears that the T2D-PrCa association may be due, in part, to common genetic variation. Further knowledge of T2D gene-PrCa mechanisms may improve understanding of PrCa etiology and may inform PrCa prevention and treatment.^