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

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

SAGE analysis highlights the putative role of underexpression of ribosomal proteins in GH-secreting pituitary adenomas

Background: Although the molecular pathogenesis of pituitary adenomas has been assessed by several different techniques, it still remains partially unclear. Ribosomal proteins (RPs) have been recently related to human tumorigenesis, but they have not yet been evaluated in pituitary tumorigenesis. Objective: The aim of this study was to introduce serial analysis of gene expression (SAGE), a high-throughput method, in pituitary research in order to compare differential gene expression. Methods: Two SAGE cDNA libraries were constructed, one using a pool of mRNA obtained from five GH-secreting pituitary tumors and another from three normal pituitaries. Genes differentially expressed between the libraries were further validated by real-time PCR in 22 GH-secreting pituitary tumors and in 15 normal pituitaries. Results: Computer-generated genomic analysis tools identified 13 722 and 14 993 exclusive genes in normal and adenoma libraries respectively. Both shared 6497 genes, 2188 were underexpressed and 4309 overexpressed in tumoral library. In adenoma library, 33 genes encoding RPs were underexpressed. Among these, RPSA, RPS3, RPS14, and RPS29 were validated by real-time PCR. Conclusion: We report the first SAGE library from normal pituitary tissue and GH-secreting pituitary tumor, which provide quantitative assessment of cellular transcriptome. We also validated some downregulated genes encoding RPs. Altogether, the present data suggest that the underexpression of the studied RP genes possibly collaborates directly or indirectly with other genes to modify cell cycle arrest, DNA repair, and apoptosis, leading to an environment that might have a putative role in the tumorigenesis, introducing new perspectives for further studies on molecular genesis of somatotrophinomas.

Serial microanalysis of renal transcriptomes

Large-scale gene expression studies can now be routinely performed on macroamounts of cells, but it is unclear to which extent current methods are valuable for analyzing complex tissues. In the present study, we used the method of serial analysis of gene expression (SAGE) for quantitative mRNA profiling in the mouse kidney. We first performed SAGE at the whole-kidney level by sequencing 12,000 mRNA tags. Most abundant tags corresponded to transcripts widely distributed or enriched in the predominant kidney epithelial cells (proximal tubular cells), whereas transcripts specific for minor cell types were barely evidenced. To better explore such cells, we set up a SAGE adaptation for downsized extracts, enabling a 1,000-fold reduction of the amount of starting material. The potential of this approach was evaluated by studying gene expression in microdissected kidney tubules (50,000 cells). Specific gene expression profiles were obtained, and known markers (e.g., uromodulin in the thick ascending limb of Henle's loop and aquaporin-2 in the collecting duct) were found appropriately enriched. In addition, several enriched tags had no databank match, suggesting that they correspond to unknown or poorly characterized transcripts with specific tissue distribution. It is concluded that SAGE adaptation for downsized extracts makes possible large-scale quantitative gene expression measurements in small biological samples and will help to study the tissue expression and function of genes not evidenced with other high-throughput methods.

Heterogeneity in polyadenylation cleavage sites in mammalian mRNA sequences: implications for SAGE analysis

The analysis of a human thyroid serial analysis of gene expression (SAGE) library shows the presence of an abundant SAGE tag corresponding to the mRNA of thyroglobulin (TG). Additional, less abundant tags are present that can not be linked to any other known gene, but show considerable homology to the wild-type TG tag. To determine whether these tags represent TG mRNA molecules with alternative cleavage, 3′-RACE clones were sequenced. The results show that the three putative TG SAGE tags can be attributed to TG transcripts and reflect the use of alternative polyadenylation cleavage sites downstream of a single polyadenylation signal in vivo. By screening more than 300 000 sequences corresponding to human, mouse and rat transcripts for this phenomenon we show that a considerable percentage of mRNA transcripts (44% human, 22% mouse and 22% rat) show cleavage site heterogeneity. When analyzing SAGE-generated expression data, this phenomenon should be considered, since, according to our calculations, 2.8% of human transcripts show two or more different SAGE tags corresponding to a single gene because of alternative cleavage site selection. Both experimental and in silico data show that the selection of the specific cleavage site for poly(A) addition using a given polyadenylation signal is more variable than was previously thought.

Poly (A)(+) Transcriptome Assessment of ERBB2-Induced Alterations in Breast Cell Lines

We report the first quantitative and qualitative analysis of the poly (A)(+) transcriptome of two human mammary cell lines, differentially expressing (human epidermal growth factor receptor) an oncogene over-expressed in approximately 25% of human breast tumors. Full-length cDNA populations from the two cell lines were digested enzymatically, individually tagged according to a customized method for library construction, and simultaneously sequenced by the use of the Titanium 454-Roche-platform. Comprehensive bioinformatics analysis followed by experimental validation confirmed novel genes, splicing variants, single nucleotide polymorphisms, and gene fusions indicated by RNA-seq data from both samples. Moreover, comparative analysis showed enrichment in alternative events, especially in the exon usage category, in ERBB2 over-expressing cells, data indicating regulation of alternative splicing mediated by the oncogene. Alterations in expression levels of genes, such as LOX, ATP5L, GALNT3, and MME revealed by large-scale sequencing were confirmed between cell lines as well as in tumor specimens with different ERBB2 backgrounds. This approach was shown to be suitable for structural, quantitative, and qualitative assessment of complex transcriptomes and revealed new events mediated by ERBB2 overexpression, in addition to potential molecular targets for breast cancer that are driven by this oncogene.

A score system for quality evaluation of RNA sequence tags: an improvement for gene expression profiling

Background: High-throughput molecular approaches for gene expression profiling, such as Serial Analysis of Gene Expression (SAGE), Massively Parallel Signature Sequencing (MPSS) or Sequencing-by-Synthesis (SBS) represent powerful techniques that provide global transcription profiles of different cell types through sequencing of short fragments of transcripts, denominated sequence tags. These techniques have improved our understanding about the relationships between these expression profiles and cellular phenotypes. Despite this, more reliable datasets are still necessary. In this work, we present a web-based tool named S3T: Score System for Sequence Tags, to index sequenced tags in accordance with their reliability. This is made through a series of evaluations based on a defined rule set. S3T allows the identification/selection of tags, considered more reliable for further gene expression analysis. Results: This methodology was applied to a public SAGE dataset. In order to compare data before and after filtering, a hierarchical clustering analysis was performed in samples from the same type of tissue, in distinct biological conditions, using these two datasets. Our results provide evidences suggesting that it is possible to find more congruous clusters after using S3T scoring system. Conclusion: These results substantiate the proposed application to generate more reliable data. This is a significant contribution for determination of global gene expression profiles. The library analysis with S3T is freely available at http://gdm.fmrp.usp.br/s3t/.S3T source code and datasets can also be downloaded from the aforementioned website.

Global gene expression profile in myelodysplastic syndromes using SAGE

The molecular pathogenesis of myelodysplastic syndromes (MDS) is poorly understood. In order to expand our knowledge of genetic defects in MDS, we determined the overall profile of genes expressed in bone marrow from patients with refractory anemia with excess blasts ( RAEB) by serial analysis of gene expression ( SAGE). The present report describes a partial transcriptome of RAEB bone marrow derived from 56,694 sequenced tags that provides information about expressed gene products. This is the first attempt to determine an overall profile of gene expression specifically in RAEB at diagnosis using SAGE, which should be useful in the understanding of the physiopathology of MDS and in identifying the genes involved.

SAGE transcript profiling of the juvenile cambial region of Eucalyptus grandis

Despite the importance of Eucalyptus spp. in the pulp and paper industry, functional genomic approaches have only recently been applied to understand wood formation in this genus. We attempted to establish a global view of gene expression in the juvenile cambial region of Eucalyptus grandis Hill ex Maiden. The expression profile was obtained from serial analysis of gene expression (SAGE) library data produced from 3- and 6-year-old trees. Fourteen-base expressed sequence tags (ESTs) were searched against public Eucalyptus ESTs and annotated with GenBank. Altogether 43,304 tags were generated producing 3066 unigenes with three or more copies each, 445 with a putative identity, 215 with unknown function and 2406 without an EST match. The expression profile of the juvenile cambial region revealed the presence of highly frequent transcripts related to general metabolism and energy metabolism, cellular processes, transport, structural components and information pathways. We made a quantitative analysis of a large number of genes involved in the biosynthesis of cellulose, pectin, hemicellulose and lignin. Our findings provide insight into the expression of functionally related genes involved in juvenile wood formation in young fast-growing E. grandis trees.

Gene expression profiling of papillary thyroid carcinoma identifies transcripts correlated with BRAF mutational status and lymph node metastasis

Purpose: To identify papillary thyroid carcinoma (PTC)-associated transcripts, we compared the gene expression profiles of three Serial Analysis of Gene Expression libraries generated from thyroid tumors and a normal thyroid tissue. Experimental Design: Selected transcripts were validated in a panel of 57 thyroid tumors using quantitative PCR (qPCR). An independent set of 71 paraffin-embedded sections was used for validation using immunohistochemical analysis. To determine if PTC-associated gene expression could predict lymph node involvement, a separate cohort of 130 primary PTC (54 metastatic and 76 nonmetastatic) was investigated. The BRAF(V600E) mutational status was compared with qPCR data to identify genes that might be regulated by abnormal BRAF/MEK/extracellular signal-regulated kinase signaling. Results: We identified and validated new PTC-associated transcripts. Three genes (CST6, CXCL14, and DHRS3) are strongly associated with PTC. Immunohistochemical analysis of CXCL14 confirmed the qPCR data and showed protein expression in PTC epithelial cells. We also observed that CST6, CXCL14, DHRS3, and SPP1 were associated with PTC lymph node metastasis, with CST6, CXCL14, and SPP1 being positively correlated with metastasis and DHRS3 being negatively correlated. Finally, we found a strong correlation between CST6 and CXCL14 expression and BRAF(V600E) mutational status, suggesting that these genes may be induced subsequently to BRAF activation and therefore may be downstream in the BRAF/MEK/extracellular signal-regulated kinase signaling pathway. Conclusion: CST6, CXCL14, DHRS3, and SPP1 may play a role in PTC pathogenesis and progression and are possible molecular targets for FTC therapy.

Transcriptome changes induced by docetaxel in human mammary cell lines expressing different levels of ERBB2

The taxane docetaxel is currently the most effective chemotherapeutic drug for the treatment of advanced breast cancer. However, a considerable proportion of breast cancer patients do not respond positively to docetaxel. The mechanisms of docetaxel resistance are poorly understood. Overexpression of ERBB2 occurs in 15-30% of breast tumors and is associated with chemoresistance to a variety of anticancer drugs. In the present study, we sought to identify genes involved in ERBB2-mediated chemoresistance to docetaxel. We generated SAGE libraries from two human mammary cell lines expressing basal (HB4a) and high (C5.2) levels of ERBB2 before and after intensive exposure to docetaxel and identified potential ERBB2 target genes implicated in a variety of cellular processes including cell proliferation, cell adhesion, apoptosis and cytoskeleton organization. Comparison of the transcriptome of the cell lines before and after docetaxel exposure revealed substantially different expression patterns. Twenty-one differentially expressed genes between HB4a and C5.2 cell lines, before and after docetaxel treatment, were further analyzed by qPCR. The alterations in the expression patterns in HB4a and C5.2 cell lines in response to docetaxel treatment observed by SAGE analysis were confirmed by qPCR for the majority of the genes analyzed. Our study provides a comprehensive view of the expression changes induced in two human mammary cells expressing different levels of ERBB2 in response to docetaxel that could contribute to the elucidation of the mechanisms involved in ERBB2-mediated chemoresistance in breast cancer.

Annexin A1 subcellular expression in laryngeal squamous cell carcinoma

Annexin A1 (ANXA1) is a soluble cytoplasmic protein, moving to membranes when calcium levels are elevated. ANXA1 has also been shown to move to the nucleus or outside the cells, depending on tyrosine-kinase signalling, thus interfering in cytoskeletal organization and cell differentiation, mostly in inflammatory and neoplastic processes. The aim was to investigate subcellular patterns of immunohistochemical expression of ANXA1 in neoplastic and non-neoplastic samples from patients with laryngeal squamous cell carcinomas (LSCC), to elucidate the role of ANXA1 in laryngeal carcinogenesis. Serial analysis of gene expression experiments detected reduced expression of ANXA1 gene in LSCC compared with the corresponding non-neoplastic margins. Quantitative polymerase chain reaction confirmed ANXA1 low expression in 15 LSCC and eight matched normal samples. Thus, we investigated subcellular patterns of immunohistochemical expression of ANXA1 in 241 paraffin-embedded samples from 95 patients with LSCC. The results showed ANXA1 down-regulation in dysplastic, tumourous and metastatic lesions and provided evidence for the progressive migration of ANXA1 from the nucleus towards the membrane during laryngeal tumorigenesis. ANXA1 dysregulation was observed early in laryngeal carcinogenesis, in intra-epithelial neoplasms; it was not found related to prognostic parameters, such as nodal metastases.

Multipotent mesenchymal stromal cells obtained from diverse human tissues share functional properties and gene-expression profile with CD146(+) perivascular cells and fibroblasts

Objective. The relationship of multipotent mesenchymal stromal cells (MSC) with pericytes and fibroblasts has not been established thus far, although they share many markers of primitive marrow stromal cells and the osteogenic, adipogenic, and chondrogenic differentiation potentials. Materials and Methods. We compared MSCs from adult or fetal tissues, MSC differentiated in vitro, fibroblasts and cultures of retinal pericytes obtained either by separation with anti-CD146 or adhesion. The characterizations included morphological, immunophenotypic, gene-expression profile, and differentiation potential. Results. Osteogenic, adipocytic, and chondrocytic differentiation was demonstrated for MSC, retinal perivascular cells, and fibroblasts. Cell morphology and the phenotypes defined by 22 markers were very similar. Analysis of the global gene expression obtained by serial analysis of gene expression for 17 libraries and by reverse transcription polymerase chain reaction of 39 selected genes from 31 different cell cultures, revealed similarities among MSC, retinal perivascular cells, and hepatic stellate cells. Despite this overall similarity, there was a heterogeneous expression of genes related to angiogenesis, in MSC derived from veins, artery, perivascular cells, and fibroblasts. Evaluation of typical pericyte and MSC transcripts, such as NG2, CD146, CD271, and CD140B on CD146 selected perivascular cells and MSC by real-time polymerase chain reaction confirm the relationship between these two cell types. Furthermore, the inverse correlation between fibroblast-specific protein-1 and CD146 transcripts observed on pericytes, MSC, and fibroblasts highlight their potential use as markers of this differentiation pathway. Conclusion. Our results indicate that human MSC and pericytes are similar cells located in the wall of the vasculature, where they function as cell sources for repair and tissue maintenance, whereas fibroblasts are more differentiated cells with more restricted differentiation potential. (C) 2008 ISEH - Society for Hematology and Stem Cells. Published by Elsevier Inc.

Role of NFKB2 on the early myeloid differentiation of CD34+ hematopoietic stem/progenitor cells

To better understand the early events regulating lineage-specific hematopoietic differentiation, we analyzed the transcriptional profiles of CD34+ human hematopoietic stem and progenitor cells (HSPCs) subjected to differentiation stimulus. CD34+ cells were cultured for 12 and 40 h in liquid cultures with supplemented media favoring myeloid or erythroid commitment. Serial analysis of gene expression (SAGE) was employed to generate four independent libraries. By analyzing the differentially expressed regulated transcripts between the un-stimulated and the stimulated CD34+ cells, we observed a set of genes that was initially up-regulated at 12 h but were then down-regulated at 40 h, exclusively after myeloid stimulus. Among those we found transcripts for NFKB2, RELB, IL1B, LTB, LTBR, TNFRSF4, TGFB1, and IKBKA. Also, the inhibitor NFKBIA (IKBA) was more expressed at 12 h. All those transcripts code for signaling proteins of the nuclear factor kappaB pathway. NFKB2 is a subunit of the NF-kappa B transcription factor that with RELB mediates the non-canonical NF-kappa B pathway. Interference RNA (RNAi) against NFKB1, NFKB2 and control RNAi were transfected into bone marrow CD34+HSPC. The percentage and the size of the myeloid colonies derived from the CD34+ cells decreased after inhibition of NFKB2. Altogether, our results indicate that NFKB2 gene has a role in the early commitment of CD34+HSPC towards the myeloid lineage. (C) 2010 International Society of Differentiation. Published by Elsevier Ltd. All rights reserved.

A destruição da muralha e do castelo de Braga no inicio do século XX: elementos para a compreensão do castelo medieval

Relatório de estágio de mestrado em História

Evaluation of the chicken transcriptome by SAGE of B cells and the DT40 cell line

Background: The understanding of whole genome sequences in higher eukaryotes depends to a large degree on the reliable definition of transcription units including exon/intron structures, translated open reading frames (ORFs) and flanking untranslated regions. The best currently available chicken transcript catalog is the Ensembl build based on the mappings of a relatively small number of full length cDNAs and ESTs to the genome as well as genome sequence derived in silico gene predictions.Results: We use Long Serial Analysis of Gene Expression (LongSAGE) in bursal lymphocytes and the DT40 cell line to verify the quality and completeness of the annotated transcripts. 53.6% of the more than 38,000 unique SAGE tags (unitags) match to full length bursal cDNAs, the Ensembl transcript build or the genome sequence. The majority of all matching unitags show single matches to the genome, but no matches to the genome derived Ensembl transcript build. Nevertheless, most of these tags map close to the 3' boundaries of annotated Ensembl transcripts.Conclusions: These results suggests that rather few genes are missing in the current Ensembl chicken transcript build, but that the 3' ends of many transcripts may not have been accurately predicted. The tags with no match in the transcript sequences can now be used to improve gene predictions, pinpoint the genomic location of entirely missed transcripts and optimize the accuracy of gene finder software.