When one depends on the other: reporting of interaction in case-control and cohort studies

BACKGROUND: Interaction refers to the situation in which the effect of 1 exposure on an outcome differs across strata of another exposure. We did a survey of epidemiologic studies published in leading journals to examine how interaction is assessed and reported. METHODS: We selected 150 case-control and 75 cohort studies published between May 2001 and May 2007 in leading general medicine, epidemiology, and clinical specialist journals. Two reviewers independently extracted data on study characteristics. RESULTS: Of the 225 studies, 138 (61%) addressed interaction. Among these, 25 (18%) presented no data or only a P value or a statement of statistical significance; 40 (29%) presented stratum-specific effect estimates but no meaningful comparison of these estimates; and 58 (42%) presented stratum-specific estimates and appropriate tests for interaction. Fifteen articles (11%) presented the individual effects of both exposures and also their joint effect or a product term, providing sufficient information to interpret interaction on an additive and multiplicative scale. Reporting was poorest in articles published in clinical specialist articles and most adequate in articles published in general medicine journals, with epidemiology journals in an intermediate position. CONCLUSIONS: A majority of articles reporting cohort and case-control studies address possible interactions between exposures. However, in about half of these, the information provided was unsatisfactory, and only 1 in 10 studies reported data that allowed readers to interpret interaction effects on an additive and multiplicative scale.

Potent and broad-spectrum antimicrobial activity of CXCL14 suggests an immediate role in skin infections

The skin is constantly exposed to commensal microflora and pathogenic microbes. The stratum corneum of the outermost skin layer employs distinct tools such as harsh growth conditions and numerous antimicrobial peptides (AMPs) to discriminate between beneficial cutaneous microflora and harmful bacteria. How the skin deals with microbes that have gained access to the live part of the skin as a result of microinjuries is ill defined. In this study, we report that the chemokine CXCL14 is a broad-spectrum AMP with killing activity for cutaneous gram-positive bacteria and Candida albicans as well as the gram-negative enterobacterium Escherichia coli. Based on two separate bacteria-killing assays, CXCL14 compares favorably with other tested AMPs, including human beta-defensin and the chemokine CCL20. Increased salt concentrations and skin-typical pH conditions did not abrogate its AMP function. This novel AMP is highly abundant in the epidermis and dermis of healthy human skin but is down-modulated under conditions of inflammation and disease. We propose that CXCL14 fights bacteria at the earliest stage of infection, well before the establishment of inflammation, and thus fulfills a unique role in antimicrobial immunity.

Long-term changes in bone metabolism, bone mineral density, quantitative ultrasound parameters, and fracture incidence after spinal cord injury: a cross-sectional observational study in 100 paraplegic men

To study the time course of demineralization and fracture incidence after spinal cord injury (SCI), 100 paraplegic men with complete motor loss were investigated in a cross-sectional study 3 months to 30 years after their traumatic SCI. Fracture history was assessed and verified using patients' files and X-rays. BMD of the lumbar spine (LS), femoral neck (FN), distal forearm (ultradistal part = UDR, 1/3 distal part = 1/3R), distal tibial diaphysis (TDIA), and distal tibial epiphysis (TEPI) was measured using DXA. Stiffness of the calcaneus (QUI.CALC), speed of sound of the tibia (SOS.TIB), and amplitude-dependent SOS across the proximal phalanges (adSOS.PHAL) were measured using QUS. Z-Scores of BMD and quantitative ultrasound (QUS) were plotted against time-since-injury and compared among four groups of paraplegics stratified according to time-since-injury (<1 year, stratum I; 1-9 years, stratum II; 10-19 years, stratum III; 20-29 years, stratum IV). Biochemical markers of bone turnover (deoxypyridinoline/creatinine (D-pyr/Cr), osteocalcin, alkaline phosphatase) and the main parameters of calcium phosphate metabolism were measured. Fifteen out of 98 paraplegics had sustained a total of 39 fragility fractures within 1,010 years of observation. All recorded fractures were fractures of the lower limbs, mean time to first fracture being 8.9 +/- 1.4 years. Fracture incidence increased with time-after-SCI, from 1% in the first 12 months to 4.6%/year in paraplegics since >20 years ( p<.01). The overall fracture incidence was 2.2%/year. Compared with nonfractured paraplegics, those with a fracture history had been injured for a longer time ( p<.01). Furthermore, they had lower Z-scores at FN, TEPI, and TDIA ( p<.01 to <.0001), the largest difference being observed at TDIA, compared with the nonfractured. At the lower limbs, BMD decreased with time at all sites ( r=.49 to.78, all p<.0001). At FN and TEPI, bone loss followed a log curve which leveled off between 1 to 3 years after injury. In contrast, Z-scores of TDIA continuously decreased even beyond 10 years after injury. LS BMD Z-score increased with time-since-SCI ( p<.05). Similarly to DXA, QUS allowed differentiation of early and rapid trabecular bone loss (QUI.CALC) vs slow and continuous cortical bone loss (SOS.TIB). Biochemical markers reflected a disproportion between highly elevated bone resorption and almost normal bone formation early after injury. Turnover declined following a log curve with time-after-SCI, however, D-pyr/Cr remained elevated in 30% of paraplegics injured >10 years. In paraplegic men early (trabecular) and persistent (cortical) bone loss occurs at the lower limbs and leads to an increasing fracture incidence with time-after-SCI.

Nuclear localization of HBD-1 in human keratinocytes.

OBJECTIVE: Human defensins and cathelicidins are a family of cationic antimicrobial peptides (AMPs), which play multiple roles in both innate and adaptive immune systems. They have direct antimicrobial activity against several microorganisms including burn pathogens. The majority of components of innate and adaptive immunity either express naturally occurring defensins or are otherwise chemoattracted or functionally affected by them. They also enhance adaptive immunity and wound healing and alter antibody production. All mechanisms to explain multiple functions of AMPs are not clearly understood. Prior studies to localize defensins in normal and burned skin using deconvolution fluorescence scanning microscopy indicate localization of defensins in the nucleus, perinuclear regions, and cytoplasm. The objective of this study is to further confirm the identification of HBD-1 in the nucleus by deconvolution microscopic studies involving image reconstruction and wire frame modeling. RESULTS: Our study demonstrated the presence of intranuclear HBD-1 in keratinocytes throughout the stratum spinosum by costaining with the nuclear probe DAPI. In addition, HBD-1 sequence does show some homology with known cationic nuclear localization signal sequences. CONCLUSION: To our knowledge, this is the first report to localize HBD-1 in the nuclear region, suggesting a role for this peptide in gene expression and providing new data that may help determine mechanisms of defensin functions.

Migrant Background and Higher Education Participation in Europe. The Effect of the Educational Systems

The main aim of this article is to shed light on the extent to which differences in higher education participation between people with and without a migrant background of low/higher social origin can be explained by two macro-level characteristics of national educational institutions: stratification of the secondary school system and provision of alternative access to higher education. General assumptions are that people with a migrant background of low social origin benefit in low-stratified secondary school systems and in systems that provide alternative access to institutions of higher education more than their native peers in the same social stratum, owing to primary and secondary effects of migrant background. Database is a pooled dataset of the five waves of the European Social Survey. Results of logistic multi-level analyses indicate that a low-stratified secondary school system improves the probability of people with a migrant background/low social origin attaining a higher education degree. On the other hand, a stratified secondary school system reduces their chances regarding this educational stage. The provision of alternative access to an institution of higher education improves their likelihood of becoming higher education graduates.

Surface Deformations as Indicators of Deep Ebullition Fluxes in a Large Northern Peatland

Peatlands deform elastically during precipitation cycles by small (+/- 3 cm) oscillations in surface elevation. In contrast, we used a Global Positioning System network to measure larger oscillations that exceeded 20 cm over periods of 4 - 12 hours during two seasonal droughts at a bog and fen site in northern Minnesota. The second summer drought also triggered 19 depressuring cycles in an overpressured stratum under the bog site. The synchronicity between the largest surface deformations and the depressuring cycles indicates that both phenomena are produced by the episodic release of large volumes of gas from deep semi-elastic compartments confined by dense wood layers. We calculate that the three largest surface deformations were associated with the release of 136 g CH4 m(-2), which exceeds by an order of magnitude the annual average chamber fluxes measured at this site. Ebullition of gas from the deep peat may therefore be a large and previously unrecognized source of radiocarbon depleted methane emissions from northern peatlands.

An index to identify stroke-related vs incidental patent foramen ovale in cryptogenic stroke

OBJECTIVE We aimed to create an index to stratify cryptogenic stroke (CS) patients with patent foramen ovale (PFO) by their likelihood that the stroke was related to their PFO. METHODS Using data from 12 component studies, we used generalized linear mixed models to predict the presence of PFO among patients with CS, and derive a simple index to stratify patients with CS. We estimated the stratum-specific PFO-attributable fraction and stratum-specific stroke/TIA recurrence rates. RESULTS Variables associated with a PFO in CS patients included younger age, the presence of a cortical stroke on neuroimaging, and the absence of these factors: diabetes, hypertension, smoking, and prior stroke or TIA. The 10-point Risk of Paradoxical Embolism score is calculated from these variables so that the youngest patients with superficial strokes and without vascular risk factors have the highest score. PFO prevalence increased from 23% (95% confidence interval [CI]: 19%-26%) in those with 0 to 3 points to 73% (95% CI: 66%-79%) in those with 9 or 10 points, corresponding to attributable fraction estimates of approximately 0% to 90%. Kaplan-Meier estimated stroke/TIA 2-year recurrence rates decreased from 20% (95% CI: 12%-28%) in the lowest Risk of Paradoxical Embolism score stratum to 2% (95% CI: 0%-4%) in the highest. CONCLUSION Clinical characteristics identify CS patients who vary markedly in PFO prevalence, reflecting clinically important variation in the probability that a discovered PFO is likely to be stroke-related vs incidental. Patients in strata more likely to have stroke-related PFOs have lower recurrence risk.

Thyroid antibody status, subclinical hypothyroidism, and the risk of coronary heart disease: an individual participant data analysis.

CONTEXT Subclinical hypothyroidism has been associated with increased risk of coronary heart disease (CHD), particularly with thyrotropin levels of 10.0 mIU/L or greater. The measurement of thyroid antibodies helps predict the progression to overt hypothyroidism, but it is unclear whether thyroid autoimmunity independently affects CHD risk. OBJECTIVE The objective of the study was to compare the CHD risk of subclinical hypothyroidism with and without thyroid peroxidase antibodies (TPOAbs). DATA SOURCES AND STUDY SELECTION A MEDLINE and EMBASE search from 1950 to 2011 was conducted for prospective cohorts, reporting baseline thyroid function, antibodies, and CHD outcomes. DATA EXTRACTION Individual data of 38 274 participants from six cohorts for CHD mortality followed up for 460 333 person-years and 33 394 participants from four cohorts for CHD events. DATA SYNTHESIS Among 38 274 adults (median age 55 y, 63% women), 1691 (4.4%) had subclinical hypothyroidism, of whom 775 (45.8%) had positive TPOAbs. During follow-up, 1436 participants died of CHD and 3285 had CHD events. Compared with euthyroid individuals, age- and gender-adjusted risks of CHD mortality in subclinical hypothyroidism were similar among individuals with and without TPOAbs [hazard ratio (HR) 1.15, 95% confidence interval (CI) 0.87-1.53 vs HR 1.26, CI 1.01-1.58, P for interaction = .62], as were risks of CHD events (HR 1.16, CI 0.87-1.56 vs HR 1.26, CI 1.02-1.56, P for interaction = .65). Risks of CHD mortality and events increased with higher thyrotropin, but within each stratum, risks did not differ by TPOAb status. CONCLUSIONS CHD risk associated with subclinical hypothyroidism did not differ by TPOAb status, suggesting that biomarkers of thyroid autoimmunity do not add independent prognostic information for CHD outcomes.

Pododermatitis in Captive and Free-Ranging Greater Flamingos (Phoenicopterus roseus)

Pododermatitis is frequent in captive flamingos worldwide, but little is known about the associated histopathologic lesions. Involvement of a papillomavirus or herpesvirus has been suspected. Histopathologic evaluation and viral assessment of biopsies from 19 live and 10 dead captive greater flamingos were performed. Selected samples were further examined by transmission electron microscopy and immunohistochemistry. Feet from 10 dead free-ranging greater flamingos were also evaluated. The histologic appearance of lesions of flamingos of increasing age was interpreted as the progression of pododermatitis. Mild histologic lesions were seen in a 3-week-old flamingo chick with no macroscopic lesions, and these were characterized by Micrococcus-like bacteria in the stratum corneum associated with exocytosis of heterophils. The inflammation associated with these bacteria may lead to further histologic changes: irregular columnar proliferations, papillary squirting, and dyskeratosis. In more chronic lesions, hydropic degeneration of keratinocytes, epidermal hyperplasia, and dyskeratosis were seen at the epidermis, as well as proliferation of new blood vessels and increased intercellular matrix in the dermis. Papillomavirus DNA was not identified in any of the samples, while herpesvirus DNA was seen only in a few cases; therefore, these viruses were not thought to be the cause of the lesions. Poor skin health through suboptimal husbandry may weaken the epidermal barrier and predispose the skin to invasion of Micrococcus-like bacteria. Histologic lesions were identified in very young flamingos with no macroscopic lesions; this is likely to be an early stage lesion that may progress to macroscopic lesions.

Bevacizumab, Pemetrexed, and Cisplatin, or Bevacizumab and Erlotinib for Patients With Advanced Non-Small-Cell Lung Cancer Stratified by Epidermal Growth Factor Receptor Mutation: Phase II Trial SAKK19/09

Treatment allocation by epidermal growth factor receptor mutation status is a new standard in patients with metastatic nonesmall-cell lung cancer. Yet, relatively few modern chemotherapy trials were conducted in patients characterized by epidermal growth factor receptor wild type. We describe the results of a multicenter phase II trial, testing in parallel 2 novel combination therapies, predefined molecular markers, and tumor rebiopsy at progression. Objective: The goal was to demonstrate that tailored therapy, according to tumor histology and epidermal growth factor receptor (EGFR) mutation status, and the introduction of novel drug combinations in the treatment of advanced nonesmall-cell lung cancer are promising for further investigation. Methods: We conducted a multicenter phase II trial with mandatory EGFR testing and 2 strata. Patients with EGFR wild type received 4 cycles of bevacizumab, pemetrexed, and cisplatin, followed by maintenance with bevacizumab and pemetrexed until progression. Patients with EGFR mutations received bevacizumab and erlotinib until progression. Patients had computed tomography scans every 6 weeks and repeat biopsy at progression. The primary end point was progression-free survival (PFS) ≥ 35% at 6 months in stratum EGFR wild type; 77 patients were required to reach a power of 90% with an alpha of 5%. Secondary end points were median PFS, overall survival, best overall response rate (ORR), and tolerability. Further biomarkers and biopsy at progression were also evaluated. Results: A total of 77 evaluable patients with EGFR wild type received an average of 9 cycles (range, 1-25). PFS at 6 months was 45.5%, median PFS was 6.9 months, overall survival was 12.1 months, and ORR was 62%. Kirsten rat sarcoma oncogene mutations and circulating vascular endothelial growth factor negatively correlated with survival, but thymidylate synthase expression did not. A total of 20 patients with EGFR mutations received an average of 16.

THE SITE AND MECHANISM OF ACTION OF BENZODIAZEPINE IN RABBIT HIPPOCAMPUS

The purpose of this study was to determine the effects of benzodiazepine in area CA1 of the hippocampus, and to explore possible mechanisms of action for these agents in this brain area. Two distinctly different benzodiazepine-induced changes in hippocampal physiology have been identified. First, benzodiazepine depresses the population spike recorded in stratum pyramidale, indicating a decrease in action potential generation. Second, benzodiazepine decreases the magnitude of post-tetanic potentiation of the population EPSP recorded in stratum radiatum, and shortens the duration. The effect of benzodiazepine on pyramidal cell excitation was reversed by the GABA antagonis bicuculline, and mimicked by GABA itself. Thus the available evidence is consistent with the hypothesis that benzodiazepine acts by enhancing the effect of GABA in this area. In stratum radiatum, on the other hand, the effect of benzodiazepine on post tetanic potentiation of the population EPSP was not altered by bicuculline although bicuculline did antagonize GABA in this area. In addition, application of GABA, while it caused profound changes in the population EPSP,p, did not cause the same changes that were induced by benzodiazepine. Thus the evidence does not support the hypothesis that benzodiazepine is acting in stratum radiatum by enhancing the effects of GABA. ^

Association between glycemic index and glycemic load and the risk of incident coronary heart disease among whites and African Americans with and without type 2 diabetes: The Atherosclerosis Risk in Communities study

Several studies have examined the association between high glycemic index (GI) and glycemic load (GL) diets and the risk for coronary heart disease (CHD). However, most of these studies were conducted primarily on white populations. The primary aim of this study was to examine whether high GI and GL diets are associated with increased risk for developing CHD in whites and African Americans, non-diabetics and diabetics, and within stratifications of body mass index (BMI) and hypertension (HTN). Baseline and 17-year follow-up data from ARIC (Atherosclerosis Risk in Communities) study was used. The study population (13,051) consisted of 74% whites, 26% African Americans, 89% non-diabetics, 11% diabetics, 43% male, 57% female aged 44 to 66 years at baseline. Data from the ARIC food frequency questionnaire at baseline were analyzed to provide GI and GL indices for each subject. Increases of 25 and 30 units for GI and GL respectively were used to describe relationships on incident CHD risk. Adjusted hazard ratios for propensity score with 95% confidence intervals (CI) were used to assess associations. During 17 years of follow-up (1987 to 2004), 1,683 cases of CHD was recorded. Glycemic index was associated with 2.12 fold (95% CI: 1.05, 4.30) increased incident CHD risk for all African Americans and GL was associated with 1.14 fold (95% CI: 1.04, 1.25) increased CHD risk for all whites. In addition, GL was also an important CHD risk factor for white non-diabetics (HR=1.59; 95% CI: 1.33, 1.90). Furthermore, within stratum of BMI 23.0 to 29.9 in non-diabetics, GI was associated with an increased hazard ratio of 11.99 (95% CI: 2.31, 62.18) for CHD in African Americans, and GL was associated with 1.23 fold (1.08, 1.39) increased CHD risk in whites. Body mass index modified the effect of GI and GL on CHD risk in all whites and white non-diabetics. For HTN, both systolic blood pressure and diastolic blood pressure modified the effect on GI and GL on CHD risk in all whites and African Americans, white and African American non-diabetics, and white diabetics. Further studies should examine other factors that could influence the effects of GI and GL on CHD risk, including dietary factors, physical activity, and diet-gene interactions. ^

Information bias and lifetime mortality risks of radiation-induced cancer

Additive and multiplicative models of relative risk were used to measure the effect of cancer misclassification and DS86 random errors on lifetime risk projections in the Life Span Study (LSS) of Hiroshima and Nagasaki atomic bomb survivors. The true number of cancer deaths in each stratum of the cancer mortality cross-classification was estimated using sufficient statistics from the EM algorithm. Average survivor doses in the strata were corrected for DS86 random error ($\sigma$ = 0.45) by use of reduction factors. Poisson regression was used to model the corrected and uncorrected mortality rates with covariates for age at-time-of-bombing, age at-time-of-death and gender. Excess risks were in good agreement with risks in RERF Report 11 (Part 2) and the BEIR-V report. Bias due to DS86 random error typically ranged from $-$15% to $-$30% for both sexes, and all sites and models. The total bias, including diagnostic misclassification, of excess risk of nonleukemia for exposure to 1 Sv from age 18 to 65 under the non-constant relative projection model was $-$37.1% for males and $-$23.3% for females. Total excess risks of leukemia under the relative projection model were biased $-$27.1% for males and $-$43.4% for females. Thus, nonleukemia risks for 1 Sv from ages 18 to 85 (DRREF = 2) increased from 1.91%/Sv to 2.68%/Sv among males and from 3.23%/Sv to 4.02%/Sv among females. Leukemia excess risks increased from 0.87%/Sv to 1.10%/Sv among males and from 0.73%/Sv to 1.04%/Sv among females. Bias was dependent on the gender, site, correction method, exposure profile and projection model considered. Future studies that use LSS data for U.S. nuclear workers may be downwardly biased if lifetime risk projections are not adjusted for random and systematic errors. (Supported by U.S. NRC Grant NRC-04-091-02.) ^

THE RELATIONSHIP BETWEEN DEGREE OF BLOOD PRESSURE REDUCTION AND MORTALITY AMONG HYPERTENSIVES IN THE HYPERTENSION DETECTION AND FOLLOW-UP PROGRAM

The relationship between degree of diastolic blood pressure (DBP) reduction and mortality was examined among hypertensives, ages 30-69, in the Hypertension Detection and Follow-up Program (HDFP). The HDFP was a multi-center community-based trial, which followed 10,940 hypertensive participants for five years. One-year survival was required for inclusion in this investigation since the one-year annual visit was the first occasion where change in blood pressure could be measured on all participants. During the subsequent four years of follow-up on 10,052 participants, 568 deaths occurred. For levels of change in DBP and for categories of variables related to mortality, the crude mortality rate was calculated. Time-dependent life tables were also calculated so as to utilize available blood pressure data over time. In addition, the Cox life table regression model, extended to take into account both time-constant and time-dependent covariates, was used to examine the relationship change in blood pressure over time and mortality.^ The results of the time-dependent life table and time-dependent Cox life table regression analyses supported the existence of a quadratic function which modeled the relationship between DBP reduction and mortality, even after adjusting for other risk factors. The minimum mortality hazard ratio, based on a particular model, occurred at a DBP reduction of 22.6 mm Hg (standard error = 10.6) in the whole population and 8.5 mm Hg (standard error = 4.6) in the baseline DBP stratum 90-104. After this reduction, there was a small increase in the risk of death. There was not evidence of the quadratic function after fitting the same model using systolic blood pressure. Methodologic issues involved in studying a particular degree of blood pressure reduction were considered. The confidence interval around the change corresponding to the minimum hazard ratio was wide and the obtained blood pressure level should not be interpreted as a goal for treatment. Blood pressure reduction was attributed, not only to pharmacologic therapy, but also to regression to the mean, and to other unknown factors unrelated to treatment. Therefore, the surprising results of this study do not provide direct implications for treatment, but strongly suggest replication in other populations. ^

REGRESSION ANALYSIS FOR STANDARDIZED MORTALITY RATIOS WITH APPLICATION TO OCCUPATIONAL FOLLOW-UP STUDIES

Traditional comparison of standardized mortality ratios (SMRs) can be misleading if the age-specific mortality ratios are not homogeneous. For this reason, a regression model has been developed which incorporates the mortality ratio as a function of age. This model is then applied to mortality data from an occupational cohort study. The nature of the occupational data necessitates the investigation of mortality ratios which increase with age. These occupational data are used primarily to illustrate and develop the statistical methodology.^ The age-specific mortality ratio (MR) for the covariates of interest can be written as MR(,ij...m) = ((mu)(,ij...m)/(theta)(,ij...m)) = r(.)exp (Z('')(,ij...m)(beta)) where (mu)(,ij...m) and (theta)(,ij...m) denote the force of mortality in the study and chosen standard populations in the ij...m('th) stratum, respectively, r is the intercept, Z(,ij...m) is the vector of covariables associated with the i('th) age interval, and (beta) is a vector of regression coefficients associated with these covariables. A Newton-Raphson iterative procedure has been used for determining the maximum likelihood estimates of the regression coefficients.^ This model provides a statistical method for a logical and easily interpretable explanation of an occupational cohort mortality experience. Since it gives a reasonable fit to the mortality data, it can also be concluded that the model is fairly realistic. The traditional statistical method for the analysis of occupational cohort mortality data is to present a summary index such as the SMR under the assumption of constant (homogeneous) age-specific mortality ratios. Since the mortality ratios for occupational groups usually increase with age, the homogeneity assumption of the age-specific mortality ratios is often untenable. The traditional method of comparing SMRs under the homogeneity assumption is a special case of this model, without age as a covariate.^ This model also provides a statistical technique to evaluate the relative risk between two SMRs or a dose-response relationship among several SMRs. The model presented has application in the medical, demographic and epidemiologic areas. The methods developed in this thesis are suitable for future analyses of mortality or morbidity data when the age-specific mortality/morbidity experience is a function of age or when there is an interaction effect between confounding variables needs to be evaluated. ^