Next Article in Journal
A Study of Polish Family with Scoliosis and Limb Contractures Expands the MYH3 Disease Spectrum
Next Article in Special Issue
Human Reproduction and Disturbed Genomic Imprinting
Previous Article in Journal
Exploring the Therapeutic Significance of microRNAs and lncRNAs in Kidney Diseases
Previous Article in Special Issue
Genetic Variants of Gonadotropins and Their Receptors Could Influence Controlled Ovarian Stimulation: IVF Data from a Prospective Multicenter Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Circular RNA Derived from the Pumilio 1 Gene Could Regulate PTEN in Human Cumulus Cells

1
Department of Biomedical and Biotechnological Sciences, Section of Biology and Genetics “G. Sichel”, University of Catania, 95123 Catania, Italy
2
Department of Physics and Astronomy “Ettore Majorana”, University of Catania, 95123 Catania, Italy
3
Department of General Surgery and Medical-Surgical Specialties, University of Catania, 95123 Catania, Italy
4
IVF Unit, Cannizzaro Hospital, 95123 Catania, Italy
5
Department of Medicine and Surgery, Kore University, 94100 Enna, Italy
6
Maternal and Child Department, Obstetrics and Gynecology Unit, Cannizzaro Hospital, 95123 Catania, Italy
*
Author to whom correspondence should be addressed.
Genes 2024, 15(1), 124; https://doi.org/10.3390/genes15010124
Submission received: 30 November 2023 / Revised: 17 January 2024 / Accepted: 17 January 2024 / Published: 19 January 2024
(This article belongs to the Special Issue Genetics and Genomics of Female Reproduction)

Abstract

:
CircRNAs are a class of non-coding RNAs able to regulate gene expression at multiple levels. Their involvement in physiological processes, as well as their altered regulation in different human diseases, both tumoral and non-tumoral, is well documented. However, little is known about their involvement in female reproduction. This study aims to identify circRNAs potentially involved in reproductive women’s health. Candidate circRNAs expressed in ovary and sponging miRNAs, already known to be expressed in the ovary, were selected by a computational approach. Using real time PCR, we verified their expression and identified circPUM1 as the most interesting candidate circRNA for further analyses. We assessed the expression of circPUM1 and its linear counterpart in all the follicle compartments and, using a computational and experimental approach, identified circPUM1 direct and indirect targets, miRNAs and mRNAs, respectively, in cumulus cells. We found that both circPUM1 and its mRNA host gene are co-expressed in all the follicle compartments and proposed circPUM1 as a potential regulator of PTEN, finding a strong positive correlation between circPUM1 and PTEN mRNA. These results suggest a possible regulation of PTEN by circPUM1 in cumulus cells and point out the important role of circRNA inside the pathways related to follicle growth and oocyte maturation.

1. Introduction

In the last few decades, there has been a spread of large-scale genome sequencing studies that have revealed that only 2% of the mammalian genome produces translated RNAs (mRNAs), while the remaining part is transcribed into non-coding protein RNAs (ncRNAs) [1]. Originally, this portion of the human genome was considered a by-product of massive transcription with less biological meaning. However, today, it is commonly thought that ncRNAs have many important regulatory functions [2]. They participate in multiple biological pathways, regulating physiological and pathological processes [3,4,5,6,7]. NcRNAs may be classified according to their function into the following: housekeeping ncRNAs (rRNAs, tRNAs, small nuclear RNAs, small nucleolar RNAs, and telomerase RNAs) and regulatory ncRNAs. The last ones, depending on their lengths, are further categorized into small non-coding RNAs (ranging from 50 to 200 nt) and long non-coding RNAs (more than 200 nt). A particular class of regulatory ncRNA is the circular RNAs (circRNAs) [8]. CircRNAs are characterized by a peculiar circular shape, lacking free terminal structures, derived from non-canonical splicing events of linear pre-mRNA (back-splicing) that generate 3′ and 5′ ends covalently joined together. This peculiar round structure promotes their stability and protects them from degradation using RNA exonucleases or RNase R [9]. Usually, circRNAs are expressed at low levels and show cell type- or tissue-specific expression, indicating their potential regulatory roles [10]. Their ability to sponge microRNAs (miRNAs), acting within the competing endogenous RNA (ceRNA) network, is relatively well known. Based on this system, a specific circRNA can impair miRNA activity, influencing the expression of mRNA miRNA targets [11,12]. Moreover, circRNAs can regulate the expression of their host genes using the interaction with RNA polymerase II; they can interact with RNA-binding proteins, and some of them can be translated into proteins. Their involvement in physiological processes as well as their altered regulation in different human diseases is well documented [13,14,15,16].
Cumulus cells (CCs) are a subset of granulosa cells (GC) surrounding oocytes. The existence of an oocyte–cumulus cell regulatory loop is necessary for follicle development and the acquisition of oocyte competence.
Various papers have shown the importance of regulatory networks based on ncRNAs in human oocytes and CCs and, in particular, the involvement of circRNAs in regulating molecular functions related to human reproduction [17,18,19,20,21]. CircRNAs have been characterized by high-throughput sequencing in human CCs, and they have been proposed to be potentially involved in endometriosis-associated infertility [22] and in polycystic ovary syndrome (PCOS) [20,23,24].
In animal models, circRNA expression has been well characterized in the ovary (GCs, CCs, and oocytes), in the testes, placenta, and embryo, and also in germ line stem cells, embryonic stem cells, and induced pluripotent stem cells [25]. However, only in a few cases has the regulatory role of circRNAs been functionally demonstrated [26]. Cao Z. and colleagues demonstrated for the first time in 2019 that circRNAs are dynamically and highly expressed in a developmental stage-specific manner in CCs and oocytes and that maternally expressed circRNAs are essential for porcine oocyte meiotic maturation and early embryo development [25].
In this study, using different approaches, we proved the expression of the little-known circPUM1 specific isoform in the different compartments of a mature human ovarian follicle (Figure 1). This circRNA arises from the PUM1 gene, a member of the Pumilio RNA-binding protein family. PUM family proteins are ubiquitously present across eukaryotes, including yeast, plants, and humans, and perform their regulatory role at post-transcriptional/translational levels [27]. The high evolutionary conservation of the RNA-binding domain and functions suggest that human proteins play a crucial role in the translational regulations of embryogenesis, cell development, and differentiation [27]. We verified the co-expression of circPUM1 with the mRNA of its host gene and, using a computational and experimental approach, identified circPUM1 as a potential regulator of PTEN, finding a strong positive correlation between circPUM1 and PTEN mRNA. These results suggest the possible regulation of PTEN by circPUM1 in CCs and point out the important role of circRNA in the female reproduction system.

2. Materials and Methods

2.1. Selection of circRNAs Expressed in the Ovary

In order to pinpoint the circRNAs performing a significant role inside the ovarian follicle and to understand their biological functions, we intersected different public databases, performing a bioinformatic analysis, as described below and in Figure 1.
Using the Tissue-Specific CircRNA Database (TSCD) “http://gb.whu.edu.cn/TSCD (accessed on 10 February 2023)”, we obtained a list of circRNAs specifically expressed in the ovary. This list of circRNAs was matched to high-throughput sequencing analysis of circRNAs in normal human ovaries by Hongcai Cai et al., 2018 [28]. We considered the common circRNAs and selected only the circRNAs annotated in the circBase database “http://www.circbase.org/ (accessed on 12 February 2023)”. In order to detect circRNA–miRNA interactions significant in the ovary, in parallel, we compared the list of circRNAs expressed in the ovary (Figure 1) with the list of circRNAs sponging miRNAs expressed in the ovary [29]. This last list was obtained by ENCORI (Encyclopedia of RNA Interactomes) “https://rnasysu.com/encori/ (accessed on 12 February 2023)”, which identified the interactions, supported by CLIP-seq data, between miRNAs already known to be expressed in the ovary [29] and all human-annotated circRNAs.

2.2. Patient Enrollment

Participants enrolled in this study were women having undergone ovarian stimulation protocol and intracytoplasmatic sperm injection (ICSI) at the IVF Center Cannizzaro Hospital Catania (Italy). Healthy women whose primary infertility was caused by male factors were included in the study; those with pathologies that could influence oocyte quality, such as endometriosis, polycystic ovary, and ovarian insufficiency, were excluded from the study, as well as heavy smokers and overweight women.
Our research followed the tenets of the Declaration of Helsinki. The patients gave their informed consent to use their CCs and follicular fluid (FF). Women undergoing infertility treatment and already having a successful pregnancy donated the supernumerary oocytes. These oocytes were vitrified after their collection and thawed for our research application.

2.3. Sample Collection

CCs, FFs, and oocytes were independently collected. CC samples were collected according to the protocol reported by Caponnetto et al., 2022 [21]. Three to seven CC samples per patient were collected from the oocytes mechanically using needles and set aside individually. The total removal of the cumulus and corona cells continued using both enzymatic (hyaluronidase) and mechanical (pipetting) methods in order to identify oocytes in metaphase II (MII). Only the cumulus cells from the MII oocytes were washed in the culture medium and frozen at −20 °C.
The collection of FF samples was performed as reported by Caponnetto et al., 2021 [5]. Briefly, in order to remove residual follicular cells and any blood traces, FF samples were centrifuged for 20 min at 1000× g at 4 °C; the supernatant was collected and stored at −20 °C until use. Only FF samples without any massive blood contaminations and in which nuclear mature oocytes (MII) were identified were included in the study.
At the end of the injection, the oocyte supernumeraries for the fertilization were vitrified. After the decision of the couples to donate surplus oocytes for research, 3 oocytes coming from 2 different patients were thawed and used for molecular applications. Commercially available total RNAs from physiological brain, heart, and liver tissues were obtained from Ambion (Austin, TX, USA).

2.4. RNA Isolation and Purification

Total RNA isolation and purification from CC samples were performed by TRIzol® (Thermo Fisher Scientific, Waltham, MA, USA), according to the manufacturer’s instructions. Total RNA from FF samples was isolated as follows: RNA was isolated separately from two aliquots of 200 μL of the same FF sample by using the Qiagen miRNeasy Mini Kit (Qiagen, GmbH, Hilden, Germany), according to Qiagen Supplementary Protocol, for the purification of the RNA (including small RNAs) from serum or plasma. The two aliquots of the same FF sample were mixed together in the final step of the procedure. RNA precipitation protocol was performed in order to increase the final yield of RNA. Briefly, RNA was eluted in 200 μL of RNAse-free water; combined with 20 μg of UltraPure Glycogen (Thermo Fisher Scientific, Waltham, MA, USA), 0.1 volumes of 3 M sodium acetate, and 2.5 volumes of ice-cold absolute ethanol; and incubated at –80 °C overnight. RNA was then washed three times in ice-cold 75% ethanol, resuspended in 7 μL of RNAse-free water, and stored at –80 °C for further analyses.
Total RNAs from CCs and FFs were quantified by a spectrophotometer. Samples showing a 260–280 ratio absorbance of approximately 2.0 were considered.
According to previously published protocols, RNA from human oocytes was extracted by thermolysis; the samples were incubated for 1 min at 100 °C in order to release nucleic acids [30,31].

2.5. Gene Expression Analysis

In order to amplify circRNAs and mRNAs, primer design was performed by using the NCBI primer blast tool “https://www.ncbi.nlm.nih.gov/tools/primer-blast/ (accessed on 30 March 2023)”. For the detection of circular isoforms, divergent primers were designed; convergent primers were used for the detection of linear mRNAs (Table 1). Primer pairs were tested in silico, and all those recognizing more than one circular isoform produced by the same host gene were discarded. Only in the case where it was not possible to obtain primer pairs for the specific circRNA isoform, two primer pairs were designed and tested. All specific circRNA primer pairs were tested in human CCs, as well as the expression of the circRNAs selected.
Total RNA was amplified in a one-step reaction by the RNA-to-Ct™ 1-Step Kit (Thermo Fisher Scientific, Waltham, MA, USA) using convergent and divergent primers already mentioned and reported in Table 1.
MiRNA gene expression analysis was performed by using single TaqMan microRNA assays. Total RNA was reverse transcribed through the TaqMan™ microRNA Reverse Transcription Kit (Thermo Fisher Scientific, Waltham, MA, USA) using miRNA-specific primers and then amplified by the TaqMan™ Universal Master Mix II (Thermo Fisher Scientific, Waltham, MA, USA) by using specific TaqMan™ assays, according to the manufacturer’s instructions.
All quantitative RT-qPCR reactions were performed on the QuantStudio™ 7 Flex Real-Time PCR system (Thermo Fisher Scientific), according to the supplier’s protocol. GAPDH was used as the endogenous control gene for the normalization of linear and circular transcripts [32,33,34,35]. RNU6 was used as the endogenous control for miRNA gene expression data normalization.

2.6. Prediction of circRNA–miRNA–mRNA Interactions

After the identification of possible miRNAs interacting with the circPUM1 isoform (as reported in Section 2.1), we identified the mRNAs targeted by these miRNAs using the miRTarbase 9.0 beta database “https://miRTarBase.cuhk.edu.cn/ (accessed on 20 March 2023)” [36]. We queried the database by selecting only the validated bindings assessed by the reporter luciferase assay experiment.

2.7. Statistical Analysis

Pearson’s correlation coefficient (r) was used to evaluate the correlation between the expression of circPUM1 and its linear counterpart and circPUM1 and its direct and indirect interactors. Statistical significance was assessed by setting a p-value ≤ 0.05 as the cut-off.

2.8. Genomics of circPUM1

The sequences of human circPUM1 were retrieved from circBase “http://www.circbase.org/ (accessed on 1 April 2023)”. The chromosome, the nucleotide position, the genomic and spliced sequence lengths, and the gene symbol were determined by the Circular RNA Interactome database “https://circinteractome.nia.nih.gov/index.html (accessed on 1 April 2023)”.

3. Results

3.1. CircRNA Selection

The comparison between the list of circRNAs specifically expressed in ovaries from TSCD and the circRNAs expressed in normal ovaries reported by Hongcai Cai et al., 2018 [28] resulted in a list of 545 common circRNAs, 106 of which were annotated in the circBase database. Among them, seven circRNAs were potentially able to interact with miRNAs highly expressed in the ovary [29]. Table 2 shows the selected circRNAs, their host genes, and the miRNAs interacting with them.

3.2. CircRNA Expression

The amplification of the seven circRNAs revealed that (i) all circRNAs were expressed in human CC samples and showed Ct values ranging from 26 to 39 (Figure S1A), and (ii) amplified products showed only a single peak and non-specific amplification or primer dimer formation (Figure S1B). Only hsa_circ_0077150, hsa_circ_0004079 (1), hsa_circ_0004079 (2), and hsa_circ_0011233 showed Ct values ≤ 33 (Figure S1A,C). circPUM1 (hsa_circ_0011233) was selected for further analyses, since it showed a lower Ct value and was the most efficiently amplified (Figure S1).

3.3. Genomics of circPUM1

circPUM1 (hsa_circ_0011233) originates from the PUM1 host gene (Pumilio RNA-binding family member 1). This gene encodes a member of the PUF family, evolutionarily conserved RNA-binding proteins related to the Pumilio proteins of Drosophila, and the fem-3 mRNA-binding factor proteins of Caenorhabditis elegans. PUM1 (chr1: 30931506-31065717 GRCh38/hg38) is composed of 22 exons, and two spliced transcript variants encoding different isoforms have been described: NM_014676.3 and NM_001020658.2. PUM1 is the host gene of 62 circRNA isoforms, as annotated in the circBase database [8,10,37,38,39]. Among these, our isoform of circPUM1 is an exonic circRNA, 527 nt long, originating from the 15th, 16th, and 17th exons of the NM_0014676.3 transcript by a back-splicing process (Table 3 and Figure S2).

3.4. CircPUM1 and Its Linear Counterpart Are Expressed in the Follicular Microenvironment

CircPUM1 and its linear counterpart are expressed in the follicular microenvironment, as well as in three different physiological tissues: the brain, heart, and liver. A similar expression trend, reported as the ratio circPUM1/PUM1, was observed in CC samples, oocytes, and the other tissues (Figure 2A,C,D); otherwise, FF samples showed a higher variability of circPUM1 and PUM1 expressions (Figure 2B). A significant and positive correlation between circPUM1 and PUM1 linear transcript in CCs (Pearson’s correlation coefficient: r = 0.84, p = 0.002) (Figure 3A) was observed. The correlation was stronger when all samples were considered together (r = 0.94, and p < 0.0001) (Figure 3B).

3.5. CircPUM1 and Its RNA Interactors in CC Samples

The computational analysis revealed that miR-16-5p and miR-214-3p were potentially sponged by circPUM1. We experimentally assessed that miR-16-5p and miR-214-3p were expressed in our CC sample cohort and identified their validated targets using the miRTarbase. As expected, the miRTarbase returned several mRNA targets for miR16-5p and miR-214-3p, among which we selected BCL2, PTEN, and TP53. In particular, BCL2 was selected, since its interaction with miR-16-5p was demonstrated by both strong and less-strong evidences and was assessed by several studies reported in literature. For the same reasons, PTEN was chosen as the interactor of miR-214-3p. Moreover, the proteins encoded by these mRNAs, as well as the TP53 protein, perform important and well-known functions in folliculogenesis [36].
Among these, we focused our attention on BCL2, PTEN, and TP53. In particular, both miRNAs shared two of them, BCL2 and TP53, while PTEN is only targeted by miR-214-3p. All the mRNA targets selected were found to be expressed in our studied CC cohort (Figure 4).

3.6. Correlation between circPUM1 and mRNA Expression Levels in CCs

In our CC samples, we found a strong positive correlation between circPUM1 and PTEN (r = 0.92, p = 0.004) (Figure 5 and Figure 6A) and also between miR-214-3p and PTEN (r = 0.86, p = 0.016). A significant positive correlation between circPUM1 and miR-16-5p was observed (r = 0.68, p = 0.04) (Figure 5 and Figure 6B). However, there was no significant correlation between circPUM1 and miR-214-3p expressions in CCs (Figure 5 and Figure 6C).

4. Discussion

In recent years, the discovery of novel ncRNAs has been greatly increased by high-throughput sequencing. Advances in sequencing technology and library preparation protocols, as well as computational biology have increased our knowledge about the presence of ncRNAs in prokaryotes and eukaryotes [40]. RNA sequencing provides researchers with insights into transcriptome in a cell-specific or tissue-specific manner and in different states of cellular differentiation [40,41]. However, in order to further investigate the biological and functional roles of RNA molecules, experimental validation using biochemical techniques and functional assays are needed, revealing the crucial roles of many ncRNAs in human health and disease [40,42].
Our study, using computational and experimental approaches, allowed us to identify a little-known circPUM1 isoform that could potentially play an important regulatory role inside ovarian follicles, controlling follicle maturation and meiosis resumption (Figure 1).
Several circPUM1 isoforms have already been identified in different cellular models, and their roles in the regulation of gene expression inside specific cellular processes have been demonstrated [43,44,45,46,47,48,49]. However, to date, little is known about the hsa_circ_0011233 isoform identified in this study.
The host gene, PUM1, is a member of the Pumilio RNA-binding protein family, evolutionarily conserved across species. In the human genome, we found two PUM genes closely related to each other, PUM1 and PUM2, spanning about 150 kb (on chromosome 1p35.2) and 80 kb (on chromosome 2p23–24), respectively. The indistinguishable RNA-binding specificities of PUM1 and PUM2 suggest potential overlapping functions [50]. PUM1 regulates several biological processes, such as embryonic development, functioning and development of the nervous system, rRNA processing, ribosome biogenesis, chemotactic cell movement, and also stem cell and germ cell maintenance [27]. The maintenance and self-renewal of stem cells might be an ancestral function of Pum, as demonstrated in Drosophila [51], Plasmodium falciparum [52], and zebrafish [53], as well as other model organisms. In mice, it has been demonstrated that maternal PUM1 is an essential post-transcriptional regulator during mammalian embryogenesis, regulating the stability of maternal mRNAs. PUM1 has been found to be expressed in oocytes of different organisms, such as Xenopus [54] and zebrafish [55], and it is involved in the control of meiosis. One of the mechanisms of action could be the post-transcriptional regulation of the mRNA of Cdk1, the kinase that dimerizes with cyclin B1, forming the metaphase-promoting factor (MPF) [56]. Moreover, knock-out experiments have revealed that PUM1 is essential in the mammalian germline, such as in the establishment of the primordial follicle pool, meiosis, and female reproductive competency, while PUM2 does not have a detectable function in these processes [57].
Here, we have demonstrated that circPUM1 is expressed in the ovarian follicle microenvironment, specifically in CCs and oocytes and also in FF samples (Figure 2). We found a strong and positive correlation between circPUM1 and its linear isoform (PUM1 mRNA), not only in the ovary but also in different tissues (Figure 3). Indeed, in FF, we found a higher variability between our circRNA and its linear counterpart, which could be expected because in biological fluids, the expression levels of RNA transcripts are affected by the secretion mechanisms and do not depend exclusively on their transcriptional levels inside the cells. However, the inclusion of FF samples in the correlation analysis, shown in the Figure 3B, did not negatively affect the result.
These data suggest that our circRNA isoform could perform similar functions as the linear counterpart and could play synergistic roles in regulating specific cellular processes involved in follicle growth and oocyte maturation. The synergic roles of circRNAs and their host genes have already been described in literature [58].
In addition, the correlation with their linear counterparts could be dynamically regulated, and the regulation of circRNAs on their corresponding host genes is a critical mechanism for their function [58]. In this regard, it is important to note that in some regulatory circuits, circRNAs and their host genes can be part of feedback loops. Changes in the expression of one component (circRNA or host gene) can influence the expression of the other, contributing to dynamic regulatory mechanisms [58]. Evaluating the expressions of PUM1 mRNA and the circRNA could be interesting under specific conditions that impair oocyte quality.
Another important characteristic supporting the potential role of circPUM1 inside the ovarian follicle is the strong positive correlation with PTEN (phosphatase and tensin homolog) mRNA (Figure 5 and Figure 6A). PTEN is a well-known tumor suppressor gene that plays a significant role in ovarian follicle development, GC function, and overall female reproductive health. PTEN is an indispensable molecule that maintains the dormancy of the primordial follicle pool, inhibiting the PI3K signaling pathway that controls primordial follicle survival and activation; regulates cyclic follicular recruitment; causes ovulation in GCs; and stimulates meiosis resumption in the oocyte [59,60]. PTEN down-regulation is an important step in the cyclical activation of primordial follicles. The depletion of the primordial follicle pool leads to the end of female reproductive life, and it has been demonstrated that in mice lacking PTEN, the entire primordial follicle pool becomes activated, causing premature ovarian failure [61,62]. Dysregulation of PTEN in the ovaries can also have implications for fertility, ovarian disorders, and even ovarian cancer [63]. Recently, it has been demonstrated that in human GCs, PTEN expression is associated with In Vitro Fertilization (IVF) outcomes [64].
The relationship between circPUM1 and PTEN implies common mechanisms of regulation that could be mediated by miRNAs. However, although we found a strong and positive correlation between circPUM1 and PTEN, we did not observe any correlation between circPUM1 and miR-214-3p (Figure 5 and Figure 6C). It is important to note that the interaction between miRNAs and circRNAs within the ceRNA regulatory mechanism decreases miRNA availability, promoting target gene expression without necessarily changing miRNA levels [65]. However, the regulation of PTEN by circPUM1 could be due to a different network of regulation not mediated by miR-214-3p. Different miRNAs or other ncRNAs could be involved in this mechanism of regulation.
Unexpectedly, we found a strong and positive correlation between miR-214-3p and PTEN (Figure 5). The positive correlation between a miRNA and its target genes has been described in several cancer models, and specifically, a positive correlation between PTEN and miR-214 has been found in breast cancer [66,67]. Recent studies indicate that miRNAs could also be implicated in the positive regulation of gene transcription by a phenomenon known as RNA activation (RNAa) [68].
In addition, we did not find any correlations between circPUM1 and the validated targets of miR-16-5p, BCL2 and TP53 (Figure 5), although we found a strong and positive correlation between circPUM1 and miR-16-5p (Figure 5 and Figure 6B). These data could suggest that circPUM1 does not play any regulatory role in the expressions of BCL2 and TP53, even if it could share regulatory mechanisms with miR-16-5p. However, miR-16-5p could exert its function inhibiting the translation of the targets without any effect on transcript degradation.
It is important to underline that our data were obtained in a physiological model and allowed us to identify a little-known cirRNA isoform, derived from the PUM1 gene and able to regulate PTEN. We propose that this study, suggesting circPUM1 as a possible regulator of PTEN, represents an important starting point for further investigations in order to understand the complex molecular mechanisms involved in reproductive women’s health. Further studies investigating the relative expressions of PUM1 and its interactors or regulators in female disorders will be needed.

5. Conclusions

The field of circRNA research in female reproduction is still relatively young, and ongoing studies continue to provide new insights into the specific roles and mechanisms of circRNAs in these processes. The complexity of circRNA functions and their interactions with other biomolecules make this area of research both exciting and challenging for scientists in the field of reproductive biology.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes15010124/s1, Figure S1: Amplification plot of the selected circRNAs (A). Report of Melting curve plot (B). Electrophoretic analysis of the circRNA products obtained by real-time PCR amplification (C). Figure S2: Representative illustration of PUM1 pre-mRNA and circular transcript (circPUM1) derived from back-splicing 15th, 16th, and 17th exons.

Author Contributions

Conceptualization, C.D.P., R.B. and A.C.; methodology, A.C., D.B., N.A. and M.S. (Marianna Scribano); validation, A.C., C.F., A.F., C.B. and M.S. (Michele Stella); formal analysis, A.C., C.F. and A.F.; investigation, A.C., C.F. and A.F.; resources, M.E.V., P.B. and P.S.; writing—original draft preparation, A.C., C.D.P. and R.B.; writing—review and editing, M.R., D.B., C.D.P., R.B. and A.C.; visualization, A.C., C.F., A.F., N.A., M.S. (Marianna Scribano), M.E.V., P.B., C.B., M.S. (Michele Stella), M.R., P.S., D.B., M.P., C.D.P. and R.B.; supervision, C.D.P. and R.B.; funding acquisition, C.D.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Acknowledgments

The authors would like to thank the Scientific Bureau of the University of Catania for language support and Fondi di Ateneo 2020–2022, University of Catania, Open Access line.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Palazzo, A.F.; Lee, E.S. Non-coding RNA: What is functional and what is junk? Front. Genet. 2015, 6, 2. [Google Scholar] [CrossRef] [PubMed]
  2. Elgar, G.; Vavouri, T. Tuning in to the signals: Noncoding sequence conservation in vertebrate genomes. Trends Genet. 2008, 24, 344–352. [Google Scholar] [CrossRef] [PubMed]
  3. Calore, F.; Lovat, F.; Garofalo, M. Non-coding RNAs and cancer. Int. J. Mol. Sci. 2013, 14, 17085–17110. [Google Scholar] [CrossRef] [PubMed]
  4. Arraiano, C.M. Regulatory noncoding RNAs: Functions and applications in health and disease. FEBS J. 2021, 288, 6308–6309. [Google Scholar] [CrossRef] [PubMed]
  5. Caponnetto, A.; Battaglia, R.; Ragusa, M.; Barbagallo, D.; Lunelio, F.; Borzi, P.; Scollo, P.; Purrello, M.; Vento, M.E.; Di Pietro, C. Molecular profiling of follicular fluid microRNAs in young women affected by Hodgkin lymphoma. Reprod. Biomed. Online 2021, 43, 1045–1056. [Google Scholar] [CrossRef] [PubMed]
  6. Battaglia, R.; Caponnetto, A.; Caringella, A.M.; Cortone, A.; Ferrara, C.; Smirni, S.; Iannitti, R.; Purrello, M.; D’Amato, G.; Fioretti, B.; et al. Resveratrol Treatment Induces Mito-miRNome Modification in Follicular Fluid from Aged Women with a Poor Prognosis for In Vitro Fertilization Cycles. Antioxidants 2022, 11, 1019. [Google Scholar] [CrossRef] [PubMed]
  7. Battaglia, R.; Alonzo, R.; Pennisi, C.; Caponnetto, A.; Ferrara, C.; Stella, M.; Barbagallo, C.; Barbagallo, D.; Ragusa, M.; Purrello, M.; et al. MicroRNA-Mediated Regulation of the Virus Cycle and Pathogenesis in the SARS-CoV-2 Disease. Int. J. Mol. Sci. 2021, 22, 13192. [Google Scholar] [CrossRef]
  8. Jeck, W.R.; Sorrentino, J.A.; Wang, K.; Slevin, M.K.; Burd, C.E.; Liu, J.; Marzluff, W.F.; Sharpless, N.E. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA 2013, 19, 141–157. [Google Scholar] [CrossRef]
  9. Kristensen, L.S.; Andersen, M.S.; Stagsted, L.V.W.; Ebbesen, K.K.; Hansen, T.B.; Kjems, J. The biogenesis, biology and characterization of circular RNAs. Nat. Rev. Genet. 2019, 20, 675–691. [Google Scholar] [CrossRef]
  10. Rybak-Wolf, A.; Stottmeister, C.; Glazar, P.; Jens, M.; Pino, N.; Giusti, S.; Hanan, M.; Behm, M.; Bartok, O.; Ashwal-Fluss, R.; et al. Circular RNAs in the Mammalian Brain Are Highly Abundant, Conserved, and Dynamically Expressed. Mol. Cell 2015, 58, 870–885. [Google Scholar] [CrossRef]
  11. Taulli, R.; Loretelli, C.; Pandolfi, P.P. From pseudo-ceRNAs to circ-ceRNAs: A tale of cross-talk and competition. Nat. Struct. Mol. Biol. 2013, 20, 541–543. [Google Scholar] [CrossRef] [PubMed]
  12. Thomson, D.W.; Dinger, M.E. Endogenous microRNA sponges: Evidence and controversy. Nat. Rev. Genet. 2016, 17, 272–283. [Google Scholar] [CrossRef] [PubMed]
  13. Barbagallo, D.; Caponnetto, A.; Brex, D.; Mirabella, F.; Barbagallo, C.; Lauretta, G.; Morrone, A.; Certo, F.; Broggi, G.; Caltabiano, R.; et al. CircSMARCA5 Regulates VEGFA mRNA Splicing and Angiogenesis in Glioblastoma Multiforme Through the Binding of SRSF1. Cancers 2019, 11, 194. [Google Scholar] [CrossRef]
  14. Chen, X.; Zhou, M.; Yant, L.; Huang, C. Circular RNA in disease: Basic properties and biomedical relevance. Wiley Interdiscip. Rev. RNA 2022, 13, e1723. [Google Scholar] [CrossRef] [PubMed]
  15. Haque, S.; Harries, L.W. Circular RNAs (circRNAs) in Health and Disease. Genes 2017, 8, 353. [Google Scholar] [CrossRef] [PubMed]
  16. Merulla, A.E.; Stella, M.; Barbagallo, C.; Battaglia, R.; Caponnetto, A.; Broggi, G.; Altieri, R.; Certo, F.; Caltabiano, R.; Ragusa, M.; et al. CircSMARCA5 Is an Upstream Regulator of the Expression of miR-126-3p, miR-515-5p, and Their mRNA Targets, Insulin-like Growth Factor Binding Protein 2 (IGFBP2) and NRAS Proto-Oncogene, GTPase (NRAS) in Glioblastoma. Int. J. Mol. Sci. 2022, 23, 13676. [Google Scholar] [CrossRef]
  17. Xu, X.; Jia, S.Z.; Dai, Y.; Zhang, J.J.; Li, X.; Shi, J.; Leng, J.; Lang, J. The Relationship of Circular RNAs With Ovarian Endometriosis. Reprod. Sci. 2018, 25, 1292–1300. [Google Scholar] [CrossRef]
  18. Liu, L.; Li, L.; Ma, X.; Yue, F.; Wang, Y.; Wang, L.; Jin, P.; Zhang, X. Altered Circular RNA Expression in Patients with Repeated Implantation Failure. Cell Physiol. Biochem. 2017, 44, 303–313. [Google Scholar] [CrossRef]
  19. Dang, Y.; Yan, L.; Hu, B.; Fan, X.; Ren, Y.; Li, R.; Lian, Y.; Yan, J.; Li, Q.; Zhang, Y.; et al. Tracing the expression of circular RNAs in human pre-implantation embryos. Genome Biol. 2016, 17, 130. [Google Scholar] [CrossRef]
  20. Che, Q.; Liu, M.; Xu, J.; Liu, Y.; Cao, X.; Dong, X.; Liu, S. Characterization of circular RNA expression profiles in cumulus cells from patients with polycystic ovary syndrome. Fertil. Steril. 2019, 111, 1243–1251.e1241. [Google Scholar] [CrossRef]
  21. Caponnetto, A.; Battaglia, R.; Ferrara, C.; Vento, M.E.; Borzi, P.; Paradiso, M.; Scollo, P.; Purrello, M.; Longobardi, S.; D’Hooghe, T.; et al. Down-regulation of long non-coding RNAs in reproductive aging and analysis of the lncRNA-miRNA-mRNA networks in human cumulus cells. J. Assist. Reprod. Genet. 2022, 39, 919–931. [Google Scholar] [CrossRef] [PubMed]
  22. Wu, R.; Li, J.; Li, J.; Yan, X.; Zhou, W.; Ling, C.; Ye, Y.; Chen, Q.; Li, Y. Circular RNA expression profiling and bioinformatic analysis of cumulus cells in endometriosis infertility patients. Epigenomics 2020, 12, 2093–2108. [Google Scholar] [CrossRef] [PubMed]
  23. Ma, Z.; Zhao, H.; Zhang, Y.; Liu, X.; Hao, C. Novel circular RNA expression in the cumulus cells of patients with polycystic ovary syndrome. Arch. Gynecol. Obstet. 2019, 299, 1715–1725. [Google Scholar] [CrossRef]
  24. Wang, L.P.; Peng, X.Y.; Lv, X.Q.; Liu, L.; Li, X.L.; He, X.; Lv, F.; Pan, Y.; Wang, L.; Liu, K.F.; et al. High throughput circRNAs sequencing profile of follicle fluid exosomes of polycystic ovary syndrome patients. J. Cell. Physiol. 2019, 234, 15537–15547. [Google Scholar] [CrossRef]
  25. Cao, Z.; Gao, D.; Xu, T.; Zhang, L.; Tong, X.; Zhang, D.; Wang, Y.; Ning, W.; Qi, X.; Ma, Y.; et al. Circular RNA profiling in the oocyte and cumulus cells reveals that circARMC4 is essential for porcine oocyte maturation. Aging 2019, 11, 8015–8034. [Google Scholar] [CrossRef] [PubMed]
  26. Yu, C.Y.; Li, T.C.; Wu, Y.Y.; Yeh, C.H.; Chiang, W.; Chuang, C.Y.; Kuo, H.C. The circular RNA circBIRC6 participates in the molecular circuitry controlling human pluripotency. Nat. Commun. 2017, 8, 1149. [Google Scholar] [CrossRef] [PubMed]
  27. Nishanth, M.J.; Simon, B. Functions, mechanisms and regulation of Pumilio/Puf family RNA binding proteins: A comprehensive review. Mol. Biol. Rep. 2020, 47, 785–807. [Google Scholar] [CrossRef]
  28. Cai, H.; Li, Y.; Li, H.; Niringiyumukiza, J.D.; Zhang, M.; Chen, L.; Chen, G.; Xiang, W. Identification and characterization of human ovary-derived circular RNAs and their potential roles in ovarian aging. Aging 2018, 10, 2511–2534. [Google Scholar] [CrossRef]
  29. Battaglia, R.; Musumeci, P.; Ragusa, M.; Barbagallo, D.; Scalia, M.; Zimbone, M.; Lo Faro, J.M.; Borzi, P.; Scollo, P.; Purrello, M.; et al. Ovarian aging increases small extracellular vesicle CD81(+) release in human follicular fluid and influences miRNA profiles. Aging 2020, 12, 12324–12341. [Google Scholar] [CrossRef]
  30. Di Pietro, C.; Vento, M.; Ragusa, M.; Barbagallo, D.; Guglielmino, M.R.; Maniscalchi, T.; Duro, L.R.; Tomasello, L.; Majorana, A.; De Palma, A.; et al. Expression analysis of TFIID in single human oocytes: New potential molecular markers of oocyte quality. Reprod. Biomed. Online 2008, 17, 338–349. [Google Scholar] [CrossRef]
  31. Battaglia, R.; Vento, M.E.; Ragusa, M.; Barbagallo, D.; La Ferlita, A.; Di Emidio, G.; Borzi, P.; Artini, P.G.; Scollo, P.; Tatone, C.; et al. MicroRNAs Are Stored in Human MII Oocyte and Their Expression Profile Changes in Reproductive Aging. Biol. Reprod. 2016, 95, 131. [Google Scholar] [CrossRef] [PubMed]
  32. Nielsen, M.E.; Rasmussen, I.A.; Kristensen, S.G.; Christensen, S.T.; Mollgard, K.; Wreford Andersen, E.; Byskov, A.G.; Yding Andersen, C. In human granulosa cells from small antral follicles, androgen receptor mRNA and androgen levels in follicular fluid correlate with FSH receptor mRNA. Mol. Hum. Reprod. 2011, 17, 63–70. [Google Scholar] [CrossRef] [PubMed]
  33. Juusola, J.; Ballantyne, J. mRNA profiling for body fluid identification by multiplex quantitative RT-PCR. J. Forensic. Sci. 2007, 52, 1252–1262. [Google Scholar] [CrossRef] [PubMed]
  34. Borgbo, T.; Jeppesen, J.V.; Lindgren, I.; Lundberg Giwercman, Y.; Hansen, L.L.; Yding Andersen, C. Effect of the FSH receptor single nucleotide polymorphisms (FSHR 307/680) on the follicular fluid hormone profile and the granulosa cell gene expression in human small antral follicles. Mol. Hum. Reprod. 2015, 21, 255–261. [Google Scholar] [CrossRef]
  35. Jeppesen, J.V.; Nielsen, M.E.; Kristensen, S.G.; Yding Andersen, C. Concentration of activin A and follistatin in follicular fluid from human small antral follicles associated to gene expression of the corresponding granulosa cells. Mol. Cell Endocrinol. 2012, 356, 48–54. [Google Scholar] [CrossRef] [PubMed]
  36. Huang, H.Y.; Lin, Y.C.; Cui, S.; Huang, Y.; Tang, Y.; Xu, J.; Bao, J.; Li, Y.; Wen, J.; Zuo, H.; et al. miRTarBase update 2022: An informative resource for experimentally validated miRNA-target interactions. Nucleic Acids Res. 2022, 50, D222–D230. [Google Scholar] [CrossRef] [PubMed]
  37. Memczak, S.; Jens, M.; Elefsinioti, A.; Torti, F.; Krueger, J.; Rybak, A.; Maier, L.; Mackowiak, S.D.; Gregersen, L.H.; Munschauer, M.; et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature 2013, 495, 333–338. [Google Scholar] [CrossRef]
  38. Salzman, J.; Chen, R.E.; Olsen, M.N.; Wang, P.L.; Brown, P.O. Cell-type specific features of circular RNA expression. PLoS Genet. 2013, 9, e1003777. [Google Scholar] [CrossRef]
  39. Maass, P.G.; Glazar, P.; Memczak, S.; Dittmar, G.; Hollfinger, I.; Schreyer, L.; Sauer, A.V.; Toka, O.; Aiuti, A.; Luft, F.C.; et al. A map of human circular RNAs in clinically relevant tissues. J. Mol. Med. 2017, 95, 1179–1189. [Google Scholar] [CrossRef]
  40. Micheel, J.; Safrastyan, A.; Wollny, D. Advances in Non-Coding RNA Sequencing. Noncoding RNA 2021, 7, 70. [Google Scholar] [CrossRef]
  41. Ozsolak, F.; Milos, P.M. RNA sequencing: Advances, challenges and opportunities. Nat. Rev. Genet. 2011, 12, 87–98. [Google Scholar] [CrossRef] [PubMed]
  42. Zhu, J.; Chen, G.; Zhu, S.; Li, S.; Wen, Z.; Bin, L.; Zheng, Y.; Shi, L. Identification of Tissue-Specific Protein-Coding and Noncoding Transcripts across 14 Human Tissues Using RNA-seq. Sci. Rep. 2016, 6, 28400. [Google Scholar] [CrossRef] [PubMed]
  43. Guan, X.; Zong, Z.H.; Liu, Y.; Chen, S.; Wang, L.L.; Zhao, Y. circPUM1 Promotes Tumorigenesis and Progression of Ovarian Cancer by Sponging miR-615-5p and miR-6753-5p. Mol. Ther. Nucleic Acids 2019, 18, 882–892. [Google Scholar] [CrossRef] [PubMed]
  44. Zhu, L.; Shi, L.; Ye, W.; Li, S.; Liu, X.; Zhu, Z. Circular RNA PUM1 (CircPUM1) attenuates trophoblast cell dysfunction and inflammation in recurrent spontaneous abortion via the MicroRNA-30a-5p (miR-30a-5p)/JUNB axis. Bioengineered 2021, 12, 6878–6890. [Google Scholar] [CrossRef] [PubMed]
  45. Chen, J.; Xu, S.; Chen, S.; Zong, Z.; Han, X.; Zhao, Y.; Shang, H. CircPUM1 promotes the malignant behavior of lung adenocarcinoma by regulating miR-326. Biochem. Biophys. Res. Commun. 2019, 508, 844–849. [Google Scholar] [CrossRef]
  46. Gong, W.; Xu, J.; Wang, Y.; Min, Q.; Chen, X.; Zhang, W.; Chen, J.; Zhan, Q. Nuclear genome-derived circular RNA circPUM1 localizes in mitochondria and regulates oxidative phosphorylation in esophageal squamous cell carcinoma. Signal Transduct. Target Ther. 2022, 7, 40. [Google Scholar] [CrossRef]
  47. Zhang, Y.; Wang, D.; Zhu, T.; Yu, J.; Wu, X.; Lin, W.; Zhu, M.; Dai, Y.; Zhu, J. CircPUM1 promotes hepatocellular carcinoma progression through the miR-1208/MAP3K2 axis. J. Cell. Mol. Med. 2021, 25, 600–612. [Google Scholar] [CrossRef]
  48. Li, M.; Wang, Q.; Zhang, X.; Yan, N.; Li, X. CircPUM1 promotes cell growth and glycolysis in NSCLC via up-regulating METTL3 expression through miR-590-5p. Cell Cycle 2021, 20, 1279–1294. [Google Scholar] [CrossRef]
  49. Chen, Y.; Liu, Y.; Tao, Q.; Fan, Y.; Ma, C.; Li, D.; Huang, F.; Tang, D. circPUM1 Activates the PI3K/AKT Signaling Pathway by Sponging to Promote the Proliferation, Invasion and Glycolysis of Pancreatic Cancer. Curr. Pharm. Biotechnol. 2022, 23, 1405–1414. [Google Scholar] [CrossRef]
  50. Spassov, D.S.; Jurecic, R. Cloning and comparative sequence analysis of PUM1 and PUM2 genes, human members of the Pumilio family of RNA-binding proteins. Gene 2002, 299, 195–204. [Google Scholar] [CrossRef]
  51. Parisi, M.; Lin, H. The Drosophila pumilio gene encodes two functional protein isoforms that play multiple roles in germline development, gonadogenesis, oogenesis and embryogenesis. Genetics 1999, 153, 235–250. [Google Scholar] [CrossRef] [PubMed]
  52. Miao, J.; Li, J.; Fan, Q.; Li, X.; Li, X.; Cui, L. The Puf-family RNA-binding protein PfPuf2 regulates sexual development and sex differentiation in the malaria parasite Plasmodium falciparum. J. Cell Sci. 2010, 123, 1039–1049. [Google Scholar] [CrossRef] [PubMed]
  53. Kuo, M.W.; Wang, S.H.; Chang, J.C.; Chang, C.H.; Huang, L.J.; Lin, H.H.; Yu, A.L.; Li, W.H.; Yu, J. A novel puf-A gene predicted from evolutionary analysis is involved in the development of eyes and primordial germ-cells. PLoS ONE 2009, 4, e4980. [Google Scholar] [CrossRef] [PubMed]
  54. Ota, R.; Kotani, T.; Yamashita, M. Biochemical characterization of Pumilio1 and Pumilio2 in Xenopus oocytes. J. Biol. Chem. 2011, 286, 2853–2863. [Google Scholar] [CrossRef] [PubMed]
  55. Saitoh, A.; Takada, Y.; Horie, M.; Kotani, T. Pumilio1 phosphorylation precedes translational activation of its target mRNA in zebrafish oocytes. Zygote 2018, 26, 372–380. [Google Scholar] [CrossRef] [PubMed]
  56. Mak, W.; Xia, J.; Cheng, E.C.; Lowther, K.; Lin, H. A role of Pumilio 1 in mammalian oocyte maturation and maternal phase of embryogenesis. Cell Biosci. 2018, 8, 54. [Google Scholar] [CrossRef] [PubMed]
  57. Mak, W.; Fang, C.; Holden, T.; Dratver, M.B.; Lin, H. An Important Role of Pumilio 1 in Regulating the Development of the Mammalian Female Germline. Biol. Reprod. 2016, 94, 134. [Google Scholar] [CrossRef] [PubMed]
  58. Wei, J.; Li, M.; Xue, C.; Chen, S.; Zheng, L.; Deng, H.; Tang, F.; Li, G.; Xiong, W.; Zeng, Z.; et al. Understanding the roles and regulation patterns of circRNA on its host gene in tumorigenesis and tumor progression. J. Exp. Clin. Cancer Res. 2023, 42, 86. [Google Scholar] [CrossRef]
  59. Zheng, W.; Nagaraju, G.; Liu, Z.; Liu, K. Functional roles of the phosphatidylinositol 3-kinases (PI3Ks) signaling in the mammalian ovary. Mol. Cell Endocrinol. 2012, 356, 24–30. [Google Scholar] [CrossRef]
  60. Di Pietro, C. Exosome-mediated communication in the ovarian follicle. J. Assist. Reprod. Genet. 2016, 33, 303–311. [Google Scholar] [CrossRef]
  61. Yao, J.; Li, M.; Lin, L.; Li, Y.; Zhuang, J.; Huang, Y.; Peng, W.; Lian, J.; Huang, R.; Yang, X. PTEN expression in human cumulus cells is associated with embryo development competence. Zygote 2022, 30, 611–618. [Google Scholar] [CrossRef] [PubMed]
  62. Reddy, P.; Liu, L.; Adhikari, D.; Jagarlamudi, K.; Rajareddy, S.; Shen, Y.; Du, C.; Tang, W.; Hamalainen, T.; Peng, S.L.; et al. Oocyte-specific deletion of Pten causes premature activation of the primordial follicle pool. Science 2008, 319, 611–613. [Google Scholar] [CrossRef] [PubMed]
  63. De Felici, M.; Klinger, F.G. PI3K/PTEN/AKT Signaling Pathways in Germ Cell Development and Their Involvement in Germ Cell Tumors and Ovarian Dysfunctions. Int. J. Mol. Sci. 2021, 22, 9838. [Google Scholar] [CrossRef] [PubMed]
  64. Yao, J.; Huang, R.; Li, M.; Jiang, Y.; Wu, P.; Li, Y.; Peng, W.; Hua, C.; Huang, Y.; You, H.; et al. PTEN Expression in Human Granulosa Cells Is Associated with Ovarian Responses and Clinical Outcomes in IVF. Reprod. Sci. 2021, 28, 1910–1921. [Google Scholar] [CrossRef] [PubMed]
  65. Singh, S.; Sinha, T.; Panda, A.C. Regulation of microRNA by circular RNA. Wiley Interdiscip. Rev. RNA 2023, e1820. [Google Scholar] [CrossRef]
  66. Vasudevan, S.; Tong, Y.; Steitz, J.A. Switching from repression to activation: MicroRNAs can up-regulate translation. Science 2007, 318, 1931–1934. [Google Scholar] [CrossRef]
  67. Sadighparvara, S.; Targhazehb, N.; Karimian, A.; Shafiei-Irannejad, V.; Farsad-Akhtar, N.; Rafieian, S.; Salamati, A.; Bastami, M.; Samadi Kafil, H.; Yousefi, M.; et al. Downregulation of microRNA-214 and PTEN in tissue samples of patients with breast cancer. Meta Gene 2020, 24, 100668. [Google Scholar] [CrossRef]
  68. Xiao, M.; Li, J.; Li, W.; Wang, Y.; Wu, F.; Xi, Y.; Zhang, L.; Ding, C.; Luo, H.; Li, Y.; et al. MicroRNAs activate gene transcription epigenetically as an enhancer trigger. RNA Biol. 2017, 14, 1326–1334. [Google Scholar] [CrossRef]
Figure 1. Workflow edited by BioRender “https://www.biorender.com/ (accessed on 20 November 2023)”.
Figure 1. Workflow edited by BioRender “https://www.biorender.com/ (accessed on 20 November 2023)”.
Genes 15 00124 g001
Figure 2. circPUM1 and PUM1 mRNA expression in ovarian follicle compartments: CCs (A), follicular fluids (FF) (B), oocytes (C), and other tissues (D) represented as the ratio of −ΔCt values.
Figure 2. circPUM1 and PUM1 mRNA expression in ovarian follicle compartments: CCs (A), follicular fluids (FF) (B), oocytes (C), and other tissues (D) represented as the ratio of −ΔCt values.
Genes 15 00124 g002
Figure 3. Expression changes in circPUM1 correlate with the linear transcript of PUM1 in CC samples (A) and all considered samples (B).
Figure 3. Expression changes in circPUM1 correlate with the linear transcript of PUM1 in CC samples (A) and all considered samples (B).
Genes 15 00124 g003
Figure 4. Expressions of PTEN, TP53, and BCL2 mRNAs in CC samples. Data are shown as −ΔCt values. In the x-axis, the number of CC samples is reported.
Figure 4. Expressions of PTEN, TP53, and BCL2 mRNAs in CC samples. Data are shown as −ΔCt values. In the x-axis, the number of CC samples is reported.
Genes 15 00124 g004
Figure 5. Correlation matrix among circPUM1, miRNAs 16-5p and 214-3p, and their target genes. The white–blue color scale bar indicates r values. Statistically significant correlations are reported as “*” and “**” when p < 0.05, and p < 0.01, respectively.
Figure 5. Correlation matrix among circPUM1, miRNAs 16-5p and 214-3p, and their target genes. The white–blue color scale bar indicates r values. Statistically significant correlations are reported as “*” and “**” when p < 0.05, and p < 0.01, respectively.
Genes 15 00124 g005
Figure 6. Expression trend of PTEN (A), miR-16-5p (B), and miR-214-3p (C) with respect to circPUM1. Data are shown as −ΔCt values. In the x-axis, the number of CC samples is reported.
Figure 6. Expression trend of PTEN (A), miR-16-5p (B), and miR-214-3p (C) with respect to circPUM1. Data are shown as −ΔCt values. In the x-axis, the number of CC samples is reported.
Genes 15 00124 g006
Table 1. Primer sequences.
Table 1. Primer sequences.
#Gene IDHost GeneForward Primer SequenceReverse Primer SequenceProduct Length (bp)
1hsa_circ_0088943NUP188TGAAACTGGCATTCTCCGTCAGGGCACTGAAGTAGCATGTTG107
2hsa_circ_0077150PHIPCATCACTACAGGTGAAGTATTACAGCTCCAGGAACACTTTGAGGAAT111
3hsa_circ_0077147PHIPTCTGGGATGCTGGAACCCTTAGTCATCAGAACCCAGCACTAA112
4hsa_circ_0004079 (1)TRIP12CACACGCCAAAAGACCACGTAGACTGCACTTTGGGTGCC125
hsa_circ_0004079 (2)GCTAGTACCAGGTCACATTTAGTTAGACTGCACTTTGGGTGC98
5hsa_circ_0068888WHSC1AGTGTCGGGTTACCCTTGGTCACAGAAAAGCAGACAGCTCG271
6hsa_circ_0011233PUM1GAGGCCACGTCCTGTCATTGTGAGTCCTCCTGCTGGTCTGA95
7hsa_circ_0043064 (1)MYO1DCCTGAAGGCAAACTGAGCATGCACAATCACTCGAGACTTTGA121
hsa_circ_0043064 (2)TACAGAGGTGACCAAGCGACCTCCTGGCTGTTGCACAATC113
8GAPDHGAPDHTGCACCACCAACTGCTTAGCGGCATGGACTGTGGTCATGAG87
9linear PUM1PUM1CCTTCAGACCAGCAGAATGAGAATTGCAAAGACTGGGGCTGT127
10BCL2BCL2CATGTGTGTGGAGAGCGTCAAGCCGGTTCAGGTACTCAGTCA83
11PTENPTENGACAATCATGTTGCAGCAATTCACCCATAGAAATCTAGGGCCTCT123
12TP53TP53GAGCACTGCCCAACAACACGTCTGGTCCTGAAGGGTGAA85
Table 2. Selected circRNAs.
Table 2. Selected circRNAs.
#circRNA IDHost GeneHost Gene IDmiRNAs
1hsa_circ_0088943Nucleoporin 188NUP188hsa-miR-16-5p
hsa-miR-214-3p
2hsa_circ_0077150Pleckstrin homology domain Interacting proteinPHIPhsa-miR-16-5p
3hsa_circ_0077147hsa-miR-16-5p
4hsa_circ_0004079 (1)Thyroid hormone receptor interactor 12TRIP12hsa-miR-372-3p
hsa-miR-18a-5p
hsa-miR-214-3p
hsa_circ_0004079 (2)
5hsa_circ_0068888Wolf–Hirschhorn syndrome candidate 1WHSC1hsa-miR-16-5p
hsa-miR-214-3p
6hsa_circ_0011233Pumilio RNA-binding family member 1PUM1hsa-miR-16-5p
hsa-miR-214-3p
7hsa_circ_0043064 (1)Myosin IDMYO1Dhsa-miR-16-5p
hsa_circ_0043064 (2)
Table 3. Information from the CircInteractome database (https://circinteractome.nia.nih.gov/index.html (accessed on 20 November 2023)).
Table 3. Information from the CircInteractome database (https://circinteractome.nia.nih.gov/index.html (accessed on 20 November 2023)).
circRNA IDHsa_circ_0011233
Locationchr1:31422979-31426828
Genomic Length3849 bp
Spliced Seq Length527 bp
Best TranscriptNM_014676
Gene SymbolPUM1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Caponnetto, A.; Ferrara, C.; Fazzio, A.; Agosta, N.; Scribano, M.; Vento, M.E.; Borzì, P.; Barbagallo, C.; Stella, M.; Ragusa, M.; et al. A Circular RNA Derived from the Pumilio 1 Gene Could Regulate PTEN in Human Cumulus Cells. Genes 2024, 15, 124. https://doi.org/10.3390/genes15010124

AMA Style

Caponnetto A, Ferrara C, Fazzio A, Agosta N, Scribano M, Vento ME, Borzì P, Barbagallo C, Stella M, Ragusa M, et al. A Circular RNA Derived from the Pumilio 1 Gene Could Regulate PTEN in Human Cumulus Cells. Genes. 2024; 15(1):124. https://doi.org/10.3390/genes15010124

Chicago/Turabian Style

Caponnetto, Angela, Carmen Ferrara, Anna Fazzio, Noemi Agosta, Marianna Scribano, Maria Elena Vento, Placido Borzì, Cristina Barbagallo, Michele Stella, Marco Ragusa, and et al. 2024. "A Circular RNA Derived from the Pumilio 1 Gene Could Regulate PTEN in Human Cumulus Cells" Genes 15, no. 1: 124. https://doi.org/10.3390/genes15010124

APA Style

Caponnetto, A., Ferrara, C., Fazzio, A., Agosta, N., Scribano, M., Vento, M. E., Borzì, P., Barbagallo, C., Stella, M., Ragusa, M., Scollo, P., Barbagallo, D., Purrello, M., Di Pietro, C., & Battaglia, R. (2024). A Circular RNA Derived from the Pumilio 1 Gene Could Regulate PTEN in Human Cumulus Cells. Genes, 15(1), 124. https://doi.org/10.3390/genes15010124

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop