Evaluation of reference-based two-color methods for measurement of gene expression ratios using spotted cDNA microarrays

Abstract Background Spotted cDNA microarrays generally employ co-hybridization of fluorescently-labeled RNA targets to produce gene expression ratios for subsequent analysis. Direct comparison of two RNA samples in the same microarray provides the highest level of accuracy; however, due to the number of combinatorial pair-wise comparisons, the direct method is impractical for studies including large number of individual samples (e.g., tumor classification studies). For such studies, indirect comparisons using a common reference standard have been the preferred method. Here we evaluated the precision and accuracy of reconstructed ratios from three indirect methods relative to ratios obtained from direct hybridizations, herein considered as the gold-standard. Results We performed hybridizations using a fixed amount of Cy3-labeled reference oligonucleotide (RefOligo) against distinct Cy5-labeled targets from prostate, breast and kidney tumor samples. Reconstructed ratios between all tissue pairs were derived from ratios between each tissue sample and RefOligo. Reconstructed ratios were compared to (i) ratios obtained in parallel from direct pair-wise hybridizations of tissue samples, and to (ii) reconstructed ratios derived from hybridization of each tissue against a reference RNA pool (RefPool). To evaluate the effect of the external references, reconstructed ratios were also calculated directly from intensity values of single-channel (One-Color) measurements derived from tissue sample data collected in the RefOligo experiments. We show that the average coefficient of variation of ratios between intra- and inter-slide replicates derived from RefOligo, RefPool and One-Color were similar and 2 to 4-fold higher than ratios obtained in direct hybridizations. Correlation coefficients calculated for all three tissue comparisons were also similar. In addition, the performance of all indirect methods in terms of their robustness to identify genes deemed as differentially expressed based on direct hybridizations, as well as false-positive and false-negative rates, were found to be comparable. Conclusion RefOligo produces ratios as precise and accurate as ratios reconstructed from a RNA pool, thus representing a reliable alternative in reference-based hybridization experiments. In addition, One-Color measurements alone can reconstruct expression ratios without loss in precision or accuracy. We conclude that both methods are adequate options in large-scale projects where the amount of a common reference RNA pool is usually restrictive.

FGFR4 Profile as a Prognostic Marker in Squamous Cell Carcinoma of the Mouth and Oropharynx

Background: Fibroblast growth factor receptor 4 (FGFR4) is a member of a receptor tyrosine kinase family of enzymes involved in cell cycle control and proliferation. A common single nucleotide polymorphism (SNP) Gly388Arg variant has been associated with increased tumor cell motility and progression of breast cancer, head and neck cancer and soft tissue sarcomas. The present study evaluated the prognostic significance of FGFR4 in oral and oropharynx carcinomas, finding an association of FGFR4 expression and Gly388Arg genotype with tumor onset and prognosis. Patients and Methods: DNA from peripheral blood of 122 patients with oral and oropharyngeal squamous cell carcinomas was used to determine FGFR4 genotype by PCR-RFLP. Protein expression was assessed by immunohistochemistry (IHC) on paraffin-embedded tissue microarrays. Results: Presence of allele Arg388 was associated with lymphatic embolization and with disease related premature death. In addition, FGFR4 low expression was related with lymph node positivity and premature relapse of disease, as well as disease related death. Conclusion: Our results propose FGFR4 profile, measured by the Gly388Arg genotype and expression, as a novel marker of prognosis in squamous cell carcinoma of the mouth and oropharynx.

Structural synaptic elements are differentially regulated in superior temporal cortex of schizophrenia patients

Inaccurate wiring and synaptic pathology appear to be major hallmarks of schizophrenia. A variety of gene products involved in synaptic neurotransmission and receptor signaling are differentially expressed in brains of schizophrenia patients. However, synaptic pathology may also develop by improper expression of intra- and extra-cellular structural elements weakening synaptic stability. Therefore, we have investigated transcription of these elements in the left superior temporal gyrus of 10 schizophrenia patients and 10 healthy controls by genome-wide microarrays (Illumina). Fourteen up-regulated and 22 downregulated genes encoding structural elements were chosen from the lists of differentially regulated genes for further qRT-PCR analysis. Almost all genes confirmed by this method were downregulated. Their gene products belonged to vesicle-associated proteins, that is, synaptotagmin 6 and syntaxin 12, to cytoskeletal proteins, like myosin 6, pleckstrin, or to proteins of the extracellular matrix, such as collagens, or laminin C3. Our results underline the pivotal roles of structural genes that control formation and stabilization of pre- and post-synaptic elements or influence axon guidance in schizophrenia. The glial origin of collagen or laminin highlights the close interrelationship between neurons and glial cells in establishment and maintenance of synaptic strength and plasticity. It is hypothesized that abnormal expression of these and related genes has a major impact on the pathophysiology of schizophrenia.

Transcriptional response to GAA deficiency (Pompe disease) in infantile-onset patients

Pompe disease is a genetic disorder resulting from a deficiency of lysosomal acid alpha-glucosidase (GAA) that manifests as a clinical spectrum with regard to symptom severity and rate of progression. In this study, we used microarrays to examine gene expression from the muscle of two cohorts of infantile-onset Pompe patients to identify transcriptional differences that may contribute to the disease phenotype. We found strong similarities among the gene expression profiles generated from biceps and quadriceps, and identified a number of signaling pathways altered in both cohorts. We also found that infantile-onset Pompe patient muscle had a gene expression pattern characteristic of immature or regenerating muscle, and exhibited many transcriptional markers of inflammation, despite having few overt signs of inflammatory infiltrate. Further, we identified genes exhibiting correlation between expression at baseline and response to therapy. This combined dataset can serve as a foundation for biological discovery and biomarker development to improve the treatment of Pompe disease. (C) 2012 Elsevier Inc. All rights reserved.

Coordinated expression of oestrogen and androgen receptors in HER2-positive breast carcinomas: impact on proliferative activity

Aims Human epidermal growth factor receptor 2 (HER2)-positive breast cancers are aggressive neoplasms associated with a variable response to systemic therapies. Therefore, the identification of biomarkers to better characterise this heterogeneity would improve treatment efficacy. The aim of this study was to evaluate the influence of androgen receptor (AR) and oestrogen receptor (ER) on clinicopathological features in a series of HER2-positive breast carcinomas. Methods A total of 104 carcinomas were selected and reviewed. Immunohistochemical studies for ER, progesterone receptor and Ki-67 were analysed on tumour whole histological sections. AR expression was analysed on samples represented on tissue microarrays. According to steroid receptor expression, cases were classified into three groups: AR positive/ER positive (48 cases), AR positive/ER negative (41 cases) and AR negative/ER negative (13 cases). Results AR-positive tumours corresponded to 89 (85.6%) of 104 carcinomas. AR-positive carcinomas were associated with a higher frequency of ER and progesterone receptor co-expression and lower proliferative activity determined by the expression of Ki-67. AR-negative carcinomas were more often high grade. The group of AR-positive/ER-negative carcinomas was associated with the highest frequency of apocrine morphological features. The group of AR-negative/ER-negative carcinomas was associated with the highest proliferative activity and the highest frequency of high histological and nuclear grade. The lowest frequency of high-grade tumours and the lowest proliferative activity were seen among tumours with expression of both receptors. Conclusions These results suggest that co-expression of AR and ER can provide a protective effect based on phenotypical presentation of HER2-positive carcinomas. Furthermore, lack of both steroid hormone receptors characterises the most aggressive phenotype.

Homeobox gene expression profile indicates HOXA5 as a candidate prognostic marker in oral squamous cell carcinoma

The search for molecular markers to improve diagnosis, individualize treatment and predict behavior of tumors has been the focus of several studies. This study aimed to analyze homeobox gene expression profile in oral squamous cell carcinoma (OSCC) as well as to investigate whether some of these genes are relevant molecular markers of prognosis and/or tumor aggressiveness. Homeobox gene expression levels were assessed by microarrays and qRT-PCR in OSCC tissues and adjacent non-cancerous matched tissues (margin), as well as in OSCC cell lines. Analysis of microarray data revealed the expression of 147 homeobox genes, including one set of six at least 2-fold up-regulated, and another set of 34 at least 2-fold down-regulated homeobox genes in OSCC. After qRT-PCR assays, the three most up-regulated homeobox genes (HOXA5, HOXD10 and HOXD11) revealed higher and statistically significant expression levels in OSCC samples when compared to margins. Patients presenting lower expression of HOXA5 had poorer prognosis compared to those with higher expression (P=0.03). Additionally, the status of HOXA5, HOXD10 and HOXD11 expression levels in OSCC cell lines also showed a significant up-regulation when compared to normal oral keratinocytes. Results confirm the presence of three significantly upregulated (>4-fold) homeobox genes (HOXA5, HOXD10 and HOXD11) in OSCC that may play a significant role in the pathogenesis of these tumors. Moreover, since lower levels of HOXA5 predict poor prognosis, this gene may be a novel candidate for development of therapeutic strategies in OSCC.

Germline DNA copy number variation in familial and early-onset breast cancer

Introduction: Genetic factors predisposing individuals to cancer remain elusive in the majority of patients with a familial or clinical history suggestive of hereditary breast cancer. Germline DNA copy number variation (CNV) has recently been implicated in predisposition to cancers such as neuroblastomas as well as prostate and colorectal cancer. We evaluated the role of germline CNVs in breast cancer susceptibility, in particular those with low population frequencies (rare CNVs), which are more likely to cause disease." Methods: Using whole-genome comparative genomic hybridization on microarrays, we screened a cohort of women fulfilling criteria for hereditary breast cancer who did not carry BRCA1/BRCA2 mutations. Results: The median numbers of total and rare CNVs per genome were not different between controls and patients. A total of 26 rare germline CNVs were identified in 68 cancer patients, however, a proportion that was significantly different (P = 0.0311) from the control group (23 rare CNVs in 100 individuals). Several of the genes affected by CNV in patients and controls had already been implicated in cancer. Conclusions: This study is the first to explore the contribution of germline CNVs to BRCA1/2-negative familial and early-onset breast cancer. The data suggest that rare CNVs may contribute to cancer predisposition in this small cohort of patients, and this trend needs to be confirmed in larger population samples.

Whole blood genome-wide gene expression profile in males after prolonged wakefulness and sleep recovery

Pellegrino R, Sunaga DY, Guindalini C, Martins RC, Mazzotti DR, Wei Z, Daye ZJ, Andersen ML, Tufik S. Whole blood genome-wide gene expression profile in males after prolonged wakefulness and sleep recovery. Physiol Genomics 44: 1003-1012, 2012. First published September 4, 2012; doi: 10.1152/physiolgenomics.00058.2012.-Although the specific functions of sleep have not been completely elucidated, the literature has suggested that sleep is essential for proper homeostasis. Sleep loss is associated with changes in behavioral, neurochemical, cellular, and metabolic function as well as impaired immune response. Using high-resolution microarrays we evaluated the gene expression profiles of healthy male volunteers who underwent 60 h of prolonged wakefulness (PW) followed by 12 h of sleep recovery (SR). Peripheral whole blood was collected at 8 am in the morning before the initiation of PW (Baseline), after the second night of PW, and one night after SR. We identified over 500 genes that were differentially expressed. Notably, these genes were related to DNA damage and repair and stress response, as well as diverse immune system responses, such as natural killer pathways including killer cell lectin-like receptors family, as well as granzymes and T-cell receptors, which play important roles in host defense. These results support the idea that sleep loss can lead to alterations in molecular processes that result in perturbation of cellular immunity, induction of inflammatory responses, and homeostatic imbalance. Moreover, expression of multiple genes was downregulated following PW and upregulated after SR compared with PW, suggesting an attempt of the body to re-establish internal homeostasis. In silico validation of alterations in the expression of CETN3, DNAJC, and CEACAM genes confirmed previous findings related to the molecular effects of sleep deprivation. Thus, the present findings confirm that the effects of sleep loss are not restricted to the brain and can occur intensely in peripheral tissues.

Transcription profiling of Prss16 (Tssp) can be used to find additional peptidase genes that are candidates for self-peptide generation in the thymus

Positive selection (PS) in the thymus involves the presentation of self-peptides that are bound to MHC class II on the surface of cortical thymus epithelial cells (cTECs). Prss16 gene corresponds to one important element regulating the PS of CD4(+) T lymphocytes, which encodes Thymus-specific serine protease (Tssp), a cTEC serine-type peptidase involved in the proteolytic generation of self-peptides. Nevertheless, additional peptidase genes participating in the generation of self-peptides need to be found. Because of its role in the mechanism of PS and its expression in cTECs, the Prss16 gene might be used as a transcriptional marker to identify new genes that share the same expression profile and that encode peptidases in the thymus. To test this hypothesis, we compared the differential thymic expression of 4,500 mRNAs of wild-type (WT) C57BL/6 mice with their respective Prss16-knockout (KO) mutants by using microarrays. From these, 223 genes were differentially expressed, of which 115 had known molecular/biological functions. Four endopeptidase genes (Casp1, Casp2, Psmb3 and Tpp2) share the same expression profile as the Prss16 gene; i.e., induced in WT and repressed in KO while one endopeptidase gene, Capns1, features opposite expression profile. The Tpp2 gene is highlighted because it encodes a serine-type endopeptidase functionally similar to the Tssp enzyme. Profiling of the KO mice featured down-regulation of Prss16, as expected, along with the genes mentioned above. Considering that the Prss16-KO mice featured impaired PS, the shared regulation of the four endopeptidase genes suggested their participation in the mechanism of self-peptide generation and PS.

Bladder expression of CD cell surface antigens and cell-type-specific transcriptomes

Many cell types have no known functional attributes. In the bladder and prostate, basal epithelial and stromal cells appear similar in cytomorphology and share several cell surface markers. Their total gene expression (transcriptome) should provide a clear measure of the extent to which they are alike functionally. Since urologic stromal cells are known to mediate organ-specific tissue formation, these cells in cancers might exhibit aberrant gene expression affecting their function. For transcriptomes, cluster designation (CD) antigens have been identified for cell sorting. The sorted cell populations can be analyzed by DNA microarrays. Various bladder cell types have unique complements of CD molecules. CD9(+) urothelial, CD104(+) basal and CD13(+) stromal cells of the lamina propria were therefore analyzed, as were CD9(+) cancer and CD13(+) cancer-associated stromal cells. The transcriptome datasets were compared by principal components analysis for relatedness between cell types; those with similarity in gene expression indicated similar function. Although bladder and prostate basal cells shared CD markers such as CD104, CD44 and CD49f, they differed in overall gene expression. Basal cells also lacked stem cell gene expression. The bladder luminal and stromal transcriptomes were distinct from their prostate counterparts. In bladder cancer, not only the urothelial but also the stromal cells showed gene expression alteration. The cancer process in both might thus involve defective stromal signaling. These cell-type transcriptomes provide a means to monitor in vitro models in which various CD-isolated cell types can be combined to study bladder differentiation and bladder tumor development based on cell-cell interaction.

Relationship between B-Cell-specific moloney murine leukemia virus integration site 1 (BMI-1) and homologous recombination regulatory genes in invasive ductal breast carcinomas

B-cell-specific Moloney murine leukemia virus integration site 1 (Bmi-1) is a Polycomb group protein that is able to induce telomerase activity, enabling the immortalization of epithelial cells. Immortalized cells are more susceptible to double-strand breaks (DSB), which are subsequently repaired by homologous recombination (HR). BRCA1 is among the HR regulatory genes involved in the response to DNA damage associated with the RAD51 protein, which accumulates in DNA damage foci after signaling H2AX, another important marker of DNA damage. Topoisomerase III beta (topoIII beta) removes HR intermediates before chromosomal segregation, preventing damage to cellular DNA structure. In breast carcinomas positive for BMI-1 the role of proteins involved in HR remains to be investigated. The aim of this study was to evaluate the association between BMI-1 and homologous recombination proteins. Using tissue microarrays containing 239 cases of primary breast tumors, the expression of Bmi-1, BRCA-1, H2AX, Rad51, p53, Ki-67, topoIII beta, estrogen receptors (ER), progesterone receptors (PR), and HER-2 was analyzed by immunohistochemistry. We observed high Bmi-1 expression in 66 cases (27.6%). Immunohistochemical overexpression of BMI-1 was related to ER (p=0.004), PR (p<0.001), Ki-67 (p<0.001), p53 (p=0.003), BRCA1 (p=0.003), H2AX (p=0.024) and topoIII beta (p<0,001). Our results show a relationship between the expression of BMI-1 and HR regulatory genes, suggesting that Bmi-1 overexpression might be an important event in HR regulation. However, further studies are necessary to understand the mechanisms in which Bmi-1 could regulate HR pathways in invasive ductal breast carcinomas.

Effect of human TGF-beta on the gene expression profile of Schistosoma mansoni adult worms

Schistosoma mansoni is responsible for schistosomiasis, a parasitic disease that affects 200 million people worldwide. Molecular mechanisms of host-parasite interaction are complex and involve a crosstalk between host signals and parasite receptors. TGF-beta signaling pathway has been shown to play an important role in S. mansoni development and embryogenesis. In particular human (h) TGF-beta has been shown to bind to a S. mansoni receptor, transduce a signal that regulates the expression of a schistosome target gene. Here we describe 381 parasite genes whose expression levels are affected by in vitro treatment with hTGF-beta. Among these differentially expressed genes we highlight genes related to morphology, development and cell cycle that could be players of cytokine effects on the parasite. We confirm by qPCR the expression changes detected with microarrays for 5 out of 7 selected genes. We also highlight a set of non-coding RNAs transcribed from the same loci of protein-coding genes that are differentially expressed upon hTCF-beta treatment. These datasets offer potential targets to be explored in order to understand the molecular mechanisms behind the possible role of hTGF-beta effects on parasite biology. (C) 2012 Elsevier B.V. All rights reserved.

Differential Gene Expression Profiles May Differentiate Responder and Nonresponder Patients with Rheumatoid Arthritis for Methotrexate (MTX) Monotherapy and MTX plus Tumor Necrosis Factor Inhibitor Combined Therapy

Objective. We aimed to evaluate whether the differential gene expression profiles of patients with rheumatoid arthritis (RA) could distinguish responders from nonresponders to methotrexate (MTX) and, in the case of MTX nonresponders, responsiveness to MTX plus anti-tumor necrosis factor-alpha (anti-TNF) combined therapy. Methods. We evaluated 25 patients with RA taking MTX 15-20 mg/week as a monotherapy (8 responders and 17 nonresponders). All MTX nonresponders received intliximab and were reassessed after 20 weeks to evaluate their anti-TNF responsiveness using the European League Against Rheumatism response criteria. A differential gene expression analysis from peripheral blood mononuclear cells was performed in terms of hierarchical gene clustering, and an evaluation of differentially expressed genes was performed using the significance analysis of microarrays program. Results. Hierarchical gene expression clustering discriminated MTX responders from nonresponders, and MTX plus anti-TNF responders from nonresponders. The evaluation of only highly modulated genes (fold change > 1.3 or < 0.7) yielded 5 induced (4 antiapoptotic and CCL4) and 4 repressed (4 proapoptotic) genes in MTX nonresponders compared to responders. In MTX plus anti-TNF nonresponders, the CCL4, CD83, and BCL2A1 genes were induced in relation to responders. Conclusion. Study of the gene expression profiles of RA peripheral blood cells permitted differentiation of responders from nonresponders to MTX and anti-TNF. Several candidate genes in MTX non-responders (CCL4, HTRA2, PRKCD, BCL2A1, CAV1, TNIP1 CASP8AP2, MXD1, and BTG2) and 3 genes in MTX plus anti-TNF nonresponders (CCL4, CD83, and BCL2A1) were identified for further study. (First Release July 1 2012; J Rheumatol 2012;39:1524-32; doi:10.3899/jrheum.120092)

Identification of protein expression signatures in gastric carcinomas using clustering analysis

Background and Aim: The identification of gastric carcinomas (GC) has traditionally been based on histomorphology. Recently, DNA microarrays have successfully been used to identify tumors through clustering of the expression profiles. Random forest clustering is widely used for tissue microarrays and other immunohistochemical data, because it handles highly-skewed tumor marker expressions well, and weighs the contribution of each marker according to its relatedness with other tumor markers. In the present study, we e identified biologically- and clinically-meaningful groups of GC by hierarchical clustering analysis of immunohistochemical protein expression. Methods: We selected 28 proteins (p16, p27, p21, cyclin D1, cyclin A, cyclin B1, pRb, p53, c-met, c-erbB-2, vascular endothelial growth factor, transforming growth factor [TGF]-beta I, TGF-beta II, MutS homolog-2, bcl-2, bax, bak, bcl-x, adenomatous polyposis coli, clathrin, E-cadherin, beta-catenin, mucin (MUC) 1, MUC2, MUC5AC, MUC6, matrix metalloproteinase [ MMP]-2, and MMP-9) to be investigated by immunohistochemistry in 482 GC. The analyses of the data were done using a random forest-clustering method. Results: Proteins related to cell cycle, growth factor, cell motility, cell adhesion, apoptosis, and matrix remodeling were highly expressed in GC. We identified protein expressions associated with poor survival in diffuse-type GC. Conclusions: Based on the expression analysis of 28 proteins, we identified two groups of GC that could not be explained by any clinicopathological variables, and a subgroup of long-surviving diffuse-type GC patients with a distinct molecular profile. These results provide not only a new molecular basis for understanding the biological properties of GC, but also better prediction of survival than the classic pathological grouping.

Evaluating different methods of microarray data normalization

Abstract Background With the development of DNA hybridization microarray technologies, nowadays it is possible to simultaneously assess the expression levels of thousands to tens of thousands of genes. Quantitative comparison of microarrays uncovers distinct patterns of gene expression, which define different cellular phenotypes or cellular responses to drugs. Due to technical biases, normalization of the intensity levels is a pre-requisite to performing further statistical analyses. Therefore, choosing a suitable approach for normalization can be critical, deserving judicious consideration. Results Here, we considered three commonly used normalization approaches, namely: Loess, Splines and Wavelets, and two non-parametric regression methods, which have yet to be used for normalization, namely, the Kernel smoothing and Support Vector Regression. The results obtained were compared using artificial microarray data and benchmark studies. The results indicate that the Support Vector Regression is the most robust to outliers and that Kernel is the worst normalization technique, while no practical differences were observed between Loess, Splines and Wavelets. Conclusion In face of our results, the Support Vector Regression is favored for microarray normalization due to its superiority when compared to the other methods for its robustness in estimating the normalization curve.