Abstract 3755: Dlx4os: a lncRNA associated with malignant melanocyte transformation

It is estimated that more than 90% of the human genome is transcribed into non-coding RNAs (ncRNA). It has been demonstrated that ncRNAs have an important role in biological processes such as proliferation, differentiation, and cell migration. Changes in their expression have been associated with several diseases, including cancer. Melanoma is one of the most aggressive types of cancer, with high probability of metastasis and an unfavorable response to therapies. Our lab has developed an in vitro melanoma progression model, which was established by exposing nontumorigenic melanocytes (melan-a) to sustained stress conditions (cycles of adhesion blockage), sequentially generating different cell lines (i.e., 4C, premalignant melanocytes; 4C11-, non-metastatic melanoma cells; and 4C11+, metastatic melanoma cells), in which each cell line represents a distinct stage of tumor progression. Previous RNA sequencing data comparing these 4 cell lines identified a total of 3404 long non-coding transcripts which displayed differential expression (DE) levels when comparing two or more cell lines. We have clustered 531 of them into three distinct signatures: malignancy, epithelial-to-mesenchymal transition (EMT), or metastasis, according to their expression profiles. We also analyzed the neighboring genes of those lncRNAs, as it is known that lncRNAs might regulate gene expression in cis. We identified the lncRNA Dlx4os, highly expressed in undifferentiated and mesenchymal-like 4C and 4C11- cells, and 6 neighboring genes whose expression correlate to the EMT signature and overall survival in human melanoma patients. Among these genes we found Dlx4, which shares the same promoter site with Dlx4os and whose expression correlates to a poor prognosis. We silenced Dlx4os in the 4C11- cells through interference RNA technology. The downregulation of Dlx4os led to expression changes in EMT-related genes, such as Snail1, Sox10 and Mitf, indicating that it might play a role in EMT. Dlx4os downregulation also inhibited cell migration in vitro and tumorigenesis in vivo. These results overall indicate a role of Dlx4os in melanocyte malignant transformation.

Citation Format: Ana Luísa P. Ayub, Hatylas Azevedo, Diogo Pessoa, Leinal Sejour, Ioannis Vlachos, Eduardo M. Reis, Frank J. Slack, Miriam G. Jasiulionis. Dlx4os: a lncRNA associated with malignant melanocyte transformation. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3755.

Necroptosis activation is associated with greater methylene blue-photodynamic therapy-induced cytotoxicity in human pancreatic ductal adenocarcinoma cells

Pancreatic ductal adenocarcinomas (PDAC) are the fourth leading cause of death due to neoplasms. In view of the urgent need of effective treatments for PDAC, photodynamic therapy (PDT) appears as a promising alternative. However, its efficacy against PDAC and the mechanisms involved in cell death induction remain unclear. In this study, we set out to evaluate PDT's cytotoxicity using methylene blue (MB) as a photosensitizer (PS) (MB-PDT) and to evaluate the contribution of necroptosis in its effect in human PDAC cells. Our results demonstrated that MB-PDT induced significant death of different human PDAC models presenting two different susceptibility profiles. This effect was independent of MB uptake or its subcellular localization. We found that the ability of triggering necroptosis was determinant to increase the treatment efficiency. Analysis of single cell RNA-seq data from normal and neoplastic human pancreatic tissues showed that specific necroptosis proteins RIPK1, RIPK3 and MLKL presented significant higher expression levels in cells displaying a transformed phenotype providing further support to the use of approaches that activate necroptosis, like MB-PDT, as useful adjunct to surgery of PDAC to tackle the problem of microscopic residual disease as well as to minimize the chance of local and metastatic recurrence.

Genome-wide promoter methylation profiling in a cellular model of melanoma progression reveals markers of malignancy and metastasis that predict melanoma survival

The epigenetic changes associated with melanoma progression to advanced and metastatic stages are still poorly understood. To shed light on the CpG methylation dynamics during melanoma development, we analyzed the methylome profiles of a four-stage cell line model of melanoma progression: non-tumorigenic melanocytes (melan-a), premalignant melanocytes (4C), non-metastatic melanoma cells (4C11−), and metastatic melanoma cells (4C11+). We identified 540 hypo- and 37 hypermethylated gene promoters that together characterized a malignancy signature, and 646 hypo- and 520 hypermethylated promoters that distinguished a metastasis signature. Differentially methylated genes from these signatures were correlated with overall survival using TCGA-SKCM methylation data. Moreover, multivariate Cox analyses with LASSO regularization identified panels of 33 and 31 CpGs, respectively, from the malignancy and metastasis signatures that predicted poor survival. We found a concordant relationship between DNA methylation and transcriptional levels for genes from the malignancy (Pyroxd2 and Ptgfrn) and metastasis (Arnt2, Igfbp4 and Ptprf) signatures, which were both also correlated with melanoma prognosis. Altogether, this study reveals novel CpGs methylation markers associated with malignancy and metastasis that collectively could improve the survival prediction of melanoma patients.

Abstract 2392: Transcriptome analysis from a linear melanoma progression model reveals lncRNAs involved with malignancy, EMT and metastasis

In the human genome, it is estimated that more than 90% of all transcripts are transcribed into non-coding RNAs (ncRNA), classified as small (such as miRNAs) and long (lncRNA). It has been demonstrated that ncRNAs have an important role in biological processes, such as proliferation, differentiation and cellular migration. Changes in their expression have been associated to several diseases, including cancer. Melanoma is one of the most aggressive types of cancer, with high probability of metastasis and showing an unfavorable response to therapies. Our lab has developed a melanoma progression model, which was established by exposing nontumorigenic melanocytes to sustained stressed conditions (cycles of adhesion blockage), sequentially generating 4 cell lines (melan-a, nontumorigenic melanocytes; 4C, premalignant melanocytes; 4C11-, non-metastatic melanoma cells; 4C11+, metastatic melanoma cells) in which different cell lines represent distinct stages of the progression. Previous data comparing the RNA sequencing results from all 4 cell lines identified a total of 3662 genes, coding and non-coding, which displayed differentiated expression (DE) levels when comparing two or more cell lines. Out of them, 492 were DE lncRNAs. We categorized them into three distinct signatures: malignancy, epithelial mesenchymal transition (EMT), or metastasis, according to their expression profile. We analyzed the neighboring genes from those lncRNAs, as it is well described that lncRNA might regulate the gene expression of its neighbors in cis. We identified the lncRNA Dlx4os, with 6 neighbor genes correlated to EMT signature and overall survival of melanoma patients using two different melanoma cohorts (TCGA and Leeds). Among them, the Dlx4 gene, which expression correlates to poor prognosis, shares the same promoter site on opposite strands in chromosome 11 with Dlx4os. We also identified HOTAIR, a well described lncRNA, altered in the EMT signature. We silenced HOTAIR, Dlx4os and Dlx4 on 4C and 4C11- cell lines through RNA interference. We studied their impact on EMT-related genes and observed a change in RNA and protein expression levels, indicating that they might play a role on EMT. Dlx4os also altered cell migration in vitro and tumorigenesis in vivo. We aim to identify lncRNAs that have specific roles in melanoma progression and understand their mode of action. The results could be used as the basis for new prognostics techniques, diagnostic and melanoma therapy, as it may contribute to a better understanding of the role of lncRNAs on the melanoma genesis and progression.

Citation Format: Ana L. Ayub, Diogo Pessoa, Lenial A. Sejour, Ioannis Vlachos, Jérémie Nsengimana, Julia Newton-Bishop, Eduardo M. Reis, Frank J. Slack, Miriam G. Jasiulionis. Transcriptome analysis from a linear melanoma progression model reveals lncRNAs involved with malignancy, EMT and metastasis [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2392.

Transcriptome of iPSC-derived neuronal cells reveals a module of co-expressed genes consistently associated with autism spectrum disorder

Evaluation of expression profile in autism spectrum disorder (ASD) patients is an important approach to understand possible similar functional consequences that may underlie disease pathophysiology regardless of its genetic heterogeneity. Induced pluripotent stem cell (iPSC)-derived neuronal models have been useful to explore this question, but larger cohorts and different ASD endophenotypes still need to be investigated. Moreover, whether changes seen in this in vitro model reflect previous findings in ASD postmortem brains and how consistent they are across the studies remain underexplored questions. We examined the transcriptome of iPSC-derived neuronal cells from a normocephalic ASD cohort composed mostly of high-functioning individuals and from non-ASD individuals. ASD patients presented expression dysregulation of a module of co-expressed genes involved in protein synthesis in neuronal progenitor cells (NPC), and a module of genes related to synapse/neurotransmission and a module related to translation in neurons. Proteomic analysis in NPC revealed potential molecular links between the modules dysregulated in NPC and in neurons. Remarkably, the comparison of our results to a series of transcriptome studies revealed that the module related to synapse has been consistently found as upregulated in iPSC-derived neurons-which has an expression profile more closely related to fetal brain-while downregulated in postmortem brain tissue, indicating a reliable association of this network to the disease and suggesting that its dysregulation might occur in different directions across development in ASD individuals. Therefore, the expression pattern of this network might be used as biomarker for ASD and should be experimentally explored as a therapeutic target.

Abstract 140: Deregulated genes by abnormal promoter methylation in early and late stages in a linear model of melanoma progression

Metastatic melanomas are extremely aggressive tumors, accounting for 70% of skin cancer deaths, which incidence is increasing globally. Despite the advances in their treatment, metastatic melanomas are still associated with a poor prognosis, with a median survival of 6 to 12 months. The study of the molecular alterations affecting key pathways during melanoma progression may contribute to a better understanding of the biology underlying the aggressiveness of the disease and lead to the identification of biomarkers for prognosis and novel targets for therapy. In addition to mutations, epigenetic alterations are important contributors to malignant transformation and tumor progression. Aberrant DNA methylation is an epigenetic hallmark of cancer, known to play important roles in melanoma development. Although different genes were described as presenting methylation alterations in their promoters in melanomas, the dynamics of these alterations along melanocyte malignant transformation and melanoma progression is still unclear. Our laboratory developed a linear model of melanoma progression in which different pre-malignant and melanoma cell lines were established after subjecting murine melanocytes to a sustained stressful condition (adhesion impediment). In this study, we used the spontaneously immortalized melanocyte lineage (melan-a); pre-malignant melanocytes (4C), obtained after subjecting melan-a cells to 4 cycles of adhesion impediment; a non-metastatic tumorigenic cell line (4C11-), obtained after a limiting dilution of spheroids formed by 4C adhesion impediment; and a metastatic cell line (4C11+), resulting from the spontaneously loss of p53 by 4C11- cells. The aim of this work was to identify genes that have their expression altered by promoter and gene body methylation that are related to the malignant transformation of melanocytes and metastasis. Through the integrative analysis of methylome (ERRBS) with transcriptome (RNAseq) data, it was possible to identify the dynamic changes in DNA methylation in promoters and gene bodies, and their relation with gene expression in both early and late stages of melanoma progression. Among those genes are Foxd1, Lrrk2, Adcy1 and Adcy3, which had their expression validated by qPCR, and their epigenetic regulation validated by the treatment of the cells with epigenetic drugs (5-aza-2'-deoxycytidine and Trichostatin A). The gene expression and promoter methylation profiles found in the murine model of melanoma progression were correlated with clinical data from melanoma patients available in cancers cohorts, revealing potential novel prognostic markers. Supported by FAPESP, CNPQ and CAPES

Citation Format: Débora Papaiz, Flávia E. Rius, Diogo Pessoa, Eduardo M. Reis, Christopher E. Mason, Gangning Liang, Miriam G. Jasiulionis. Deregulated genes by abnormal promoter methylation in early and late stages in a linear model of melanoma progression [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 140.

Where do we aspire to publish? A position paper on scientific communication in biochemistry and molecular biology

The scientific publication landscape is changing quickly, with an enormous increase in options and models. Articles can be published in a complex variety of journals that differ in their presentation format (online-only or in-print), editorial organizations that maintain them (commercial and/or society-based), editorial handling (academic or professional editors), editorial board composition (academic or professional), payment options to cover editorial costs (open access or pay-to-read), indexation, visibility, branding, and other aspects. Additionally, online submissions of non-revised versions of manuscripts prior to seeking publication in a peer-reviewed journal (a practice known as pre-printing) are a growing trend in biological sciences. In this changing landscape, researchers in biochemistry and molecular biology must re-think their priorities in terms of scientific output dissemination. The evaluation processes and institutional funding for scientific publications should also be revised accordingly. This article presents the results of discussions within the Department of Biochemistry, University of São Paulo, on this subject.

Chromatin Landscape Distinguishes the Genomic Loci of Hundreds of Androgen-Receptor-Associated LincRNAs From the Loci of Non-associated LincRNAs

Cell signaling events triggered by androgen hormone in prostate cells is dependent on activation of the androgen receptor (AR) transcription factor. Androgen hormone binding to AR promotes its displacement from the cytoplasm to the nucleus and AR binding to DNA motifs, thus inducing activatory and inhibitory transcriptional programs through a complex regulatory mechanism not yet fully understood. In this work, we performed RNA-seq deep-sequencing of LNCaP prostate cancer cells and found over 7000 expressed long intergenic non-coding RNAs (lincRNAs), of which ∼4000 are novel lincRNAs, and 258 lincRNAs have their expression activated by androgen. Immunoprecipitation of AR, followed by large-scale sequencing of co-immunoprecipitated RNAs (RIP-Seq) has identified in the LNCaP cell line a total of 619 lincRNAs that were significantly enriched (FDR < 10%, DESeq2) in the anti-Androgen Receptor (antiAR) fraction in relation to the control fraction (non-specific IgG), and we named them Androgen-Receptor-Associated lincRNAs (ARA-lincRNAs). A genome-wide analysis showed that protein-coding gene neighbors to ARA-lincRNAs had a significantly higher androgen-induced change in expression than protein-coding genes neighboring lincRNAs not associated to AR. To find relevant epigenetic signatures enriched at the ARA-lincRNAs’ transcription start sites (TSSs) we used a machine learning approach and identified that the ARA-lincRNA genomic loci in LNCaP cells are significantly enriched with epigenetic marks that are characteristic of in cis enhancer RNA regulators, and that the H3K27ac mark of active enhancers is conspicuously enriched at the TSS of ARA-lincRNAs adjacent to androgen-activated protein-coding genes. In addition, LNCaP topologically associating domains (TADs) that comprise chromatin regions with ARA-lincRNAs exhibit transcription factor contents, epigenetic marks and gene transcriptional activities that are significantly different from TADs not containing ARA-lincRNAs. This work highlights the possible involvement of hundreds of lincRNAs working in synergy with the AR on the genome-wide androgen-induced gene regulatory program in prostate cells.

Evaluating the Stability of mRNAs and Noncoding RNAs

Changes in RNA stability have an important impact in the gene expression regulation. Different methods based on the transcription blockage with RNA polymerase inhibitors or metabolic labeling of newly synthesized RNAs have been developed to evaluate RNA decay rates in cultured cell. Combined with techniques to measure transcript abundance genome-wide, these methods have been used to reveal novel features of the eukaryotic transcriptome. The stability of protein-coding mRNAs is in general closely associated to the physiological function of their encoded proteins, with short-lived mRNAs being significantly enriched among regulatory genes whereas genes associated with housekeeping functions are predominantly stable. Likewise, the stability of noncoding RNAs (ncRNAs) seems to reflect their functional role in the cell. Thus, investigating RNA stability can provide insights regarding the function of yet uncharacterized regulatory ncRNAs. In this chapter, we discuss the methodologies currently used to estimate RNA decay and outline an experimental protocol for genome-wide estimation of RNA stability of protein-coding and lncRNAs. This protocol details the transcriptional blockage of cultured cells with actinomycin D, followed by RNA isolation at different time points, the determination of transcript abundance by qPCR/DNA oligoarray hybridization, and the calculation of individual transcript half-lives.

Identification of novel biomarkers associated with poor patient outcomes in invasive breast carcinoma

Breast carcinoma (BC) corresponds to 23 % of all cancers in women, with 1.38 million new cases and 460,000 deaths worldwide annually. Despite the significant advances in the identification of molecular markers and different modalities of treatment for primary BC, the ability to predict its metastatic behavior is still limited. The purpose of this study was to identify novel molecular markers associated with distinct clinical outcomes in a Brazilian cohort of BC patients. We generated global gene expression profiles using tumor samples from 24 patients with invasive ductal BC who were followed for at least 5 years, including a group of 15 patients with favorable outcomes and another with nine patients who developed metastasis. We identified a set of 58 differentially expressed genes (p ≤ 0.01) between the two groups. The prognostic value of this metastasis signature was corroborated by its ability to stratify independent BC patient datasets according to disease-free survival and overall survival. The upregulation of B3GNT7, PPM1D, TNKS2, PHB, and GTSE1 in patients with poor outcomes was confirmed by quantitative reverse transcription polymerase chain reaction (RT-qPCR) in an independent sample of patients with BC (47 with good outcomes and eight that presented metastasis). The expression of BCL2-associated agonist of cell death (BAD) protein was determined in 1276 BC tissue samples by immunohistochemistry and was consistent with the reduced BAD mRNA expression levels in metastatic cases, as observed in the oligoarray data. These findings point to novel prognostic markers that can distinguish breast carcinomas with metastatic potential from those with favorable outcomes.

Abstract 105: Bioinformatic analysis of RNA-Seq data to search for novel diagnostic/prognostic biomarkers of pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is one of the most deadly human malignancies. The only curative treatment available is the surgical removal of the tumor in early stages of the disease. Current methods for early detection and treatment are poor, justifying more studies in this field. We aim to generate a high-resolution catalog of the PDAC transcriptome, to reveal transcriptional alterations associated with the malignant phenotype of PDAC. To that end 28 patient matched samples of PDAC and nontumor adjacent pancreatic tissue were processed for the production of rRNA subtracted cDNAs libraries and subsequently sequenced in the Illumina HiSeq platform producing 17 million 100 nt paired-end reads per sample on average. We are implementing an informatics pipeline to detect and evaluate the differential expression in PDAC of well-annotated genes, including long noncoding RNAs (lncRNAs). In addition, we will search for novel lncRNAs and alternative splicing isoforms of protein-coding genes expressed in pancreatic tissues. Following the initial quality control (FastQC) and trimming of low quality regions (Trimmomatic), reads from each sample were aligned to the human genome (TopHat2) and a count table with the number of reads per gene for each sample was created (HTSeq). For evidence of changes across experimental conditions we use DESeq2, identifying 310 genes up-regulated and 354 down-regulated (p-adjusted&lt;0.001 and lfc&lt;|3.3|), including 156 lncRNAs. The differentially expressed genes detected are involved in pathways related to cancer, regulation of cell cycle process, membrane and secreted proteins. De novo transcript assembly (Trinity) revealed a number of yet unannotated transcripts and isoforms expressed in pancreatic tissues. Altogether, the transcript catalog generated herein provides a valuable resource for the identification of novel putative candidate biomarkers to improve the early detection or evaluation of therapeutic response in PDAC.

Work supported by FAPESP, CNPq and CAPES.

Citation Format: Omar J. Sosa, Vinícius F. Paixão, João C. Setubal, Eduardo M. Reis. Bioinformatic analysis of RNA-Seq data to search for novel diagnostic/prognostic biomarkers of pancreatic ductal adenocarcinoma. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 105.

Global analysis of biogenesis, stability and sub-cellular localization of lncRNAs mapping to intragenic regions of the human genome

Long noncoding RNAs (lncRNAs) that map to intragenic regions of the human genome with the same (intronic lncRNAs) or opposite orientation (antisense lncRNAs) relative to protein-coding mRNAs have been largely dismissed from biochemical and functional characterization due to the belief that they are mRNA precursors, byproducts of RNA splicing or simply transcriptional noise. In this work, we used a custom microarray to investigate aspects of the biogenesis, processing, stability, evolutionary conservation, and cellular localization of ∼ 6,000 intronic lncRNAs and ∼ 10,000 antisense lncRNAs. Most intronic (2,903 of 3,427, 85%) and antisense lncRNAs (4,945 of 5,214, 95%) expressed in HeLa cells showed evidence of 5' cap modification, compatible with their transcription by RNAP II. Antisense lncRNAs (median t1/2 = 3.9 h) were significantly (p < 0.0001) more stable than mRNAs (median t1/2 = 3.2 h), whereas intronic lncRNAs (median t1/2 = 2.1 h) comprised a more heterogeneous class that included both stable (t1/2 > 3 h) and unstable (t1/2 < 1 h) transcripts. Intragenic lncRNAs display evidence of evolutionary conservation, have little/no coding potential and were ubiquitously detected in the cytoplasm. Notably, a fraction of the intronic and antisense lncRNAs (13 and 15%, respectively) were expressed from loci at which the corresponding host mRNA was not detected. The abundances of a subset of intronic/antisense lncRNAs were correlated (r ≥ |0.8|) with those of genes encoding proteins involved in cell division and DNA replication. Taken together, the findings of this study contribute novel biochemical and genomic information regarding intronic and antisense lncRNAs, supporting the notion that these classes include independently transcribed RNAs with potentials for exerting regulatory functions in the cell.

CIP4 promotes metastasis in triple-negative breast cancer and is associated with poor patient prognosis

Signaling via epidermal growth factor receptor (EGFR) and Src kinase pathways promote triple-negative breast cancer (TNBC) cell invasion and tumor metastasis. Here, we address the role of Cdc42-interacting protein-4 (CIP4) in TNBC metastasis in vivo, and profile CIP4 expression in human breast cancer patients. In human TNBC cells, CIP4 knock-down (KD) led to less sustained activation of Erk kinase and impaired cell motility compared to control cells. This correlated with significant defects in 3D invasion of surrounding extracellular matrix by CIP4 KD TNBC cells when grown as spheroid colonies. In mammary orthotopic xenograft assays using both human TNBC cells (MDA-MB-231, HCC 1806) and rat MTLn3 cells, CIP4 silencing had no overt effect on tumor growth, but significantly reduced the incidence of lung metastases in each tumor model. In human invasive breast cancers, high CIP4 levels was significantly associated with high tumor stage, TNBC and HER2 subtypes, and risk of progression to metastatic disease. Together, these results implicate CIP4 in promoting metastasis in TNBCs.

Long non-coding RNA INXS is a critical mediator of BCL-XS induced apoptosis

BCL-X mRNA alternative splicing generates pro-apoptotic BCL-XS or anti-apoptotic BCL-XL gene products and the mechanism that regulates splice shifting is incompletely understood. We identified and characterized a long non-coding RNA (lncRNA) named INXS, transcribed from the opposite genomic strand of BCL-X, that was 5- to 9-fold less abundant in tumor cell lines from kidney, liver, breast and prostate and in kidney tumor tissues compared with non-tumors. INXS is an unspliced 1903 nt-long RNA, is transcribed by RNA polymerase II, 5′-capped, nuclear enriched and binds Sam68 splicing-modulator. Three apoptosis-inducing agents increased INXS lncRNA endogenous expression in the 786-O kidney tumor cell line, increased BCL-XS/BCL-XL mRNA ratio and activated caspases 3, 7 and 9. These effects were abrogated in the presence of INXS knockdown. Similarly, ectopic INXS overexpression caused a shift in splicing toward BCL-XS and activation of caspases, thus leading to apoptosis. BCL-XS protein accumulation was detected upon INXS overexpression. In a mouse xenograft model, intra-tumor injections of an INXS-expressing plasmid caused a marked reduction in tumor weight, and an increase in BCL-XS isoform, as determined in the excised tumors. We revealed an endogenous lncRNA that induces apoptosis, suggesting that INXS is a possible target to be explored in cancer therapies.

Expression analysis and in silico characterization of intronic long noncoding RNAs in renal cell carcinoma: emerging functional associations

BACKGROUND

Intronic and intergenic long noncoding RNAs (lncRNAs) are emerging gene expression regulators. The molecular pathogenesis of renal cell carcinoma (RCC) is still poorly understood, and in particular, limited studies are available for intronic lncRNAs expressed in RCC.

METHODS

Microarray experiments were performed with custom-designed arrays enriched with probes for lncRNAs mapping to intronic genomic regions. Samples from 18 primary RCC tumors and 11 nontumor adjacent matched tissues were analyzed. Meta-analyses were performed with microarray expression data from three additional human tissues (normal liver, prostate tumor and kidney nontumor samples), and with large-scale public data for epigenetic regulatory marks and for evolutionarily conserved sequences.

RESULTS

A signature of 29 intronic lncRNAs differentially expressed between RCC and nontumor samples was obtained (false discovery rate (FDR) < 5%). A signature of 26 intronic lncRNAs significantly correlated with the RCC five-year patient survival outcome was identified (FDR < 5%, p-value ≤ 0.01). We identified 4303 intronic antisense lncRNAs expressed in RCC, of which 22% were significantly (p < 0.05) cis correlated with the expression of the mRNA in the same locus across RCC and three other human tissues. Gene Ontology (GO) analysis of those loci pointed to 'regulation of biological processes' as the main enriched category. A module map analysis of the protein-coding genes significantly (p < 0.05) trans correlated with the 20% most abundant lncRNAs, identified 51 enriched GO terms (p < 0.05). We determined that 60% of the expressed lncRNAs are evolutionarily conserved. At the genomic loci containing the intronic RCC-expressed lncRNAs, a strong association (p < 0.001) was found between their transcription start sites and genomic marks such as CpG islands, RNA Pol II binding and histones methylation and acetylation.

CONCLUSION

Intronic antisense lncRNAs are widely expressed in RCC tumors. Some of them are significantly altered in RCC in comparison with nontumor samples. The majority of these lncRNAs is evolutionarily conserved and possibly modulated by epigenetic modifications. Our data suggest that these RCC lncRNAs may contribute to the complex network of regulatory RNAs playing a role in renal cell malignant transformation.

Mycobacterium leprae DNA in peripheral blood may indicate a bacilli migration route and high-risk for leprosy onset

Leprosy epidemiological studies have been restricted to Mycobacterium leprae DNA detection in nasal and oral mucosa samples with scarce literature on peripheral blood. We present the largest study applying quantitative real-time PCR (qPCR) for the detection of M. leprae DNA in peripheral blood samples of 200 untreated leprosy patients and 826 household contacts, with results associated with clinical and laboratory parameters. To detect M. leprae DNA a TaqMan qPCR assay targeting the M. leprae ML0024 genomic region was performed. The ML0024 qPCR in blood samples detected the presence of bacillus DNA in 22.0% (44/200) of the leprosy patients: 23.2% (16/69) in paucibacillary (PB), and 21.4% (28/131) in multibacillary (MB) patients. Overall positivity among contacts was 1.2% (10/826), with similar percentages regardless of whether the index case was PB or MB. After a follow-up period of 7 years, 26 contacts have developed leprosy. Comparing the results of healthy contacts with those that become ill, ML0024 qPCR positivity at the time of diagnosis of their index case represented an impressive 14.78-fold greater risk for leprosy onset (95% CI 3.6-60.8; p <0.0001). In brief, contacts with positive PCR in blood at diagnosis of index cases are at higher risk of later leprosy onset and this marker might be combined with other prognostic markers for management of contacts, which requires further studies.

The intronic long noncoding RNA ANRASSF1 recruits PRC2 to the RASSF1A promoter, reducing the expression of RASSF1A and increasing cell proliferation

The down-regulation of the tumor-suppressor gene RASSF1A has been shown to increase cell proliferation in several tumors. RASSF1A expression is regulated through epigenetic events involving the polycomb repressive complex 2 (PRC2); however, the molecular mechanisms modulating the recruitment of this epigenetic modifier to the RASSF1 locus remain largely unknown. Here, we identify and characterize ANRASSF1, an endogenous unspliced long noncoding RNA (lncRNA) that is transcribed from the opposite strand on the RASSF1 gene locus in several cell lines and tissues and binds PRC2. ANRASSF1 is transcribed through RNA polymerase II and is 5′-capped and polyadenylated; it exhibits nuclear localization and has a shorter half-life compared with other lncRNAs that bind PRC2. ANRASSF1 endogenous expression is higher in breast and prostate tumor cell lines compared with non-tumor, and an opposite pattern is observed for RASSF1A. ANRASSF1 ectopic overexpression reduces RASSF1A abundance and increases the proliferation of HeLa cells, whereas ANRASSF1 silencing causes the opposite effects. These changes in ANRASSF1 levels do not affect the RASSF1C isoform abundance. ANRASSF1 overexpression causes a marked increase in both PRC2 occupancy and histone H3K27me3 repressive marks, specifically at the RASSF1A promoter region. No effect of ANRASSF1 overexpression was detected on PRC2 occupancy and histone H3K27me3 at the promoter regions of RASSF1C and the four other neighboring genes, including two well-characterized tumor suppressor genes. Additionally, we demonstrated that ANRASSF1 forms an RNA/DNA hybrid and recruits PRC2 to the RASSF1A promoter. Together, these results demonstrate a novel mechanism of epigenetic repression of the RASSF1A tumor suppressor gene involving antisense unspliced lncRNA, in which ANRASSF1 selectively represses the expression of the RASSF1 isoform overlapping the antisense transcript in a location-specific manner. In a broader perspective, our findings suggest that other non-characterized unspliced intronic lncRNAs transcribed in the human genome might contribute to a location-specific epigenetic modulation of genes.

RASSF1A is a tumor suppressor gene whose expression is repressed through epigenetic events in a wide range of different cancers. Repression is effected by DNA hypermethylation of the RASSF1A promoter, which in turn is triggered through histone H3K9/H3K27 trimethylation repressive marks. The addition of the H3K27me3 mark at the RASSF1A promoter locus involves the polycomb repressive complex 2 (PRC2). The molecular mechanisms that control the recruitment of PRC2 to the promoter to initiate H3K27 trimethylation and repress RASSF1A expression have not been described. Here, we identified a long noncoding RNA (lncRNA), termed ANRASSF1 for antisense noncoding RASSF1, that is transcribed from the opposite strand of the RASSF1A gene and is responsible for recruiting PRC2 to the RASSF1A promoter region in a highly location-specific manner. No effect of ANRASSF1 was detected on the promoter of the RASSF1C isoform or the promoters of the four other genes within the vicinity of RASSF1, including two other well-characterized tumor suppressor genes. This work provides evidence that the epigenetic modulation of the tumor suppressor gene RASSF1A is dependent on the lncRNA ANRASSF1 and highlights the importance of further studies on the involvement of ANRASSF1 in tumorigenesis.

Loss of caspase 7 expression is associated with poor prognosis in renal cell carcinoma clear cell subtype

OBJECTIVE

To investigate the expression of CASP7 protein in renal cell carcinoma clear cell subtype (ccRCC) and its value to predict cancer-specific survival (CSS).

METHODS

A tissue microarray containing 120 samples of ccRCC, 45 non-ccRCC, and 66 nontumor paired samples from patients who underwent partial or radical nephrectomy was hybridized with anti-CASP7 antibody. Tissue sections were scored according to intensity and the percentage of stained cells. CASP7 immunostaining scores were used to estimate the association with clinicopathologic parameters and calculate Kaplan-Meier survival curves.

RESULTS

Reduced CASP7 expression was observed in ccRCC and non-ccRCC subtypes in comparison with nontumor renal tissues (P <.0001). CASP7 immunostaining was associated (P <.05) with clinicopathologic parameters (size, incidental tumor, clinical stage, renal vein invasion, and tumor necrosis) and correlated with CSS (P = .032) and global survival (P = .046) of patients with ccRCC. In addition, CASP7 expression was able to substratify patients with ccRCC with favorable prognosis according to low clinical stage, in which negative CASP7 staining was associated with patients with lower CSS (P = .045). Finally, CASP7 staining was able to provide significant stratification according to CSS (P = .018) among patients with ccRCC with disease relapse.

CONCLUSION

Our results implicate the loss of CASP7 expression in the aggressiveness of ccRCC and indicate its potential use as a clinical prognostic marker of the disease.

Perspectives of Long Non-Coding RNAs in Cancer Diagnostics

Long non-coding RNAs (lncRNAs) transcribed from intergenic and intronic regions of the human genome constitute a broad class of cellular transcripts that are under intensive investigation. While only a handful of lncRNAs have been characterized, their involvement in fundamental cellular processes that control gene expression highlights a central role in cell homeostasis. Not surprisingly, aberrant expression of regulatory lncRNAs has been increasingly documented in different types of cancer, where they can mediate both oncogenic or tumor suppressor effects. Interaction with chromatin remodeling complexes that promote silencing of specific genes or modulation of splicing factor proteins seem to be two general modes of lncRNA regulation, but it is conceivable that additional mechanisms of action are yet to be unveiled. LncRNAs show greater tissue specificity compared to protein-coding mRNAs making them attractive in the search of novel diagnostics/prognostics cancer biomarkers in body fluid samples. In fact, lncRNA prostate cancer antigen 3 can be detected in urine samples and has been shown to improve diagnosis of prostate cancer. We suggest that an unbiased screening of the presence of RNAs in easily accessible body fluids such as serum and urine might reveal novel circulating lncRNAs as potential biomarkers in many types of cancer. Annotation and functional characterization of the lncRNA complement of the cancer transcriptome will conceivably provide new venues for early diagnosis and treatment of the disease.

Long noncoding intronic RNAs are differentially expressed in primary and metastatic pancreatic cancer

Background

Pancreatic ductal adenocarcinoma (PDAC) is known by its aggressiveness and lack of effective therapeutic options. Thus, improvement in current knowledge of molecular changes associated with pancreatic cancer is urgently needed to explore novel venues of diagnostics and treatment of this dismal disease. While there is mounting evidence that long noncoding RNAs (lncRNAs) transcribed from intronic and intergenic regions of the human genome may play different roles in the regulation of gene expression in normal and cancer cells, their expression pattern and biological relevance in pancreatic cancer is currently unknown. In the present work we investigated the relative abundance of a collection of lncRNAs in patients' pancreatic tissue samples aiming at identifying gene expression profiles correlated to pancreatic cancer and metastasis.

Methods

Custom 3,355-element spotted cDNA microarray interrogating protein-coding genes and putative lncRNA were used to obtain expression profiles from 38 clinical samples of tumor and non-tumor pancreatic tissues. Bioinformatics analyses were performed to characterize structure and conservation of lncRNAs expressed in pancreatic tissues, as well as to identify expression signatures correlated to tissue histology. Strand-specific reverse transcription followed by PCR and qRT-PCR were employed to determine strandedness of lncRNAs and to validate microarray results, respectively.

Results

We show that subsets of intronic/intergenic lncRNAs are expressed across tumor and non-tumor pancreatic tissue samples. Enrichment of promoter-associated chromatin marks and over-representation of conserved DNA elements and stable secondary structure predictions suggest that these transcripts are generated from independent transcriptional units and that at least a fraction is under evolutionary selection, and thus potentially functional.

Statistically significant expression signatures comprising protein-coding mRNAs and lncRNAs that correlate to PDAC or to pancreatic cancer metastasis were identified. Interestingly, loci harboring intronic lncRNAs differentially expressed in PDAC metastases were enriched in genes associated to the MAPK pathway. Orientation-specific RT-PCR documented that intronic transcripts are expressed in sense, antisense or both orientations relative to protein-coding mRNAs. Differential expression of a subset of intronic lncRNAs (PPP3CB, MAP3K14 and DAPK1 loci) in metastatic samples was confirmed by Real-Time PCR.

Conclusion

Our findings reveal sets of intronic lncRNAs expressed in pancreatic tissues whose abundance is correlated to PDAC or metastasis, thus pointing to the potential relevance of this class of transcripts in biological processes related to malignant transformation and metastasis in pancreatic cancer.

Gene expression profiling reveals molecular marker candidates of laryngeal squamous cell carcinoma

Laryngeal squamous cell carcinoma is very common in head and neck cancer, with high mortality rates and poor prognosis. In this study, we compared expression profiles of clinical samples from 13 larynx tumors and 10 non-neoplastic larynx tissues using a custom-built cDNA microarray containing 331 probes for 284 genes previously identified by informatics analysis of EST databases as markers of head and neck tumors. Thirty-five genes showed statistically significant differences (SNR > or = | 1.0 |, p< or =0.001) in the expression between tumor and non-tumor larynx tissue samples. Functional annotation indicated that these genes are involved in cellular processes relevant to the cancer phenotype, such as apoptosis, cell cycle, DNA repair, proteolysis, protease inhibition, signal transduction and transcriptional regulation. Six of the identified transcripts map to intronic regions of protein-coding genes and may comprise non-annotated exons or as yet uncharacterized long ncRNAs with a regulatory role in the gene expression program of larynx tissue. The differential expression of 10 of these genes (ADCY6, AES, AL2SCR3, CRR9, CSTB, DUSP1, MAP3K5, PLAT, UBL1 and ZNF706) was independently confirmed by quantitative real-time RT-PCR. Among these, the CSTB gene product has cysteine protease inhibitor activity that has been associated with an antimetastatic function. Interestingly, CSTB showed a low expression in the tumor samples analyzed (p<0.0001). The set of genes identified here contribute to a better understanding of the molecular basis of larynx cancer, and provide candidate markers for improving diagnosis, prognosis and treatment of this carcinoma.

Gene expression profiling of cultured cells from brainstem of newborn spontaneously hypertensive and Wistar Kyoto rats

The spontaneously hypertensive rat (SHR) is a good model to study several diseases such as the attention-deficit hyperactivity disorder, cardiopulmonary impairment, nephropathy, as well as hypertension, which is a multifactor disease that possibly involves alterations in gene expression in hypertensive relative to normotensive subjects. In this study, we used high-density oligoarrays to compare gene expression profiles in cultured neurons and glia from brainstem of newborn normotensive Wistar Kyoto (WKY) and SHR rats. We found 376 genes differentially expressed between SHR and WKY brainstem cells that preferentially map to 17 metabolic/signaling pathways. Some of the pathways and regulated genes identified herein are obviously related to cardiovascular regulation; in addition there are several genes differentially expressed in SHR not yet associated to hypertension, which may be attributed to other differences between SHR and WKY strains. This constitute a rich resource for the identification and characterization of novel genes associated to phenotypic differences observed in SHR relative to WKY, including hypertension. In conclusion, this study describes for the first time the gene profiling pattern of brainstem cells from SHR and WKY rats, which opens up new possibilities and strategies of investigation and possible therapeutics to hypertension, as well as for the understanding of the brain contribution to phenotypic differences between SHR and WKY rats.

Identification of protein-coding and intronic noncoding RNAs down-regulated in clear cell renal carcinoma

The clear cell subtype of renal cell carcinoma (RCC) is the most lethal and prevalent cancer of the urinary system. To investigate the molecular changes associated with malignant transformation in clear cell RCC, the gene expression profiles of matched samples of tumor and adjacent non-neoplastic tissue were obtained from six patients. A custom-built cDNA microarray platform was used, comprising 2292 probes that map to exons of genes and 822 probes for noncoding RNAs mapping to intronic regions. Intronic transcription was detected in all normal and neoplastic renal tissues. A subset of 55 transcripts was significantly down-regulated in clear cell RCC relative to the matched nontumor tissue as determined by a combination of two statistical tests and leave-one-out patient cross-validation. Among the down-regulated transcripts, 49 mapped to untranslated or coding exons and 6 were intronic relative to known exons of protein-coding genes. Lower levels of expression of SIN3B, TRIP3, SYNJ2BP and NDE1 (P < 0.02), and of intronic transcripts derived from SND1 and ACTN4 loci (P < 0.05), were confirmed in clear cell RCC by Real-time RT-PCR. A subset of 25 transcripts was deregulated in additional six nonclear cell RCC samples, pointing to common transcriptional alterations in RCC irrespective of the histological subtype or differentiation state of the tumor. Our results indicate a novel set of tumor suppressor gene candidates, including noncoding intronic RNAs, which may play a significant role in malignant transformations of normal renal cells.

Detailed molecular characterization of cord blood-derived endothelial progenitors

OBJECTIVE

Given their involvement in pathological and physiological angiogenesis, there has been growing interest in understanding and manipulating endothelial progenitor cells (EPC) for therapeutic purposes. However, detailed molecular analysis of EPC before and during endothelial differentiation is lacking and is the subject of the present study.

MATERIALS AND METHODS

We report a detailed microarray gene-expression profile of freshly isolated (day 0) human cord blood (CB)-derived EPC (CD133+KDR+ or CD34+KDR+), and at different time points during in vitro differentiation (early: day 13; late: day 27).

RESULTS

Data obtained reflect an EPC transcriptome enriched in genes related to stem/progenitor cells properties (chromatin remodeling, self-renewal, signaling, cytoskeleton organization and biogenesis, recruitment, and adhesion). Using a complementary DNA microarray enriched in intronic transcribed sequences, we observed, as well, that naturally transcribed intronic noncoding RNAs were specifically expressed at the EPC stage.

CONCLUSION

Taken together, we have defined the global gene-expression profile of CB-derived EPC during the process of endothelial differentiation, which can be used to identify genes involved in different vascular pathologies.

Splice variants of TLE family genes and up-regulation of a TLE3 isoform in prostate tumors

The TLE genes constitute a family of important transcriptional co-repressors involved in many cellular processes. We found evidence of alternatively spliced mRNAs for human TLE1-4 containing premature stop codons, thus encoding putative shortened proteins. Microarray experiments and Real-time RT-PCR assays showed that alternatively spliced isoforms of TLE1, TLE2 an d TLE3 were preferentially expressed in prostate in comparison to liver and kidney tissues. We identified by orientation-specific R T-PCR an antisense partially intronic non-coding RNA that overlaps a novel exon of the TLE3 gene, raising the possibility of regulation of alternative splicing by this non-coding transcript. The alternatively spliced isoform of TLE3 was up-regulated (6- to 17-fo ld) in prostate tumors in comparison to matched non-tumor adjacent tissue from 7 out of 11 (64%) patients and in four prostate tumor cell lines in comparison to a normal prostate cell line. These results demonstrate that different isoforms of TLE genes are commonly transcribed in human tissues and suggest that TLE3 could be involved in prostate cancer development.

Apert p.Ser252Trp mutation in FGFR2 alters osteogenic potential and gene expression of cranial periosteal cells

Apert syndrome (AS), a severe form of craniosynostosis, is caused by dominant gain-of-function mutations in FGFR2. Because the periosteum contribution to AS cranial pathophysiology is unknown, we tested the osteogenic potential of AS periosteal cells (p.Ser252Trp mutation) and observed that these cells are more committed toward the osteoblast lineage. To delineate the gene expression profile involved in this abnormal behavior, we performed a global gene expression analysis of coronal suture periosteal cells from seven AS patients (p.Ser252Trp), and matched controls. We identified 263 genes with significantly altered expression in AS samples (118 upregulated, 145 downregulated; SNR >or= |0.4|, P <or= 0.05). Several upregulated genes are involved in positive regulation of cell proliferation and nucleotide metabolism, whereas several downregulated genes are involved in inhibition of cell proliferation, gene expression regulation, cell adhesion, and extracellular matrix organization, and in PIK3-MAPK cascades. AS expression profile was confirmed through real-time PCR of a selected set of genes using RNAs from AS and control cells as well as from control cells treated with high FGF2 concentration, and through the analysis of genes involved in FGF-FGFR signaling. Our results allowed us to: (a) suggest that AS periosteal cells present enhanced osteogenic potential, (b) unravel a specific gene expression signature characteristic of AS periosteal cells which may be associated with their osteogenic commitment, (c) identify a set of novel genes involved in the pathophysiology of AS or other craniosynostotic conditions, and (d) suggest for the first time that the periosteum might be involved in the pathophysiology of AS.

Apert p.Ser252Trp mutation in FGFR2 alters osteogenic potential and gene expression of cranial periosteal cells

Apert syndrome (AS), a severe form of craniosynostosis, is caused by dominant gain-of-function mutations in FGFR2. Because the periosteum contribution to AS cranial pathophysiology is unknown, we tested the osteogenic potential of AS periosteal cells (p.Ser252Trp mutation) and observed that these cells are more committed toward the osteoblast lineage. To delineate the gene expression profile involved in this abnormal behavior, we performed a global gene expression analysis of coronal suture periosteal cells from seven AS patients (p.Ser252Trp), and matched controls. We identified 263 genes with significantly altered expression in AS samples (118 upregulated, 145 downregulated; SNR >or= |0.4|, P <or= 0.05). Several upregulated genes are involved in positive regulation of cell proliferation and nucleotide metabolism, whereas several downregulated genes are involved in inhibition of cell proliferation, gene expression regulation, cell adhesion, and extracellular matrix organization, and in PIK3-MAPK cascades. AS expression profile was confirmed through real-time PCR of a selected set of genes using RNAs from AS and control cells as well as from control cells treated with high FGF2 concentration, and through the analysis of genes involved in FGF-FGFR signaling. Our results allowed us to: (a) suggest that AS periosteal cells present enhanced osteogenic potential, (b) unravel a specific gene expression signature characteristic of AS periosteal cells which may be associated with their osteogenic commitment, (c) identify a set of novel genes involved in the pathophysiology of AS or other craniosynostotic conditions, and (d) suggest for the first time that the periosteum might be involved in the pathophysiology of AS.

Androgen responsive intronic non-coding RNAs

BACKGROUND

Transcription of large numbers of non-coding RNAs originating from intronic regions of human genes has been recently reported, but mechanisms governing their biosynthesis and biological functions are largely unknown. In this work, we evaluated the existence of a common mechanism of transcription regulation shared by protein-coding mRNAs and intronic RNAs by measuring the effect of androgen on the transcriptional profile of a prostate cancer cell line.

RESULTS

Using a custom-built cDNA microarray enriched in intronic transcribed sequences, we found 39 intronic non-coding RNAs for which levels were significantly regulated by androgen exposure. Orientation-specific reverse transcription-PCR indicated that 10 of the 13 were transcribed in the antisense direction. These transcripts are long (0.5-5 kb), unspliced and apparently do not code for proteins. Interestingly, we found that the relative levels of androgen-regulated intronic transcripts could be correlated with the levels of the corresponding protein-coding gene (asGAS6 and asDNAJC3) or with the alternative usage of exons (asKDELR2 and asITGA6) in the corresponding protein-coding transcripts. Binding of the androgen receptor to a putative regulatory region upstream from asMYO5A, an androgen-regulated antisense intronic transcript, was confirmed by chromatin immunoprecipitation.

CONCLUSION

Altogether, these results indicate that at least a fraction of naturally transcribed intronic non-coding RNAs may be regulated by common physiological signals such as hormones, and further corroborate the notion that the intronic complement of the transcriptome play functional roles in the human gene-expression program.

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

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.

Auto-antibodies in prostate cancer: humoral immune response to antigenic determinants coded by the differentially expressed transcripts FLJ23438 and VAMP3

BACKGROUND

Here we evaluate auto-antibody response against two potential antigenic determinants of genes highly expressed in low Gleason Score prostate cancer (PC) tumor samples, namely FLJ23438 and VAMP3.

METHODS

RT-PCR assays were used to analyze mRNA expression profiles of FLJ23438 and VAMP3 transcripts. The auto-antibody response against FLJ23438 and VAMP3 recombinant proteins was tested by immunoblot assays using PC, benign prostate hyperplasia (BPH), healthy donors (HD), and other human cancers plasma samples.

RESULTS

Our data showed that 37% (10/27) and 7.4% (2/27) of PC plasma samples presented auto-antibodies against FLJ23438 and VAMP3, respectively. Only 8.3% (1/12) of BPH plasma samples were reactive for both auto-antibodies, while none (0/12) of HD plasma samples tested were reactive.

CONCLUSIONS

The prevalence of 37% of positive PC plasma samples for anti-FLJ23438 antibodies suggests that humoral immune response against this antigenic determinant could be a potential serum marker for this cancer.

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

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.

As antisense RNA gets intronic

Recent work describing the transcriptional output of the human genome points to the existence of a significant number of non-coding RNA transcripts coming from intronic regions, with a fraction of these being oriented antisense relative to the protein-coding mRNA of the known gene. In this article, we survey the main findings of the large-scale expression analysis projects that led to the identification of antisense intronic messages and which demonstrate their ubiquitous expression in the human genome. We review the current knowledge on long, unspliced, intronic antisense transcripts, a new class of non-coding RNAs, recently described by our group to be correlated with the degree of tumor differentiation in prostate cancer, which we postulate is involved in the fine tuning of gene expression in eukaryotes. Possible mechanisms of antisense intronic transcript biogenesis and function in gene expression regulation are discussed, as is their involvement in human diseases. While there is still no conclusive evidence demonstrating a functional role for these long, intronic antisense messages, the far-reaching implications of their existence for the mechanisms regulating gene expression certainly warrant further experimentation.

Gene expression arrays in cancer research: methods and applications

During the last 5 years, the number of papers describing data obtained by microarray technology increased exponentially with about 3000 papers in 2003. Undoubtedly, cancer is by far the disease that received most of the attention as far as the amount of data generated. As array technology is rather new and highly dependent on bioinformatics, mathematics and statistics, a clear understanding of the knowledge and information derived from array-based experiments is not widely appreciated. We shall review herein some of the issues related to the construction of DNA arrays, quantities and heterogeneity of probes and targets, the consequences of the physical characteristics of the probes, data extraction and data analysis as well as the applications of array technology. Our goal is to bring to the general audience, some of the basics of array technology and its possible application in oncology. By discussing some of the basic aspects of the methodology, we hope to stimulate criticism concerning the conclusions proposed by authors, especially in the light of the very low degree of reproducibility already proven when commercially available platforms were compared . Regardless of its pitfalls, it is unquestionable that array technology will have a great impact in the management of cancer and its applications will range from the discovery of new drug targets, new molecular tools for diagnosis and prognosis as well as for a tailored treatment that will take into account the molecular determinants of a given tumor. Hence, we shall also highlight some of the already available and promising applications of array technology on the day-to-day practice of oncology.

Large-scale Transcriptome Analyses Reveal New Genetic Marker Candidates of Head, Neck, and Thyroid Cancer

A detailed genome mapping analysis of 213,636 expressed sequence tags (EST) derived from nontumor and tumor tissues of the oral cavity, larynx, pharynx, and thyroid was done. Transcripts matching known human genes were identified; potential new splice variants were flagged and subjected to manual curation, pointing to 788 putatively new alternative splicing isoforms, the majority (75%) being insertion events. A subset of 34 new splicing isoforms (5% of 788 events) was selected and 23 (68%) were confirmed by reverse transcription-PCR and DNA sequencing. Putative new genes were revealed, including six transcripts mapped to well-studied chromosomes such as 22, as well as transcripts that mapped to 253 intergenic regions. In addition, 2,251 noncoding intronic RNAs, eventually involved in transcriptional regulation, were found. A set of 250 candidate markers for loss of heterozygosis or gene amplification was selected by identifying transcripts that mapped to genomic regions previously known to be frequently amplified or deleted in head, neck, and thyroid tumors. Three of these markers were evaluated by quantitative reverse transcription-PCR in an independent set of individual samples. Along with detailed clinical data about tumor origin, the information reported here is now publicly available on a dedicated Web site as a resource for further biological investigation. This first in silico reconstruction of the head, neck, and thyroid transcriptomes points to a wealth of new candidate markers that can be used for future studies on the molecular basis of these tumors. Similar analysis is warranted for a number of other tumors for which large EST data sets are available.

Antisense intronic non-coding RNA levels correlate to the degree of tumor differentiation in prostate cancer

A large fraction of transcripts are expressed antisense to introns of known genes in the human genome. Here we show the construction and use of a cDNA microarray platform enriched in intronic transcripts to assess their biological relevance in pathological conditions. To validate the approach, prostate cancer was used as a model, and 27 patient tumor samples with Gleason scores ranging from 5 to 10 were analyzed. We find that a considerably higher fraction (6.6%, [23/346]) of intronic transcripts are significantly correlated (P< or =0.001) to the degree of prostate tumor differentiation (Gleason score) when compared to transcripts from unannotated genomic regions (1%, [6/539]) or from exons of known genes (2%, [27/1369]). Among the top twelve transcripts most correlated to tumor differentiation, six are antisense intronic messages as shown by orientation-specific RT-PCR or Northern blot analysis with strand-specific riboprobe. Orientation-specific real-time RT-PCR with six tumor samples, confirmed the correlation (P=0.024) between the low/high degrees of tumor differentiation and antisense intronic RASSF1 transcript levels. The need to use intron arrays to reveal the transcriptome profile of antisense intronic RNA in cancer has clearly emerged.

RASL11A, member of a novel small monomeric GTPase gene family, is down-regulated in prostate tumors

We performed a genome-wide search for novel loci encoding for Ras-related proteins based on the genome mapping coordinates of the cancer-derived EST dataset at GenBank. Partial sequences from two novel human genes were identified and subsequently used for full length transcript cloning. RASL11A and ARL9 belong to two novel subfamilies coding for small GTPases that we found to be highly conserved among eukaryotes. The Arl9/Arl10 subfamily displays a conserved interswitch toggle that places it evolutionarily closer to the Arf family. Rasl11 proteins are more closely related to the Ras branch of GTPases. All orthologues newly identified here exhibit an Asn residue in place of the highly conserved Thr35 of the G domain, suggesting that the universal switch mechanism of small GTPases may be structurally different in this subfamily. We determined by Northern blot that RASL11A is transcribed in several human tissues and that it is down-regulated in prostate tumors as measured by quantitative real-time PCR. These results highlight a previously uncharacterized subfamily of Ras-related genes that may have a tumor suppressor role in prostate cancer.

Transcriptome analysis of the acoelomate human parasite Schistosoma mansoni

Schistosoma mansoni is the primary causative agent of schistosomiasis, which affects 200 million individuals in 74 countries. We generated 163,000 expressed-sequence tags (ESTs) from normalized cDNA libraries from six selected developmental stages of the parasite, resulting in 31,000 assembled sequences and 92% sampling of an estimated 14,000 gene complement. By analyzing automated Gene Ontology assignments, we provide a detailed view of important S. mansoni biological systems, including characterization of metazoa-specific and eukarya-conserved genes. Phylogenetic analysis suggests an early divergence from other metazoa. The data set provides insights into the molecular mechanisms of tissue organization, development, signaling, sexual dimorphism, host interactions and immune evasion and identifies novel proteins to be investigated as vaccine candidates and potential drug targets.

ESTWeb: bioinformatics services for EST sequencing projects

ESTWeb is an internet based software package designed for uniform data processing and storage for large-scale EST sequencing projects. The package provides for: (a) reception of sequencing chromatograms; (b) sequence processing such as base-calling, vector screening, comparison with public databases; (c) storage of data and analysis in a relational database, (d) generation of a graphical report of individual sequence quality; and (e) issuing of reports with statistics of productivity and redundancy. The software facilitates real-time monitoring and evaluation of EST sequence acquisition progress along an EST sequencing project.

Zerg: a very fast BLAST parser library

SUMMARY

Zerg is a library of sub-routines that parses the output from all NCBI BLAST programs (Blastn, Blastp, Blastx, Tblastn and Tblastx) and returns the attributes of a BLAST report to the user. It is optimized for speed, being especially useful for large-scale genomic analysis. Benchmark tests show that Zerg is over two orders of magnitude faster than some widely used BLAST parsers.

AVAILABILITY

http://bioinfo.iq.usp.br/zerg

The contribution of 700,000 ORF sequence tags to the definition of the human transcriptome

Open reading frame expressed sequences tags (ORESTES) differ from conventional ESTs by providing sequence data from the central protein coding portion of transcripts. We generated a total of 696,745 ORESTES sequences from 24 human tissues and used a subset of the data that correspond to a set of 15,095 full-length mRNAs as a means of assessing the efficiency of the strategy and its potential contribution to the definition of the human transcriptome. We estimate that ORESTES sampled over 80% of all highly and moderately expressed, and between 40% and 50% of rarely expressed, human genes. In our most thoroughly sequenced tissue, the breast, the 130,000 ORESTES generated are derived from transcripts from an estimated 70% of all genes expressed in that tissue, with an equally efficient representation of both highly and poorly expressed genes. In this respect, we find that the capacity of the ORESTES strategy both for gene discovery and shotgun transcript sequence generation significantly exceeds that of conventional ESTs. The distribution of ORESTES is such that many human transcripts are now represented by a scaffold of partial sequences distributed along the length of each gene product. The experimental joining of the scaffold components, by reverse transcription-PCR, represents a direct route to transcript finishing that may represent a useful alternative to full-length cDNA cloning.

Vibrio cholerae resistant to 2,4-diamino-6,7-diisopropylpteridine (O/129) isolated from patients with enteritis in Ceará, Brazil

This paper reports the characterization of clinical Vibrio cholerae resistant to vibriostatic agent O/129, using classical and plasmid analysis. In a study conducted during December 1991-September 1993, two of 7,058 V. cholerae strains, obtained from patients suspected to have cholera in the State of Ceará, northeast Brazil, were resistant to 150 micrograms of the vibriostatic agent O/129 (2,4-diamino-6,7-diisopropylpteridine). One strain was identified as V. cholerae O1 El Tor Inaba and the other one as serogroup O22. Only one O1 strain harboured a plasmid of 147 kb transferable to Escherichia coli K12, and five strains of V. cholerae O1 and non-O1 were sensitive to O/129 and plasmid-negative at a frequency between 8 x 10(-2) and 3.6 x 10(-5). Additionally, O/129-resistant strains of V. cholerae O1 and O22 were resistant to trimethoprim/sulphamethoxazole.

Biological treatment of wastewater from the cassava meal industry

Biological treatments of cassava meal processing wastewaters were investigated by aerobic and combined anaerobic/aerobic reactors. As a pretreatment, flocculation and sedimentation with aluminum salts and natural polyelectrolytes were employed, in order to change the effluent concentration of organics from 14,000 to 2000 mg/L in the bench scale reactor, with a hydraulic retention time of 37 min and an influent rate of 0.56 L/cm. Biological degradation by an activated sludge process was carried out using a reactor volume of 20 L. The observed organic matter removal rates were 89 to 93%, the cyanide removal was 95 to 99%, and the food to microorganisms ratio was found to be in the range of 0.166 to 0.242/day, with a hydraulic residence time of 1.4 to 4.2 days. As a posttreatment, flocculation and coagulation were employed, resulting in an effluent of good quality, as shown by turbidity measurements and negative tests for fecal and total coliforms.

N-terminal chimeric constructs improve the expression of sarcoplasmic reticulum Ca(2+)-ATPase in yeast

Wild-type and chimeric constructs comprising rabbit sarcoplasmic reticulum (SR) Ca(2+)-ATPase and the N-terminal cytoplasmic portion of yeast plasma membrane H(+)-ATPase were expressed in yeast under control of a heat-shock regulated promoter. The wild-type ATPase was found predominantly in endoplasmic reticulum (ER) membranes. Addition of the first 88 residues of H(+)-ATPase to the Ca(2+)-ATPase N-terminal end promoted a marked shift in the localization of chimeric H(+)/Ca(2+)-ATPase which accumulated in a light membrane fraction associated with yeast smooth ER. Furthermore, there was a three-fold increase in the overall level of expression of chimeric H(+)/Ca(2+)-ATPase. Similar results were obtained for a chimeric Ca(2+)-ATPase containing a hexahistidine sequence added to its N-terminal end. Both H(+)/Ca(2+)-ATPase and 6xHis-Ca(2+)-ATPase were functional as demonstrated by their ability to form a phosphorylated intermediate and undergo fast turnover. Conversely, a replacement chimera in which the N-terminal end of SR Ca(2+)-ATPase was replaced by the corresponding segment of H(+)-ATPase was not stably expressed in yeast membranes. These results indicate that the N-terminal segment of Ca(2+)-ATPase plays an important role in enzyme assembly and contains structural determinants necessary for ER retention of the ATPase.

Structural and phylogenetic relationships among plant and animal cystatins

The plant cystatins or phytocystatins (PhyCys) are cysteine proteinase inhibitors containing the QxVxG motif and have been placed in the cystatin superfamily of proteins. The primary sequences of PhyCys have a high degree of homology with the members of the cystatin family, but they resemble stefins by the absence of disulfide bonds and cysteine residues. A multialignment and a phylogenetic analysis of 63 cystatins, 32 of which are PhyCys, demonstrate that all PhyCys cluster in a major evolutionary tree branch and support the classification of PhyCys as a new cystatin family. The PhyCys also possess a specific consensus sequence [LVI]-[AGT]-[RKE]-[FY]-[AS]-[VI]-x-[EDQV]-[HYFQ] -N placed on the region corresponding to a predictable amino-terminal alpha-helix. This sequence can be used to specifically identify PhyCys on protein data banks and to differentiate them from the other members of the superfamily.

Colimetry of marine waters off Fortaleza (Ceará State, Brazil) and detection of enteropathogenic Escherichia coli strains

Bacteriological analyses of seawater from three main beaches in Fortaleza, Brazil were performed during 1997. Thirty-six samples per beach were collected for a total of 108 samples. For Meireles Beach, 44% of the samples had MPN total coliforms values of at least 1100 or over 2400/100 ml, followed by Formosa and Diários beaches showing lower counts. For fecal coliforms the highest numbers were demonstrated for Formosa, followed by Meireles and Diários beaches in this descending order: 13.0%, 11.1% and 8.3%, respectively. Escherichia coli strains were identified in 76.8% of the 108 samples. Among 295 strains of E. coli, 21 belonged to serogroups O25, O26, O91, O112, O119, O158 and O164. Strains from serogroup O26 were tested using PCR, ELISA and Vero cells to detect Verotoxins VT1 and VT2 and all strains were negative. No LT and ST, as determined by ELISA and suckling mice assays, were detected among the 295 strains. All strains of E. coli were sensitive to ampicillin, cephalothin, gentamicin, tetracycline, sulfametox-trimethoprim, chloramphenicol and ciprofloxacin. Although the E. coli strains were not toxigenic, their presence in high numbers could be of public health significance.

Heterologous expression of sarcoplasmic reticulum Ca(2+)-ATPase

In recent years, expression of rabbit sarcoplasmic reticulum (SR) Ca(2+)-ATPase in heterologous systems has been a widely used strategy to study altered enzymes generated by site-directed mutagenesis. Various eukaryotic expression systems have been tested, all of them yielding comparable amounts of recombinant protein. However, the relatively low yield of recombinant protein obtained so far suggests that novel purification techniques will be required to allow further characterization of this enzyme based on direct ligand-binding measurements.

Helicobacter pylori in dyspeptic children and adults: endoscopic, bacteriologic and histologic correlations

Using different bacteriological (urease test, Gram staining and culture) and histological (Steiner staining and modified Giemsa staining) techniques, we searched for the presence of Helicobacter pylori in the gastric antrum of 200 dyspeptic Brazilian patients (106 females and 94 males aged 19 days to 81 years). The presence of bacteria was then correlated with the endoscopic and histological findings. H. pylori was present in 59.5% of the population studied. In Brazil, colonization occurs early, involving 37% of the dyspeptic population by 20 years of age. The presence of H. pylori in the gastric antrum was strongly associated with duodenal ulcer (P < 0.001) and a normal endoscopic examination did not exclude the possibility of colonization of the gastric antrum by H. pylori. The most sensitive test was the preformed urease test (89%). We conclude that more than one diagnostic method should preferably be used for the detection of H. pylori and that the presence of H. pylori is closely correlated with active chronic gastritis (P < 0.001).