Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2018 Jan 10;13(1):e0190485.
doi: 10.1371/journal.pone.0190485. eCollection 2018.

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

Garima Singh  1 , Jyoti Roy  1 , Pratiti Rout  1 , Bibekanand Mallick  1
Affiliations

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

Garima Singh et al. PLoS One. .

Abstract

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

PubMed Disclaimer

Conflict of interest statement

Competing Interests: The authors have declared that no competing interests exist.

Figures

Fig 1
Fig 1. Expression profile of PIWIL genes and proteins in human normal ovarian and cancer tissues.
A. The relative expression of thee PIWIL mRNAs analyzed by qRT-PCR analysis. The amount of PIWIL mRNAs was normalized to the endogenous control, β-actin mRNA. The fold-change was calculated based on the ratio of the normalized values of the ENOCa and SOCa to that of normal ovary. B. Western blot of three PIWIL proteins in normal ovary, ENOCa and SOCa.
Fig 2
Fig 2. QC analysis of tissue samples of ENOCa, SOCa and normal ovary for small RNA sequencing and generation of reads.
Agilent RNA Bioanalyzer profile of A. ENOCa; B. SOCa; C. normal ovary; D. RIN values revealing small RNA intactness optimal for sequencing; E. Number of trimmed reads of 16–40 nts generated from each sample type.
Fig 3
Fig 3. Characteristic properties of piRNAs identified in ENOCa, SOCa and normal ovary.
A. Length and B. nucleotide bias observed among the piRNAs identified in each samples.
Fig 4
Fig 4. The chromosomal origin and genomic contexts of piRNAs in human genome.
A. Origin of piRNAs on all human chromosomes in three samples; B. piRNAs mapped to various genomic contexts of human genome.
Fig 5
Fig 5. The biogenetic sequence signatures of the intron derived piRNAs (piRtrons).
The preference of 1U & 10A are observed among piRtrons in ENOCa, SOCa and normal ovary.
Fig 6
Fig 6. The chromosomal distribution of piRNA clusters.
Location of piRNA clusters on human chromosomes in A. Normal ovary, B. ENOCa and C. SOCa samples.
Fig 7
Fig 7. Box plot analysis showing distribution of conservation scores of piRNAs originated from various genomic contexts in three samples (Normal ovary, ENOCa, SOCa).
Center lines show the medians; box limits indicate the 25th and 75th percentiles as determined by R program; whiskers extend to minimum and maximum values; crosses represent sample means; bars indicate 90% confidence intervals of the means; width of the boxes is proportional to the square root of the sample size.
Fig 8
Fig 8. piRNA-target duplexes in ENOCa and SOCa.
A. The binding sites of piR-52207 with its targets NUDT4, MTR, EIF2S3 and MPHOSPH8 in ENOCa; B. The binding sites of piR-33733 and piR-52207 with its targets LIAS and ACTR10, PLEKHA5 respectively in SOCa.
Fig 9
Fig 9. Validation of expression of piRNAs in ovarian cancer subtypes with respect to normal ovarian tissue by qRT-PCR.
Fig 10
Fig 10. Expression profile of piRNA targets by qRT-PCR in EOCa subtypes.
A. The relative expression of genes targeted by piR-52207 in ENOCa; B. The relative expression of genes targeted by piR-33733 and piR-52207) in SOCa.
Fig 11
Fig 11. The possible effects of piR-52207 on target genes and subsequent pathophysiological consequences in ENOCa.
Fig 12
Fig 12. The possible effects of piR-33733 and piR-52207 on target genes and subsequent pathophysiological consequences in SOCa.
See this image and copyright information in PMC

References

    1. Iorio MV, Visone R, Di Leva G, Donati V, Petrocca F, Casalini P, et al. MicroRNA signatures in human ovarian cancer. Cancer Res. 2007;67(18):8699–707. doi: 10.1158/0008-5472.CAN-07-1936 . - DOI - PubMed
    1. Devouassoux-Shisheboran M, Genestie C. Pathobiology of ovarian carcinomas. Chin J Cancer. 2015;34(1):50–5. doi: 10.5732/cjc.014.10273 - DOI - PMC - PubMed
    1. Sung PL, Chang YH, Chao KC, Chuang CM, Task Force on Systematic R, Meta-analysis of Ovarian C. Global distribution pattern of histological subtypes of epithelial ovarian cancer: a database analysis and systematic review. Gynecol Oncol. 2014;133(2):147–54. doi: 10.1016/j.ygyno.2014.02.016 . - DOI - PubMed
    1. Cunningham JM, Cicek MS, Larson NB, Davila J, Wang C, Larson MC, et al. Clinical characteristics of ovarian cancer classified by BRCA1, BRCA2, and RAD51C status. Scientific reports. 2014;4:4026 doi: 10.1038/srep04026 - DOI - PMC - PubMed
    1. Mosig RA, Lin L, Senturk E, Shah H, Huang F, Schlosshauer P, et al. Application of RNA-Seq transcriptome analysis: CD151 is an Invasion/Migration target in all stages of epithelial ovarian cancer. Journal of ovarian research. 2012;5:4 doi: 10.1186/1757-2215-5-4 - DOI - PMC - PubMed

Publication types

MeSH terms