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.

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.

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.

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.

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.

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.

Transcriptional profiling of Neurospora crassa Δmak-2 reveals that mitogen-activated protein kinase MAK-2 participates in the phosphate signaling pathway.

The filamentous fungus Neurospora crassa is an excellent model system for examining molecular responses to ambient signals in eukaryotic microorganisms. Inorganic phosphate (Pi) is an essential growth-limiting nutrient in nature and is crucial for the synthesis of nucleic acids and the flow of genetic information. The genetic and molecular mechanisms controlling the response to Pi starvation in N. crassa include at least four genes (nuc-2, preg, pogv, and nuc-1), which are involved in a hierarchical regulatory activation network. In a previous work, we identified a number of genes modulated by NUC-2 protein, including the mak-2 gene, which codes for a mitogen-activated protein kinase (MAPK), suggesting its participation in the phosphate signaling pathway. Thus, to identify other genes involved in metabolic responses to exogenous phosphate sensing and the functioning of the MAPK MAK-2, we performed microarray experiments using a mak-2 knockout strain (Δmak-2) grown under phosphate-shortage conditions by comparing its transcription profile to that of a control strain grown in low- and high-phosphate cultures. These experiments revealed 912 unique differentially expressed genes involved in a number of physiological processes related to phosphate transport, metabolism, and regulation as well as posttranslational modification of proteins, and MAPK signaling pathways. Quantitative Real-time PCR gene expression analysis of 18 selected genes, using independent RNA samples, validated our microarray results. A high Pearson correlation between microarray and quantitative Real-time PCR data was observed. The analysis of these differentially expressed genes in the Δmak-2 strain provide evidence that the mak-2 gene participates in the hierarchical phosphate-signaling pathway in N. crassa in addition to its involvement in other metabolic routes such as the isoprenylation pathway, thus revealing novel aspects of the N. crassa phosphorus-sensing network.

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.

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.

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 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.

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.

Conserved tissue expression signatures of intronic noncoding RNAs transcribed from human and mouse loci.

It has been postulated that noncoding RNAs (ncRNAs) are involved in the posttranscriptional control of gene expression, and may have contributed to the emergence of the complex attributes observed in mammalians. We show here that the complement of ncRNAs expressed from intronic regions of the human and mouse genomes comprises at least 78,147 and 39,660 transcriptional units, respectively. To identify conserved intronic sequences expressed in both humans and mice, we used custom-designed human cDNA microarrays to separately interrogate RNA from mouse and human liver, kidney, and prostate tissues. An overlapping tissue expression signature was detected for both species, comprising 198 transcripts; among these, 22 RNAs map to intronic regions with evidence of evolutionary conservation in humans and mice. Transcription of selected human-mouse intronic ncRNAs was confirmed using strand-specific RT-PCR. Altogether, these results support an evolutionarily conserved role of intronic ncRNAs in human and mouse, which are likely to be involved in the fine tuning of gene expression regulation in different mammalian tissues.

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.

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.

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.

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.

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.

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.