Abstract
The threshold of 200 nucleotides (nt) conventionally divides non-coding RNAs (ncRNA) into long ncRNA (lincRNA, that have more than 200 nt in length) and the remaining ones which are grouped as "small" RNAs (microRNAs, small nucleolar RNAs and piwiRNAs). Promoter-associated RNAs (paRNAs) are generally 200-500 nt long and are transcribed from sequences positioned in the promoter regions of genes. Growing evidence suggests that paRNAs play a crucial role in controlling gene transcription. Here, we used deep sequencing to identify paRNA sequences that show altered expression in a melanoma cell line compared to normal melanocytes. Thousands of reads were mapped to transcription start site (TSS) regions. We limited our search to paRNAs adjacent to genes with an expression that differed between melanoma and normal melanocytes and a length of 200-500 nt that did not overlap the gene mRNA by more than 300 nt, ultimately leaving us with 11 such transcripts. Using quantitative real-time PCR (qRT-PCR), we found a significant correlation between the expression of the mRNA and its corresponding paRNA for two studied genes: TYR and HSPC152. Ectopic overexpression of the paRNA of HSPC152 (designated paHSPC) enhanced the expression of the HSPC152 mRNA, and an siRNA targeting the paHSPC152 decreased the expression of the HSPC152 mRNA. Overexpression of paHSPC also affected the epigenetic structure of its putative promoter region along with effects on several biologic features of melanoma cells. The ectopic expression of the paRNA to TYR did not have any effect. Overall, our work indicates that paRNAs may serve as an additional layer in the regulation of gene expression in melanoma, thus meriting further investigation.