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.

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.

The full Bayesian significance test for mixture models: results in gene expression clustering

Gene clustering is a useful exploratory technique to group together genes with similar expression levels under distinct cell cycle phases or distinct conditions. It helps the biologist to identify potentially meaningful relationships between genes. In this study, we propose a clustering method based on multivariate normal mixture models, where the number of clusters is predicted via sequential hypothesis tests: at each step, the method considers a mixture model of m components (m = 2 in the first step) and tests if in fact it should be m - 1. If the hypothesis is rejected, m is increased and a new test is carried out. The method continues (increasing m) until the hypothesis is accepted. The theoretical core of the method is the full Bayesian significance test, an intuitive Bayesian approach, which needs no model complexity penalization nor positive probabilities for sharp hypotheses. Numerical experiments were based on a cDNA microarray dataset consisting of expression levels of 205 genes belonging to four functional categories, for 10 distinct strains of Saccharomyces cerevisiae. To analyze the method's sensitivity to data dimension, we performed principal components analysis on the original dataset and predicted the number of classes using 2 to 10 principal components. Compared to Mclust (model-based clustering), our method shows more consistent results.

Searching for molecular markers in head and neck squamous cell carcinomas (HNSCC) by statistical and bioinformatic analysis of larynx-derived SAGE libraries

Background: Head and neck squamous cell carcinoma (HNSCC) is one of the most common malignancies in humans. The average 5-year survival rate is one of the lowest among aggressive cancers, showing no significant improvement in recent years. When detected early, HNSCC has a good prognosis, but most patients present metastatic disease at the time of diagnosis, which significantly reduces survival rate. Despite extensive research, no molecular markers are currently available for diagnostic or prognostic purposes. Methods: Aiming to identify differentially-expressed genes involved in laryngeal squamous cell carcinoma (LSCC) development and progression, we generated individual Serial Analysis of Gene Expression (SAGE) libraries from a metastatic and non-metastatic larynx carcinoma, as well as from a normal larynx mucosa sample. Approximately 54,000 unique tags were sequenced in three libraries. Results: Statistical data analysis identified a subset of 1,216 differentially expressed tags between tumor and normal libraries, and 894 differentially expressed tags between metastatic and non-metastatic carcinomas. Three genes displaying differential regulation, one down-regulated (KRT31) and two up-regulated (BST2, MFAP2), as well as one with a non-significant differential expression pattern (GNA15) in our SAGE data were selected for real-time polymerase chain reaction (PCR) in a set of HNSCC samples. Consistent with our statistical analysis, quantitative PCR confirmed the upregulation of BST2 and MFAP2 and the downregulation of KRT31 when samples of HNSCC were compared to tumor-free surgical margins. As expected, GNA15 presented a non-significant differential expression pattern when tumor samples were compared to normal tissues. Conclusion: To the best of our knowledge, this is the first study reporting SAGE data in head and neck squamous cell tumors. Statistical analysis was effective in identifying differentially expressed genes reportedly involved in cancer development. The differential expression of a subset of genes was confirmed in additional larynx carcinoma samples and in carcinomas from a distinct head and neck subsite. This result suggests the existence of potential common biomarkers for prognosis and targeted-therapy development in this heterogeneous type of tumor.

Gene Expression Profiles in Radiation Workers Occupationally Exposed to Ionizing Radiation

Ionizing radiation OR) imposes risks to human health and the environment. IR at low doses and low (lose rates has the potency to initiate carcinogenesis. Genotoxic environmental agents such as IR trigger a cascade of signal transduction pathways for cellular protection. In this study, using cDNA microarray technique, we monitored the gene expression profiles in lymphocytes derived from radiation-ex posed individuals (radiation workers). Physical dosimetry records on these patients indicated that the absorbed dose ranged from 0.696 to 39.088 mSv. Gene expression analysis revealed statistically significant transcriptional changes in a total of 78 genes (21 up-regulated and 57 clown-regulated) involved in several biological processes such as ubiquitin cycle (UHRF2 and PIAS1), DNA repair (LIG3, XPA, ERCC5, RAD52, DCLRE1C), cell cycle regulation/proliferation (RHOA, CABLES2, TGFB2, IL16), and stress response (GSTP1, PPP2R5A, DUSP22). Some of the genes that showed altered expression profiles in this study call be used as biomarkers for monitoring the chronic low level exposure in humans. Additionally, alterations in gene expression patterns observed in chronically exposed radiation workers reinforces the need for defining the effective radiation dose that causes immediate genetic damage as well as the long-term effects on genomic instability, including cancer.

A role for mammalian target of rapamycin (mTOR) pathway in non alcoholic steatohepatitis related-cirrhosis

Non-alcoholic fatty liver disease (NAFLD) encompasses the whole spectrum of steatosis, nonalcoholic steatohepatitis (NASH), and NASH-related cirrhosis (NASH/Cir). Although molecular advances have been made in this field, the pathogenesis of NAFLD is not completely understood. The gene expression profiling associated to NASH/Cir was assessed, in an attempt to better characterize the pathways involved in its etiopathogenesis. Methods: In the first step, we used cDNA microarray to evaluate the gene expression profiles in normal liver (n=3) and NASH/Cir samples (n=3) by GeneSifter (TM) analysis to identify differentially expressed genes and biological pathways. Second, tissue microarray was used to determine immunohistochemical expression of phosphorylated mTOR and 4E-BP1 in 11 normal liver samples, 10 NASH/Cir samples and in 37 samples of cirrhosis of other etiologies to further explore the involvement of the mTOR pathway evidenced by the gene expression analysis. Results: 138 and 106 genes were, respectively, up and down regulated in NASH/Cir in comparison to normal liver. Among the 9 pathways identified as significantly modulated in NASH/Cir, the participation of the mTOR pathway was confirmed, since expression of cytoplasmic and membrane phospho-mTOR were higher in NASH/Cir in comparison to cirrhosis of other etiologies and to normal liver. Conclusions: Recent findings have suggested a role for the cellular ""nutrient sensor"" mTOR in NAFLD and the present study corroborates the participation of this pathway in NASH/Cir. Phospho-mTOR evaluation might be of clinical utility as a potential marker for identification of NASH/Cir in cases mistakenly considered as cryptogenic cirrhosis owing to paucity of clinical data.

No-match ORESTES explored as tumor markers

Sequencing technologies and new bioinformatics tools have led to the complete sequencing of various genomes. However, information regarding the human transcriptome and its annotation is yet to be completed. The Human Cancer Genome Project, using ORESTES (open reading frame EST sequences) methodology, contributed to this objective by generating data from about 1.2 million expressed sequence tags. Approximately 30 of these sequences did not align to ESTs in the public databases and were considered no-match ORESTES. On the basis that a set of these ESTs could represent new transcripts, we constructed a cDNA microarray. This platform was used to hybridize against 12 different normal or tumor tissues. We identified 3421 transcribed regions not associated with annotated transcripts, representing 83.3 of the platform. The total number of differentially expressed sequences was 1007. Also, 28 of analyzed sequences could represent noncoding RNAs. Our data reinforces the knowledge of the human genome being pervasively transcribed, and point out molecular marker candidates for different cancers. To reinforce our data, we confirmed, by real-time PCR, the differential expression of three out of eight potentially tumor markers in prostate tissues. Lists of 1007 differentially expressed sequences, and the 291 potentially noncoding tumor markers were provided.

Differential gene expression between the biotrophic-like and saprotrophic mycelia of the witches` broom pathogen Moniliophthora perniciosa

Moniliophthora perniciosa is a hemibiotrophic fungus that causes witches` broom disease (WBD) in cacao. Marked dimorphism characterizes this fungus, showing a monokaryotic or biotrophic phase that causes disease symptoms and a later dikaryotic or saprotrophic phase. A combined strategy of DNA microarray, expressed sequence tag, and real-time reverse-transcriptase polymerase chain reaction analyses was employed to analyze differences between these two fungal stages in vitro. In all, 1,131 putative genes were hybridized with cDNA from different phases, resulting in 189 differentially expressed genes, and 4,595 reads were clusterized, producing 1,534 unigenes. The analysis of these genes, which represent approximately 21% of the total genes, indicates that the biotrophic-like phase undergoes carbon and nitrogen catabollite repression that correlates to the expression of phytopathogenicity genes. Moreover, downregulation of mitochondrial oxidative phosphorylation and the presence of a putative ngr1 of Saccharomyces cerevisiae could help explain its lower growth rate. In contrast, the saprotrophic mycelium expresses genes related to the metabolism of hexoses, ammonia, and oxidative phosphorylation, which could explain its faster growth. Antifungal toxins were upregulated and could prevent the colonization by competing fungi. This work significantly contributes to our understanding of the molecular mechanisms of WBD and, to our knowledge, is the first to analyze differential gene expression of the different phases of a hemibiotrophic fungus.

Molecular profiling of isolated histological components of Wilms tumor implicates a common role for the Wnt signaling pathway in kidney and tumor development

Wilms tumor (WT), a tumor composed of three histological components - blastema (BL), epithelia and stroma - is considered an appropriate model system to study the biological relationship between differentiation and tumorigenesis. To investigate molecular associations between nephrogenesis and WT, the gene expression pattern of individual cellular components was analyzed, using a customized platform containing 4,608 genes. WT gene expression patterns were compared to genes regulated during kidney differentiation. BL had a closer gene expression pattern to the earliest stage of normal renal development. The BL gene expression pattern was compared to that of fetal kidney (FK) and also between FK and mature kidney, identifying 25 common de-regulated genes supposedly involved in the earliest events of WT onset. Quantitative RT-PCR was performed, confirming the difference in expression levels for 13 of 16 genes (81.2%) in the initial set and 8 of 13 (61.5%) in an independent set of samples. An overrepresentation of genes belonging to the Wnt signaling pathway was identified, namely PLCG2, ROCK2 and adenomatous polyposis coli (APC). Activation of the Wnt pathway was confirmed in WT, using APC at protein level and PLCG2 at mRNA and protein level. APC showed positive nuclear immunostaining for an independent set of WT samples, similarly to the FK in week 11. Lack of PLCG2 expression was confirmed in WT and in FK until week 18. Taken together, these results provided molecular evidence of the recapitulation of the embryonic kidney by WT as well as involvement of the Wnt pathway in the earliest events of WT onset. Copyright (C) 2008 S. Karger AG, Basel.

Influence of the interaction between nodal fibroblast and breast cancer cells on gene expression

Our aim was to evaluate the interaction between breast cancer cells and nodal fibroblasts, by means of their gene expression profile. Fibroblast primary cultures were established from negative and positive lymph nodes from breast cancer patients and a similar gene expression pattern was identified, following cell culture. Fibroblasts and breast cancer cells (MDA-MB231, MDA-MB435, and MCF7) were cultured alone or co-cultured separated by a porous membrane (which allows passage of soluble factors) for comparison. Each breast cancer lineage exerted a particular effect on fibroblasts viability and transcriptional profile. However, fibroblasts from positive and negative nodes had a parallel transcriptional behavior when co-cultured with a specific breast cancer cell line. The effects of nodal fibroblasts on breast cancer cells were also investigated. MDA MB-231 cells viability and migration were enhanced by the presence of fibroblasts and accordingly, MDA-MB435 and MCF7 cells viability followed a similar pattern. MDA-MB231 gene expression profile, as evaluated by cDNA microarray, was influenced by the fibroblasts presence, and HNMT, COMT, FN3K, and SOD2 were confirmed downregulated in MDA-MB231 co-cultured cells with fibroblasts from both negative and positive nodes, in a new series of RT-PCR assays. In summary, transcriptional changes induced in breast cancer cells by fibroblasts from positive as well as negative nodes are very much alike in a specific lineage. However, fibroblasts effects are distinct in each one of the breast cancer lineages, suggesting that the inter-relationships between stromal and malignant cells are dependent on the intrinsic subtype of the tumor.

Gene expression profile of residual breast cancer after doxorubicin and cyclophosphamide neoadjuvant chemotherapy

In breast cancer patients, primary chemotherapy is associated with the same survival benefits as adjuvant chemotherapy. Residual tumors represent a clinical challenge, Lis they may be resistant to additional cycles of the same drugs. Our aim was to identify differential transcripts expressed in residual tumors, after neoadjuvant chemotherapy, that might be related with tumor resistance. Hence, 16 patients with paired tumor samples, collected before and after treatment (4 cycles doxorubicin/cyclophosphamide, AC) had their gene expression evaluated on cDNA microarray slides containing 4,608 genes. Three hundred and eighty-nine genes were differentially expressed (paired Student`s t-test, pFDR<0.01) between pre- and post-chemotherapy samples and among the regulated functions were the JNK cascade and cell death. Unsupervised hierarchical clustering identified one branch comprising exclusively, eight pre-chemotherapy samples and another branch, including the former correspondent eight post-chemotherapy samples and other 16 paired pre/post-chemotherapy samples. No differences in clinical and tumor parameters could explain this clustering. Another group of I I patients with paired samples had expression of selected genes determined by real-time RT-PCR and CTGF and DUSP1 were confirmed more expressed in post- as compared to pre-chemotherapy samples. After neoadjuvant chemotherapy some residual samples may retain their molecular signature while others present significant changes in their gene expression, probably induced by the treatment. CTGF and DUSP1 overexpression in residual samples may be a reflection of resistance to further administration of AC regimen.

Maternal embryonic leucine zipper kinase transcript abundance correlates with malignancy grade in human astrocytomas

We have performed cDNA microarray analyses to identify gene expression differences between highly invasive glioblastoma multiforme (GBM) and typically benign pilocytic astrocytomas (PA). Despite the significant clinical and pathological differences between the 2 tumor types, only 63 genes were found to exhibit 2-fold or greater overexpression in GBM as compared to PA. Forty percent of these genes are related to the regulation of the cell cycle and mitosis. QT-PCR validation of 6 overexpressed genes: MELK, AUKB, ASPM, PRC1, IL13RA2 and KIAA0101 confirmed at least a 5-fold increase in the average expression levels in GBM. Maternal embryonic leucine zipper kinase (MELK) exhibited the most statistically significant difference. A more detailed investigation of MELK expression was undertaken to study its oncogenic relevance. In the examination of more than 100 tumors of the central nervous system, we found progressively higher expression of MELK with astrocytoma grade and a noteworthy uniformity of high level expression in GBM. Similar level of overexpression was also observed in medulloblastoma. We found neither gene promoter hypomethylation nor amplification to be a factor in MELK expression, but were able to demonstrate that MELK knockdown in malignant astrocytoma cell lines caused a reduction in proliferation and anchorage-independent growth in in vitro assays. Our results indicate that GBM and PA differ by the expression of surprisingly few genes. Among them, MELK correlated with malignancy grade in astrocytomas and represents a therapeutic target for the management of the most frequent brain tumors in adult and children. (C) 2007 Wiley-Liss, Inc.

Transcriptional Response of Peripheral Lymphocytes to Early Fibrosarcoma: A Model System for Cancer Detection Based on Hybridization Signatures

Since circulating leukocytes, mainly B and T cells, continuously maintain vigilant and comprehensive immune surveillance, these cells could be used as reporters for signs of infection or other pathologies, including cancer. Activated lymphocyte clones trigger a sensitive transcriptional response, which could be identified by gene expression profiling. To assess this hypothesis, we conducted microarray analysis of the gene expression profile of lymphocytes isolated from immunocompetent BALB/c mice subcutaneously injected with different numbers of tumorigenic B61 fibrosarcoma cells. Flow cytometry demonstrated that the number of circulating T (CD3(+)CD4(+) or CD3(+)CD8(+)) or B (CD19(+)) cells did not change. However, the lymphocytes isolated from tumor cell-injected animals expressed a unique transcriptional profile that was identifiable before the development of a palpable tumor mass. This finding demonstrates that the transcriptional response appears before alterations in the main lymphocyte subsets and that the gene expression profile of peripheral lymphocytes can serve as a sensitive and accurate method for the early detection of cancer. Exp Biol Med 234:802-812, 2009

Shared and Unique Gene Expression in Systemic Lupus Erythematosus Depending on Disease Activity

Patients presenting with active Systemic lupus erythematosus (SLE) manifestations may exhibit distinct pathogenetic features in relation to inactive SLE. Also, cDNA microarrays may potentially discriminate the gene expression profile of a disease or disease variant. Therefore, we evaluated the expression profile of 4500 genes in peripheral blood lymphocytes (PBL) of SLE patients. We studied 11 patients with SLE (seven with active SLE and four with inactive SLE) and eight healthy controls. Total RNA was isolated from PBL, reverse transcribed into cDNA, and postlabeled with Cy3 fluorochrome. These probes were then hybridized to a glass slide cDNA microarray containing 4500 human IMAGE cDNA target sequences. An equimolar amount of total RNA from human cell lines served as reference. The microarray images were quantified, normalized, and analyzed using the R environment (ANOVA, significant analysis of microarrays, and cluster-tree view algorithms). Disease activity was assessed by the SLE disease activity index. Compared to the healthy controls, 104 genes in active SLE patients (80 repressed and 24 induced) and 52 genes in nonactive SLE patients (31 induced and 21 repressed) were differentially expressed. The modulation of 12 genes, either induced or repressed, was found in both disease variants; however, each disease variant had differential expression of different genes. Taken together, these results indicate that the two lupus variants studied have common and unique differentially expressed genes. Although the biological significance of the differentially expressed genes discussed above has not been completely understood, they may serve as a platform to further explore the molecular basis of immune deregulation in SLE.

Affected and non-affected skin fibroblasts from systemic sclerosis patients share a gene expression profile deviated from the one observed in healthy individuals

Objectives To evaluate the gene expression profile of fibroblasts from affected and non-affected skin of systemic sclerosis (SSc) patients and from controls. Materials and methods Labeled cDNA from fibroblast cultures from forearm (affected) and axillary (non-affected) skin from six diffuse SSc patients, from three normal controls, and from MOLT-4/HEp-2/normal fibroblasts (reference pool) was probed in microarrays generated with 4193 human cDNAs from the IMAGE Consortium. Microarray images were converted into numerical data and gene expression was calculated as the ratio between fibroblast cDNA (Cy5) and reference pool cDNA (Cy3) data and analyzed by R environment/Aroma, Cluster, Tree View, and SAM softwares. Differential expression was confirmed by real time PCR for a set of selected genes. Results Eighty-eight genes were up- and 241 genes down-regulated in SSc fibroblasts. Gene expression correlation was strong between affected and non-affected fibroblast samples from the same patient (r>0.8), moderate among fibroblasts from all patients (r=0.72) and among fibroblasts from all controls (r=0.70), and modest among fibroblasts from patients and controls (r=0.55). The differential expression was confirmed by real time PCR for all selected genes. Conclusions Fibroblasts from affected and non-affected skin of SSc patients shared a similar abnormal gene expression profile, suggesting that the widespread molecular disturbance in SSc fibroblasts is more sensitive than histological and clinical alterations. Novel molecular elements potentially involved in SSc pathogenesis were identified.