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.

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.

Insulin-like growth factor 2 cDNA cloning and ontogeny of gene expression in the liver of the marsupial brushtail possum (Trichosurus vulpecula)

The cDNA sequence for insulin-like growth factor 2 (IGF-2) was determined from the liver of the marsupial brushtail possum (Trichosurus vulpecula) using reverse transcription followed by polymerase chain reaction (RT-PCR) with gene-specific primers. The 359 bp of possum sequence encompassed the mature peptide, 27 bp of the signal peptide, and 125 bp of the E-peptide. Alignment of the deduced amino acid sequence with those from other species indicated that the mature peptide was 71 amino acids in length, 4 amino acids longer than most other mammals. At both the nucleotide and amino acid levels there was a high degree of sequence identity with IGF-2 from other mammalian and nonmammalian species. Amino acid identity ranged from 94.4% with a variant form of human IGF-2 to 80.3% with zebrafinch IGF-2. Northern analysis revealed that radiolabeled possum IGF-2, cDNA hybridized to multiple transcripts in the liver of both adult possums and 150-day-old pouch young and that the overall level of expression was greater in pouch young. Semiquantitative RT-PCR with total RNA from liver samples of pouch young aged 12 to 150 days postpartum and adults confirmed that IGF-2 gene expression was two to three times more abundant in pouch young than in adults but there was no significant change in the level of expression during pouch life. Unlike other mammalian species, in which there is a decline in levels of liver IGF-2 gene expression around the time of birth, levels in the marsupial brushtail possum remain elevated for at least 150 days after birth. This suggests that the decline in liver IGF-2 expression in marsupials and eutherians occurs at a similar stage of development and may reflect a role for this growth factor during the postnatal growth and development of the marsupial, (C) 2001 Academic Press.

Alterations in gene expression and activity during squamous cell carcinoma development

This study focuses on characterizing the genetic and biological alterations associated with squamous cell carcinoma development. Normal human epidermal keratinocytes (HEKs), cells isolated from a preneoplastic lesion (IEC-1), and two neoplastic cell lines, SCC-25 and COLD-16, were grown as raft cultures, and their gene expression profiles were screened using cDNA arrays. Our data indicated that the expression levels of at least 37 genes were significantly (P less than or equal to 0.05; 1.9% of genes screened) altered in neoplastic cells compared with normal cells. Of these genes, 10 genes were up-regulated and 27 genes were down-regulated in the neoplastic cells. In addition, 51% of the genes altered in the neoplastic cells were already altered in the preneoplastic IEC-1 cells. Immunohistochemical staining of patient tumors was used to verify the cDNA array analysis. Our analysis indicated that alterations in genes associated with extracellular matrix production and apoptosis are disrupted in preneoplastic cells, whereas later stages of neoplasia are associated with alterations in gene expression for genes involved in DNA repair or epidermal growth factor (EGF) receptor/mitogen-activated protein kinase kinase (MAPKK)/MAPK/activator protein-1 (AP-1) signaling. Subsequent functional analysis of the alterations in expression of the EGF receptor/MAPKK/MAPK/AP-1 genes suggested they did not contribute to the neoplastic phenotype.

Patterns of gene expression are altered in the frontal and motor cortices of human alcoholics

Alcoholism is a major health problem in Western countries, yet relatively little is known about the mechanisms by which chronic alcohol abuse causes the pathologic changes associated with the disease. It is likely that chronic alcoholism affects a number of signaling cascades and transcription factors, which in turn result in distinct gene expression patterns. These patterns are difficult to detect by traditional experiments measuring a few mRNAs at a time, but are well suited to microarray analyses. We used cDNA microarrays to analyze expression of approximately 10 000 genes in the frontal and motor cortices of three groups of chronic alcoholic and matched control cases. A functional hierarchy was devised for classification of brain genes and the resulting groups were compared based on differential expression. Comparison of gene expression patterns in these brain regions revealed a selective reprogramming of gene expression in distinct functional groups. The most pronounced differences were found in myelin-related genes and genes involved in protein trafficking. Significant changes in the expression of known alcohol-responsive genes, and genes involved in calcium, cAMP, and thyroid signaling pathways were also identified. These results suggest that multiple pathways may be important for neuropathology and altered neuronal function observed in alcoholism.

Improved detection of lacZ reporter gene expression in transgenic epithelia by immunofluorescence microscopy

The bacterial lacZ gene is commonly used as a reporter for the in vivo analysis of gene regulation in transgenic mice. However, several laboratories have reported poor detection of beta-galactosidase (the lacZ gene product) using histochemical techniques, particularly in skin. Here we report the difficulties we encountered in assessing lacZ expression in transgenic keratinocytes using classic X-gal histochemical protocols in tissues shown to express the transgene by mRNA in situ hybridization. We found that lacZ reporter gene expression could be reliably detected in frozen tissue sections by immunofluorescence analysis using a beta-galactosidase-specific antibody. Moreover, we were able to localize both transgene and endogenous gene products simultaneously using double-label immunofluorescence. Our results suggest that antibody detection of beta-galactosidase should be used to verify other assays of lacZ expression, particularly where low expression levels are suspected or patchy expression is observed.

Expression analysis of delta-catenin and prostate-specific membrane antigen: Their potential as diagnostic markers for prostate cancer

The current approach to prostate cancer diagnosis has major limitations including the inability of prostate-specific antigen (PSA) assays to accurately differentiate between prostate cancer and benign prostate hyperplasia (BPH) and the imprecision of transrectal ultrasound (TRUS) biopsy sampling. We have employed cDNA microarray screening to compare gene expression patterns in BPH and tumour samples to identify expression markers that may be useful in discriminating between these conditions. Screening of 3 individual cDNA arrays identified 8 genes with expression 3-fold greater in 6 tumour tissues than in 1 nontumour sample and I BPH sample. Real-time PCR was used to confirm the overexpression of these 8 genes and 12 genes selected from the literature against a panel of 17 tumours and I 1 BPH samples. Two genes, delta-catenin (delta-catenin; CTNND2) and prostate-specific membrane antigen (PSMA; FOLH1), were significantly overexpressed in prostate cancer compared to BPH. Prostate epithelial cells stained positively for S-catenin and PSMA in our prostate cancer tissues, whereas the majority of our BPH tissues were negative for both markers. Thus we have identified delta-catenin (not previously associated with prostatic adenocarcinoma) and confirmed the potential of PSMA as potential candidates for the diagnosis and management of prostate cancer. (C) 2002 Wiley-Liss. Inc.

Gene Expression in Brain: A Window on Ethanol Dependence, Neuroadaptation, and Preference

This article represents the proceedings of a symposium at the 2002 joint RSA/ISBRA Conference in San Francisco, California. The organizer was Paula L. Hoffman and the co-chairs were Paula L. Hoffman and Michael Miles. The presentations were (1) Introduction and overview of the use of DNA microarrays, by Michael Miles; (2) DNA microarray analysis of gene expression in brains of P and NP rats, by Howard J. Edenberg; (3) Gene expression patterns in brain regions of AA and ANA rats, by Wolfgang Sommer; (4) Patterns of gene expression in brains of selected lines of mice that differ in ethanol tolerance, by Boris Tabakoff; (5) Gene expression profiling related to initial sensitivity and tolerance in gamma-protein kinase C mutants, by Jeanne Wehner; and (6) Gene expression patterns in human alcoholic brain: from microarrays to protein profiles, by Joanne Lewohl.

Model-based clustering in gene expression microarrays: An application to breast cancer data

In microarray studies, the application of clustering techniques is often used to derive meaningful insights into the data. In the past, hierarchical methods have been the primary clustering tool employed to perform this task. The hierarchical algorithms have been mainly applied heuristically to these cluster analysis problems. Further, a major limitation of these methods is their inability to determine the number of clusters. Thus there is a need for a model-based approach to these. clustering problems. To this end, McLachlan et al. [7] developed a mixture model-based algorithm (EMMIX-GENE) for the clustering of tissue samples. To further investigate the EMMIX-GENE procedure as a model-based -approach, we present a case study involving the application of EMMIX-GENE to the breast cancer data as studied recently in van 't Veer et al. [10]. Our analysis considers the problem of clustering the tissue samples on the basis of the genes which is a non-standard problem because the number of genes greatly exceed the number of tissue samples. We demonstrate how EMMIX-GENE can be useful in reducing the initial set of genes down to a more computationally manageable size. The results from this analysis also emphasise the difficulty associated with the task of separating two tissue groups on the basis of a particular subset of genes. These results also shed light on why supervised methods have such a high misallocation error rate for the breast cancer data.

Influence of ageing, heat shock treatment and in vivo total antioxidant status on gene-expression profile and protein synthesis in human peripheral lymphocytes

Ageing results in a progressive, intrinsic and generalised imbalance of the control of regulatory systems. A key manifestation of this complex biological process includes the attenuation of the universal stress response. Here we provide the first global assessment of the ageing process as it affects the heat shock response, utilising human peripheral lymphocytes and cDNA microarray analysis. The genomic approach employed in our preliminary study was supplemented with a proteomic approach. In addition, the current study correlates the in vivo total antioxidant status with the age-related differential gene expression as well as the translational kinetics of heat shock proteins (hsps). Most of the genes encoding stress response proteins on the 4224 element microarray used in this study were significantly elevated after heat shock treatment of lymphocytes obtained from both young and old individuals albeit to a greater extent in the young. Cell signaling and signal transduction genes as well as some oxidoreductases showed varied response. Results from translational kinetics of induction of major hsps, from 0 to 24 It recovery period were broadly consistent with the differential expression of HSC 70 and HSP 40 genes. Total antioxidant levels in plasma from old individuals were found to be significantly lower by comparison with young, in agreement with the widely acknowledged role of oxidant homeostasis in the ageing process. (C) 2002 Elsevier Science Ireland Ltd. All rights reserved.

Expression analysis of four Pinus radiata male cone promoters in the heterologous host Arabidopsis

Four male cone-specific promoters were isolated from the genome of Pinus radiata D. Don, fused to the beta-glucuronidase (GUS) reporter gene and analysed in the heterologous host Arabidopsis thaliana (L.) Heynh. The temporal and spatial activities of the promoters PrCHS1, PrLTP2, PrMC2 and PrMALE1 during seven anther developmental stages are described in detail. The two promoters PrMC2 and PrMALE1 confer an identical GUS expression pattern on Arabidopsis anthers. DNA sequence analysis of the PrMC2 and PrMALE1 promoters revealed an 88% sequence identity over 276 bp and divergence further upstream (

Polymorphisms in the 5 '-untranslated region of the human serotonin receptor 1B (HTR1B) gene affect gene expression

We present evidence of complex balancing regulation of HTR1B transcription by common polymorphisms in its promoter. Computational analysis of the HTR1B gene predicted that a 50 segment, spanning common DNA sequence variations, T-261G, A-161T, and -182INS/DEL-181, contained a putative functional promoter. Using a secreted alkaline phosphatase (SEAP) reporter gene system, we found that the haplotype -261G_-182INS-181_A-161 enhanced transcriptional activity 2.3-fold compared with the haplotype T-261_-182INS-181_A-161. Conversely, -161T reversed this, and the net effect when -261G and -161T were in the same haplotype (-261G_-182INS-181_-161T) was equivalent to the major haplotype (T-261_-182INS-181_A-161). Electrophoretic mobility shift experiments showed that -261G and -161T modify the binding of transcription factors (TFs): -261G generates a new AP2 binding site, while alleles A-161 and -161T exhibit different binding characteristics to AP1. T-261G and A-161T were found to be in linkage disequilibrium (LD) with G861C in a European ancestry population. Interestingly, G861C has been reported to be associated with several psychiatric disorders. Our results indicate that HTR1B is the target of substantial transcriptional genetic regulation by common haplotypes, which are in LD with the HTR1B single-nucleotide polymorphism (SNP) most commonly used in association studies.

Systemic gene expression in arabidopsis during an incompatible interaction with Alternaria brassicicola

Pathogen challenge can trigger an integrated set of signal transduction pathways, which ultimately leads to a state of high alert, otherwise known as systemic or induced resistance in tissue remote to the initial infection. Although large-scale gene expression during systemic acquired resistance, which is induced by salicylic acid or necrotizing pathogens has been previously reported using a bacterial pathogen, the nature of systemic defense responses triggered by an incompatible necrotrophic fungal pathogen is not known. We examined transcriptional changes that occur during systemic defense responses in Arabidopsis plants inoculated with the incompatible fungal pathogen Alternaria brassicicola. Substantial changes (2.00-fold and statistically significant) were demonstrated in distal tissue of inoculated plants for 35 genes (25 up-regulated and 10 down-regulated), and expression of a selected subset of systemically expressed genes was confirmed using real-time quantitative polymerase chain reaction. Genes with altered expression in distal tissue included those with putative functions in cellular housekeeping, indicating that plants modify these vital processes to facilitate a coordinated response to pathogen attack. Transcriptional up-regulation of genes encoding enzymes functioning in the beta-oxidation pathway of fatty acids was particularly interesting. Transcriptional up-regulation was also observed for genes involved in cell wall synthesis and modification and genes putatively involved in signal transduction. The results of this study, therefore, confirm the notion that distal tissue of a pathogen-challenged plant has a heightened preparedness for subsequent pathogen attacks.

Characterizing embryonic gene expression patterns in the mouse using nonredundant sequence-based selection

This article investigates the expression patterns of 160 genes that are expressed during early mouse development. The cDNAs were isolated from 7.5 d postcoitum (dpc) encloderm, a region that comprises visceral encloderm (VE), definitive encloderm, and the node-tissues that are required for the initial steps of axial specification and tissue patterning in the mouse. To avoid examining the same gene more than once, and to exclude potentially ubiquitously expressed housekeeping genes, cDNA sequence was derived from 1978 clones of the Endoderm library. These yielded 1440 distinct cDNAs, of which 123 proved to be novel in the mouse. In situ hybridization analysis was carried out on 160 of the cDNAs, and of these, 29 (18%) proved to have restricted expression patterns.