Selection bias in gene extraction on the basis of microarray gene-expression data

In the context of cancer diagnosis and treatment, we consider the problem of constructing an accurate prediction rule on the basis of a relatively small number of tumor tissue samples of known type containing the expression data on very many (possibly thousands) genes. Recently, results have been presented in the literature suggesting that it is possible to construct a prediction rule from only a few genes such that it has a negligible prediction error rate. However, in these results the test error or the leave-one-out cross-validated error is calculated without allowance for the selection bias. There is no allowance because the rule is either tested on tissue samples that were used in the first instance to select the genes being used in the rule or because the cross-validation of the rule is not external to the selection process; that is, gene selection is not performed in training the rule at each stage of the cross-validation process. We describe how in practice the selection bias can be assessed and corrected for by either performing a cross-validation or applying the bootstrap external to the selection process. We recommend using 10-fold rather than leave-one-out cross-validation, and concerning the bootstrap, we suggest using the so-called. 632+ bootstrap error estimate designed to handle overfitted prediction rules. Using two published data sets, we demonstrate that when correction is made for the selection bias, the cross-validated error is no longer zero for a subset of only a few genes.

Assessment of transient gene expression in plant tissues using the green fluorescent protein as a reference

Four different promoters (35S and enhanced 35S of the cauliflower mosaic virus, polyubiquitin of maize and actin1 of rice) were compared in a transient assay using maize leaves and particle bombardment. A gene encoding the jellyfish green fluorescent protein (GFP) driven by the 358 promoter was used as an internal standard to monitor the effectiveness of each bombardment. Normalisation of the transient expression assay using the GFP reference significantly reduced the variability between separate bombardments and allowed for a rapid and accurate evaluation of different promoters in microprojectile-bombarded leaves.

Sustained gene expression in transplanted skin fibroblasts in rats

Retrovirus-mediated gene transfer into adult skin fibroblasts has provided measurable amounts of therapeutic proteins in animal models. However, the major problem emerging from these experiments was a limited time of vector encoded gene expression once transduced cells were engrafted We hypothesized that sustained transduced gene expression in quiescent fibroblasts in vivo might be obtained by using a fibronectin (Fn) promoter. Fibronectin plays a key role in cell adhesion, migration and wound healing and is up-regulated in quiescent fibroblasts. Retroviral vectors containing human adenosine deaminase (ADA) cDNA linked to rat fibronectin promoter (LNFnA) or viral LTR promoter (LASN) were compared for their ability to express ADA from transduced primary rat skin fibroblasts in vivo. Skin grafts formed from fibroblasts transduced with LNFnA showed strong human ADA enzyme activity from 1 week to 3 months. In contrast, skin grafts containing LASN-transduced fibroblasts tested positive for human ADA for weeks 1 and 2, were faintly positive at week 3 and showed no human ADA expression at 1, 2 and 3 months. Thus, a fibronectin promoter provided sustained transduced gene expression at high levels for at least 3 months in transplanted rat skin fibroblasts, perhaps permitting the targeting of this tissue for human gene therapy.

Regulation of the Hsp90-binding immunophilin, cyclophilin 40, is mediated by multiple sites for CA-binding protein (GABP)

Within steroid receptor heterocomplexes the large tetraticopeptide repeat-containing immunophilins, cyclophilin 40 (CyP40), FKBP51, and FKBP52, target a common interaction site in heat shock protein 90 (HspSO) and act coordinately with HspSO to modulate receptor activity. The reversible nature of the interaction between the immunophilins and HspSO suggests that relative cellular abundance might be a key determinant of the immunophilin component within steroid receptor complexes. To investigate CyP40 gene regulation, we have isolated a fi-kilobase (kb) 5 ' -flanking region of the human gene and demonstrated that a similar to 50 base pair (bp) sequence adjacent to the transcription start site is essential for CyP40 basal expression. Three tandemly arranged Ets sites within this critical region were identified as binding elements for the multimeric Ets-related transcription factor, GA binding protein (GABP). Functional studies of this proximal promoter sequence, in combination with mutational analysis, confirmed these sites to be crucial for basal promoter function. Furthermore, overexpression of both GABP alpha and GABP beta subunits in Cos1 cells resulted in increased endogenous CyP40 mRNA levels. Significantly, a parallel increase in FKBP52 mRNA expression was not observed, highlighting an important difference in the mode of regulation of the CyP40 and FKBP52 genes. Our results identify GABP as a key regulator of CyP40 expression. GAFF is a common target of mitogen and stress-activated pathways and may integrate these diverse extracellular signals to regulate CyP40 gene expression.

Using biplots to interpret gene expression patterns in plants

Sum: Plant biologists in fields of ecology, evolution, genetics and breeding frequently use multivariate methods. This paper illustrates Principal Component Analysis (PCA) and Gabriel's biplot as applied to microarray expression data from plant pathology experiments. Availability: An example program in the publicly distributed statistical language R is available from the web site (www.tpp.uq.edu.au) and by e-mail from the contact. Contact: scott.chapman@csiro.au.

Myb-binding protein 1a augments AhR-dependent gene expression

We have studied the mechanism by which an acidic domain (amino acids 515-583) of the aromatic hydrocarbon receptor (AhR) transactivates a target gene. Studies with glutathione S-transferase fusion proteins demonstrate that the wild-type acidic domain associates in vitro with Myb-binding protein la, whereas a mutant domain (F542A, 1569A) does not. AhR-defective cells reconstituted with an AhR containing the wild-type acidic domain exhibit normal AhR function; however, cells reconstituted with an AhR containing the mutant acidic domain do not function normally. Transient transfection of Myb-binding protein la into mouse hepatoma cells is associated with augmentation of AhR-dependent gene expression. Such augmentation does not occur when Myb-binding protein la is transfected into AhR-defective cells that have been reconstituted with an AhR that lacks the acidic domain. We infer that 1) Myb-binding protein la associates with AhR, thereby enhancing transactivation, and 2) the presence of AhR's acidic domain is both necessary and sufficient for Myb-binding protein la to increase AhR-dependent gene expression.

Gene expression profiling of papillary thyroid carcinoma identifies transcripts correlated with BRAF mutational status and lymph node metastasis

Purpose: To identify papillary thyroid carcinoma (PTC)-associated transcripts, we compared the gene expression profiles of three Serial Analysis of Gene Expression libraries generated from thyroid tumors and a normal thyroid tissue. Experimental Design: Selected transcripts were validated in a panel of 57 thyroid tumors using quantitative PCR (qPCR). An independent set of 71 paraffin-embedded sections was used for validation using immunohistochemical analysis. To determine if PTC-associated gene expression could predict lymph node involvement, a separate cohort of 130 primary PTC (54 metastatic and 76 nonmetastatic) was investigated. The BRAF(V600E) mutational status was compared with qPCR data to identify genes that might be regulated by abnormal BRAF/MEK/extracellular signal-regulated kinase signaling. Results: We identified and validated new PTC-associated transcripts. Three genes (CST6, CXCL14, and DHRS3) are strongly associated with PTC. Immunohistochemical analysis of CXCL14 confirmed the qPCR data and showed protein expression in PTC epithelial cells. We also observed that CST6, CXCL14, DHRS3, and SPP1 were associated with PTC lymph node metastasis, with CST6, CXCL14, and SPP1 being positively correlated with metastasis and DHRS3 being negatively correlated. Finally, we found a strong correlation between CST6 and CXCL14 expression and BRAF(V600E) mutational status, suggesting that these genes may be induced subsequently to BRAF activation and therefore may be downstream in the BRAF/MEK/extracellular signal-regulated kinase signaling pathway. Conclusion: CST6, CXCL14, DHRS3, and SPP1 may play a role in PTC pathogenesis and progression and are possible molecular targets for FTC therapy.

PVALB, a New Hurthle Adenoma Diagnostic Marker Identified through Gene Expression

Context: A better means to accurately identify malignant thyroid nodules and to distinguish them from benign tumors is needed. We previously identified markers for detecting thyroid malignancy, with sensitivity estimated at or close to 100%. One lingering problem with these markers was that false positives occurred with Hurthle cell adenomas (HCA) which lowered test specificity. Methods: To locate accurate diagnostic markers, we profiled in depth the transcripts of a HCA and a Hurthle cell carcinoma (HCC). From 1146 differentially expressed genes, 18 transcripts specifically expressed in HCA were tested by quantitative PCR in a wide range of thyroid tumors (n = 76). Sensibility and specificity were calculated using receiver operating characteristic (ROC). Selected markers were further validated in an independent set of thyroid tumors (n = 82) by immunohistochemistry. To define the panel that would yield best diagnostic accuracy, these markers were tested in combination with our previous identified markers. Results: Seventeen of the 18 genes showed statistical significance based on a mean relative level of expression (P < 0.05). KLK1 (sensitivity = 0.97) and PVALB (sensitivity = 0.94) were the best candidate markers. The combination of PVALB and C1orf24 increased specificity to > 97% and maintained sensitivity for detection of carcinoma. Conclusion: We identified tumor markers that can be used in combination for a more accurate preoperative diagnosis of thyroid nodules and for postoperative diagnosis of thyroid carcinoma in tumor sections. This improved test would help physicians rapidly focus treatment on true malignancies and avoid unnecessary treatment of benign tumors, simultaneously improving medical care and reducing costs. (J Clin Endocrinol Metab 96: E151-E160, 2011)