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.

The Oligocene-Neogene deep-sea Columbia Channel system in the South Brazilian Basin: Seismic stratigraphy and environmental changes

The Columbia Channel (CCS) system is a depositional system located in the South Brazilian Basin, south of the Vitoria-Trindade volcanic chain. It lies in a WNW-ESE direction on the continental rise and abyssal plain, at a depth of between 4200 and 5200 m. It is formed by two depocenters elongated respectively south and north of the channel that show different sediment patterns. The area is swept by a deep western boundary current formed by AABW. The system has been previously interpreted has a mixed turbidite-contourite system. More detailed study of seismic data permits a more precise definition of the modern channel morphology, the system stratigraphy as well as the sedimentary processes and control. The modern CCS presents active erosion and/or transport along the channel. The ancient Oligo-Neogene system overlies a ""upper Cretaceous-Paleogene"" sedimentary substratum (Unit U1) bounded at the top by a major erosive ""late Eocene-early Oligocene"" discordance (D2). This ancient system is subdivided into 2 seismic units (U2 and U3). The thick basal U2 unit constitutes the larger part of the system. It consists of three subunits bounded by unconformities: D3 (""Oligocene-Miocene boundary""), D4 (""late Miocene"") and D5 (""late Pliocene""). The subunits have a fairly tabular geometry in the shallow NW depocenter associated with predominant turbidite deposits. They present a mounded shape in the deep NE depocenter, and are interpreted as forming a contourite drift. South of the channel, the deposits are interpreted as a contourite sheet drift. The surficial U3 unit forms a thin carpet of deposits. The beginning of the channel occurs at the end of U1 and during the formation of D2. Its location seems to have been determined by active faults. The channel has been active throughout the late Oligocene and Neogene and its depth increased continuously as a consequence of erosion of the channel floor and deposit aggradation along its margins. Such a mixed turbidite-contourite system (or fan drift) is characterized by frequent, rapid lateral facies variations and by unconformities that cross the whole system and are associated with increased AABW circulation. (C) 2009 Elsevier B.V. All rights reserved.

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.

Gene expression profiles displayed by peripheral blood mononuclear cells from patients with type 2 diabetes mellitus focusing on biological processes implicated on the pathogenesis of the disease

Patients with type 2 diabetes mellitus (T2DM) exhibit insulin resistance associated with obesity and inflammatory response, besides an increased level of oxidative DNA damage as a consequence of the hyperglycemic condition and the generation of reactive oxygen species (ROS). In order to provide information on the mechanisms involved in the pathophysiology of T2DM, we analyzed the transcriptional expression patterns exhibited by peripheral blood mononuclear cells (PBMCs) from patients with T2DM compared to non-diabetic subjects, by investigating several biological processes: inflammatory and immune responses, responses to oxidative stress and hypoxia, fatty acid processing, and DNA repair. PBMCs were obtained from 20 T2DM patients and eight non-diabetic subjects. Total RNA was hybridized to Agilent whole human genome 4x44K one-color oligo-microarray. Microarray data were analyzed using the GeneSpring GX 11.0 software (Agilent). We used BRB-ArrayTools software (gene set analysis - GSA) to investigate significant gene sets and the Genomica tool to study a possible influence of clinical features on gene expression profiles. We showed that PBMCs from T2DM patients presented significant changes in gene expression, exhibiting 1320 differentially expressed genes compared to the control group. A great number of genes were involved in biological processes implicated in the pathogenesis of T2DM. Among the genes with high fold-change values, the up-regulated ones were associated with fatty acid metabolism and protection against lipid-induced oxidative stress, while the down-regulated ones were implicated in the suppression of pro-inflammatory cytokines production and DNA repair. Moreover, we identified two significant signaling pathways: adipocytokine, related to insulin resistance; and ceramide, related to oxidative stress and induction of apoptosis. In addition, expression profiles were not influenced by patient features, such as age, gender, obesity, pre/post-menopause age, neuropathy, glycemia, and HbA(1c) percentage. Hence, by studying expression profiles of PBMCs, we provided quantitative and qualitative differences and similarities between T2DM patients and non-diabetic individuals, contributing with new perspectives for a better understanding of the disease. (C) 2012 Elsevier B.V. All rights reserved.

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.

Primary antiphospholipid syndrome: morphofunctional penile abnormalities with normal sperm analysis

Objective: To perform a global gonadal and sexual functions assessment in primary antiphospholipid syndrome (PAPS) patients. Methods: A cross-sectional study was conducted in 12 male PAPS patients and 20 healthy controls. They were assessed by demographic data, clinical features, systematic urological examination, sexual function, testicular ultrasound, seminal parameters according to the World Health Organization (WHO), seminal sperm antibodies, and hormone profile, including follicle stimulating hormone (FSH), luteinizing hormone (LH), morning total testosterone, and thyroid hormones. Results: The median of current age and age of spermarche were similar in PAPS patients and controls (37.5 vs. 32.4 years, p = 0.270, and 13.1 vs. 12.85 years, p = 0.224, respectively), with a higher frequency of erectile dysfunction in the former group (25% vs. 0%, p = 0.044). Further analysis of PAPS patients with and without previous arterial thrombosis demonstrated that the median penis circumference was significantly lower in PAPS with arterial thrombosis than in PAPS without this complication (8.1 [6-10] vs. 10.2 [10-11] cm, p = 0.007). In addition, the median penis circumference was significantly lower in PAPS patients with erectile dysfunction than in patients without this complication (7.5 [6-9.5] vs. 9.5 [7.5-11] cm, p = 0.039). Regarding seminal analysis, the median sperm concentration, sperm motility, and normal sperm forms by WHO guidelines were comparable in PAPS patients and controls (141.5 [33-575] vs. 120.06 [34.5-329] x 106/ml, p = 0.65; 61.29 [25-80] vs. 65.42 [43-82]%, p = 0.4; 21.12 [10-42.5] vs. 23.95 [10-45]%, p = 0.45, respectively), and none of them had oligo/azoospermia. No differences were observed between PAPS patients and controls regarding the frequency of antisperm antibodies, testicular volume by ultrasound, or hormone profile (FSH, LH, morning total testosterone, and thyroid hormone) (p > 0.05). Conclusions: Normal testicular function has been identified in PAPS patients, in spite of morphofunctional penile abnormalities. Previous arterial thrombosis may underlie penile anthropometry alteration. Lupus (2012) 21, 251-256.

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.