Genome-wide profiling of the PIWI-interacting RNA-mRNA regulatory networks in epithelial ovarian cancers

PIWI-interacting (piRNAs), ~23-36 nucleotide-long small non-coding RNAs (sncRNAs), earlier believed to be germline-specific, have now been identified in somatic cells, including cancer cells. These sncRNAs impact critical biological processes by fine-tuning gene expression at post-transcriptional and epigenetic levels. The expression of piRNAs in ovarian cancer, the most lethal gynecologic cancer is largely uncharted. In this study, we investigated the expression of PIWILs by qRT-PCR and western blotting and then identified piRNA transcriptomes in tissues of normal ovary and two most prevalent epithelial ovarian cancer subtypes, serous and endometrioid by small RNA sequencing. We detected 219, 256 and 234 piRNAs in normal ovary, endometrioid and serous ovarian cancer samples respectively. We observed piRNAs are encoded from various genomic regions, among which introns harbor the majority of them. Surprisingly, piRNAs originated from different genomic contexts showed the varied level of conservations across vertebrates. The functional analysis of predicted targets of differentially expressed piRNAs revealed these could modulate key processes and pathways involved in ovarian oncogenesis. Our study provides the first comprehensive piRNA landscape in these samples and a useful resource for further functional studies to decipher new mechanistic views of piRNA-mediated gene regulatory networks affecting ovarian oncogenesis. The RNA-seq data is submitted to GEO database (GSE83794).