Expression a La Bimode

The recently published article by Wang, and commented on by Ertel in this journal, describes development of a method, named the Bimodality Index, to objectively identify and rank meaningful and reliable bimodal patterns from large-scale gene expression datasets. This is an important step because genes with a bimodal distribution may be prime candidates to develop as molecular diagnostic tests for distinguishing clinically important groups. In their introduction, they make the point that a bimodal expression pattern may be observed when two distinct subgroups of samples are measured as one group, with each mode representing the mean expression of a gene in one of the sub-groups. Recently, it was reported that bimodal gene expression may represent another important phenomenon. Specifically, two clinically important groups may be distinguished on the basis of a gene being unimodal in one group and bimodal in the other group. For example, certain key antioxidant, DNA repair, and transcription factor genes each display a unimodal pattern in non-lung cancer subjects but a bimodal pattern in lung cancer subjects. Importantly, when two groups are distinguished on the basis of unimodal vs bimodal distribution of an analyte (in this case gene expression value), biomarker development requires use of two cut-points, one on either side of the unimodal distribution. It is timely then that an improved method for separating groups on the basis of two cut-points was described recently. Thus, the analytical method to discover genes with bimodal expression distribution described by Wang and recent statistical methods that enable separation of groups on the basis of two-cut points are likely to accelerate both diagnostic test discovery, as well as mechanistic understanding of disease and disease risk. It is important to recognize that transcript abundance methods vary in their linear dynamic range, and that the ability to identify bimodal distribution of genes will depend to a significant degree on the linear dynamic range of the method used

Quality Control Methods for Optimal BCR-ABL1 Clinical Testing in Human Whole Blood Samples

Reliable breakpoint cluster region (BCR)–Abelson (ABL) 1 measurement is essential for optimal management of chronic myelogenous leukemia. There is a need to optimize quality control, sensitivity, and reliability of methods used to measure a major molecular response and/or treatment failure. The effects of room temperature storage time, different primers, and RNA input in the reverse transcription (RT) reaction on BCR-ABL1 and β-glucuronidase (GUSB) cDNA yield were assessed in whole blood samples mixed with K562 cells. BCR-ABL1 was measured relative to GUSB to control for sample loading, and each gene was measured relative to known numbers of respective internal standard molecules to control for variation in quality and quantity of reagents, thermal cycler conditions, and presence of PCR inhibitors. Clinical sample and reference material measurements with this test were concordant with results reported by other laboratories. BCR-ABL1 per 103 GUSB values were significantly reduced (P = 0.004) after 48-hour storage. Gene-specific primers yielded more BCR-ABL1 cDNA than random hexamers at each RNA input. In addition, increasing RNA inhibited the RT reaction with random hexamers but not with gene-specific primers. Consequently, the yield of BCR-ABL1 was higher with gene-specific RT primers at all RNA inputs tested, increasing to as much as 158-fold. We conclude that optimal measurement of BCR-ABL1 per 103 GUSB in whole blood is obtained when gene-specific primers are used in RT and samples are analyzed within 24 hours after blood collection

Stable low-level expression of p21(WAF1/CIP1 )in A549 human bronchogenic carcinoma cell line-derived clones down-regulates E2F1 mRNA and restores cell proliferation control

BACKGROUND: Deregulated cell cycle progression and loss of proliferation control are key properties of malignant cells. In previous studies, an interactive transcript abundance index (ITAI) comprising three cell cycle control genes, [MYC × E2F1]/p21 accurately distinguished normal from malignant bronchial epithelial cells (BEC), using a cut-off threshold of 7,000. This cut-off is represented by a line with a slope of 7,000 on a bivariate plot of p21 versus [MYC × E2F1], with malignant BEC above the line and normal BEC below the line. This study was an effort to better quantify, at the transcript abundance level, the difference between normal and malignant BEC. The hypothesis was tested that experimental elevation of p21 in a malignant BEC line would decrease the value of the [MYC × E2F1]/p21 ITAI to a level below this line, resulting in loss of immortality and limited cell population doubling capacity. In order to test the hypothesis, a p21 expression vector was transfected into the A549 human bronchogenic carcinoma cell line, which has low constitutive p21 TA expression relative to normal BEC. RESULTS: Following transfection of p21, four A549/p21 clones with stable two-fold up-regulated p21 expression were isolated and expanded. For each clone, the increase in p21 transcript abundance (TA) was associated with increased total p21 protein level, more than 5-fold reduction in E2F1 TA, and 10-fold reduction in the [MYC × E2F1]/p21 ITAI to a value below the cut-off threshold. These changes in regulation of cell cycle control genes were associated with restoration of cell proliferation control. Specifically, each transfectant was capable of only 15 population doublings compared with unlimited population doublings for parental A549. This change was associated with an approximate 2-fold increase in population doubling time to 38.4 hours (from 22.3 hrs), resumption of contact-inhibition, and reduced dividing cell fraction as measured by flow cytometric DNA analysis. CONCLUSION: These results, likely due to increased p21-mediated down-regulation of E2F1 TA at the G1/S phase transition, are consistent with our hypothesis. Specifically, they provide experimental confirmation that a line with slope of 7,000 on the p21 versus [MYC × E2F1] bivariate plot quantifies the difference between normal and malignant BEC at the level of transcript abundance

ABCC5, ERCC2, XPA and XRCC1 transcript abundance levels correlate with cisplatin chemoresistance in non-small cell lung cancer cell lines

BACKGROUND: Although 40–50% of non-small cell lung cancer (NSCLC) tumors respond to cisplatin chemotherapy, there currently is no way to prospectively identify potential responders. The purpose of this study was to determine whether transcript abundance (TA) levels of twelve selected DNA repair or multi-drug resistance genes (LIG1, ERCC2, ERCC3, DDIT3, ABCC1, ABCC4, ABCC5, ABCC10, GTF2H2, XPA, XPC and XRCC1) were associated with cisplatin chemoresistance and could therefore contribute to the development of a predictive marker. Standardized RT (StaRT)-PCR, was employed to assess these genes in a set of NSCLC cell lines with a previously published range of sensitivity to cisplatin. Data were obtained in the form of target gene molecules relative to 10(6 )β-actin (ACTB) molecules. To cancel the effect of ACTB variation among the different cell lines individual gene expression values were incorporated into ratios of one gene to another. Each two-gene ratio was compared as a single variable to chemoresistance for each of eight NSCLC cell lines using multiple regression. In an effort to validate these results, six additional lines then were evaluated. RESULTS: Following validation, single variable models best correlated with chemoresistance (p < 0.001), were ERCC2/XPC, ABCC5/GTF2H2, ERCC2/GTF2H2, XPA/XPC and XRCC1/XPC. All single variable models were examined hierarchically to achieve two variable models. The two variable model with the highest correlation was (ABCC5/GTF2H2, ERCC2/GTF2H2) with an R(2 )value of 0.96 (p < 0.001). CONCLUSION: These results provide markers suitable for assessment of small fine needle aspirate biopsies in an effort to prospectively identify cisplatin resistant tumors

Variation in transcriptional regulation of cyclin dependent kinase inhibitor p21(waf1/cip1 )among human bronchogenic carcinomas

BACKGROUND: Cell proliferation control depends in part on the carefully ordered regulation of transcription factors. The p53 homolog p73, contributes to this control by directly upregulating the cyclin dependent kinase inhibitor, p21(waf1/cip1). E2F1, an inducer of cell proliferation, directly upregulates p73 and in some systems upregulates p21 directly. Because of its central role in controlling cell proliferation, upregulation of p21 has been explored as a modality for treating bronchogenic carcinoma (BC). Improved understanding of p21 transcriptional regulation will facilitate identification of BC tissues that are responsive to p21-directed therapies. Toward this goal, we investigated the role that E2F1 and p73 each play in the transcriptional regulation of p21. RESULTS: Among BC samples (N = 21) p21 transcript abundance (TA) levels varied over two orders of magnitude with values ranging from 400 to 120,000 (in units of molecules/10(6 )molecules β-actin). The p21 values in many BC were high compared to those observed in normal bronchial epithelial cells (BEC) (N = 18). Among all BC samples, there was no correlation between E2F1 and p21 TA but there was positive correlation between E2F1 and p73α (p < 0.001) TA. Among BC cell lines with inactivated p53 and wild type p73 (N = 7) there was positive correlation between p73α and p21 TA (p < 0.05). Additionally, in a BC cell line in which both p53 and p73 were inactivated (H1155), E2F1 TA level was high (50,000), but p21 TA level was low (470). Transiently expressed exogenous p73α in the BC cell line Calu-1, was associated with a significant (p < 0.05) 90% increase in p21 TA and a 20% reduction in E2F1 TA. siRNA mediated reduction of p73 TA in the N417 BC cell line was associated with a significant reduction in p21 TA level (p < 0.01). CONCLUSION: p21 TA levels vary considerably among BC patients which may be attributable to 1) genetic alterations in Rb and p53 and 2) variation in TA levels of upstream transcription factors E2F1 and p73. Here we provide evidence that p73 upregulates p21 TA in BC tissues and upregulated p21 TA may result from E2F1 upregulation of p73 but not from E2F1 directly

MARCKS regulates growth and radiation sensitivity and is a novel prognostic factor for glioma

Purpose This study assessed whether Myristoylated Alanine Rich C-Kinase Substrate (MARCKS) can regulate glioblastoma (GBM) growth, radiation sensitivity and clinical outcome. Experimental Design MARCKS protein levels were analyzed in five GBM explant cell lines and eight patient-derived xenograft tumors by immunoblot, and these levels were correlated to proliferation rates and intracranial growth rates, respectively. Manipulation of MARCKS protein levels was assessed by lentiviral-mediated shRNA knockdown in the U251 cell line and MARCKS over-expression in the U87 cell line. The effect of manipulation of MARCKS on proliferation, radiation sensitivity and senescence was assessed. MARCKS gene expression was correlated with survival outcomes in the Repository of Molecular Brain Neoplasia Data (REMBRANDT) Database and The Cancer Genome Atlas (TCGA). Results MARCKS protein expression was inversely correlated with GBM proliferation and intracranial xenograft growth rates. Genetic silencing of MARCKS promoted GBM proliferation and radiation resistance, while MARCKS overexpression greatly reduced GBM growth potential and induced senescence. We found MARCKS gene expression to be directly correlated with survival in both the REMBRANDT and TCGA databases. Specifically, patients with high MARCKS expressing tumors of the Proneural molecular subtype had significantly increased survival rates. This effect was most pronounced in tumors with unmethylated O6-methylguanine DNA methyltransferase (MGMT) promoters, a traditionally poor prognostic factor. Conclusions MARCKS levels impact GBM growth and radiation sensitivity. High MARCKS expressing GBM tumors are associated with improved survival, particularly with unmethylated MGMT promoters. These findings suggest the use of MARCKS as a novel target and biomarker for prognosis in the Proneural subtype of GBM

Kinomic exploration of temozolomide and radiation resistance in Glioblastoma multiforme xenolines

Glioblastoma multiforme (GBM) represents the most common and deadly primary brain malignancy, particularly due to temozolomide (TMZ) and radiation (RT) resistance. To better understand resistance mechanisms, we examined global kinase activity (kinomic profiling) in both treatment sensitive and resistant human GBM patient-derived xenografts (PDX or “xenolines”)

Kinomic Profiling of Electromagnetic Navigational Bronchoscopy Specimens: A New Approach for Personalized Medicine

Purpose Researchers are currently seeking relevant lung cancer biomarkers in order to make informed decisions regarding therapeutic selection for patients in so-called “precision medicine.” However, there are challenges to obtaining adequate lung cancer tissue for molecular analyses. Furthermore, current molecular testing of tumors at the genomic or transcriptomic level are very indirect measures of biological response to a drug, particularly for small molecule inhibitors that target kinases. Kinase activity profiling is therefore theorized to be more reflective of in vivo biology than many current molecular analysis techniques. As a result, this study seeks to prove the feasibility of combining a novel minimally invasive biopsy technique that expands the number of lesions amenable for biopsy with subsequent ex vivo kinase activity analysis. Methods Eight patients with lung lesions of varying location and size were biopsied using the novel electromagnetic navigational bronchoscopy (ENB) technique. Basal kinase activity (kinomic) profiles and ex vivo interrogation of samples in combination with tyrosine kinase inhibitors erlotinib, crizotinib, and lapatinib were performed by PamStation 12 microarray analysis. Results Kinomic profiling qualitatively identified patient specific kinase activity profiles as well as patient and drug specific changes in kinase activity profiles following exposure to inhibitor. Thus, the study has verified the feasibility of ENB as a method for obtaining tissue in adequate quantities for kinomic analysis and has demonstrated the possible use of this tissue acquisition and analysis technique as a method for future study of lung cancer biomarkers. Conclusions We demonstrate the feasibility of using ENB-derived biopsies to perform kinase activity assessment in lung cancer patients

The Reproducibility of Lists of Differentially Expressed Genes in Microarray Studies

Reproducibility is a fundamental requirement in scientific experiments and clinical contexts. Recent publications raise concerns about the reliability of microarray technology because of the apparent lack of agreement between lists of differentially expressed genes (DEGs). In this study we demonstrate that (1) such discordance may stem from ranking and selecting DEGs solely by statistical significance (P) derived from widely used simple t-tests; (2) when fold change (FC) is used as the ranking criterion, the lists become much more reproducible, especially when fewer genes are selected; and (3) the instability of short DEG lists based on P cutoffs is an expected mathematical consequence of the high variability of the t-values. We recommend the use of FC ranking plus a non-stringent P cutoff as a baseline practice in order to generate more reproducible DEG lists. The FC criterion enhances reproducibility while the P criterion balances sensitivity and specificity

Targeting mitochondrial responses to intra-articular fracture to prevent posttraumatic osteoarthritis

We tested whether inhibiting mechanically responsive articular chondrocyte mitochondria after severe traumatic injury and preventing oxidative damage represent a viable paradigm for posttraumatic osteoarthritis (PTOA) prevention. We used a porcine hock intra-articular fracture (IAF) model well suited to human-like surgical techniques and with excellent anatomic similarities to human ankles. After IAF, amobarbital or N-acetylcysteine (NAC) was injected to inhibit chondrocyte electron transport or downstream oxidative stress, respectively. Effects were confirmed via spectrophotometric enzyme assays or glutathione/glutathione disulfide assays and immunohistochemical measures of oxidative stress. Amobarbital or NAC delivered after IAF provided substantial protection against PTOA at 6 months, including maintenance of proteoglycan content, decreased histological disease scores, and normalized chondrocyte metabolic function. These data support the therapeutic potential of targeting chondrocyte metabolism after injury and suggest a strong role for mitochondria in mediating PTOA