PubMed Health⌕ Search

Biomedical subjects

Lang Li

Publications and source records attributed to Lang Li.

At least 19 recordsLinked to original sources

ToxiVerse: chemical bioprofiling, toxicity data sharing and customizable predictive modeling.

MOTIVATION: Chemical toxicity assessment is critical for drug development and environmental safety. Computational models have emerged as a promising alternative to animal testing and now play a significant role in efficiently evaluating new chemicals. To address the urgent need for user-friendly machine learning tools in computational toxicology, we developed ToxiVerse, a public web-based platform. RESULTS: ToxiVerse provides automatic chemical bioprofiling, curated toxicity datasets, and a predictive modeling interface designed for researchers who lack programming expertise. The platform comprises three integrated modules: (i) Bioprofiler, which provides chemical descriptors by combining chemical-bioactivity data from PubChem assays with a machine learning-based data gap-filling procedure; (ii) Database, which hosts ∼50 000 curated chemicals covering diverse toxicity endpoints; and (iii) Cheminformatics, which enables dataset upload, chemical curation, and automatic generation of quantitative structure-activity relationship models for toxicity prediction. AVAILABILITY: The tool is accessible at www.toxiverse.com, and source code is available at https://github.com/zhu-research-group/toxiverse.

Quantitative Structure-Activity Relationship↗

Diverse gene expression and DNA methylation profiles correlate with differential adaptation of breast cancer cells to the antiestrogens tamoxifen and fulvestrant.

The development of targeted therapies for antiestrogen-resistant breast cancer requires a detailed understanding of its molecular characteristics. To further elucidate the molecular events underlying acquired resistance to the antiestrogens tamoxifen and fulvestrant, we established drug-resistant sublines from a single colony of hormone-dependent breast cancer MCF7 cells. These model systems allowed us to examine the cellular and molecular changes induced by antiestrogens in the context of a uniform clonal background. Global changes in both basal and estrogen-induced gene expression profiles were determined in hormone-sensitive and hormonal-resistant sublines using Affymetrix Human Genome U133 Plus 2.0 Arrays. Changes in DNA methylation were assessed by differential methylation hybridization, a high-throughput promoter CpG island microarray analysis. By comparative studies, we found distinct gene expression and promoter DNA methylation profiles associated with acquired resistance to fulvestrant versus tamoxifen. Fulvestrant resistance was characterized by pronounced up-regulation of multiple growth-stimulatory pathways, resulting in estrogen receptor alpha (ERalpha)-independent, autocrine-regulated proliferation. Conversely, acquired resistance to tamoxifen correlated with maintenance of the ERalpha-positive phenotype, although receptor-mediated gene regulation was altered. Activation of growth-promoting genes, due to promoter hypomethylation, was more frequently observed in antiestrogen-resistant cells compared with gene inactivation by promoter hypermethylation, revealing an unexpected insight into the molecular changes associated with endocrine resistance. In summary, this study provides an in-depth understanding of the molecular changes specific to acquired resistance to clinically important antiestrogens. Such knowledge of resistance-associated mechanisms could allow for identification of therapy targets and strategies for resensitization to these well-established antihormonal agents.

Breast Neoplasms↗

Hoosier Oncology Group randomized phase II study of docetaxel, vinorelbine, and estramustine in combination in hormone-refractory prostate cancer with pharmacogenetic survival analysis.

PURPOSE: To determine the safety and efficacy of two docetaxel doublets in hormone-refractory prostate cancer (HRPC) patients and to examine the prognostic role of polymorphisms in host genes important to docetaxel metabolism and transport. EXPERIMENTAL DESIGN: Sixty-four chemotherapy-naive patients with HRPC were randomized to docetaxel and vinorelbine (D, 20 mg/m2 i.v. days 1 and 8; V, 25 mg/m2 i.v. days 1 and 8) or docetaxel and estramustine phosphate (D, 60-70 mg/m2 i.v. day 1; E, 280 mg oral thrice daily days 1-5) administered q21d. Primary end point was clinically significant toxicity. A pharmacogenetic analysis of host genes was done in patients who received at least one cycle of docetaxel therapy. RESULTS: Grade 3/4 toxicity occurred in 15.6% of DV patients and in 28.6% DE patients. Neither arm exceeded the threshold of clinically significant toxicity. In the DV arm, objective response rate was 33%, prostate-specific antigen response rate was 20%, and median survival was 16.2 months. In the DE arm, objective response rate was 67%, prostate-specific antigen response rate was 43%, and median survival was 19.7 months. Pharmacogenetic analyses showed a significant association between survival beyond 15 months and the ABCG2 421 C > A (Q141K) polymorphism compared with the wild-type (C/C) genotype (66% versus 27%; P = 0.05). CONCLUSIONS: DV and DE doublets are active with a tolerable toxicity profile in patients with HRPC; however, efficacy does not seem superior to standard single-agent docetaxel. The ABCG2 421 C > A (Q141K) polymorphism may be an important predictor of response and survival in HRPC patients treated with docetaxel-based chemotherapy.

Adult↗

The amino terminus of the human multidrug resistance transporter ABCC1 has a U-shaped folding with a gating function.

Multidrug resistance is a serious problem in successful cancer chemotherapy. Studies using model cell lines have demonstrated that overexpression of some members of the ATP-binding cassette (ABC) transporter superfamily, such as ABCC1, causes enhanced efflux and, thus, decreased accumulation of multiple anticancer drugs, which leads to increased cell survival. Unlike most other ABC transporters, ABCC1 has an additional membrane-spanning domain (MSD0) with a putative extracellular amino terminus of 32 amino acids. However, the function of MSD0 and the role of the extracellular amino terminus are largely unknown. In this study, we examined the structural folding and the function of the amino terminus. We found that it has a U-shaped folding with the bottom of the U-structure facing cytoplasm and both ends in extracellular space. We also found that this U-shaped amino terminus probably functions as a gate to regulate the drug transport activity of human ABCC1.

Amino Acid Sequence↗

A statistical framework to discover true associations from multiprotein complex pull-down proteomics data sets.

Experimental processes to collect and process proteomics data are increasingly complex, and the computational methods to assess the quality and significance of these data remain unsophisticated. These challenges have led to many biological oversights and computational misconceptions. We developed an empirical Bayes model to analyze multiprotein complex (MPC) proteomics data derived from peptide mass spectrometry detections of purified protein complex pull-down experiments. Using our model and two yeast proteomics data sets, we estimated that there should be an average of about 20 true associations per MPC, almost 10 times as high as was previously estimated. For data sets generated to mimic a real proteome, our model achieved on average 80% sensitivity in detecting true associations, as compared with the 3% sensitivity in previous work, while maintaining a comparable false discovery rate of 0.3%. Cross-examination of our results with protein complexes confirmed by various experimental techniques demonstrates that many true associations that cannot be identified by previous approach are identified by our method.

Algorithms↗

A mixture model-based discriminate analysis for identifying ordered transcription factor binding site pairs in gene promoters directly regulated by estrogen receptor-alpha.

MOTIVATION: To detect and select patterns of transcription factor binding sites (TFBSs) which distinguish genes directly regulated by estrogen receptor-alpha (ERalpha), we developed an innovative mixture model-based discriminate analysis for identifying ordered TFBS pairs. RESULTS: Biologically, our proposed new algorithm clearly suggests that TFBSs are not randomly distributed within ERalpha target promoters (P-value < 0.001). The up-regulated targets significantly (P-value < 0.01) possess TFBS pairs, (DBP, MYC), (DBP, MYC/MAX heterodimer), (DBP, USF2) and (DBP, MYOGENIN); and down-regulated ERalpha target genes significantly (P-value < 0.01) possess TFBS pairs, such as (DBP, c-ETS1-68), (DBP, USF2) and (DBP, MYOGENIN). Statistically, our proposed mixture model-based discriminate analysis can simultaneously perform TFBS pattern recognition, TFBS pattern selection, and target class prediction; such integrative power cannot be achieved by current methods. AVAILABILITY: The software is available on request from the authors. CONTACT: lali@iupui.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Algorithms↗

A phase I, pharmacokinetic and pharmacodynamic dose escalation trial of weekly paclitaxel with interferon-alpha2b in patients with solid tumors.

PURPOSE: Paclitaxel and interferon have demonstrated anti-angiogenic activity in vitro and in vivo. The toxicity, pharmacokinetics, and pharmacodynamics of paclitaxel with interferon-alpha2b (IFN-alpha2b) were assessed in patients with solid tumors to assess the feasibility of this novel anti-angiogenic regimen. METHODS: IFN-alpha2b (1 million units) was administered twice daily by subcutaneous injection. Paclitaxel was given weekly over 1 h starting at 30 mg/m2 and increased to 50 mg/m2. Cycles were repeated every 4 weeks. RESULTS: Nineteen patients with a variety of solid tumors were enrolled. Dose-limiting toxicity in cycle 1 was observed at 50 mg/m2. Eleven patients were treated at 40 mg/m2 with no undue toxicity. Pharmacokinetic parameter comparison studies were completed in 11 patients who received days 1 and 29 paclitaxel. Mean paclitaxel clearance and area under the curve (0-infinity) were not statistically different from days 1 to 29. There was a 50% increase in the average Cmax from days 1 to 29. There was also a 73% decrease of matrix metalloproteinase-9 (MMP-9) levels in these 11 patients from days 1 to 29 (p < 0.0005). All three patients with cutaneous angiosarcomas experienced clinically meaningful remissions. In addition, minor responses were observed in one patient with heavily pretreated ovarian cancer and another with adrenocortical carcinoma. CONCLUSION: This trial details the inability to dose escalate to the maximum tolerated dose of weekly paclitaxel when combined with low-dose interferon. However, this low-dose regimen caused a significant decrease in MMP-9 and demonstrated anti-cancer activity in cutaneous angiosarcomas.

Adult↗

Endoxifen, a secondary metabolite of tamoxifen, and 4-OH-tamoxifen induce similar changes in global gene expression patterns in MCF-7 breast cancer cells.

We recently demonstrated that endoxifen (4-hydroxy-N-desmethyl-tamoxifen), a pharmacogenetically regulated metabolite of tamoxifen, is equipotent to 4-hydroxy-tamoxifen (4-OH-Tam) with respect to estrogen receptor binding and inhibition of 17beta-estradiol (E2)-induced cell proliferation. Endoxifen was also found to be more abundant in human plasma than 4-OH-Tam, and its formation has been shown to be primarily catalyzed by cytochrome P450 2D6 (CYP2D6). Here, we report studies evaluating the effects of endoxifen, 4-OH-Tam, and E2 on gene expression in MCF-7 cells using Affymetrix U133A GeneChip Arrays (Santa Clara, CA). We detected 4062 genes that were E2-regulated (1924 induced; 2138 suppressed), and the ratio of E2-induced versus E2-suppressed genes was consistent regardless of the cutoff value. In the presence of E2, 2444 and 2390 genes were affected by 4-OH-Tam and endoxifen, respectively, when no minimal -fold change cutoff was implemented. The majority of genes regulated by the tamoxifen metabolites were also E2-responsive (74.4 and 73.3%, respectively). Endoxifen and 4-OH-Tam had overlapping effects on 1365 E2-sensitive genes, whose -fold effects between these metabolites were highly correlated (R2 = 0.99). A significant correlation was also found between the -fold effects of 249 E2-insensitive genes coregulated by both metabolites (R2 = 0.99). Hierarchical clustering analysis demonstrated similar gene regulation patterns between these metabolites, which were distinct from E2 or vehicle treatment patterns. Using real time-polymerase chain reaction, we validated the gene expression patterns of five genes that were differentially regulated by endoxifen and 4-OH-Tam. We conclude that endoxifen and 4-OH-Tam have similar effects on global gene expression patterns in MCF-7 cells and that the majority of the affected genes are estrogen-regulated genes.

Breast Neoplasms↗

Prognostic DNA methylation biomarkers in ovarian cancer.

PURPOSE: Aberrant DNA methylation, now recognized as a contributing factor to neoplasia, often shows definitive gene/sequence preferences unique to specific cancer types. Correspondingly, distinct combinations of methylated loci can function as biomarkers for numerous clinical correlates of ovarian and other cancers. EXPERIMENTAL DESIGN: We used a microarray approach to identify methylated loci prognostic for reduced progression-free survival (PFS) in advanced ovarian cancer patients. Two data set classification algorithms, Significance Analysis of Microarray and Prediction Analysis of Microarray, successfully identified 220 candidate PFS-discriminatory methylated loci. Of those, 112 were found capable of predicting PFS with 95% accuracy, by Prediction Analysis of Microarray, using an independent set of 40 advanced ovarian tumors (from 20 short-PFS and 20 long-PFS patients, respectively). Additionally, we showed the use of these predictive loci using two bioinformatics machine-learning algorithms, Support Vector Machine and Multilayer Perceptron. CONCLUSION: In this report, we show that highly prognostic DNA methylation biomarkers can be successfully identified and characterized, using previously unused, rigorous classifying algorithms. Such ovarian cancer biomarkers represent a promising approach for the assessment and management of this devastating disease.

Adenocarcinoma↗

Effects of HIV protease inhibitor ritonavir on Akt-regulated cell proliferation in breast cancer.

PURPOSE: These studies were designed to determine whether ritonavir inhibits breast cancer in vitro and in vivo and, if so, how. EXPERIMENTAL DESIGN: Ritonavir effects on breast cancer cell growth were studied in the estrogen receptor (ER)-positive lines MCF7 and T47D and in the ER-negative lines MDA-MB-436 and MDA-MB-231. Effects of ritonavir on Rb-regulated and Akt-mediated cell proliferation were studied. Ritonavir was tested for inhibition of a mammary carcinoma xenograft. RESULTS: ER-positive estradiol-dependent lines (IC50, 12-24 micromol/L) and ER-negative (IC50, 45 micromol/L) lines exhibit ritonavir sensitivity. Ritonavir depletes ER-alpha levels notably in ER-positive lines. Ritonavir causes G1 arrest, depletes cyclin-dependent kinases 2, 4, and 6 and cyclin D1 but not cyclin E, and depletes phosphorylated Rb and Ser473 Akt. Ritonavir induces apoptosis independent of G1 arrest, inhibiting growth of cells that have passed the G1 checkpoint. Myristoyl-Akt, but not activated K-Ras, rescues ritonavir inhibition. Ritonavir inhibited a MDA-MB-231 xenograft and intratumoral Akt activity at a clinically attainable serum Cmax of 22 +/- 8 micromol/L. Because heat shock protein 90 (Hsp90) substrates are depleted by ritonavir, ritonavir effects on Hsp90 were tested. Ritonavir binds Hsp90 (K(D), 7.8 micromol/L) and partially inhibits its chaperone function. Ritonavir blocks association of Hsp90 with Akt and, with sustained exposure, notably depletes Hsp90. Stably expressed Hsp90alpha short hairpin RNA also depletes Hsp90, inhibiting proliferation and sensitizing breast cancer cells to low ritonavir concentrations. CONCLUSIONS: Ritonavir inhibits breast cancer growth in part by inhibiting Hsp90 substrates, including Akt. Ritonavir may be of interest for breast cancer therapeutics and its efficacy may be increased by sustained exposure or Hsp90 RNA interference.

Animals↗

Analysis of angiogenesis genes from paraffin-embedded breast tumor and lymph nodes.

BACKGROUND: Angiogenesis is important in tumor growth and metastasis. Germ-line polymorphisms critical to the angiogenesis pathway have been shown to confer prognostic information in multiple tumor types. These genes include vascular endothelial growth factor (VEGF) and endothelial nitric oxide synthase (eNOS). EXPERIMENTAL DESIGN: We extracted DNA from 53 specimens obtained from 21 patients including a primary breast tumor, and/or a histologically involved lymph node, and/or a histologically normal lymph node. We subsequently genotyped all specimens to evaluate two polymorphisms in the eNOS gene and one polymorphism in the VEGF gene. RESULTS: Chromatographs were generated in 145/159 (91%) samples. When assessing all polymorphisms by site, chromatographs were generated in 42/51 (83%) samples obtained from the primary tumor and 103/108 (95%) from lymph nodes. Chromatographs were generated in 46/53 (87%) samples from the T(-786)C polymorphism in the 5'-flanking region in the eNOS gene, 49/53 (92%) when assessing the Glu298Asp polymorphism in exon 7 in the eNOS gene and 50/53 samples (94%) for the C(936)T polymorphism in the VEGF gene. There was 100% concordance between analyses from the primary tumor, uninvolved lymph node, and involved lymph node from the same case. CONCLUSION: We successfully extracted DNA and genotyped several polymorphisms in two genes important in angiogenesis. These genotypes were determined in breast tumors, but also in involved and uninvolved lymph nodes. There was concordance between the genotypes of germline DNA obtained from uninvolved lymph nodes and those determined in tumor samples, implying that the host angiogenic genotype imprints the tumor genotype.

Breast Neoplasms↗

Quantitative effect of CYP2D6 genotype and inhibitors on tamoxifen metabolism: implication for optimization of breast cancer treatment.

BACKGROUND AND OBJECTIVES: N-Desmethyltamoxifen (NDM), a major primary metabolite of tamoxifen, is hydroxylated by cytochrome P450 (CYP) 2D6 to yield endoxifen. Because of its high antiestrogenic potency, endoxifen may play an important role in the clinical activity of tamoxifen. We conducted a prospective trial in 158 patients with breast cancer who were taking tamoxifen to further understand the effect of CYP2D6 genotype and concomitant medications on endoxifen plasma concentrations. METHODS: Medication history, genotype for 33 CYP2D6 alleles, and plasma concentrations of tamoxifen and its metabolites were determined at the fourth month of tamoxifen treatment. RESULTS: By use of a mixture model approach, endoxifen plasma concentration identified 2 phenotypic groups, whereas 4 were defined by the endoxifen/NDM plasma concentration ratio. Three distinct genotype groups were identified in the distribution of endoxifen/NDM ratio: (1) low ratios composed of patients lacking any functional allele (mean, 0.04 +/- 0.02); (2) intermediate ratios represented by patients with 1 active allele (mean, 0.08 +/- 0.04); and (3) high ratios composed of patients with 2 or more functional alleles (mean, 0.15 +/- 0.09). Endoxifen/NDM plasma ratios were significantly different between these groups (P < .001). The mean endoxifen plasma concentration was significantly lower in CYP2D6 extensive metabolizers who were taking potent CYP2D6 inhibitors than in those who were not taking CYP2D6 inhibitors (23.5 +/- 9.5 nmol/L versus 84.1 +/- 39.4 nmol/L, P < .001). CONCLUSION: CYP2D6 genotype and concomitant potent CYP2D6 inhibitors are highly associated with endoxifen plasma concentration and may have an impact on the response to tamoxifen therapy. These iterative approaches may be valuable in the study of other complex genotype-phenotype relationships.

Adult↗

Differential gene induction by type I and type II interferons and their combination.

Type I and type II interferons (IFNs) bind to different cell surface receptors but activate overlapping signal transduction pathways. We examined the effects of a type I IFN (IFN-alphacon1) and a type II IFN (IFN-gamma1b) on gene expression in A549 cells and demonstrate that there is a common set of genes modulated by both IFNs as well as a set of genes specifically regulated by each, reflecting the activation of different signaling pathways. In particular, IFN-gamma induced many more genes of the signaling pathways, apoptosis, and cytokine interactions than did IFN-alpha. Even with genes induced by both IFNs there were distinctive quantitative differences in expression. IFN-gamma1b plays a major role in the induction and regulation of the complement pathway. Previous work has shown a synergistic antiviral and antiproliferative effect of type I and type II IFNs in cell culture and in the treatment of tumors in mice. We demonstrate that a majority of genes showed an additive effect of IFN-alphacon1 and IFN-gamma1b, but a subset of genes is synergistically induced; these include ISG20, MX2, OAS2, and other genes known to be involved in the antiviral response, TRAIL (TNFSF10) and caspases involved in apoptosis and chemokine genes RANTES, CXCL10, and CXCL11. Greater than additive transcription of some of these genes in the presence of both IFNs was confirmed by real-time kinetic RT-PCR. Elevated induction of many of these genes may be sufficient to explain the synergistic antiviral and antitumor effects of this combination of IFNs in vivo.

Antineoplastic Agents↗

Dextromethorphan to dextrorphan urinary metabolic ratio does not reflect dextromethorphan oral clearance.

Dextromethorphan urinary metabolic ratio is widely used to determine the CYP2D6 phenotype, but its utility to reflect subtle differences in catalytic activity is unclear. We evaluated the capability of dextromethorphan urinary metabolic ratio to predict dextromethorphan oral clearance as a measure of CYP2D6 activity. Data from 10 healthy extensive metabolizers of CYP2D6 were given 30 mg of dextromethorphan hydrobromide orally on two occasions. Blood and urine samples were collected for 72 h. Dextromethorphan and dextrorphan were determined in urine by high-performance liquid chromatography with fluorescence detection and in serum by liquid chromatography-mass spectrometry. The urinary metabolic ratio was very weakly correlated with dextromethorphan oral clearance (r = 0.24; p = 0.04). In contrast, the dextromethorphan oral clearance was highly correlated with the dextromethorphan to dextrorphan area under the concentration-time curve ratio (r = 0.84; p = 0.005) and the 3-h (r = 0.60; p = 0.003), 4-h (r = 0.72, p < 0.001), 6-h (r = 0.67; p < 0.001), and 8-h (r = 0.74; p < 0.001) dextromethorphan to dextrorphan serum ratios. Assuming an effect size of 30%, the number of volunteers required for crossover and cross-sectional studies using the urinary metabolic ratio as the CYP2D6 index was calculated to be 56 and 524, respectively, whereas 14 and 60 subjects are needed if oral clearance is used. Considering the required sample size and the low correlation with oral clearance, urinary metabolic ratio is not recommended as the primary outcome variable in studies requiring the detection of modest changes in CYP2D6 activity.

Administration, Oral↗

Effect of cirrhosis and liver transplantation on the gender difference in QT interval.

The effect of gender on the QT interval in patients with cirrhosis before and after liver transplantation and the relation between the QT interval and changes in gender hormones after liver transplantation were assessed. The study showed that (1) physiologic gender difference in the QTc interval is abolished in cirrhosis, and it is not restored after liver transplantation, and (2) gender hormone concentrations had no effect on the QTc interval.

Adolescent↗

CYP2D6 genotype, antidepressant use, and tamoxifen metabolism during adjuvant breast cancer treatment.

BACKGROUND: The efficacy of tamoxifen therapy for the treatment of breast cancer varies widely among individuals. Plasma concentrations of the active tamoxifen metabolite endoxifen are associated with the cytochrome P450 (CYP) 2D6 genotype. We examined the effects of concomitant use of selective serotonin reuptake inhibitor antidepressants, which are CYP2D6 enzyme inhibitors commonly prescribed to treat hot flashes in women who take tamoxifen, and genotypes for genes that encode tamoxifen-metabolizing enzymes on plasma concentrations of tamoxifen and its metabolites. METHODS: Eighty patients with newly diagnosed with breast cancer who were beginning tamoxifen therapy (20 mg/day orally), 24 of whom were taking CYP2D6 inhibitors, were genotyped for common alleles of the CYP2D6, CYP2C9, CYP3A5, and sulfotransferase (SULT) 1A1 genes. Plasma concentrations of tamoxifen and its metabolites were measured after 1 and 4 months of tamoxifen therapy. Differences in plasma concentrations of tamoxifen and its metabolites between genotype groups were analyzed by the Wilcoxon rank sum test. All statistical tests were two-sided. RESULTS: Among all women, plasma endoxifen concentrations after 4 months of tamoxifen therapy were statistically significantly lower in subjects with a CYP2D6 homozygous variant genotype (20.0 nM, 95% confidence interval [CI] = 11.1 to 28.9 nM) or a heterozygous genotype (43.1 nM, 95% CI = 33.3 to 52.9 nM) than in those with a homozygous wild-type genotype (78.0 nM, 95%CI = 65.9 to 90.1 nM) (both P = .003). Among subjects who carried a homozygous wild-type genotype, the mean plasma endoxifen concentration for those who were using CYP2D6 inhibitors was 58% lower than that for those who were not (38.6 nM versus 91.4 nM, difference = -52.8 nM, 95% CI = -86.1 to -19.5 nM, P = .0025). The plasma endoxifen concentration was slightly reduced in women taking venlafaxine, a weak inhibitor of CYP2D6, whereas the plasma endoxifen concentration was reduced substantially in subjects who took paroxetine (a potent inhibitor of CYP2D6). Genetic variations of CYP2C9, CYP3A5, or SULT1A1 had no statistically significant associations with plasma concentrations of tamoxifen or its metabolites. CONCLUSION: Interactions between CYP2D6 polymorphisms and coadministered antidepressants and other drugs that are CYP2D6 inhibitors may be associated with altered tamoxifen activity.

Adult↗

Epigenetic hypothesis tests for methylation and acetylation in a triple microarray system.

To fully elucidate the functional relationship between DNA methylation and histone hypoacetylation in gene silencing, we have developed an integrated "triple" microarray system that allows us to begin to decipher the influence of epigenetic hierarchies on the regulation of gene expression in cancer cells. Our hypothesis is that in the promoter region of a silenced gene, reversal of two epigenetic factors (i.e., DNA demethylation and/or histone hyperacetylation) is highly correlated with gene reexpression after treatment of the human epithelial ovarian cancer cell line CP70 with the drug combination 5-aza-2'-deoxycytidine (DAC), a demethylating agent, and trichostatin A (TSA), an inhibitor of histone deacetylases. To estimate the posterior probabilities for genes with altered expression, DNA methylation and histone acetylation status measured with a triple-microarray system, we have employed an established empirical Bayes model. Two methods have been proposed to test our hypothesis that DNA demethylation and histone hyperacetylation are highly correlated among those up-regulated genes. One method follows a weighted least squares regression, while the other is derived from a chi-square statistic. The data derived by these approaches, which have been further verified through bootstrap analyses, support the proposed epigenetic correlation (p-values are less than 0.001). Further simulations suggest that even if the constant variance and normality assumptions do not hold, the power of those two tests is robust.

Acetylation↗