Imputation in data fusion of heterogeneous data sets a model-based numerical experiment

Given the very large amount of data obtained everyday through population surveys, much of the new research again could use this information instead of collecting new samples. Unfortunately, relevant data are often disseminated into different files obtained through different sampling designs. Data fusion is a set of methods used to combine information from different sources into a single dataset. In this article, we are interested in a specific problem: the fusion of two data files, one of which being quite small. We propose a model-based procedure combining a logistic regression with an Expectation-Maximization algorithm. Results show that despite the lack of data, this procedure can perform better than standard matching procedures.

25-hydroxyvitamin D concentration is inversely associated with serum MMP-9 in a cross-sectional study of African American ESRD patients.

BACKGROUND: Circulating 25-hydroxyvitamin D [25(OH)D] concentration is inversely associated with peripheral arterial disease and hypertension. Vascular remodeling may play a role in this association, however, data relating vitamin D level to specific remodeling biomarkers among ESRD patients is sparse. We tested whether 25(OH)D concentration is associated with markers of vascular remodeling and inflammation in African American ESRD patients.METHODS: We conducted a cross-sectional study among ESRD patients receiving maintenance hemodialysis within Emory University-affiliated outpatient hemodialysis units. Demographic, clinical and dialysis treatment data were collected via direct patient interview and review of patients records at the time of enrollment, and each patient gave blood samples. Associations between 25(OH)D and biomarker concentrations were estimated in univariate analyses using Pearson's correlation coefficients and in multivariate analyses using linear regression models. 25(OH) D concentration was entered in multivariate linear regression models as a continuous variable and binary variable (<15 ng/ml and =15 ng/ml). Adjusted estimate concentrations of biomarkers were compared between 25(OH) D groups using analysis of variance (ANOVA). Finally, results were stratified by vascular access type.RESULTS: Among 91 patients, mean (standard deviation) 25(OH)D concentration was 18.8 (9.6) ng/ml, and was low (<15 ng/ml) in 43% of patients. In univariate analyses, low 25(OH) D was associated with lower serum calcium, higher serum phosphorus, and higher LDL concentrations. 25(OH) D concentration was inversely correlated with MMP-9 concentration (r = -0.29, p = 0.004). In multivariate analyses, MMP-9 concentration remained negatively associated with 25(OH) D concentration (P = 0.03) and anti-inflammatory IL-10 concentration positively correlated with 25(OH) D concentration (P = 0.04).CONCLUSIONS: Plasma MMP-9 and circulating 25(OH) D concentrations are significantly and inversely associated among ESRD patients. This finding may suggest a potential mechanism by which low circulating 25(OH) D functions as a cardiovascular risk factor.

Direct Effect of Birth Weight on Childhood Blood Pressure: a Causal Mediation Analysis

Background: Numerous studies have shown a negative association between birth weight (BW) and blood pressure (BP) later in life. To estimate the direct effect of BW on BP, it is conventional to condition on current weight (CW). However, such conditioning can induce collider stratification bias in the estimate of the direct effect. Objective: To bound the potential bias due to U, an unmeasured common cause of CW and BP, on the estimate of the (controlled) direct effect of BW on BP. Methods: Data from a school based study in Switzerland were used (N = 4,005; 2,010 B/1,995 G; mean age: 12.3 yr [range: 10.1-14.9]). Measured common causes of BW-BP (SES, smoking, body weight, and hypertension status of the mother) and CW-BP (breastfeeding and child's physical activity and diet) were identified with DAGs. Linear regression models were fitted to estimate the association between BW and BP. Sensitivity analyses were conducted to assess the potential effect of U on the association between BW and BP. U was assumed 1) to be a binary variable that affected BP by the same magnitude in low BWand in normal BW children and 2) to have a different prevalence in low BW children and in normal BW children for a given CW. Results: A small negative association was observed between BW and BP [beta: -0.3 mmHg/kg (95% CI: -0.9 to 0.3)]. The association was strengthened upon conditioning for CW [beta: -1.5 mmHg/kg (95% CI: -2.1 to -0.9)]. Upon further conditioning on common causes of BW-BP and CW-BP, the association did not change substantially [beta: -1.4 mmHg/kg (95% CI: -2.0 to -0.8)]. The negative association could be explained by U only if U was strongly associated with BP and if there was a large difference in the prevalence of U between low BWand normal BW children. Conclusion: The observed negative association between BW and BP upon adjustment for CW was not easily explained by an unmeasured common cause of CWand BP.

A comparison of three methods of Mendelian randomization when the genetic instrument, the risk factor and the outcome are all binary.

The method of instrumental variable (referred to as Mendelian randomization when the instrument is a genetic variant) has been initially developed to infer on a causal effect of a risk factor on some outcome of interest in a linear model. Adapting this method to nonlinear models, however, is known to be problematic. In this paper, we consider the simple case when the genetic instrument, the risk factor, and the outcome are all binary. We compare via simulations the usual two-stages estimate of a causal odds-ratio and its adjusted version with a recently proposed estimate in the context of a clinical trial with noncompliance. In contrast to the former two, we confirm that the latter is (under some conditions) a valid estimate of a causal odds-ratio defined in the subpopulation of compliers, and we propose its use in the context of Mendelian randomization. By analogy with a clinical trial with noncompliance, compliers are those individuals for whom the presence/absence of the risk factor X is determined by the presence/absence of the genetic variant Z (i.e., for whom we would observe X = Z whatever the alleles randomly received at conception). We also recall and illustrate the huge variability of instrumental variable estimates when the instrument is weak (i.e., with a low percentage of compliers, as is typically the case with genetic instruments for which this proportion is frequently smaller than 10%) where the inter-quartile range of our simulated estimates was up to 18 times higher compared to a conventional (e.g., intention-to-treat) approach. We thus conclude that the need to find stronger instruments is probably as important as the need to develop a methodology allowing to consistently estimate a causal odds-ratio.

Variable viral clearance despite adequate ganciclovir plasma levels during valganciclovir treatment for cytomegalovirus disease in D+/R- transplant recipients.

ABSTRACT: BACKGROUND: Valganciclovir, the oral prodrug of ganciclovir, has been demonstrated equivalent to iv ganciclovir for CMV disease treatment in solid organ transplant recipients. Variability in ganciclovir exposure achieved with valganciclovir could be implicated as a contributing factor for explaining variations in the therapeutic response. This prospective observational study aimed to correlate clinical and cytomegalovirus (CMV) viral load response (DNAemia) with ganciclovir plasma concentrations in patients treated with valganciclovir for CMV infection/disease. METHODS: Seven CMV D+/R- transplant recipients (4 kidney, 2 liver and 1 heart) were treated with valganciclovir (initial dose was 900-1800 mg/day for 3-6.5 weeks, followed by 450-900 mg/day for 2-9 weeks). DNAemia was monitored by real time quantitative PCR and ganciclovir plasma concentration was measured at trough (Ctrough) and 3 h after drug administration (C3h) by HPLC. RESULTS: Four patients presented with CMV syndrome, two had CMV tissue-invasive disease after prophylaxis discontinuation, and one liver recipient was treated pre-emptively for asymptomatic rising CMV viral load 5 weeks post-transplantation in the absence of prophylaxis. CMV DNAemia decreased during the first week of treatment in all recipients except in one patient (median decrease: -1.2 log copies/mL, range: -1.8 to 0) despite satisfactory ganciclovir exposure (AUC0-12 = 48 mg.h/L, range for the 7 patients: 40-118 mg.h/L). Viral clearance was obtained in five patients after a median of time of 34 days (range: 28-82 days). Two patients had recurrent CMV disease despite adequate ganciclovir exposure (65 mg.h/L, range: 44-118 mg.h/L). CONCLUSIONS: Valganciclovir treatment for CMV infection/disease in D+/R- transplant recipients can thus result in variable viral clearance despite adequate ganciclovir plasma concentrations, probably correlating inversely with anti-CMV immune responses after primary infection.

Variable phenotypic expressivity in a Swiss family with autosomal dominant retinitis pigmentosa due to mutation T494m in the Prpf3 gene

Purpose:to describe the clinical features in a five generations family segregating autosomal dominant retinitis pigmentosa and to identify the causative gene Patient and Methods:Twenty five individuals of a large five-generation family originating from Western Switzerland were ascertained for phenotypic and genotypic characterization. Ophthalmologic evaluations included color vision testing, Goldman perimetry and digital fundus photography. Some patients had autofluorescence (AF) imaging, ocular coherence tomography (OCT) and ISCEV-standard full-field electroretinography (ERG). Blood samples were collected from 10 affected (4 to 70 years of age) and 15 unaffected members after informed consent. DNA was isolated and exons and intron-exons junctions of known adRP genes were sequenced using a Big Dye sequencing kit 1.1. Results:Age of onset of nightblindness and severity of progression of the disease was variable between members of the family. Some patients had early onset of nightblindess aged 3, others at mid-twenties. Most patients had visual acuity above 0.6 for the first 4 decades. Two older patients still had good vision (0.4) in their seventies. Myopia (range: -2 to -5) was noticed in most affected subjects. Fundus findings showed areas of atrophy along the arcades. The AF imaging showed a large high density ring bilaterally. A T494M change was found in exon 11 of PRPF3 gene. The change segregates with the disease in the family. Conclusion: A mutation in the PRPF3 gene is rare compared with other genes causing ADRP. Although a T494M change has been reported, our family is the first one with a variable expressivity. Mutations in PRPF3 gene can cause a variable phenotype of ADRP unlike the previously described Danish and English families. Our report gives a better understanding as to the phenotype/genotype description of ADRP due to PRPF3 mutation.

Two monoclonal rat antibodies with specificity for the beta-chain variable region V beta 6 of the murine T-cell receptor.

Two rat monoclonal antibodies (mAbs), 44-22-1 and 46-6B5, which recognize an alloreactive cytotoxic clone, 3F9, have been further tested on a panel of T hybridomas and cytotoxic T-cell clones for binding and functional activities. The mAbs recognized only those cells sharing the expression of the T-cell receptor beta-chain variable region gene V beta 6 with 3F9. All V beta 6+ cells were activated by these mAbs under cross-linking conditions and their antigen-specific activation was blocked by soluble mAb. Furthermore, depletion of 46-6B5+ normal lymph node T cells eliminated all cells expressing the epitope recognized by 44-22-1 and V beta 6 mRNA.

Predictors of nonresponse in a questionnaire-based outcome study of vocational rehabilitation patients.

OBJECTIVE: To identify predictors of nonresponse to a self-report study of patients with orthopedic trauma hospitalized for vocational rehabilitation between November 15, 2003, and December 31, 2005. The role of biopsychosocial complexity, assessed using the INTERMED, was of particular interest. DESIGN: Cohort study. Questionnaires with quality of life, sociodemographic, and job-related questions were given to patients at hospitalization and 1 year after discharge. Sociodemographic data, biopsychosocial complexity, and presence of comorbidity were available at hospitalization (baseline) for all eligible patients. Logistic regression models were used to test a number of baseline variables as potential predictors of nonresponse to the questionnaires at each of the 2 time points. SETTING: Rehabilitation clinic. PARTICIPANTS: Patients (N=990) hospitalized for vocational rehabilitation over a period of 2 years. INTERVENTIONS: Not applicable. MAIN OUTCOME MEASURE: Nonresponse to the questionnaires was the binary dependent variable. RESULTS: Patients with high biopsychosocial complexity, foreign native language, or low educational level were less likely to respond at both time points. Younger patients were less likely to respond at 1 year. Those living in a stable partnership were less likely than singles to respond at hospitalization. Sex, psychiatric, and somatic comorbidity and alcoholism were never associated with nonresponse. CONCLUSIONS: We stress the importance of assessing biopsychosocial complexity to predict nonresponse. Furthermore, the factors we found to be predictive of nonresponse are also known to influence treatment outcome and vocational rehabilitation. Therefore, it is important to increase the response rate of the groups of concern in order to reduce selection bias in epidemiologic investigations.

Immunodéficience commune variable: ce qu'il faut savoir [Common variable immune deficiency: what you need to know].

Common variable immune deficiency is the most frequent primary immune deficiency, characterized mainly by a disorder of B lymphocytes differentiation and a deficit in immunoglobulins. The clinical manifestations include recurrent infections, non-infectious lung and digestive involvements, autoimmune diseases, and an increased susceptibility to cancers. Recent breakthroughs have been made in the understanding of some genetic mechanisms of the disease. Replacement therapy with intravenous immunoglobulins remains the treatment of choice, which allows significant improvement in the survival and quality of life. However progress should be made in the understanding of the pathophysiology and in the early detection of this disease, since a delay in the diagnosis may have harmful consequences in terms of morbidity and mortality.

Assessing the Relationship between the Baseline Value of a Continuous Variable and Subsequent Change over Time

Analyzing the relationship between the baseline value and subsequent change of a continuous variable is a frequent matter of inquiry in cohort studies. These analyses are surprisingly complex, particularly if only two waves of data are available. It is unclear for non-biostatisticians where the complexity of this analysis lies and which statistical method is adequate.With the help of simulated longitudinal data of body mass index in children,we review statistical methods for the analysis of the association between the baseline value and subsequent change, assuming linear growth with time. Key issues in such analyses are mathematical coupling, measurement error, variability of change between individuals, and regression to the mean. Ideally, it is better to rely on multiple repeated measurements at different times and a linear random effects model is a standard approach if more than two waves of data are available. If only two waves of data are available, our simulations show that Blomqvist's method - which consists in adjusting for measurement error variance the estimated regression coefficient of observed change on baseline value - provides accurate estimates. The adequacy of the methods to assess the relationship between the baseline value and subsequent change depends on the number of data waves, the availability of information on measurement error, and the variability of change between individuals.

Binary classification of dyslipidemia from the waist-to-hip ratio and body mass index: a comparison of linear, logistic, and CART models.

BACKGROUND: We sought to improve upon previously published statistical modeling strategies for binary classification of dyslipidemia for general population screening purposes based on the waist-to-hip circumference ratio and body mass index anthropometric measurements. METHODS: Study subjects were participants in WHO-MONICA population-based surveys conducted in two Swiss regions. Outcome variables were based on the total serum cholesterol to high density lipoprotein cholesterol ratio. The other potential predictor variables were gender, age, current cigarette smoking, and hypertension. The models investigated were: (i) linear regression; (ii) logistic classification; (iii) regression trees; (iv) classification trees (iii and iv are collectively known as "CART"). Binary classification performance of the region-specific models was externally validated by classifying the subjects from the other region. RESULTS: Waist-to-hip circumference ratio and body mass index remained modest predictors of dyslipidemia. Correct classification rates for all models were 60-80%, with marked gender differences. Gender-specific models provided only small gains in classification. The external validations provided assurance about the stability of the models. CONCLUSIONS: There were no striking differences between either the algebraic (i, ii) vs. non-algebraic (iii, iv), or the regression (i, iii) vs. classification (ii, iv) modeling approaches. Anticipated advantages of the CART vs. simple additive linear and logistic models were less than expected in this particular application with a relatively small set of predictor variables. CART models may be more useful when considering main effects and interactions between larger sets of predictor variables.

Characteristics of eyes with secondary loss of visual acuity receiving variable dosing ranibizumab for neovascular age-related macular degeneration.

PURPOSE: The aim of this work is to investigate the characteristics of eyes failing to maintain visual acuity (VA) receiving variable dosing ranibizumab for neovascular age-related macular degeneration (nAMD) after three initial loading doses. METHODS: A consecutive series of patients with nAMD, who, after three loading doses of intravitreal ranibizumab (0.5 mg each), were re-treated for fluid seen on optical coherence tomography. After exclusion of eyes with previous treatment, follow-up less than 12 months, or missed visits, 99 patients were included in the analysis. The influence of baseline characteristics, initial VA response, and central retinal thickness (CRT) fluctuations on the VA stability from month 3 to month 24 were analyzed using subgroups and multiple regression analyses. RESULTS: Mean follow-up duration was 21.3 months (range 12-40 months, 32 patients followed-up for ≥24 months). Secondary loss of VA (loss of five letters or more) after month 3 was seen in 30 patients (mean VA improvement from baseline +5.8 letters at month 3, mean loss from baseline -5.3 letters at month 12 and -9.7 at final visit up to month 24), while 69 patients maintained vision (mean gain +8.9 letters at month 3, +10.4 letters at month 12, and +12.8 letters at final visit up to month 24). Secondary loss of VA was associated with the presence of pigment epithelial detachment (PED) at baseline (p 0.01), but not with baseline fibrosis/atrophy/hemorrhage, CRT fluctuations, or initial VA response. Chart analysis revealed additional individual explanations for the secondary loss of VA, including retinal pigment epithelial tears, progressive fibrosis, and atrophy. CONCLUSIONS: Tissue damage due to degeneration of PED, retinal pigment epithelial tears, progressive fibrosis, progressive atrophy, or massive hemorrhage, appears to be relevant in causing secondary loss of VA despite vascular endothelial growth factor suppression. PED at baseline may represent a risk factor.

Polymeric IgA is superior to monomeric IgA and IgG carrying the same variable domain in preventing Clostridium difficile toxin A damaging of T84 monolayers.

The two exotoxins A and B produced by Clostridium difficile are responsible for antibiotic-associated enterocolitis in human and animals. When added apically to human colonic carcinoma-derived T84 cell monolayers, toxin A, but not toxin B, abolished the transepithelial electrical resistance and altered the morphological integrity. Apical addition of suboptimal concentration of toxin A made the cell monolayer sensitive to toxin B. Both toxins induced drastic and rapid epithelial alterations when applied basolaterally with a complete disorganization of tight junctions and vacuolization of the cells. Toxin A-specific IgG2a from hybridoma PCG-4 added apically with toxin A alone or in combination with toxin B abolished the toxin-induced epithelial alterations for up to 8 h. The Ab neutralized basolateral toxin A for 4 h, but not the mixture of the two toxins. Using an identical Ab:Ag ratio, we found that recombinant polymeric IgA (IgAd/p) with the same Fv fragments extended protection against toxin A for at least 24 h in both compartments. In contrast, the recombinant monomeric IgA counterpart behaved as the PCG-4 IgG2a Ab. The direct comparison between different Ig isotype and molecular forms, but of unique specificity, demonstrates that IgAd/p Ab is more efficient in neutralizing toxin A than monomeric IgG and IgA. We conclude that immune protection against C. difficile toxins requires toxin A-specific secretory Abs in the intestinal lumen and IgAd/p specific for both toxins in the lamina propria.

Highly variable social organisation of colonies in the ant Formica cinerea.

Social organisation of colonies was examined in the ant Formica cinerea by estimating the coefficient of genetic relatedness among worker nest mates. The estimates based on microsatellite genotypes at three loci ranged from values close to zero to 0.61 across the populations studied in Finland. These results showed that a fundamental feature of colonies, the number of reproductive queens, varied greatly among the populations. Colonies in some populations had a single queen, whereas the nests could have a high number number of queens in other populations. There was a weak but non-significant correlation between the genetic and metric distance of nests within two populations with intermediate level of relatedness. Differentiation among nearby populations (within the dispersal distance of individuals) in one locality indicated limited dispersal or founder effects. This could occur when females are philopatric and stay in the natal polygynous colony which expands by building a network of nest galleries within a single habitat patch.