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.

Hedgehog signaling and osteoblast gene expression are regulated by purmorphamine in human mesenchymal stem cells

Several biological events are controlled by Hedgehog (Hh) signaling, including osteoblast phenotype development. This study aimed at evaluating the gene expression profile of human mesenchymal stem cells (hMSCs) treated with the Hh agonist, purmorphamine, focusing on Hh signaling and osteoblast differentiation. hMSCs from bone marrow were cultured in non-osteogenic medium with or without purmorphamine (2 mu M) for periods of up to 14 days. Purmorphamine up-regulated gene expression of the mediators of Hh pathway, SMO, PTCH1, GLI1, and GLI2. The activation of Hh pathway by purmorphamine increased the expression of several genes (e.g., RUNX2 and BMPs) related to osteogenesis. Our results indicated that purmorphamine triggers Hh signaling pathway in hMSCs, inducing an increase in the expression of a set of genes involved in the osteoblast differentiation program. Thus, we conclude that Hh is a crucial pathway in the commitment of undifferentiated cells to the osteoblast lineage. J. Cell. Biochem. 113: 204208, 2012. (C) 2011 Wiley Periodicals, Inc.

Four-gene expression model predictive of lymph node metastases in oral squamous cell carcinoma

Background. Previous knowledge of cervical lymph node compromise may be crucial to choose the best treatment strategy in oral squamous cell carcinoma (OSCC). Here we propose a set four genes, whose mRNA expression in the primary tumor predicts nodal status in OSCC, excluding tongue. Material and methods. We identified differentially expressed genes in OSCC with and without compromised lymph nodes using Differential Display RT-PCR. Known genes were chosen to be validated by means of Northern blotting or real time RT-PCR (qRT-PCR). Thereafter we constructed a Nodal Index (NI) using discriminant analysis in a learning set of 35 patients, which was further validated in a second independent group of 20 patients. Results. Of the 63 differentially expressed known genes identified comparing three lymph node positive (pN+) and three negative (pN0) primary tumors, 23 were analyzed by Northern analysis or RT-PCR in 49 primary tumors. Six genes confirmed as differentially expressed were used to construct a NI, as the best set predictive of lymph nodal status, with the final result including four genes. The NI was able to correctly classify 32 of 35 patients comprising the learning group (88.6%; p = 0.009). Casein kinase 1alpha1 and scavenger receptor class B, member 2 were found to be up regulated in pN + group in contrast to small proline-rich protein 2B and Ras-GTPase activating protein SH3 domain-binding protein 2 which were upregulated in the pN0 group. We validated further our NI in an independent set of 20 primary tumors, 11 of them pN0 and nine pN+ with an accuracy of 80.0% (p = 0.012). Conclusions. The NI was an independent predictor of compromised lymph nodes, taking into the consideration tumor size and histological grade. The genes identified here that integrate our "Nodal Index" model are predictive of lymph node metastasis in OSCC.

Analysis of MLL2 gene in the first Brazilian family with Kabuki syndrome

Most patients with Kabuki syndrome (KS) are the only person in their family with the condition. However, familial cases of KS have been described showing evidence that this syndrome can be inherited as a dominant trait with variable expressivity. We report on two related individuals with facial findings characteristic of KS. The proposita had arched eyebrows, long and upward slanting palpebral fissures, cleft lip and palate, retromicrognathia, brachydactyly of hands and feet, stubby fingers, nail hypoplasia, and prominent finger pads. Her mother had eyebrows with dispersed lateral half, long and upward slanting palpebral fissures, retrognathia, abnormal and posteriorly rotated ears, prominent finger pads, brachydactyly of feet, learning difficulties, and psychomotor development delay. DNA sequencing revealed a novel missense mutation in the MLL2 gene in both the proposita and her mother. The mutation (p.R5432Q) was found in the exon 51, within the SET domain of the gene, which confers methyltransferase activity on the protein. Therefore, the epigenetic and transcriptional regulatory properties of this protein may be altered and this suggests that the mutation is the cause of phenotype observed in both the patient and her mother. The clinical signs and the molecular evidence in this family further support the notion that KS is an autosomal dominant condition with variable expressivity. To our knowledge this is the first report of a Brazilian family with recurrence of this syndrome. (C) 2012 Wiley Periodicals, Inc.

Hox Gene Expression Leads to Differential Hind Leg Development between Honeybee Castes

Beyond the physiological and behavioural, differences in appendage morphology between the workers and queens of Apis mellifera are pre-eminent. The hind legs of workers, which are highly specialized pollinators, deserve special attention. The hind tibia of worker has an expanded bristle-free region used for carrying pollen and propolis, the corbicula. In queens this structure is absent. Although the morphological differences are well characterized, the genetic inputs driving the development of this alternative morphology remain unknown. Leg phenotype determination takes place between the fourth and fifth larval instar and herein we show that the morphogenesis is completed at brown-eyed pupa. Using results from the hybridization of whole genome-based oligonucleotide arrays with RNA samples from hind leg imaginal discs of pre-pupal honeybees of both castes we present a list of 200 differentially expressed genes. Notably, there are castes preferentially expressed cuticular protein genes and members of the P450 family. We also provide results of qPCR analyses determining the developmental transcription profiles of eight selected genes, including abdominal-A, distal-less and ultrabithorax (Ubx), whose roles in leg development have been previously demonstrated in other insect models. Ubx expression in workers hind leg is approximately 25 times higher than in queens. Finally, immunohistochemistry assays show that Ubx localization during hind leg development resembles the bristles localization in the tibia/basitarsus of the adult legs in both castes. Our data strongly indicate that the development of the hind legs diphenism characteristic of this corbiculate species is driven by a set of caste-preferentially expressed genes, such as those encoding cuticular protein genes, P450 and Hox proteins, in response to the naturally different diets offered to honeybees during the larval period.

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.

Sugarwin: A Sugarcane Insect-Induced Gene with Antipathogenic Activity

In sugarcane fields, colonization of the stalk by opportunistic fungi usually occurs after the caterpillar Diatraea saccharalis attacks the sugarcane plant. Plants respond to insect attack by inducing and accumulating a large set of defense proteins. Two homologues of a barley wound-inducible protein (BARWIN), sugarcane wound-inducible proteins SUGARWIN1 and SUGARWIN2, have been identified in sugarcane by an in silico analysis. Antifungal properties have been described for a number of BARWIN homologues. We report that a SUGARWIN:green fluorescent protein fusion protein is located in the endoplasmic reticulum and in the extracellular space of sugarcane plants. The induction of sugarwin transcripts occurs in response to mechanical wounding, D. saccharalis damage, and methyl jasmonate treatment. The accumulation of transcripts is late induced and is restricted to the site of the wound. Although the transcripts of sugarwin genes were strongly increased following insect attack, the protein itself did not show any effect on insect development; rather, it altered fungal morphology, leading to the apoptosis of the germlings. These results suggest that, in the course of evolution, sugarwin-encoding genes were recruited by sugarcane due to their antipathogenic activity. We rationalize that sugarcane is able to induce sugarwin gene expression in response to D. saccharalis feeding as a concerted plant response to the anticipated invasion by the fungi that typically penetrate the plant stalk after insect damage.

Gene Expression Profile in Response to Doxorubicin-Rapamycin Combined Treatment of HER-2-Overexpressing Human Mammary Epithelial Cell Lines

HER-2-positive breast cancers frequently sustain elevated AKT/mTOR signaling, which has been associated with resistance to doxorubicin treatment. Here, we investigated whether rapamycin, an mTOR inhibitor, increased the sensitivity to doxorubicin therapy in two HER-2-overexpressing cell lines: C5.2, which was derived from the parental HB4a by transfection with HER-2 and SKBR3, which exhibits HER-2 amplification. The epithelial mammary cell line HB4a was also analyzed. The combined treatment using 20 nmol/L of rapamycin and 30 nmol/L of doxorubicin arrested HB4a and C5.2 cells in S to G(2)-M, whereas SKBR3 cells showed an increase in the G(0)-G(1) phase. Rapamycin increased the sensitivity to doxorubicin in HER-2-overexpressing cells by approximately 2-fold, suggesting that the combination displayed a more effective antiproliferative action. Gene expression profiling showed that these results might reflect alterations in genes involved in canonical pathways related to purine metabolism, oxidative phosphorylation, protein ubiquitination, and mitochondrial dysfunction. A set of 122 genes modulated by the combined treatment and specifically related to HER-2 overexpression was determined by finding genes commonly regulated in both C5.2 and SKBR3 that were not affected in HB4a cells. Network analysis of this particular set showed a smaller subgroup of genes in which coexpression pattern in HB4a cells was disrupted in C5.2 and SKBR3. Altogether, our data showed a subset of genes that might be more robust than individual markers in predicting the response of HER-2-overexpressing breast cancers to doxorubicin and rapamycin combination. Mol Cancer Ther; 11(2); 464-74. (C) 2011 AACR.

A new enrichment space for the treatment of discontinuous pressures in multi-fluid flows

In this work, a new enrichment space to accommodate jumps in the pressure field at immersed interfaces in finite element formulations, is proposed. The new enrichment adds two degrees of freedom per element that can be eliminated by means of static condensation. The new space is tested and compared with the classical P1 space and to the space proposed by Ausas et al (Comp. Meth. Appl. Mech. Eng., Vol. 199, 10191031, 2010) in several problems involving jumps in the viscosity and/or the presence of singular forces at interfaces not conforming with the element edges. The combination of this enrichment space with another enrichment that accommodates discontinuities in the pressure gradient has also been explored, exhibiting excellent results in problems involving jumps in the density or the volume forces. Copyright (c) 2011 John Wiley & Sons, Ltd.

SET overexpression decreases cell detoxification efficiency: ALDH2 and GSTP1 are downregulated, DDR is impaired and DNA damage accumulates

Alcohol and tobacco consumption are risk factors for head and neck squamous cell carcinoma (HNSCC). Aldehyde dehydrogenase 2 (ALDH2) and glutathione Stransferase pi 1 (GSTP1) are important enzymes for cellular detoxification and low efficiencies are implicated in cancer. We assessed the potential role of SET protein overexpression, a histone acetylation modulator accumulated in HNSCC, in gene regulation and protein activity of ALDH2 and GSTP1. SET was knocked down in HN13, HN12 and Cal27, and overexpressed in HEK293 cells; ethanol and cisplatin were the chemical agents. Cells with SET overexpression (HEK293/SET, HN13 and HN12) showed lower ALDH2 and GSTP1 mRNA levels and trichostatin A increased them (real-time PCR). Ethanol upregulated GSTP1 and ALDH2 mRNAs, whereas cisplatin upregulated GSTP1 in HEK293 cells. SET-chromatin binding revealed SET interaction with ALDH2 and GSTP1 promoters, specifically via SET NAP domain; ethanol and cisplatin abolished SET binding. ALDH2 and GSTP1 efficiency was assessed by enzymatic and comet assay. A lower ALDH2 activity was associated with greater DNA damage (tail intensity) in HEK293/SET compared with HEK293 cells, whereas HN13/siSET showed ALDH2 activity higher than HN13 cells. HN13/siSET cells showed increased tail intensity. Cisplatin-induced DNA damage response showed negative relationship between SET overexpression and BRCA2 recruitment. SET downregulated repair genes ATM, BRCA1 and CHEK2 and upregulated TP53. Cisplatin-induced cell-cycle arrest occurred in G0/G1 and S in HEK293 cells, whereas HEK293/SET showed G2/M stalling. Overall, cisplatin was more cytotoxic for HN13 than HN13/siSET cells. Our data suggest a role for SET in cellular detoxification, DNA damage response and genome integrity.

BayGO: Bayesian analysis of ontology term enrichment in microarray data

Abstract Background The search for enriched (aka over-represented or enhanced) ontology terms in a list of genes obtained from microarray experiments is becoming a standard procedure for a system-level analysis. This procedure tries to summarize the information focussing on classification designs such as Gene Ontology, KEGG pathways, and so on, instead of focussing on individual genes. Although it is well known in statistics that association and significance are distinct concepts, only the former approach has been used to deal with the ontology term enrichment problem. Results BayGO implements a Bayesian approach to search for enriched terms from microarray data. The R source-code is freely available at http://blasto.iq.usp.br/~tkoide/BayGO in three versions: Linux, which can be easily incorporated into pre-existent pipelines; Windows, to be controlled interactively; and as a web-tool. The software was validated using a bacterial heat shock response dataset, since this stress triggers known system-level responses. Conclusion The Bayesian model accounts for the fact that, eventually, not all the genes from a given category are observable in microarray data due to low intensity signal, quality filters, genes that were not spotted and so on. Moreover, BayGO allows one to measure the statistical association between generic ontology terms and differential expression, instead of working only with the common significance analysis.

Identification of unannotated exons of low abundance transcripts in Drosophila melanogaster and cloning of a new serine protease gene upregulated upon injury

Abstract Background The sequencing of the D.melanogaster genome revealed an unexpected small number of genes (~ 14,000) indicating that mechanisms acting on generation of transcript diversity must have played a major role in the evolution of complex metazoans. Among the most extensively used mechanisms that accounts for this diversity is alternative splicing. It is estimated that over 40% of Drosophila protein-coding genes contain one or more alternative exons. A recent transcription map of the Drosophila embryogenesis indicates that 30% of the transcribed regions are unannotated, and that 1/3 of this is estimated as missed or alternative exons of previously characterized protein-coding genes. Therefore, the identification of the variety of expressed transcripts depends on experimental data for its final validation and is continuously being performed using different approaches. We applied the Open Reading Frame Expressed Sequence Tags (ORESTES) methodology, which is capable of generating cDNA data from the central portion of rare transcripts, in order to investigate the presence of hitherto unnanotated regions of Drosophila transcriptome. Results Bioinformatic analysis of 1,303 Drosophila ORESTES clusters identified 68 sequences derived from unannotated regions in the current Drosophila genome version (4.3). Of these, a set of 38 was analysed by polyA+ northern blot hybridization, validating 17 (50%) new exons of low abundance transcripts. For one of these ESTs, we obtained the cDNA encompassing the complete coding sequence of a new serine protease, named SP212. The SP212 gene is part of a serine protease gene cluster located in the chromosome region 88A12-B1. This cluster includes the predicted genes CG9631, CG9649 and CG31326, which were previously identified as up-regulated after immune challenges in genomic-scale microarray analysis. In agreement with the proposal that this locus is co-regulated in response to microorganisms infection, we show here that SP212 is also up-regulated upon injury. Conclusion Using the ORESTES methodology we identified 17 novel exons from low abundance Drosophila transcripts, and through a PCR approach the complete CDS of one of these transcripts was defined. Our results show that the computational identification and manual inspection are not sufficient to annotate a genome in the absence of experimentally derived data.

Multivariate gene expression analysis reveals functional connectivity changes between normal/tumoral prostates

Abstract Background Prostate cancer is a leading cause of death in the male population, therefore, a comprehensive study about the genes and the molecular networks involved in the tumoral prostate process becomes necessary. In order to understand the biological process behind potential biomarkers, we have analyzed a set of 57 cDNA microarrays containing ~25,000 genes. Results Principal Component Analysis (PCA) combined with the Maximum-entropy Linear Discriminant Analysis (MLDA) were applied in order to identify genes with the most discriminative information between normal and tumoral prostatic tissues. Data analysis was carried out using three different approaches, namely: (i) differences in gene expression levels between normal and tumoral conditions from an univariate point of view; (ii) in a multivariate fashion using MLDA; and (iii) with a dependence network approach. Our results show that malignant transformation in the prostatic tissue is more related to functional connectivity changes in their dependence networks than to differential gene expression. The MYLK, KLK2, KLK3, HAN11, LTF, CSRP1 and TGM4 genes presented significant changes in their functional connectivity between normal and tumoral conditions and were also classified as the top seven most informative genes for the prostate cancer genesis process by our discriminant analysis. Moreover, among the identified genes we found classically known biomarkers and genes which are closely related to tumoral prostate, such as KLK3 and KLK2 and several other potential ones. Conclusion We have demonstrated that changes in functional connectivity may be implicit in the biological process which renders some genes more informative to discriminate between normal and tumoral conditions. Using the proposed method, namely, MLDA, in order to analyze the multivariate characteristic of genes, it was possible to capture the changes in dependence networks which are related to cell transformation.

Constraint-based analysis of gene interactions using restricted boolean networks and time-series data

Abstract Background A popular model for gene regulatory networks is the Boolean network model. In this paper, we propose an algorithm to perform an analysis of gene regulatory interactions using the Boolean network model and time-series data. Actually, the Boolean network is restricted in the sense that only a subset of all possible Boolean functions are considered. We explore some mathematical properties of the restricted Boolean networks in order to avoid the full search approach. The problem is modeled as a Constraint Satisfaction Problem (CSP) and CSP techniques are used to solve it. Results We applied the proposed algorithm in two data sets. First, we used an artificial dataset obtained from a model for the budding yeast cell cycle. The second data set is derived from experiments performed using HeLa cells. The results show that some interactions can be fully or, at least, partially determined under the Boolean model considered. Conclusions The algorithm proposed can be used as a first step for detection of gene/protein interactions. It is able to infer gene relationships from time-series data of gene expression, and this inference process can be aided by a priori knowledge available.

Genome mapping and expression analyses of human intronic noncoding RNAs reveal tissue-specific patterns and enrichment in genes related to regulation of transcription

Abstract Background RNAs transcribed from intronic regions of genes are involved in a number of processes related to post-transcriptional control of gene expression. However, the complement of human genes in which introns are transcribed, and the number of intronic transcriptional units and their tissue expression patterns are not known. Results A survey of mRNA and EST public databases revealed more than 55,000 totally intronic noncoding (TIN) RNAs transcribed from the introns of 74% of all unique RefSeq genes. Guided by this information, we designed an oligoarray platform containing sense and antisense probes for each of 7,135 randomly selected TIN transcripts plus the corresponding protein-coding genes. We identified exonic and intronic tissue-specific expression signatures for human liver, prostate and kidney. The most highly expressed antisense TIN RNAs were transcribed from introns of protein-coding genes significantly enriched (p = 0.002 to 0.022) in the 'Regulation of transcription' Gene Ontology category. RNA polymerase II inhibition resulted in increased expression of a fraction of intronic RNAs in cell cultures, suggesting that other RNA polymerases may be involved in their biosynthesis. Members of a subset of intronic and protein-coding signatures transcribed from the same genomic loci have correlated expression patterns, suggesting that intronic RNAs regulate the abundance or the pattern of exon usage in protein-coding messages. Conclusion We have identified diverse intronic RNA expression patterns, pointing to distinct regulatory roles. This gene-oriented approach, using a combined intron-exon oligoarray, should permit further comparative analysis of intronic transcription under various physiological and pathological conditions, thus advancing current knowledge about the biological functions of these noncoding RNAs.