Biomedical subjects

Alan E Hubbard

Publications and source records attributed to Alan E Hubbard.

12 recordsLinked to original sources

PubMed2006

Microarray analysis of gene expression in peripheral blood mononuclear cells from dioxin-exposed human subjects.

Tetrachlorodibenzo-p-dioxin (TCDD) is classified as a human carcinogen and exerts toxic effects on the skin (chloracne). Effects on reproductive, immunological, and endocrine systems have also been observed in animal models. TCDD acts through the aryl hydrocarbon receptor (AhR) pathway influencing largely unknown gene networks. An industrial accident in Seveso, Italy in 1976 exposed thousands of people to substantial quantities of TCDD. Twenty years after the exposure, this study examines global gene expression in the mononuclear cells of 26 Seveso female never smokers, with similar age, alcohol consumption, use of medications, and background plasma levels of 22 dioxin congeners unrelated to the Seveso accident. Plasma dioxin levels were still elevated in the exposed subjects. We performed analyses in two different comparison groups. The first included high-exposed study subjects compared with individuals with background TCDD levels (average plasma levels 99.4 and 6.7ppt, respectively); the second compared subjects who developed chloracne after the accident, and those who did not develop this disease. Overall, we observed a modest alteration of gene expression based on dioxin levels or on chloracne status. In the comparison between high levels and background levels of TCDD, four histone genes were up-regulated and modified expression of HIST1H3H was confirmed by real-time PCR. In the comparison between chloracne case-control subjects, five hemoglobin genes were up-regulated. Pathway analysis revealed two major networks for each comparison, involving cell proliferation, apoptosis, immunological and hematological disease, and other pathways. Further examination of the role of these genes in dioxin induced-toxicity is warranted.

Acneiform Eruptions↗

PubMed2006

The histone deacetylase inhibitor trichostatin a has genotoxic effects in human lymphoblasts in vitro.

Histone deacetylase inhibitors (HDACi) are a class of putative chemotherapeutic agents for which the mechanism of toxicity has not been fully identified. To explore the possibility that HDACi are genotoxic, human TK6 lymphoblastoid cells were exposed to trichostatin A (TSA) and genetic damage was measured. TSA caused a dose-dependent increase of G1-arrested cells at 24 h that correlated with increasing levels of p21 and apoptosis. Significantly elevated frequencies of structural chromosomal aberrations in cells exposed to TSA were observed using both the kinetochore-antibody micronucleus assay and nonbanding metaphase chromosome analysis. Increased tail intensities, indicative of elevated levels of DNA damage, were observed using the alkaline comet assay. Elevated levels of phosphorylated histone gammaH2AX protein were observed as early as 3 h following TSA exposure and peaked at 12 h for 200nM TSA. Significant levels of aneuploidy at the 200nM TSA dose were observed using metaphase analysis, but interestingly, kinetochore-positive micronuclei were not detected at any dose using the kinetochore micronucleus assay, suggesting that TSA induces aneuploidy via a nondisjunction event rather than chromosome lagging. Increases in chromosomal loss and breakage were observed using simultaneous FISH metaphase analysis of chromosomes 5, 7, 8, and 21, consistent with data obtained from the micronucleus and metaphase chromosome analyses. We conclude that TSA is both a clastogen and aneugen in the TK6 cell line and propose that the observed cytostatic and apoptotic properties of TSA may partially be due to this genotoxicity.

Aneuploidy↗

PubMed2006

Quantile-function based null distribution in resampling based multiple testing.

Simultaneously testing a collection of null hypotheses about a data generating distribution based on a sample of independent and identically distributed observations is a fundamental and important statistical problem involving many applications. Methods based on marginal null distributions (i.e., marginal p-values) are attractive since the marginal p-values can be based on a user supplied choice of marginal null distributions and they are computationally trivial, but they, by necessity, are known to either be conservative or to rely on assumptions about the dependence structure between the test-statistics. Re-sampling based multiple testing (Westfall and Young, 1993) involves sampling from a joint null distribution of the test-statistics, and controlling (possibly in a, for example, step-down fashion) the user supplied type-I error rate under this joint null distribution for the test-statistics. A generally asymptotically valid null distribution avoiding the need for the subset pivotality condition for the vector of test-statistics was proposed in Pollard, van der Laan (2003) for null hypotheses about general real valued parameters. This null distribution was generalized in Dudoit, vanderLaan, Pollard (2004) to general null hypotheses and test-statistics. In ongoing recent work van der Laan, Hubbard (2005), we propose a new generally asymptotically valid null distribution for the test-statistics and a corresponding bootstrap estimate, whose marginal distributions are user supplied, and can thus be set equal to the (most powerful) marginal null distributions one would use in univariate testing to obtain a p-value. Previous proposed null distributions either relied on a restrictive subset pivotality condition (Westfall and Young) or did not guarantee this latter property (Dudoit, vanderLaan, Pollard, 2004). It is argued and illustrated that the resulting new re-sampling based multiple testing methods provide more accurate control of the wished Type-I error in finite samples and are more powerful. We establish formal results and investigate the practical performance of this methodology in a simulation and data analysis.

Chi-Square Distribution↗

PubMed2006

Issues of processing and multiple testing of SELDI-TOF MS proteomic data.

A new data filtering method for SELDI-TOF MS proteomic spectra data is described. We examined technical repeats (2 per subject) of intensity versus m/z (mass/charge) of bone marrow cell lysate for two groups of childhood leukemia patients: acute myeloid leukemia (AML) and acute lymphoblastic leukemia (ALL). As others have noted, the type of data processing as well as experimental variability can have a disproportionate impact on the list of "interesting'' proteins (see Baggerly et al. (2004)). We propose a list of processing and multiple testing techniques to correct for 1) background drift; 2) filtering using smooth regression and cross-validated bandwidth selection; 3) peak finding; and 4) methods to correct for multiple testing (van der Laan et al. (2005)). The result is a list of proteins (indexed by m/z) where average expression is significantly different among disease (or treatment, etc.) groups. The procedures are intended to provide a sensible and statistically driven algorithm, which we argue provides a list of proteins that have a significant difference in expression. Given no sources of unmeasured bias (such as confounding of experimental conditions with disease status), proteins found to be statistically significant using this technique have a low probability of being false positives.

Acute Disease↗

PubMed2005

Empirical Bayes and resampling based multiple testing procedure controlling tail probability of the proportion of false positives.

Simultaneously testing a collection of null hypotheses about a data generating distribution based on a sample of independent and identically distributed observations is a fundamental and important statistical problem involving many applications. In this article we propose a new re-sampling based multiple testing procedure asymptotically controlling the probability that the proportion of false positives among the set of rejections exceeds q at level alpha, where q and alpha are user supplied numbers. The procedure involves 1) specifying a conditional distribution for a guessed set of true null hypotheses, given the data, which asymptotically is degenerate at the true set of null hypotheses, and 2) specifying a generally valid null distribution for the vector of test-statistics proposed in Pollard & van der Laan (2003), and generalized in our subsequent article Dudoit, van der Laan, & Pollard (2004), van der Laan, Dudoit, & Pollard (2004), and van der Laan, Dudoit, & Pollard (2004b). Ingredient 1) is established by fitting the empirical Bayes two component mixture model (Efron (2001b)) to the data to obtain an upper bound for marginal posterior probabilities of the null being true, given the data. We establish the finite sample rational behind our proposal, and prove that this new multiple testing procedure asymptotically controls the wished tail probability for the proportion of false positives under general data generating distributions. In addition, we provide simulation studies establishing that this method is generally more powerful in finite samples than our previously proposed augmentation multiple testing procedure (van der Laan, Dudoit, & Pollard (2004b)) and competing procedures from the literature. Finally, we illustrate our methodology with a data analysis.

Journal Article↗

PubMed2005

Use of OctoChrome fluorescence in situ hybridization to detect specific aneuploidy among all 24 chromosomes in benzene-exposed workers.

Benzene is an established human leukemogen. The mechanism of benzene-induced leukemogenesis, however, remains unclear, but chromosomal damage is thought to play a critical role. We previously reported that the loss of chromosomes 5 and 7 (monosomy 5 and 7) and the gain of chromosomes 8 and 21 (trisomy 8 and 21) are significantly increased in benzene-exposed workers in comparison to matched controls. To determine if selective effects of benzene can occur, we employed three-color painting on an 8-square slide to screen numerical changes in all 24 human chromosomes (OctoChrome FISH) in a pilot study of 11 subjects (6 exposed to >5 ppm benzene and 5 age- and sex-matched controls). The effects of benzene on each chromosome were assessed as the incidence rate ratio (IRR) from a Poisson regression model with the strongest effects being reflected by the highest IRR values. Monosomy of chromosomes 5, 6, 7 and 10 had the highest IRRs and statistical significance in this preliminary study (IRR>2.5, p<0.01). On the other hand, the monosomy levels of six other chromosomes (1, 4, 9, 11, 22 and Y) were unchanged in the exposed workers with IRRs close to 1.0. Similarly, selective effects were also observed on trisomy induction with chromosomes 8, 9, 17, 21 and 22 (IRR>2.5, p<0.01). These results suggest that benzene has the capability of producing selective effects on certain chromosomes, which is supported by our in vitro findings showing that chromosomes 5 and 7 are more sensitive to loss than other chromosomes following exposure to benzene metabolites. We are currently investigating potential mechanisms for this induction of selective aneuploidy.

Adult↗

PubMed2005

Use of 'Omic' technologies to study humans exposed to benzene.

'Omic' technologies include genomics, transcriptomics (gene expression profiling), proteomics and metabolomics. We are utilizing these new technologies in an effort to develop novel biomarkers of exposure, susceptibility and response to benzene. Advances in genomics allow one to study hundreds to thousands of single nucleotide polymorphisms simultaneously on small quantities of DNA using array-based technologies. We are currently utilizing these technologies to examine genetic variation in pathways relating to biotransformation, DNA repair, folate metabolism and immune response with the goal of finding biomarkers of susceptibility to benzene hematotoxicity. Transcriptomics is used to measure the full complement of activated genes, mRNAs or transcripts in a particular tissue at a particular time typically using microarray technology. We have applied microarrays to the study of global gene expression in the peripheral blood cells of benzene-exposed workers. More than 100 genes were identified as being potentially differentially expressed, with genes related to apoptosis and immune function being the most significantly affected. Initial studies employing proteomics have also shown that several proteins are altered in the serum of exposed compared to control subjects and these proteins are potential biomarkers of benzene exposure. Omic technologies therefore have significant potential in generating novel biomarkers of exposure, susceptibility and response to benzene.

Air Pollutants, Occupational↗

PubMed2005

Nonrandom aneuploidy of chromosomes 1, 5, 6, 7, 8, 9, 11, 12, and 21 induced by the benzene metabolites hydroquinone and benzenetriol.

The loss and gain of whole chromosomes (aneuploidy) is common in the development of leukemia and other cancers. In acute myeloid leukemia, the loss (monosomy) of chromosomes 5 and 7 and the gain (trisomy) of chromosome 8 are common clonal chromosomal abnormalities. Here, we have tested the hypothesis that metabolites of the human leukemogen benzene cause a higher rate of gain and loss among the chromosomes involved in leukemogenesis and, as such, are nonrandom and selective in their effects. Human peripheral blood was exposed to two metabolites of benzene, namely, hydroquinone (HQ) and benzenetriol (BT), and the ploidy status of nine different chromosomes (1, 5, 6, 7, 8, 9, 11, 12, and 21) was examined using fluorescence in situ hybridization of metaphase spreads. Poisson regression was used to provide interpretable incidence rate ratios and corresponding P values for all nine chromosomes. Statistically significant differences were found between the sensitivity of the nine chromosomes to gain or loss. Chromosomes 5 and 7 were highly sensitive to loss following HQ and BT exposure, whereas chromosomes 7, 8, and 21 were highly sensitive to gain in comparison to other chromosomes. Significant support for the a priori hypothesis that chromosomes 5 and 7 are more sensitive to loss induced by HQ and BT than the other seven chromosomes was also obtained. These data support the notion that benzene metabolites affect the ploidy status of specific chromosomes more than others and may initiate or promote leukemia induction through these specific effects.

Acute Disease↗

PubMed2005

Polymorphisms in ghrelin and neuropeptide Y genes are associated with non-Hodgkin lymphoma.

We previously reported a positive association among body mass index, single nucleotide polymorphisms (SNP) in the leptin and leptin receptor genes that are involved in body weight regulation, and non-Hodgkin lymphoma (NHL). Polymorphisms in the ghrelin (GHRL) and neuropeptide Y (NPY) genes were examined in the same population-based case-control study of NHL to further explore the role of genes involved in energy homeostasis and obesity in susceptibility to NHL. Ghrelin is an orexigenic hormone that induces NPY release and inhibits proinflammatory cytokines via its antagonistic relationship with leptin. NPY is a potent appetite stimulator controlled by ghrelin and leptin and also acts as a mediator of immune function. DNA from 458 cases and 812 controls was genotyped. Among genotyped GHRL SNPs, the variant allele for GHRL -4427G>A was inversely associated with all NHL [odds ratios (OR), 0.78; 95% confidence interval (95% CI), 0.59-1.0] and more specifically with diffuse large cell lymphoma (DLCL; homozygous variant: OR, 0.31; 95% CI, 0.13-0.74). Another SNP, GHRL 5179A>G, decreased the risk of DLCL (homozygous variant: OR, 0.35; 95% CI, 0.10-1.2). NPY -485T>C, 1258G>A, and 5671C>T were in total linkage disequilibrium (D' = 0.99) and the homozygous variants were associated with an increased risk of NHL in NPY SNPs -485T>C (OR, 1.7; 95% CI, 1.1-2.5), 1258G>A (OR, 1.7; 95% CI, 1.1-2.5), and 5671C>T (OR, 1.9; 95% CI, 1.3-2.8). When stratified by subtype, the variant allele for NPY 1128T>C was positively associated with follicular lymphoma (OR, 2.3; 95% CI, 1.1-4.9) as were homozygous variants for NPY SNPs -485T>C (OR, 2.4; 95% CI, 1.3-4.4), 1258G>A (OR, 2.0; 95% CI, 1.1-3.5), and 5671C>T (OR, 1.8; 95% CI, 1.1-3.0). These results add further support for the hypothesis that SNPs in energy-regulating genes affect risk of NHL.

Adult↗

PubMed2005

Discovery of novel biomarkers by microarray analysis of peripheral blood mononuclear cell gene expression in benzene-exposed workers.

Benzene is an industrial chemical and component of gasoline that is an established cause of leukemia. To better understand the risk benzene poses, we examined the effect of benzene exposure on peripheral blood mononuclear cell (PBMC) gene expression in a population of shoe-factory workers with well-characterized occupational exposures using microarrays and real-time polymerase chain reaction (PCR). PBMC RNA was stabilized in the field and analyzed using a comprehensive human array, the U133A/B Affymetrix GeneChip set. A matched analysis of six exposed-control pairs was performed. A combination of robust multiarray analysis and ordering of genes using paired t-statistics, along with bootstrapping to control for a 5% familywise error rate, was used to identify differentially expressed genes in a global analysis. This resulted in a set of 29 known genes being identified that were highly likely to be differentially expressed. We also repeated these analyses on a smaller subset of 508 cytokine probe sets and found that the expression of 19 known cytokine genes was significantly different between the exposed and the control subjects. Six genes were selected for confirmation by real-time PCR, and of these, CXCL16, ZNF331, JUN, and PF4 were the most significantly affected by benzene exposure, a finding that was confirmed in a larger data set from 28 subjects. The altered expression was not caused by changes in the makeup of the PBMC fraction. Thus, microarray analysis along with real-time PCR confirmation reveals that altered expressions of CXCL16, ZNF331, JUN, and PF4 are potential biomarkers of benzene exposure.

Adult↗

PubMed2002

Statistical estimation of parameters in a disease transmission model: analysis of a Cryptosporidium outbreak.

Population dynamic models, commonly used tools in the study of epidemics and other complex population processes, are implicit non-linear mathematical equations. Inference based on such models can be difficult due to the problems associated with high dimensional parameters that may be non-identified and complex likelihood functions that are difficult to maximize. To address a problem of non-identifiability due to collinearity of parameter estimates in a mathematical model of the 1993 Milwaukee Cryptosporidium parvum outbreak, we examined the utility of a constrained profile likelihood approach. This method was used to study two parameters of interest from the mathematical model: (i). the rate of secondary transmission; (ii). the proportional increase in primary transmission due to water treatment failure. The estimated values of these parameters were shown to depend strongly on poorly understood aspects of Cryptosporidium epidemiology such as asymptomatic proportion and the population immune status. Our analysis demonstrated that the combination of a disease transmission model and a constrained profile likelihood procedure provides an effective approach for inference and estimation of important parameters regulating infectious disease outbreaks.

Animals↗

PubMed2002

Post-abortion depot medroxyprogesterone acetate continuation rates: a randomized trial of cyclic estradiol.

Irregular menstrual bleeding is a common reason for discontinuation of depot medroyprogesterone acetate (DMPA) contraception. We conducted a double-blind, randomized placebo controlled trial to estimate if transdermal estradiol used in a cyclic fashion could make DMPA bleeding patterns more acceptable and improve DMPA continuation rates. Women initiating DMPA contraception immediately post-abortion (N = 132) were randomized to receive either a low sustained dose of transdermal 17-beta estradiol (0.10 mg/day, Climara) or an identical placebo to be used in a cyclic manner for 3 months. Bleeding patterns, contraceptive use, and pregnancy rates were measured for 1 year after enrollment. Cyclic transdermal estrogen did not create regular bleeding patterns in new users of DMPA and had no effect on DMPA continuation rates (relative hazard = 0.93, 95%CI 0.65-1.33). By 12 months after enrollment, 19.4% of those randomized to estrogen and 23.9% of those randomized to placebo were still using DMPA for contraception. There was no difference in the frequency of DMPA discontinuation primarily for abnormal uterine bleeding between groups (relative risk ratio = 0.95, 95%CI 0.40-2.23). Few women who discontinued DMPA switched to highly effective methods of contraception and within 1 year of an abortion, there were 18 repeat pregnancies (21.9%). Compliance with the estrogen patches was poor resulting in early study termination and inadequate statistical power to detect a clinically significant difference in DMPA continuation rates. Continuation rates of DMPA were low and repeat pregnancy rates were high among women initiating DMPA at the time of abortion. Cyclic transdermal estradiol use is an unacceptable regimen for preventing the irregular bleeding associated with DMPA initiation.

Abortion, Induced↗