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.

Coordinated plant defense responses in Arabidopsis revealed by microarray analysis

Disease resistance is associated with a plant defense response that involves an integrated set of signal transduction pathways. Changes in the expression patterns of 2.375 selected genes were examined simultaneously by cDNA microarray analysis in Arabidopsis thaliana after inoculation with an incompatible fungal pathogen Alternaria brassicicola or treatment with the defense-related signaling molecules salicylic acid (SA), methyl jasmonate (MJ), or ethylene, Substantial changes (up- and down-regulation) in the steady-state abundance of 705 mRNAs were observed in response to one or more of the treatments, including known and putative defense-related genes and 106 genes with no previously described function or homology, In leaf tissue inoculated with A. brassicicola, the abundance of 168 mRNAs was increased more than 2.5-fold, whereas that of 39 mRNAs was reduced. Similarly, the abundance of 192, 221, and 55 mRNAs was highly (>2.5-fold) increased after treatment with SA, MJ, and ethylene, respectively. Data analysis revealed a surprising level of coordinated defense responses, including 169 mRNAs regulated by multiple treatments/defense pathways. The largest number of genes coinduced (one of four induced genes) and corepressed was found after treatments with SA and MJ. In addition, 50% of the genes induced by ethylene treatment were also induced by MJ treatment. These results indicated the existence of a substantial network of regulatory interactions and coordination occurring during plant defense among the different defense signaling pathways, notably between the salicylate and jasmonate pathways that were previously thought to act in an antagonistic fashion.

Construction and validation of a Bovine Innate Immune Microarray

Background: Microarray transcript profiling has the potential to illuminate the molecular processes that are involved in the responses of cattle to disease challenges. This knowledge may allow the development of strategies that exploit these genes to enhance resistance to disease in an individual or animal population. Results: The Bovine Innate Immune Microarray developed in this study consists of 1480 characterised genes identified by literature searches, 31 positive and negative control elements and 5376 cDNAs derived from subtracted and normalised libraries. The cDNA libraries were produced from 'challenged' bovine epithelial and leukocyte cells. The microarray was found to have a limit of detection of 1 pg/mu g of total RNA and a mean slide-to-slide correlation co-efficient of 0.88. The profiles of differentially expressed genes from Concanavalin A ( ConA) stimulated bovine peripheral blood lymphocytes were determined. Three distinct profiles highlighted 19 genes that were rapidly up-regulated within 30 minutes and returned to basal levels by 24 h; 76 genes that were upregulated between 2 - 8 hours and sustained high levels of expression until 24 h and 10 genes that were down-regulated. Quantitative real-time RT-PCR on selected genes was used to confirm the results from the microarray analysis. The results indicate that there is a dynamic process involving gene activation and regulatory mechanisms re-establishing homeostasis in the ConA activated lymphocytes. The Bovine Innate Immune Microarray was also used to determine the cross-species hybridisation capabilities of an ovine PBL sample. Conclusion: The Bovine Innate Immune Microarray has been developed which contains a set of well-characterised genes and anonymous cDNAs from a number of different bovine cell types. The microarray can be used to determine the gene expression profiles underlying innate immune responses in cattle and sheep.

cDNA microarray analysis of differentially expressed genes in penile tissue after treatment with tadalafil

To evaluate differential gene expression in penile tissue after treatment with the phosphodiesterase 5 (PDE5) inhibitor tadalafil, as of the three clinically available PDE5 inhibitors (sildenafil, tadalafil, and vardenafil) used for the treatment of erectile dysfunction (ED), tadalafil has a long half-life and low incidence of side-effects. In all, 32 adult rats were divided into two groups. The control group received 0.5 mL of drinking water alone, while the tadalafil group was treated with tadalafil at a dose of 0.27 mg/kg. At 4 h after treatment with water or tadalafil the rats were killed and the penile tissue was removed. The total RNA was isolated from the penile tissue from both groups and differentially expressed genes were identified by cDNA microarray analysis. To validate the expression data from the microarray analysis, quantitative real-time polymerase chain reaction (PCR) and immunohistochemistry were used. In all, 153 genes were differentially expressed between the control group and the tadalafil group. We validated the microarray results by quantitative PCR for the insulin-like growth factor binding protein 6 (IGFBP-6) gene and the neuronal calcium sensor 1 (NCS-1) gene, both of which were up-regulated in the tadalafil group, and for the natriuretic peptide receptor 1 (NPR-1) gene that was down-regulated in this group. Immunohistochemistry showed localization of the NCS-1 protein in sinusoid trabeculae of the corpus cavernosum in control and tadalafil-treated rats. There was differential expression in 153 genes after tadalafil treatment. Some of these genes such as IGFBP-6, NPR-1 and NCS-1, might result in new targets in the treatment of ED.

Protein profile of the luteal phase endometrium by tissue microarray assessment

To investigate the luteal phase endometrial expression of leukemia inhibitor factor (LIF), insulin-like growth factor 1 (IGF-1), progesterone receptor (PR), claudin 4 (CLDN4), vascular-endothelial growth factor receptor 3 (VEGFR-3), bone morphogenetic protein 4 (BMP-4) and citokeratin 7 (CK-7), we obtained luteal phase endometrial samples from 52 women. Samples were dated and integrated using a tissue microarray (TMA). Samples were immunostained for LIF, IGF-1, PR, CLDN4, VEGFR-3, BMP-4 and CK-7. Frequencies of positive expressions at the early, mid and late luteal phases were compared by two proportions test. Concomitant expression of these proteins was assessed with Chi-square or Fischer`s test. The frequency of LIF was positively correlated to the frequency of IGF-1 (r = 0.99; p < 0.05) and PR (r = 0.99; p < 0.05), and the correlation between IGF-1 and PR tended to be significant (r = 0.98; p < 0.1). The expression of PR was associated with the absence of CLDN4 (p < 0.001). Thus, expression of LIF, IGF-1 and PR are correlated during the luteal phase, and immunohistochemistry for these proteins might be used to assist in the assessment of endometrial maturation. In addition, the expression of CLDN4 and PR was not concomitant, warranting further investigation on the relationship of their endometrial expression.

Using a combination of MLPA kits to detect chromosomal imbalances in patients with multiple congenital anomalies and mental retardation is a valuable choice for developing countries

Conventional karyotyping detects anomalies in 3-15% of patients with multiple congenital anomalies and mental retardation (MCA/MR). Whole-genome array screening (WGAS) has been consistently suggested as the first choice diagnostic test for this group of patients, but it is very costly for large-scale use in developing countries. We evaluated the use of a combination of Multiplex Ligation-dependent Probe Amplification (MLPA) kits to increase the detection rate of chromosomal abnormalities in MCA/MR patients. We screened 261 MCA/MR patients with two subtelomeric and one microdeletion kits. This would theoretically detect up to 70% of all submicroscopic abnormalities. Additionally we scored the de Vries score for 209 patients in an effort to find a suitable cut-off for MLPA screening. Our results reveal that chromosomal abnormalities were present in 87 (33.3%) patients, but only 57 (21.8%) were considered causative. Karyotyping detected 15 abnormalities (6.9%), while MLPA identified 54 (20.7%). Our combined MLPA screening raised the total detection number of pathogenic imbalances more than three times when compared to conventional karyotyping. We also show that using the de Vries score as a cutoff for this screening would only be suitable under financial restrictions. A decision analytic model was constructed with three possible strategies: karyotype, karyotype + MLPA and karyotype + WGAS. Karyotype + MLPA strategy detected anomalies in 19.8% of cases which account for 76.45% of the expected yield for karyotype + WGAS. Incremental Cost Effectiveness Ratio (ICER) of MLPA is three times lower than that of WGAS, which means that, for the same costs, we have three additional diagnoses with MLPA but only one with WGAS. We list all causative alterations found, including rare findings, such as reciprocal duplications of regions deleted in Sotos and Williams-Beuren syndromes. We also describe imbalances that were considered polymorphisms or rare variants, such as the new SNP that confounded the analysis of the 22q13.3 deletion syndrome. (C) 2011 Elsevier Masson SAS. All rights reserved.

Gene expression profile analysis of primary glioblastomas and non-neoplastic brain tissue: identification of potential target genes by oligonucleotide microarray and real-time quantitative PCR

The prognosis of glioblastomas is still extremely poor and the discovery of novel molecular therapeutic targets can be important to optimize treatment strategies. Gene expression analyses comparing normal and neoplastic tissues have been used to identify genes associated with tumorigenesis and potential therapeutic targets. We have used this approach to identify differentially expressed genes between primary glioblastomas and non-neoplastic brain tissues. We selected 20 overexpressed genes related to cell cycle, cellular movement and growth, proliferation and cell-to-cell signaling and analyzed their expression levels by real time quantitative PCR in cDNA obtained from microdissected fresh tumor tissue from 20 patients with primary glioblastomas and from 10 samples of non-neoplastic white matter tissue. The gene expression levels were significantly higher in glioblastomas than in non-neoplastic white matter in 18 out of 20 genes analyzed: P < 0.00001 for CDKN2C, CKS2, EEF1A1, EMP3, PDPN, BNIP2, CA12, CD34, CDC42EP4, PPIE, SNAI2, GDF15 and MMP23b; and NFIA (P: 0.0001), GPS1 (P: 0.0003), LAMA1 (P: 0.002), STIM1 (P: 0.006), and TASP1 (P: 0.01). Five of these genes are located in contiguous loci at 1p31-36 and 2 at 17q24-25 and 8 of them encode surface membrane proteins. PDPN and CD34 protein expression were evaluated by immunohistochemistry and they showed concordance with the PCR results. The present results indicate the presence of 18 overexpressed genes in human primary glioblastomas that may play a significant role in the pathogenesis of these tumors and that deserve further functional investigation as attractive candidates for new therapeutic targets.

Genomic deletions at 1p and 14q are associated with an abnormal cDNA microarray gene expression pattern in meningiomas but not in schwannomas

The molecular pathology of meningiomas and shwannomas involve the inactivation of the NF2 gene to generate grade I tumors. Genomic losses at 1p and 14q are observed in both neoplasms, although more frequently in meningiomas. The inactivation of unidentified genes located in these regions appears associated with tumor progression in meningiomas, but no clues to its molecular/clinical meaning are available in schwannomas. Recent microarray gene expression studies have demonstrated the existence of molecular subgroups in both entities. In the present study, we correlated the presence of genomic deletions at 1p, 14q, and 22q with the expression patterns of 96 tumor-related genes obtained by cDNA low-density microarrays in a series of 65 tumors including 42 meningiomas and 23 schwannomas. Two expression pattern groups were identified by cDNA mycroarray analysis when compared to the expression pattern in normal control RNA in both meningiomas and schwannomas, each one with patterns similar and different from the normal control. Meningioma and schwannoma subgroups differed in the expression of 38 and 16 genes, respectively. Using MLPA and microsatellites, we identified genomic losses at 1p, 14q, and 22q at nonrandom frequencies (12.5-69%) in meningiomas and schwannomas. Losses at 22q were almost equally frequent in both molecular expression subgroups in both neoplasms. However, deletions at 1p and 14q accumulated in meningiomas with a gene expression pattern different from the normal pattern, whereas the inverse situation occurred in schwannomas. Those anomalies characterized the schwannomas with expression pattern similar to the normal control. These findings suggest that deletions at 1p and 14q enhance the development of an abnormal tumor-related gene expression pattern in meningiomas, but this fact is not corroborated in schwannomas. (C) 2010 Elsevier Inc. All rights reserved.

Loss of samples in the tissue microarray technique: comparison between slides using adhesive tape and silanized slides

The tissue microarray (TMA) technique allows multiple tissue samples in a single block. Commercial adhesive tape is used to avoid the loss of tissue samples during the immunostaining process. Few reports exist in the literature comparing the use of these adhesive tapes to other adhesive techniques. The objective of this study was to compare loss of sections adhered to slides using commercial adhesive tapes versus using silanized only slides. TMA was constructed with varying tissues using a fixed-base device (Beecher Instruments), placing 108 cylinders of 1 mm diameter in duplicate, spaced 1.2 mm apart. Section of 4 mu m were cut from the TMA block and adhered to 30 silanized slides and 30 commercial glass slides using adhesive tape, according to manufacturer`s recommendations. Vimentin immunoexpression was evaluated by immunohistochemistry. Antigenic recovery was realized in citrate buffer using a microwave oven. Cylinder loss in the immunohistochemical process was quantified and expressed as: total (>80%), almost complete (75-79%), or partial (50-74%). The commercial adhesive tape group presented lesser total loss (1.1 versus 6.4%), almost complete loss (2.2 versus 3.5%), and partial loss (2.1 versus 3.8%) than the silanized slide group (ANOVA, P < 0.05). The sum of total and almost complete losses in the silanized slide group was 9.9%, greater than the losses in slides using commercial adhesive tapes (3.3%) and less than reported and considered acceptable in the literature (10-30%). In conclusion, the use of silanized only slides presents very satisfactory results, requires less training, and reduces costs significantly, thus justifying their use in research.

Nucleophosmin, p53, and Ki-67 expression patterns on an oral squamous cell carcinoma tissue microarray

Oral cancer is the eighth most prevalent cancer worldwide. It causes significant mortality and morbidity rates, which have motivated the search for prognostic factors to better tailor the individual management of oral squamous cell carcinoma patients. Nucleophosmin is a multifunctional protein that is involved in many cellular activities, such as, regulation of the tumor suppressor genes TP53 and p14(ARF). and is associated with proliferative and growth suppressive roles in the cell. Nucleophosmin is overexpressed in many solid tumors in human, including tumors of the colon, liver, stomach, ovary, and prostate. In this study, we analyzed the expression of nucleophosmin, Ki-67, and p53 by immunohistochemistry in oral squamous cell carcinomas. Less than 10% of nuclear staining was observed in 90.3%, 50.6%, and 65.3% of the cases for nucleophosmin, p53, and Ki-67, respectively. Expression of p53 was not significantly associated with any of the clinicopathologic parameters analyzed. Increased expression of Ki-67 was associated with the presence of lymph node metastasis (P < .0001), advanced stages of disease (P = .0030), tumors occurring in the floor of mouth (P = .0018), and moderately/well-differentiated tumors (P = .0287). Local recurrence was associated with higher expression of nucleophosmin (P = .0233), and disease-free survival rate was significantly better in patients with low expression of nucleophosmin. Multivariate analysis suggested that expression of nucleophosmin could be an independent prognostic factor for oral squamous cell carcinoma patients. (C) 2010 Elsevier Inc. All rights reserved.

Oral squamous cell carcinoma: status of tight junction claudins in the different histopathological patterns and relationship with clinical parameters. A tissue-microarray-based study of 136 cases

Aims Claudins are integral transmembrane proteins of the tight junctions, critical for maintaining cell adhesion and polarity. Alterations in the expression of individual claudins have been detected in carcinomas and appear to correlate with tumour progression. Methods In this study, a panel of anti-claudin antibodies (anti-claudins 1, 2, 3, 4, 5 and 7) was employed to map claudin expression in 136 cases of oral squamous cell carcinoma (OSCC) organised in a tissue microarray. Results Claudins were expressed in a reticular pattern up to the prickle layer in normal mucosal epithelium. In OSCC, claudins were strongly present in well-differentiated tumours, they presented mild and low expression in moderately differentiated OSCC, and were negative in poorly differentiated OSCC; the absences of claudin 1 (p = 0.002) and claudin 4 (p<0.001) were associated with moderately/poorly differentiated tumours. Strong expression of claudin 4 was associated with decreased perineural infiltration (p = 0.024). Claudins 5 and 7 were mostly negative or weakly expressed in all cases studied. Expression of claudin 7 was associated with the early clinical stages of the disease, whereas loss of claudin 7 tended to be more frequent in advanced stages of OSCC (p = 0.054). Absence of claudin 7 was also associated with absent vascular infiltration (p = 0.045) and with presence of recurrence (p = 0.052). Conclusions Claudin expression patterns showed a strong correlation with histological type of OSCC; claudin expression was decreased in areas of invasion, and negative in poorly differentiated tumours. This pattern may be related to evolution and prognosis of these tumours, especially in the case of claudin 7, which seems to be associated with a poor prognosis in OSCC.

A mixture model-based approach to the clustering of microarray expression data

Motivation: This paper introduces the software EMMIX-GENE that has been developed for the specific purpose of a model-based approach to the clustering of microarray expression data, in particular, of tissue samples on a very large number of genes. The latter is a nonstandard problem in parametric cluster analysis because the dimension of the feature space (the number of genes) is typically much greater than the number of tissues. A feasible approach is provided by first selecting a subset of the genes relevant for the clustering of the tissue samples by fitting mixtures of t distributions to rank the genes in order of increasing size of the likelihood ratio statistic for the test of one versus two components in the mixture model. The imposition of a threshold on the likelihood ratio statistic used in conjunction with a threshold on the size of a cluster allows the selection of a relevant set of genes. However, even this reduced set of genes will usually be too large for a normal mixture model to be fitted directly to the tissues, and so the use of mixtures of factor analyzers is exploited to reduce effectively the dimension of the feature space of genes. Results: The usefulness of the EMMIX-GENE approach for the clustering of tissue samples is demonstrated on two well-known data sets on colon and leukaemia tissues. For both data sets, relevant subsets of the genes are able to be selected that reveal interesting clusterings of the tissues that are either consistent with the external classification of the tissues or with background and biological knowledge of these sets.

Development of polyaniline microarray electrodes for cadmium analysis

Disposable screen-printed electrodes (SPCE) were modified using a cosmetic product to partially block the electrode surface in order to obtain a microelectrode array. The microarrays formed were electropolymerized with aniline. Scanning electron microscopy was used to evaluate the modified and polymerized electrode surface. Electrochemical characteristics of the constructed sensor for cadmium analysis were evaluated by cyclic and square-wave voltammetry. Optimized stripping procedure in which the preconcentration of cadmium was achieved by depositing at –1.20 V (vs. Ag/AgCl) resulted in a well defined anodic peak at approximately –0.7 V at pH 4.6. The achieved limit of detection was 4 × 10−9 mol dm−3. Spray modified and polymerized microarray electrodes were successfully applied to quantify cadmium in fish sample digests.