- 1Unit of Basic and Applied Biosciences, University of Teramo, Teramo, Italy
- 2National Research Council, Institute of Biochemistry and Cell Biology, Rome, Italy
- 3Department of Veterinary Medicine, University of Perugia, Perugia, Italy
- 4Centre for Perinatal and Reproductive Medicine, University of Perugia, Perugia, Italy
In vitro folliculogenesis (ivF) has been proposed as an emerging technology to support follicle growth and oocyte development. It holds a great deal of attraction from preserving human fertility to improving animal reproductive biotechnology. Despite the mice model, where live offspring have been achieved,in medium-sized mammals, ivF has not been validated yet. Thus, the employment of a network theory approach has been proposed for interpreting the large amount of ivF information collected to date in different mammalian models in order to identify the controllers of the in vitro system. The WoS-derived data generated a scale-free network, easily navigable including 641 nodes and 2089 links. A limited number of controllers (7.2%) are responsible for network robustness by preserving it against random damage. The network nodes were stratified in a coherent biological manner on three layers: the input was composed of systemic hormones and somatic-oocyte paracrine factors; the intermediate one recognized mainly key signaling molecules such as PI3K, KL, JAK-STAT, SMAD4, and cAMP; and the output layer molecules were related to functional ivF endpoints such as the FSH receptor and steroidogenesis. Notably, the phenotypes of knock-out mice previously developed for hub.BN indirectly corroborate their biological relevance in early folliculogenesis. Finally, taking advantage of the STRING analysis approach, further controllers belonging to the metabolic axis backbone were identified, such as mTOR/FOXO, FOXO3/SIRT1, and VEGF, which have been poorly considered in ivF to date. Overall, this in silico study identifies new metabolic sensor molecules controlling ivF serving as a basis for designing innovative diagnostic and treatment methods to preserve female fertility.
Introduction
Assisted reproductive technologies (ARTs) represent a consolidated clinical practice, which have resulted in several million births since 1978 (Adamson et al., 2018) (International Federation of Fertility Societies’ Surveillance, 2019). However, the current ART protocols allow to use only a limited number of oocytes derived from antral follicles, while the large ovarian reserve, represented by the pool of early-stage follicles, remains a genetic patrimony that can be preserved by cryopreservation procedures but not managed (Martinez et al., 2017). In this context, the development of protocols aiming to obtain and fertilize mature oocytes from immature follicles grown outside the body could represent a useful strategy to recover the largest pool of ovarian gametes by promoting their in vitro growth and differentiation.
This follicle rescue approach could be applied as a possible future clinical strategy to preserve the ovarian reserve, as in the case of adult and prepubertal patients with cancer, where ovarian transplantation may expose to the risk of reintroducing malignant cells (Herta et al., 2018) (De Vos et al., 2014). Besides, it would benefit the veterinary medicine, providing a way to improve reproductive performance of species of zootechnical interest as well as to conserve the genetic inheritance of endangered animals.
Since the first attempt in the field of in vitro folliculogenesis (ivF) in 1996 (Eppig, 1996) (O’Brien et al., 2003), many efforts are made to set up new culture systems able to support in vitro growth of early-stage follicles toward competent oocytes (Laronda et al., 2017) (Xiao et al., 2015) (Xiao et al., 2017) (McLaughlin et al., 2018). The use of animal models to recapitulate ivF steps, driving to the production of fertilizable oocytes, has proven to be decisive for providing a knowledge basis and for developing validated methods with high translational potential for humans. Based on similarities in physiology and anatomy of the ovaries, folliculogenesis timing, and the follicle size (Barboni et al., 2011) (Bähr and Wolf, 2012) (Telfer and Zelinski, 2013), medium-sized mono-ovulatory mammals are commonly accepted as a translational model, and they are increasingly considered as being very relevant for human preimplantation reproductive research.
However, apart from the murine model where in vitro production of fertilizable oocytes has reached high levels of efficiency in terms of embryo development (Xu et al., 2006) (Hornick et al., 2013) (Hikabe et al., 2016), in medium-sized mammals, ivF remains still experimental. Indeed, a very low number of embryos produced from in vitro grown preantral follicles were reported in these models (de Figueiredo et al., 2018) even if several groups are working on bovine (Gupta et al., 2008) (Gupta and Nandi, 2012) (Antonino et al., 2019), porcine (Wu et al., 2001), caprine (Magalhães et al., 2011), ovine (Barboni et al., 2011) (Luz et al., 2012) (Arunakumari et al., 2010), and non-human primates (Xu et al., 2009) (Xu et al., 2011).
The difficulty to recapitulate in vitro the process of folliculogenesis in non-rodent animal models appeared to be related to the longer period required for follicle/oocyte growth, the greater dimension of antral follicles with competent oocytes, and the difficulty to mimic the environmentally favorable conditions to guarantee a synergic oocyte and somatic compartment development by preserving the tissue architecture (Rossetto et al., 2016).
Considering the large amount of data collected in vitro on the molecular mechanisms involved in the folliculogenesis among different species and the advances in in vitro follicle culture models, involving in vitro 2D and 3D culture approaches (Jones and Shikanov, 2019), the adoption of mathematical models might represent a valuable tool to organize the evidence collected to date by offering predictive models.
This study supports the use of a computational method based on network theory to identify the molecular events and the main factors sustaining ivF steps in mammals and to discover new molecular players in the ivF process to be targeted and/or exploited for therapeutic purposes, improving female fertility.
Materials and Methods
Data Collection: Web of Science-Mammals-Made ivF Database (WoS_MMivF)
Scientific literature published in the peer-reviewed international indexes such as the Advanced Search of Web of Science (v.5.35) “Core collection” archive (https://apps.webofknowledge.com/WOS_AdvancedSearch) of the past 30 years was considered (Bernabò et al., 2010) (Bernabò et al., 2014) using the following key words: “in vitro culture”, “Follicle culture”, “in vitro folliculogenesis”, “Oocyte”, and “Ovary”. “AND” and “NOT” were used as Boolean operators, and “TS” was used as a field tag.
Each list obtained from manual data mining was matched to create a unique database “WoS_MMivF”, including exclusively mammalian-related manuscripts, which accounts totally for 1,111 papers. The quality control of manually collected data was carried out according to Bernabò et al., (Bernabò et al., 2016).
The final WoS_MMivF database contains 513 selected ivF-related manuscripts classified as original primary research articles (444/513; 87%) and reviews (69/;513; 13%), according to Taraschi et al. (Taraschi et al., 2020). The WoS_MMivF database was enriched in Microsoft Excel 365 with the following fields (Supplementary File 1):
a) Source molecule: The molecule working as the source of interaction.
b) Interaction: The interaction the molecules carry out.
c) Target molecule: Molecules or molecular events that are the target of interaction.
d) Species: Different species of mammals in which molecular interactions occur.
e) Reference: PubMed IDentifiers (PMID).
Additional details related to the database set up can be found below:
The freely available and diffusible molecules such as H2O, CO2, Pi, and H+O2 were mainly omitted. In case the target where a single molecular determinant of the phenomenon is unknown as a target, the related ovarian function was indicated (i.e., “preantral follicle growth” and “follicle activation”).
ivF Network Creation, Visualization, and Analysis
The data, extracted from the database, were used to build the ivF network using the Cytoscape 3.6.0 software (http://www.cytoscape.org) (Paul et al., 1971). The network was analyzed with the specific plug-in Network Analyzer by computing the topological parameters described in Supplementary Data Sheet S2. The hubs, defined as hyperconnected nods, were identified as previously described (Bernabò et al., 2015b) (Bernabò et al., 2015a) by using the following equation:
where
γ = number of links per node (connectivity).
μ = mean node degree
σ = node degree standard deviation.
Closeness Centrality
Closeness centrality is a measure of how fast information spreads from one node to another reachable node (https://med.bioinf.mpi-inf.mpg.de/netanalyzer/help/2.7/index.html#refNewman2003).
This parameter is defined as the reciprocal of the average shortest path length and is computed as follows: Cc(n) = 1 / avg(L(n,m).
Here, L(n,m) is the length of the shortest path between two nodes n and m. The closeness centrality of each node is a number between 0 and 1. Network Analyzer computes the closeness centrality of all nodes and plots it against the number of neighbors. The closeness centrality of isolated nodes is equal to 0.
Betweenness Centrality
Cb(n) of a node n is computed as follows: Cb(n) = ∑s≠n≠t (σst (n) / σst), where s and t are nodes in the network different from n, σst denotes the number of shortest paths from s to t, and σst (n) is the number of shortest paths from s to t that n lies on. The betweenness centrality is computed only for networks that do not contain multiple edges. The betweenness value for each node n is normalized by dividing by the number of node pairs excluding n: (N-1) (N-2)/2, where N is the total number of nodes in the connected component that n belongs to. Thus, the betweenness centrality of each node is a number between 0 and 1 (https://med.bioinf.mpi-inf.mpg.de/netanalyzer/help/2.7/index.html#nodeBetween).
Identification of Bottlenecks (CytoHubba)
Bottlenecks (BN) were identified as follows: let Ts be a shortest path tree rooted at node s, BN(v) = Σs∈V ps(v). In detail, e ps(v) = 1 if more than |V(Ts)|/ 4 paths from node s to other nodes in Ts meet at the vertex v; otherwise ps(v) = 0 (Chin et al., 2014) (Taraschi et al., 2020) (Ordinelli et al., 2018).
In/Out Degree Ratio
In/out degree ratio (DRIO) was computed as follows:
It expresses the ratio between the number of links in the input and in the output for each node; consequently, it has been used to layer the nodes in input, processing, and output strata of the network.
Network Topology Transition
To assess the relevance of network controllers in the maintenance of network stability, the removal of the most connected nodes on the network topology (targeted attack theory) was performed by means of two cycles of attacks, removing 2.5% of hubs each time. At the end of each attack, the network topology was examined.
Enrichment Analysis
A Search Tool for the Retrieval of Interacting Genes/Proteins (STRING, http://string-db.org/newstring_cgi/show_input_page.pl?UserId=eNOo92_OQ_LS&sessionId=Cfz4mDP5ayne) (Szklarczyk et al., 2015) was used in order to enrich the database by including known and predicted protein interactions. They could be either direct (physical) or indirect (functional) associations and are derived from different sources: genomic context, high-throughput experiments, conserved coexpression, and previous knowledge. A new network was obtained (STRING_MMivF) by adopting a medium confidence score (0.400). For the enrichment procedures, the false discovery rate (FDR) value was set to be <0.05, and 4 cycles of enrichment were performed.
Gene Ontology
Gene Ontology (GO; https://www.geneontology.org) was carried out to identify the main functions, processes, and cellular compartments of the hub.BN and interactors according to their GO terms. The FRD value was set for p < 0.05.
Results
The ivf process is described by a scale-free, non-clustered, and non-hierarchical system network.
Up to half of the links within the network are referred to rodents (39%) and ruminants (ovine 17.2, caprine 3.9, and bovine 13.7%), whereas less than 9% comes from humans (Table 1).
Based on data collected in the WoS_MMivF database, it was possible to generate a biological network whose topological parameters were computed as reported in Supplementary File 3.
The network displayed a scale-free topology, according to the Barabási-Albert (BA) model. The statistical analysis of its topology (Table 2) demonstrated that it recognizes 641 nodes and 2089 links, with a very low clustering coefficient (0.076).
The main controllers of the network were identified: hub, bottleneck, and hub-bottleneck nodes.
The ivF network recognized as hyperconnected nodes (hubs) 5.9% of nodes (38 out of a total of 641) (Supplementary File 4), which were ranked based on node degree (Supplementary File 4). They include activating systemic hormones and ovarian activating factors (8 and 11 out of 38, respectively; overall 50% of hubs) more than transduction terminal events (4, 10, and 5%).
A Kernel density estimation (KDE) based on their clustering coefficient was then performed (see Supplementary File 4 for clustering coefficient value) in order to explore their distribution, showing that no node subpopulation is present (as reported in Supplementary File 5).
Furthermore, to assess the hub role in the ivF network stability control, a computational experiment was performed, carrying out subsequent cycles of hub removal (Supplementary File 6). The network topology was deeply affected by hub removal and collapsed upon two cycles of network attack. As comparison, the removal of the same number of nodes randomly identified [=random.between(min;max)] did not have detectable effects on network stability (data not shown).
In addition, seeking for controllers of the information flow within the network, 38 bottleneck nodes were identified (BN: Supplementary File 4) (Ordinelli et al., 2018) (Taraschi et al., 2020).
Finally, by intersecting the two subsets of nodes categorized as hubs and BN, 30 main nodes were identified and named hub.BN (Supplementary File 7).
Role of hub.BN Nodes in the Scale-free ivf Network
More in detail, the hub.BN enclosed 10 functional events and 20 molecules.
Among the events, four stages of follicle development (primordial follicle activation, primordial to primary follicle transition, preantral follicle growth, and preantral to antral follicle transition), three key outcomes of the early stage of follicle development (antrum differentiation, steroidogenesis, and meiotic competence), and three cell functions controlling tissue homeostasis (cell survival, proliferation, and apoptosis) were identified.
The 20 molecules belonging to hub.BN summarized in Supplementary File 7 based on their role in the ovarian folliculogenesis (hormone or paracrine/autocrine factors and driven follicular events) were classified for their biological function, process, and cellular localization (Figure 1) by using Gene Ontology (GO) (http://geneontology.org/) and setting the selection on the FDR value <0.05 (Supplementary File 8). The top 10 most abundant categories in terms of GO for hub.BN recognized hormone/paracrine binding (5 out of 10) and activities (4 out of 10) as key biological functions addressed to regulate the key cellular outcome processes of signaling (4 out of top 10 GO biological processes category) and cell to cell communication (1 out of top 10) inside the reproductive systems (4 out of 10). The main extracellular (3 out of 10) and secretory vesicle (3 out of 10) localization of hub.BN molecule actions was also highly consistent with the component GO category results.
FIGURE 1. GO enrichment analysis. Representative scheme of the top 10 most abundant GO terms identified for the hub.BN molecules in the three GO categories: biological function (red), cellular localization (green), and biological processes (blue). The x-axis indicates the hub.BN molecules in a specific category, while the y-axis indicates different GO terms.
Interestingly, knock-out mouse model data retrieved from the Mouse Genome International database (MGI, http://www.informatics.jax.org/) and WoS database (https://apps.webofknowledge.com) for the identified hub.BN genes showed altered phenotypes in the reproductive system (Supplementary File 9), supporting their key role as controllers of the network.
A 2D KDE analysis was then carried out to identify eventual subpopulation in the hub.BN population based on node degree and BN scores. Four isolated nodes and three subpopulations were identified (Figure 2). Four isolated hub.BN were characterized by higher values of node degree and BN score. Among them, it was possible to distinguish the main target of oogenesis related to the early stage of follicle development (acquisition of oocyte meiotic competence) and two follicular stages reproducible in vitro (primordial to primary follicle transition and preantral follicle growth). In addition, the 2D shape KDE computed analysis was assigned to this main subpopulation of hub.BN and also to the pituitary hormone FSH.
FIGURE 2. KDE 2D analysis of the network hubs. Centrality parameters of BN and node degree parameters of the hub.BN were considered for the 2D KDE analysis.
Subpopulations 3 and 2 recognized seven key endocrine/paracrine controllers of folliculogenesis (activin A and LH, GDF9, EGF, E2, Insulin, and IGF1) as well as the event of follicle specialization that occurs in vitro by applying the current protocols (antrum differentiation). The remaining hub.BN belongs to subpopulation 1.
To organize the network depending on in/out connectivity, the in/out degree ratio of each node was computed (Figure 3A). This analysis allowed to stratify nodes in three layers (input, processing, and output layers) by classifying the network nodes on the basis of the in/out degree (Figure 3B):
FIGURE 3. Nodes with Hub, bottleneck, and hub.BN role stratification into the ivF network. (A) Diagram showing the signal stratification of the ivF network. The color gradient varies depending on the direction of links characterizing each node, computed as in/out degree ratio, from purple (higher) to yellow (lower). The spatial network arrangement was obtained by using the Cytoscape prefuse force-directed layout. (B) Node classification depending on their role in the signal propagation in the input layer (green scale color), processing layer (yellow scale color), and output layer (red scale color).
DRio 0–35 = input layer.
DRio 36–67 = processing layer.
DRio 68–100 = output layer.
The analysis displayed that the hub.BN was mostly abundant in the input layer (15 out of 30). BN mainly operated as processing (4 out of 8), whereas hubs were distributed in either input or output layers (3 out of 8, respectively).
The analysis of the ivF signaling network designed a coherent stratification of nodes by positioning in the input layer 16 out of 22 nodes belonging to the systemic endocrine controller released by pituitary (FSH and LH), chorion (hCG), and enteric-related endocrine glands (insulin, DHT, and IGF1) as well as reproductive controlling hormone/factors secreted either from follicles (FST, activin A, AMH, EGF, bFGF, TGFbeta, BMP4, and BMP15) or from oocytes (GDF9 and cumulin).
Analogously, the processing layer recognized among the eight nodes the intracellular second messenger cAMP and components of four signaling pathways: PI3K/AKT, JAK/STAT, TGFbeta (SMAD4), and KL (Figure 3B). Finally, 14 out of 17 nodes of the output layer were follicular events, and further three nodes were key components of steroidogenesis starting from the FSH receptor, the endpoint enzyme, and the hormone of follicular steroidogenesis (CYP19A1 and P4).
Identification of New Molecular Players in the ivf Network Flow
To identify and predict new molecules involved in the ivF process, a functional protein association network (STRING_MMivF) was created by STRING using as input hub.BN genes, where possible (Figure 4A). Totally, 30 new predicted interactor molecules were identified and summarized in Supplementary File 10.
FIGURE 4. STRING-MMivF as the interaction network tuning the ivF process. (A) Known and predicted protein–protein interactions related to the hub.BN molecules of the network. The interactions include direct (physical) and indirect (functional) associations; they stem from computational prediction, from knowledge transfer between organisms, and from interactions aggregated from other (primary) databases. (B) Representative scheme of the top 10 most abundant GO terms identified for the new interactor molecules in the three GO categories: biological function (red), cellular localization (green), and biological processes (blue). The x-axis indicates the number of new interactor molecules in a specific category, while the y-axis indicates different GO terms.
Functional enrichment analysis was then performed (Supplementary File 11). The 10 most significantly enriched terms (p < 0.05) in each category are presented in Figure 4B. In the molecular function category, GO recognized the growth factor and intracellular signaling protein binding (7 out of 10) as the most abundant biological function regulating cellular processes such as metabolism control (2 out of 10), signal transduction (6 out of 10), and cell to cell communication (2 out of 10), sustaining angiogenesis. For the cellular component category, intracellular (3 out of 10), extracellular (2 out of 10), and mTOR complex (3 out of 10) localizations were reported accordingly.
Interestingly, involvement in follicular functions was found for most of the newly identified interactors, as reported in Supplementary File 10. Conversely, for four of them (KRAS, MLST8, MAPKAP1, and FGFR1), a specific role in the ovarian folliculogenesis has not been described to date. Of note, an abnormal female reproductive phenotype has been reported for KRAS and FGFR1 knock-in and knock-out mouse models, respectively (Table 3). Moreover, in the context of the STRING-MMivF network, KRAS was predicted to interact with a high confidence score (score >0.9) with hub.BN AKT1, FGF2, and EGF but FGFR1 with hub.BN FGF2 (Table 3).
Discussion
The ivF network (Albert and Barabási, 2002) is characterized by a power-law distribution of node degree and by the absence of correlation between the node degree itself and the clustering coefficient in keeping with the Barabasi-Albert model, conferring to the network the following biologically relevant properties:
1) Robustness against random damage: when a random perturbation affects the network, it is very likely that only a scarcely linked node (i.e., a node belonging to the most frequent class of nodes) will be affected. Thus, the probability that a hub of the graph is affected remains very low, and in this model, it can be estimated to be about 5.9%.
2) Controllability: the small number of highly linked nodes implies that the whole system can be modulated with high efficiency by acting just on a few molecules, thus reducing the energetic cost and facilitating/accelerating the cell response.
3) Easy navigability: the virtual absence of clustering, together with the low values of the characteristic path length and of the average number of neighbors, confers to the network a typical structure of signaling networks. Interestingly, this low value of clustering coefficient also implies a little redundancy.
After having defined the topology, the computational analysis allowed to dissect the network to identify the ivF network flow (hub.BN). Indeed, the 2D KDE approach allowed to select those with a major modulatory role on the network information flow by identifying either cultural functional endpoints or in vitro follicle controllers such as systemic hormones and local factors (Matzuk et al., 2002).
More in detail, the priority of hub.BN was assigned by 2D KDE analysis to four events related to both follicle (primordial to primary follicle transition and preantral follicle growth) and oocyte development (oocyte meiotic competence). Furthermore, the analysis identified FSH as a key hub.BN; thus, reinforcing the idea of its central role as a controller of either the gonadotropin-dependent or independent phase of folliculogenesis (Holesh and Lord, 2018) (Casarini and Crépieux, 2019). Several pieces of scientific evidence collected mainly using in vitro studies, indeed, suggest the need to supplement FSH to also stimulate the early stage phases (from primary to later stages) (Hsueh et al., 2015) of folliculogenesis by stimulating either follicle growth or differentiation (Cecconi et al., 1999) (Barboni et al., 2011) (Cadoret et al., 2017) (Candelaria et al., 2020) (Cortvrindt et al., 1997) (Kreeger et al., 2005) (Wright et al., 1999). In addition, 7 hub.BN have been identified as main controllers of the network (activin A and LH, GDF-9, EGF, E2, Insulin, and IGF-1) (Orisaka et al., 2006) (Dong et al., 1996) (Edson et al., 2009).
The results of the network stratification analysis are in agreement with current knowledge on the physiological cross-talk between molecules regulating the inter- and intra-follicular communication. Indeed, in the input layer were identified hub.BN recapitulating the main modulatory factors, either hormones exerting a remote control on the ovary (LH, FSH, hCG, IGF1, and EGF) or molecules involved in intra-ovarian control of both somatic and germinal compartments (TGF-beta superfamily members and growth factors).
The processing layer is composed of hub.BN molecules belonging to PI3K/AKT, JAK/STAT, TGFbeta, and KL signaling pathways and the second messenger cAMP which transduce and amplify the actions of the previous endocrine and paracrine controllers (Li et al., 2008) (Fujihara et al., 2014) (John et al., 2009) (Imai et al., 2014).
Finally, the output layer recognized steroidogenic molecules (P4 and CYP19A1) which may be considered as endpoints of the pathways controlling ivF. The emerging picture defined by the computational stratification analysis thus supports the biological strength of the network recapitulating and organizing in three layers the physiological feedback orchestrating ovarian follicle development.
The biological strength of the network has been further confirmed by analyzing the reproductive phenotype of the available knock-out mice generated by using some ivF selected hub.BN. The analysis of the scientific evidence collected to date not only confirmed that the in vivo silencing of these controllers always promoted a negative impact on female fertility but also confirmed the interference specifically on pathways controlling the early phase of in vivo folliculogenesis, despite the different underlying mechanisms. STRING analysis was then performed to find new molecules and define their functional connections. Taking advantage of such an approach, molecules belonging to three effector categories have been identified as potential controllers of the early stage of the in vitro follicle development. Metabolic controllers involved in the signaling leading to follicle recruitment, apoptosis, and differentiation were identified. Indeed, STRING enriched the ivF network with mTOR and IGF1 signaling pathway molecules.
mTOR represents a molecular sensor for diverse environmental inputs including nutrients and growth factors and regulates various fundamental processes including cell growth, the metabolism, differentiation, and autophagy (Long et al., 2003). In folliculogenesis, some lines of evidence have shown that the IGF1-dependent activation (hub.BN role in the ivF network) of the PI3K/AKT signaling (hub role in the ivF network) and its downstream cognate mTOR complex sustains ovarian primordial follicle dormancy and activation, oocyte maintenance and activation, and GC proliferation and differentiation (Makker et al., 2014) (Guo and Yu, 2019) (Herta et al., 2018).
More in detail, mTORC1 plays a documented role in primordial follicle-oocyte bidirectional signaling. It was shown to activate the KIT receptor in oocytes through the KIT ligand (hub role in the ivF network), which triggers a PI3K/PTEN/AKT/forkhead box 3 (FOXO3) cascade and awakens the dormant oocytes (Zhao et al., 2018) (Liu et al., 2014). In the awakened oocyte, secretion of oocyte-specific growth factors such as BMP15 (hub.BN role in the ivF network) and GDF9 (hub.BN in the ivF network) further activates receptor serine kinases and downstream SMAD (mothers against DPP homolog 1 Drosophila) proteins including SMAD4 (BN role in the ivF network) in surrounding GC, leading to their growth and proliferation (Sanfins et al., 2018). A targeted deletion of the mTORC1 negative regulator TSC2 (STRING-enriched partner) in mouse oocytes results in prematurely follicular activation due to elevated mTORC1 activity in oocytes which in turn cause depletion of follicles in the early adulthood (Adhikari and Liu, 2009). Of note, a compensatory elevation of PI3K signaling was proposed to the reason for the unaffected follicular development observed in RPTOR (STRING-enriched partner) conditional KO in primordial and all subsequent oocyte stages (Gorre et al., 2014). Conversely, conditional KO mice for RPTOR in primordial follicle GC prevent the cell differentiation, and this arrests the dormant oocytes in their quiescent states, leading to the oocyte death age (Liu et al., 2014), indicating that the KIT/PI3K cascade in oocytes is indispensable for primordial follicle survival.
The key role of mTOR signaling was also supported by the computational identification of the downstream target of mTOR, such as FOXO1, FOXO3, and SIRT1 (STRING-enriched partners). Furthermore, transgenic mice for these molecules have provided evidence on their cross-talk in regulating the dynamics of the primordial follicle pool (Adhikari and Liu, 2009) (Adhikari et al., 2013) (Chen et al., 2015).
FOXO3 functions at the earliest stages of follicular growth as a suppressor of follicular overactivation, increasing the follicle reserves in the ovary in order to extend the reproductive period of females (Castrillon et al., 2003) (Shah et al., 2018) (Lee and Chang, 2019). Accordingly, FOXO3 has been found to be highly expressed in the nuclei of oocytes of primordial follicles, and its expression is downregulated in oocytes of primary and later-growing follicles, indicating that its downregulation in oocytes could be a prerequisite for the initiation of oocyte growth during follicular activation (Liu et al., 2007).
Furthermore, FOXO3−/− female mice exhibit a distinctive ovarian phenotype of global follicular activation, leading to oocyte death, early depletion of functional ovarian follicles, and secondary infertility (Castrillon et al., 2003) (Hosaka et al., 2004) (Lin et al., 2004).
Of note, selective depletion of FOXO1 and FOXO3 in mouse GC leads to an infertile phenotype characterized by metabolic changes and the production of factors that exerts potent negative feedback to prevent gene expression of pituitary FSH (hub.BN role in the ivF network). Decreased levels of serum FSH further restrict follicle growth and development, ultimately preventing ovulation. Besides, FOXO1/3 depletion alters the expression of genes involved in follicle growth and apoptosis, disrupting cell regulatory signals associated with the granulosa cell metabolism and follicle growth. These results reveal a novel ovarian-pituitary endocrine feedback loop preventing uncontrolled proliferation and/or premature differentiation of GC in follicles where apoptosis is impaired (Liu Z. et al., 2013).
Accordingly, the STRING-enriched molecule FOXO1 was shown to act as a silent guardian of follicle development. Indeed, it is a critical factor in promoting GC apoptosis to counteract the stimulatory FSH role, one of the selected ivF network hub.BN (Liu et al., 2009) (Fan et al., 2010) (Shen et al., 2014). More in detail, FSH was found to promote the expression of GC genes required for proliferation, survival, and estrogen synthesis by decreasing FOXO1, which negatively regulates proliferation and steroidogenesis (Liu et al., 2009) (Fan et al., 2010) (Shen et al., 2014) (Rosairo et al., 2008).
Moreover, the identification of SIRT1 as a potential partner of STRING enhances the role of the mTOR pathway in the control of in vitro early mammalian folliculogenesis. Indeed, it has been proposed as the molecular mediator of the calorie restriction-dependent prevention of follicular activation, leading to ovarian reserve preservation (Bordone et al., 2007) (Long et al., 2019). At a molecular level, the effect of SIRT1 on the ovary occurs through the deacetylation dependent-activation of FOXO3 and the suppression of the mTOR signaling (Long et al., 2019). The latter could be mediated by TSC2-SIRT1 interactions (Ghosh et al., 2010). These results are in good agreement with data showing that the SIRT1 activator (SRT1720) improves the follicle reserve and prolongs the ovarian lifespan of diet-induced obesity in female mice via activating SIRT1 and suppressing mTOR signaling (Zhou et al., 2014). Conversely, a recent work shows how SIRT1 can sustain the activation of mouse primordial follicles independent of its deacetylase activity. Specifically, SIRT1 was shown to trigger primordial follicle awakening by activating the PI3K/AKT-AKT and TSC1/2-mTOR signaling pathways. Indeed, its pharmacologically (resveratrol) and/or genetic-induced activation in mouse ovary cultures leads to SIRT1 to work as a transcription cofactor which modulates the expression levels of AKT, mTOR, and genes related to classic primordial follicle activation (Zhang et al., 2019). Interestingly, SIRT1 overactivation was shown to inhibit FOXO3 activity by promoting its exclusion form the nucleus through the AKT-dependent phosphorylation (Zhang et al., 2019), suggesting that SIRT1 activators might be used to efficiently activate the primordial follicles to be applied in in vitro activation protocols and/or to avoid uncontrolled follicular atresia characterizing POF.
The evidence of mTOR and SIRT1 pathways on follicular fate has been achieved by taking advantage of their in vivo and in vitro pharmacological modulation.
In vivo, the SIRT1 activation has been achieved by caloric restriction or the application of specific sirtuin activators or mTOR blockers. Caloric restriction was shown to induce the accumulation of SIRT1 in murine ovaries, which is associated with inhibition of primordial follicle activation and impairment of ovarian follicle development (Tatone et al., 2015) (Liu et al., 2015) (Long et al., 2019) mechanisms probably involved in mice delay of puberty, prolongation of reproductive lifespan, and prevention of age-associated infertility (Liu M. et al., 2013).
Various effects on folliculogenesis have been observed in vitro by using the SIRT1 activator, mTOR inhibitor, and resveratrol.
More in detail, it promoted growth in human ovarian follicles (Hao et al., 2018) and in ovine primordial follicles (Bezerra et al., 2018). The resveratrol action is addressed to regulate GC cell proliferation and survival (Han et al., 2017) (Ortega et al., 2012) (Schube et al., 2014) (Wong et al., 2010) as well as steroidogenesis (Qasem, 2020) (Morita et al., 2012) even if with contradictory results.
In late folliculogenesis, mTOR signaling is dependent on gonadotropin FSH modulation to manage processes such as oocyte maturation (Guo et al., 2018), ovarian somatic cell proliferation, and steroidogenesis (Palaniappan and Menon, 2012) and to intensify EGFR/RAS signaling (Fan et al., 2009) (Gloaguen et al., 2011) (Fan et al., 2012) (Wayne et al., 2007).
Interestingly, the STRING network enrichment identified the isoform KRAS as a potential interactor of the AKT1 and EGF hub.BN molecules with a high confidence score, suggesting its key role as a metabolic effector of gonadotropin stimulation. Its contribution in follicular development is corroborated using conditional knock-in mouse models in which the GC expresses a constitutively active form of KRAS (KrasG12D). These models show alterations in GC differentiation, proliferation, and apoptosis at early stage of folliculogenesis, thus impairing cell responses to gonadotropins and leading to premature ovarian failure (Fan et al., 2008) (Fan and Richards, 2010) (Bulun et al., 2019). This phenotype is similar to that of the mutant mouse model with EGFR signaling defects (Hsieh et al., 2007). The altered response of KRASG12D-expressing GC to gonadotropins appears to be related to low levels of FSHR and the inability of FSH to induce expression of LHCGR mRNA and therefore to the loss of the crucial LH-MAPK3/1 signaling pathways. This conclusion is supported by the reduced expression of specific genes known to be essential for COC expansion and ovulation (Richards, 2005). Moreover, mutations causing KRAS overactivity have been reported as markers of epithelial and endometrioid ovarian cancer (Ramalingam, 2016) (Dinulescu et al., 2005).
Although the effect of these energy sensors in folliculogenesis remains to be investigated in depth, some data have been reported in the literature.
The third category of molecules enriched from the STRING network was represented by several angiogenetic factors. Indeed, it must be considered that as the follicle growth, the mTOR-dependent metabolic signaling needs to be implemented with components of the vascular system that guarantees the correct trophic supply and spreading of precursors to properly complete folliculogenesis (Martelli et al., 2017).
The angiogenesis process assumes a key role in the contest of the reproductive function (Robinson et al., 2009) (Wulff et al., 2002). Accordingly, angiogenesis inhibition leads to the reduction of follicular growth, ovulation disruption, and drastic effects on the corpus luteum activity and development (Robinson et al., 2009). The large amount of factors is able to control angiogenesis and their spatio-temporal regulatory activities, suggesting that more than one factor might be useful for angiogenesis well-functioning associated with ovulation.
Specifically, functional interactions from STRING were unveiled for VEGF members (VEGFA, VEGFB, VEGFC, and VEGFD), VEGF receptors (NRP1, KDR, and FLT), and angiogenic factors (FGF1, FGF9, PGF, and FGFR).
Several pieces of evidence confirmed that both somatic and germinal compartments contribute to modulate the expression of VEGFA by modulating the blood vessel network during preantral follicle development. More in detail, small and middle pig preantral follicles seem to behave as autonomous recruited units, where the growth of the somatic compartment is always accompanied by the simultaneous activation of the endothelial cells (Martelli et al., 2017).
Differently, once the preantral follicles reach the late stage enclosing an almost fully growth oocyte, it starts to express high and stable levels of VEGFA, essential to maintain follicle angiogenesis in a steady status of activation (Martelli et al., 2017). Furthermore, high levels of VEGFA also characterize the preantral follicular structures at the stage of antrum formation. In this context, VEGFA has been supposed to increase microvessel permeability by stimulating plasma extravasations, thus allowing the accumulation of fluids within the differentiating follicular cavity (Isobe et al., 2005).
The strict correlation between somatic and vascular parameters may represent an indirect biological validation of the key role of VEGF members identified by ivF network STRING analysis depicting a synergic action between follicle compartments and blood vessel system components as a key event driving preantral follicle activation first and then sustaining the process of transition from preantral to early antral follicles (Martelli et al., 2017) (Danforth et al., 2003) (Abramovich et al., 2006).
These molecules exert their function through their receptors (Stouffer et al., 2001) and are tightly regulated by paracrine factors such as nitric oxide. Nitric oxide has been reported to mediate positive effects on follicle development and selection related to angiogenic events and play a modulatory role in the ovarian steroidogenesis (Basini and Grasselli, 2015).
Nevertheless, proangiogenetic factors such as FGF1, FGF9, and PGF have also been described as they can potentiate the effect of VEGF signaling by increasing the VEGF expression in theca cells of cattle (Gray et al., 1987) (Gospodarowicz et al., 1987) (Berisha et al., 2004).
In conclusion, the computational analysis of the ivF network has enabled us to select among several molecules adopted to date those that are the main spatio-temporal controllers of the early stage of follicle development. Taking advantage of the biological robustness of such a network, new molecules slightly deepened or unexplored for ivF purposes were identified by taking advantage of the STRING approach. Altogether, this evidence suggests the ivF network as a sounded system biology tool to be exploited for research, technological, and innovation aims with the final goal of designing new diagnostic and therapeutic strategies for female fertility.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, and further inquiries can be directed to the corresponding author.
Author Contributions
BB and NB contributed to conception and design of the study. CD, GC, and GB organized the database. NB performed the data analysis. BB, NB, AP, CD, GC, GB, MM, MC, and UT contributed to the interpretation of the results. BB, NB, and AP wrote the first draft of the manuscript. CD, GC, and GB wrote sections of the manuscript. AP, BB, and NB contributed to figure manuscript editing. BB revised the whole manuscript. BB contributed to funding acquisition. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
This research was funded by MUR (Ministero dell’Universita e della Ricerca), Action I.2 “Attraction and International Mobility” PON (Programma Operativo Nazionale) “R&I” 2014–2020, grant number 894039_1.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at https://www.frontiersin.org/articles/10.3389/fmolb.2021.737912/full#supplementary-material
References
Abramovich, D., Parborell, F., and Tesone, M. (2006). Effect of a Vascular Endothelial Growth Factor (VEGF) Inhibitory Treatment on the Folliculogenesis and Ovarian Apoptosis in Gonadotropin-Treated Prepubertal Rats1. Biol. Reprod. 75, 434–441. doi:10.1095/biolreprod.106.051052
Adamson, G. D., de Mouzon, J., Chambers, G. M., Zegers-Hochschild, F., Mansour, R., Ishihara, O., et al. (2018). International Committee for Monitoring Assisted Reproductive Technology: World Report on Assisted Reproductive Technology, 2011. Fertil. Sterility 110, 1067–1080. doi:10.1016/j.fertnstert.2018.06.039
Adhikari, D., and Liu, K. (2009). Molecular Mechanisms Underlying the Activation of Mammalian Primordial Follicles. Endocr. Rev. 30, 438–464. doi:10.1210/er.2008-0048
Adhikari, D., Risal, S., Liu, K., and Shen, Y. (2013). Pharmacological Inhibition of mTORC1 Prevents Over-activation of the Primordial Follicle Pool in Response to Elevated PI3K Signaling. PLoS One 8, e53810–5. doi:10.1371/journal.pone.0053810
Albert, R., and Barabási, A.-L. (2002). Statistical Mechanics of Complex Networks. Rev. Mod. Phys. 74, 47–97. doi:10.1103/RevModPhys.74.47
Antonino, D. d. C., Soares, M. M., Júnior, J. d. M., de Alvarenga, P. B., Mohallem, R. d. F. F., Rocha, C. D., et al. (2019). Three-dimensional Levitation Culture Improves Iin-Vvitro Growth of Secondary Follicles in Bovine Model. Reprod. BioMedicine Online 38, 300–311. doi:10.1016/j.rbmo.2018.11.013
Arunakumari, G., Shanmugasundaram, N., and Rao, V. H. (2010). Development of Morulae from the Oocytes of Cultured Sheep Preantral Follicles. Theriogenology 74, 884–894. doi:10.1016/j.theriogenology.2010.04.013
Bähr, A., and Wolf, E. (2012). Domestic Animal Models for Biomedical Research. Reprod. Domest. Anim. 47, 59–71. doi:10.1111/j.1439-0531.2012.02056.x
Barboni, B., Russo, V., Cecconi, S., Curini, V., Colosimo, A., Garofalo, M. L. A., et al. (2011). In Vitro grown Sheep Preantral Follicles Yield Oocytes with normal Nuclear-Epigenetic Maturation. PLoS One 6, e27550. doi:10.1371/journal.pone.0027550
Basini, G., and Grasselli, F. (2015). Nitric Oxide in Follicle Development and Oocyte Competence. Reproduction 150, R1–R9. doi:10.1530/REP-14-0524
Berisha, B., Sinowatz, F., and Schams, D. (2004). Expression and Localization of Fibroblast Growth Factor (FGF) Family Members during the Final Growth of Bovine Ovarian Follicles. Mol. Reprod. Dev. 67, 162–171. doi:10.1002/mrd.10386
Bernabò, N., Barboni, B., and Maccarrone, M. (2014). The Biological Networks in Studying Cell Signal Transduction Complexity: The Examples of Sperm Capacitation and of Endocannabinoid System. Comput. Struct. Biotechnol. J. 11, 11–21. doi:10.1016/j.csbj.2014.09.002
Bernabò, N., Greco, L., Ordinelli, A., Mattioli, M., and Barboni, B. (2015a). Capacitation-Related Lipid Remodeling of Mammalian Spermatozoa Membrane Determines the Final Fate of Male Gametes: A Computational Biology Study. OMICS: A J. Integr. Biol. 19, 712–721. doi:10.1089/omi.2015.0114
Bernabò, N., Mattioli, M., and Barboni, B. (2015b). Signal Transduction in the Activation of Spermatozoa Compared to Other Signalling Pathways: A Biological Networks Study. Ijdmb 12, 59–69. doi:10.1504/IJDMB.2015.068953
Bernabò, N., Mattioli, M., and Barboni, B. (2010). The Spermatozoa Caught in the Net: The Biological Networks to Study the Male Gametes post-ejaculatory Life. BMC Syst. Biol. 4, 87. doi:10.1186/1752-0509-4-87
Bernabò, N., Ordinelli, A., Ramal Sanchez, M., Mattioli, M., and Barboni, B. (2016). Networks Models of Actin Dynamics during Spermatozoa Postejaculatory Life: A Comparison Among Human-Made and Text Mining-Based Models. Biomed. Res. Int. 2016, 1–8. doi:10.1155/2016/9795409
Bezerra, M. É. S., Gouveia, B. B., Barberino, R. S., Menezes, V. G., Macedo, T. J. S., Cavalcante, A. Y. P., et al. (2018). Resveratrol Promotes In Vitro Activation of Ovine Primordial Follicles by Reducing DNA Damage and Enhancing Granulosa Cell Proliferation via Phosphatidylinositol 3‐kinase Pathway. Reprod. Dom Anim. 53, 1298–1305. doi:10.1111/rda.13274
Bordone, L., Cohen, D., Robinson, A., Motta, M. C., Van Veen, E., Czopik, A., et al. (2007). SIRT1 Transgenic Mice Show Phenotypes Resembling Calorie Restriction. Aging Cell 6, 759–767. doi:10.1111/j.1474-9726.2007.00335.x
Bulun, S. E., Wan, Y., and Matei, D. (2019). Epithelial Mutations in Endometriosis: Link to Ovarian Cancer. Endocrinology 160, 626–638. doi:10.1210/en.2018-00794
Cadoret, V., Frapsauce, C., Jarrier, P., Maillard, V., Bonnet, A., Locatelli, Y., et al. (2017). Molecular Evidence that Follicle Development Is Accelerated In Vitro Compared to In Vivo. Reproduction 153, 493–508. doi:10.1530/REP-16-0627
Candelaria, J. I., Rabaglino, M. B., and Denicol, A. C. (2020). Ovarian Preantral Follicles Are Responsive to FSH as Early as the Primary Stage of Development. J. Endocrinol. 247, 153–168. doi:10.1530/JOE-20-0126
Casarini, L., and Crépieux, P. (2019). Molecular Mechanisms of Action of FSH. Front. Endocrinol. 10. doi:10.3389/fendo.2019.00305
Castrillon, D. H., Miao, L., Kollipara, R., Horner, J. W., and DePinho, R. A. (2003). Suppression of Ovarian Follicle Activation in Mice by the Transcription Factor Foxo3a. Science 301, 215–218. doi:10.1126/science.1086336
Cecconi, S., Barboni, B., Coccia, M., and Mattioli, M. (1999). In Vitro Development of Sheep Preantral Follicles1. Biol. Reprod. 60, 594–601. doi:10.1095/biolreprod60.3.594
Chen, Z., Kang, X., Wang, L., Dong, H., Wang, C., Xiong, Z., et al. (2015). Rictor/mTORC2 Pathway in Oocytes Regulates Folliculogenesis, and its Inactivation Causes Premature Ovarian Failure. J. Biol. Chem. 290, 6387–6396. doi:10.1074/jbc.M114.605261
Chin, C.-H., Chen, S.-H., Wu, H.-H., Ho, C.-W., Ko, M.-T., and Lin, C.-Y. (2014). cytoHubba: Identifying Hub Objects and Sub-networks from Complex Interactome. BMC Syst. Biol. 8, 1–7. doi:10.1186/1752-0509-8-S4-S11
Cortvrindt, R., Smitz, J., and Van Steirteghem, A. C. (1997). Assessment of the Need for Follicle Stimulating Hormone in Early Preantral Mouse Follicle Culture In Vitro. Hum. Reprod. 12, 759–768. doi:10.1093/humrep/12.4.759
Danforth, D. R., Arbogast, L. K., Ghosh, S., Dickerman, A., Rofagha, R., and Friedman, C. I. (2003). Vascular Endothelial Growth Factor Stimulates Preantral Follicle Growth in the Rat Ovary1. Biol. Reprod. 68, 1736–1741. doi:10.1095/biolreprod.101.000679
De Vos, M., Smitz, J., and Woodruff, T. K. (2014). Fertility Preservation in Women with Cancer. The Lancet 384, 1302–1310. doi:10.1016/S0140-6736(14)60834-5
Dinulescu, D. M., Ince, T. A., Quade, B. J., Shafer, S. A., Crowley, D., and Jacks, T. (2005). Role of K-Ras and Pten in the Development of Mouse Models of Endometriosis and Endometrioid Ovarian Cancer. Nat. Med. 11, 63–70. doi:10.1038/nm1173
Dong, J., Albertini, D. F., Nishimori, K., Kumar, T. R., Lu, N., and Matzuk, M. M. (1996). Growth Differentiation Factor-9 Is Required during Early Ovarian Folliculogenesis. Nature 383, 531–535. doi:10.1038/383531a0
Edson, M. A., Nagaraja, A. K., and Matzuk, M. M. (2009). The Mammalian Ovary from Genesis to Revelation. Endocr. Rev. 30, 624–712. doi:10.1210/er.2009-0012
Eppig, J. J. (1996). Development In Vitro of Mouse Oocytes from Primordial Follicles. Biol. Reprod. 54, 197–207. doi:10.1095/biolreprod54.1.197
Fan, H.-Y., Liu, Z., Mullany, L. K., and Richards, J. S. (2012). Consequences of RAS and MAPK Activation in the Ovary: The Good, the Bad and the Ugly. Mol. Cell Endocrinol. 356, 74–79. doi:10.1016/j.mce.2011.12.005
Fan, H.-Y., Liu, Z., Shimada, M., Sterneck, E., Johnson, P. F., Hedrick, S. M., et al. (2009). MAPK3/1 (ERK1/2) in Ovarian Granulosa Cells Are Essential for Female Fertility. Science 324, 938–941. doi:10.1126/science.1171396
Fan, H.-Y., O'Connor, A., Shitanaka, M., Shimada, M., Liu, Z., and Richards, J. S. (2010). β-Catenin (CTNNB1) Promotes Preovulatory Follicular Development but Represses LH-Mediated Ovulation and Luteinization. Mol. Endocrinol. 24, 1529–1542. doi:10.1210/me.2010-0141
Fan, H.-Y., and Richards, J. S. (2010). Minireview: Physiological and Pathological Actions of RAS in the Ovary. Mol. Endocrinol. 24, 286–298. doi:10.1210/me.2009-0251
Fan, H.-Y., Shimada, M., Liu, Z., Cahill, N., Noma, N., Wu, Y., et al. (2008). Selective Expression ofKrasG12Din Granulosa Cells of the Mouse Ovary Causes Defects in Follicle Development and Ovulation. Development 135, 2127–2137. doi:10.1242/dev.020560
Figueiredo, J. R. d., Cadenas, J., Lima, L. F. d., and Santos, R. R. (2019). Advances in In Vitro Folliculogenesis in Domestic Ruminants. Anim. Reprod. 16, 52–65. doi:10.21451/1984-3143-AR2018-0123
Fujihara, M., Comizzoli, P., Keefer, C. L., Wildt, D. E., and Songsasen, N. (2014). Epidermal Growth Factor (EGF) Sustains In Vitro Primordial Follicle Viability by Enhancing Stromal Cell Proliferation via MAPK and PI3K Pathways in the Prepubertal, but Not Adult, Cat Ovary1. Biol. Reprod. 90. doi:10.1095/biolreprod.113.115089
Ghosh, H. S., McBurney, M., and Robbins, P. D. (2010). SIRT1 Negatively Regulates the Mammalian Target of Rapamycin. PLoS One 5, e9199–8. doi:10.1371/journal.pone.0009199
Gloaguen, P., Crépieux, P., Heitzler, D., Poupon, A., and Reiter, E. (2011). Mapping the Follicle-Stimulating Hormone-Induced Signaling Networks. Front. Endocrin. 2, 1–13. doi:10.3389/fendo.2011.00045
Gorre, N., Adhikari, D., Lindkvist, R., Brännström, M., Liu, K., and Shen, Y. (2014). mTORC1 Signaling in Oocytes Is Dispensable for the Survival of Primordial Follicles and for Female Fertility. PLoS One 9, e110491. doi:10.1371/journal.pone.0110491
Gospodarowicz, D., Ferrara, N., Schweigerer, L., and Neufeld, G. (1987). Structural Characterization and Biological Functions of Fibroblast Growth Factor. Endocr. Rev. 8, 95–114. doi:10.1210/edrv-8-2-95
Gray, A. J., Redmond, A. D., and Martin, M. A. (1987). Use of the Automatic External Defibrillator-Pacemaker by Ambulance Personnel: the Stockport Experience. Br. Med. J. (Clin Res. Ed. 294, 1133. doi:10.1136/bmj.294.6580.1133
Guo, J., Zhang, T., Guo, Y., Sun, T., Li, H., Zhang, X., et al. (2018). Oocyte Stage-specific Effects of MTOR Determine Granulosa Cell Fate and Oocyte Quality in Mice. Proc. Natl. Acad. Sci. USA 115, E5326–E5333. doi:10.1073/pnas.1800352115
Guo, Z., and Yu, Q. (2019). Role of mTOR Signaling in Female Reproduction. Front. Endocrinol. 10, 692. doi:10.3389/fendo.2019.00692
Gupta, P., and Nandi, S. (2012). Isolation and Culture of Preantral Follicles for Retrieving Oocytes for the Embryo Production: Present Status in Domestic Animals. Reprod. Domest. Anim. 47, 513–519. doi:10.1111/j.1439-0531.2011.01904.x
Gupta, P. S. P., Ramesh, H. S., Manjunatha, B. M., Nandi, S., and Ravindra, J. P. (2008). Production of buffalo Embryos Using Oocytes from In Vitro Grown Preantral Follicles. Zygote 16, 57–63. doi:10.1017/S096719940700442X
Han, Y., Luo, H., Wang, H., Cai, J., and Zhang, Y. (2017). SIRT1 Induces Resistance to Apoptosis in Human Granulosa Cells by Activating the ERK Pathway and Inhibiting NF-Κb Signaling with Anti-inflammatory Functions. Apoptosis 22, 1260–1272. doi:10.1007/s10495-017-1386-y
Hao, J., Tuck, A. R., Sjödin, M. O. D., Lindberg, J., Sand, A., Niklasson, B., et al. (2018). Resveratrol Supports and Alpha-Naphthoflavone Disrupts Growth of Human Ovarian Follicles in an In Vitro Tissue Culture Model. Toxicol. Appl. Pharmacol. 338, 73–82. doi:10.1016/j.taap.2017.11.009
Herta, A. C., Lolicato, F., and Smitz, J. E. J. (2018). In Vitro follicle Culture in the Context of IVF. Reproduction 156, F59–F73. doi:10.1530/REP-18-0173
Hikabe, O., Hamazaki, N., Nagamatsu, G., Obata, Y., Hirao, Y., Hamada, N., et al. (2016). Reconstitution In Vitro of the Entire Cycle of the Mouse Female Germ Line. Nature 539, 299–303. doi:10.1038/nature20104
Hornick, J. E., Duncan, F. E., Shea, L. D., and Woodruff, T. K. (2013). Multiple Follicle Culture Supports Primary Follicle Growth through Paracrine-Acting Signals. Reproduction 145, 19–32. doi:10.1530/REP-12-0233
Hosaka, T., Biggs, W. H., Tieu, D., Boyer, A. D., Varki, N. M., Cavenee, W. K., et al. (2004). Disruption of Forkhead Transcription Factor (FOXO) Family Members in Mice Reveals Their Functional Diversification. Proc. Natl. Acad. Sci. 101, 2975–2980. doi:10.1073/pnas.0400093101
Hsieh, M., Lee, D., Panigone, S., Horner, K., Chen, R., Theologis, A., et al. (2007). Luteinizing Hormone-dependent Activation of the Epidermal Growth Factor Network Is Essential for Ovulation. Mol. Cel. Biol. 27, 1914–1924. doi:10.1128/mcb.01919-06
Hsueh, A. J. W., Kawamura, K., Cheng, Y., and Fauser, B. C. J. M. (2015). Intraovarian Control of Early Folliculogenesis. Endocr. Rev. 36, 1–24. doi:10.1210/er.2014-1020
Imai, F., Kishi, H., Nakao, K., Nishimura, T., and Minegishi, T. (2014). IL-6 Up-Regulates the Expression of Rat LH Receptors during Granulosa Cell Differentiation. Endocrinology 155, 1436–1444. doi:10.1210/en.2013-1821
International Federation of Fertility Societies’ Surveillance (2019). “International Federation of Fertility Societies’ Surveillance (IFFS) 2019,” in Global Trends in Reproductive Policy and Practice. 8th Edition. Global Reproductive Health, 4 (1, e29. doi:10.1097/GRH.0000000000000029
Isobe, N., Kitabayashi, M., and Yoshimura, Y. (2005). Microvascular Distribution and Vascular Endothelial Growth Factor Expression in Bovine Cystic Follicles. Domest. Anim. Endocrinol. 29, 634–645. doi:10.1016/j.domaniend.2005.04.007
John, G. B., Shidler, M. J., Besmer, P., and Castrillon, D. H. (2009). Kit Signaling via PI3K Promotes Ovarian Follicle Maturation but Is Dispensable for Primordial Follicle Activation. Dev. Biol. 331, 292–299. doi:10.1016/j.ydbio.2009.05.546
Jones, A. S. K., and Shikanov, A. (2019). Follicle Development as an Orchestrated Signaling Network in a 3D Organoid. J. Biol. Eng. 13, 2–12. doi:10.1186/s13036-018-0134-3
Kreeger, P. K., Fernandes, N. N., Woodruff, T. K., and Shea, L. D. (2005). Regulation of Mouse Follicle Development by Follicle-Stimulating Hormone in a Three-Dimensional In Vitro Culture System Is Dependent on Follicle Stage and Dose1. Biol. Reprod. 73, 942–950. doi:10.1095/biolreprod.105.042390
Laronda, M. M., Rutz, A. L., Xiao, S., Whelan, K. A., Duncan, F. E., Roth, E. W., et al. (2017). A Bioprosthetic Ovary Created Using 3D Printed Microporous Scaffolds Restores Ovarian Function in Sterilized Mice. Nat. Commun. 8, 1–10. doi:10.1038/ncomms15261
Lee, H. N., and Chang, E. M. (2019). Primordial Follicle Activation as New Treatment for Primary Ovarian Insufficiency. Clin. Exp. Reprod. Med. 46, 43–49. doi:10.5653/cerm.2019.46.2.43
Li, M., Liang, C.-G., Xiong, B., Xu, B.-Z., Lin, S.-L., Hou, Y., et al. (2008). PI3-kinase and Mitogen-Activated Protein Kinase in Cumulus Cells Mediate EGF-Induced Meiotic Resumption of Porcine Oocyte. Domest. Anim. Endocrinol. 34, 360–371. doi:10.1016/j.domaniend.2007.10.001
Lin, L., Hron, J. D., and Peng, S. L. (2004). Regulation of NF-Κb, Th Activation, and Autoinflammation by the Forkhead Transcription Factor Foxo3a. Immunity 21, 203–213. doi:10.1016/j.immuni.2004.06.016
Liu, L., Rajareddy, S., Reddy, P., Du, C., Jagarlamudi, K., Shen, Y., et al. (2007, Infertility Caused by Retardation of Follicular Development in Mice with Oocyte-specific Expression of Foxo3a, 134, 199, 209 doi:10.1242/dev.02667
Liu, M., Yin, Y., Ye, X., Zeng, M., Zhao, Q., Keefe, D. L., et al. (2013a). Resveratrol Protects against Age-Associated Infertility in Mice. Hum. Reprod. 28, 707–717. doi:10.1093/humrep/des437
Liu, W.-J., Zhang, X.-M., Wang, N., Zhou, X.-L., Fu, Y.-C., and Luo, L.-L. (2015). Calorie Restriction Inhibits Ovarian Follicle Development and Follicle Loss through Activating SIRT1 Signaling in Mice. Eur. J. Med. Res. 20, 22. doi:10.1186/s40001-015-0114-8
Liu, Z., Castrillon, D. H., Zhou, W., and Richards, J. S. (2013b). FOXO1/3 Depletion in Granulosa Cells Alters Follicle Growth, Death and Regulation of Pituitary FSH. Mol. Endocrinol. 27, 238–252. doi:10.1210/me.2012-1296
Liu, Z., Rudd, M. D., Hernandez-Gonzalez, I., Gonzalez-Robayna, I., Fan, H.-Y., Zeleznik, A. J., et al. (2009). FSH and FOXO1 Regulate Genes in the Sterol/steroid and Lipid Biosynthetic Pathways in Granulosa Cells. Mol. Endocrinol. 23, 649–661. doi:10.1210/me.2008-0412
Long, G.-Y., Yang, J.-Y., Xu, J.-J., Ni, Y.-H., Zhou, X.-L., Ma, J.-Y., et al. (2019). SIRT1 Knock-In Mice Preserve Ovarian reserve Resembling Caloric Restriction. Gene 686, 194–202. doi:10.1016/j.gene.2018.10.040
Long, X., Müller, F., and Avruch, J. (2004). TOR Action in Mammalian Cells and in Caenorhabditis elegans. Curr. Top. Microbiol. Immunol. 279, 115–138. doi:10.1007/978-3-642-18930-2_8
Luz, V. B., Araújo, V. R., Duarte, A. B. G., Celestino, J. J. H., Silva, T. F. P., Magalhães-Padilha, D. M., et al. (2012). Eight-cell Parthenotes Originated from In Vitro Grown Sheep Preantral Follicles. Reprod. Sci. 19, 1219–1225. doi:10.1177/1933719112446072
Magalhães, D. M., Duarte, A. B. G., Araújo, V. R., Brito, I. R., Soares, T. G., Lima, I. M. T., et al. (2011). In Vitro production of a Caprine Embryo from a Preantral Follicle Cultured in media Supplemented with Growth Hormone. Theriogenology 75, 182–188. doi:10.1016/j.theriogenology.2010.08.004
Makker, A., Goel, M. M., and Mahdi, A. A. (2014). PI3K/PTEN/Akt and TSC/mTOR Signaling Pathways, Ovarian Dysfunction, and Infertility: An Update. J. Mol. Endocrinol. 53, R103–R118. doi:10.1530/JME-14-0220
Martelli, A., Russo, V., Mauro, A., O, D. G., Nardinocchi, D., and Mattioli, M. (2017). Insights into Ovarian Follicle Angiogenesis : Morphological and Chronological Vascular Remodeling from Primordial to Ovulating Follicles. SM Vasc. Med., 2 (1), 1009.
Martinez, F., Andersen, C. Y., Barri, P. N., Brannigan, R., Cobo, A., Donnez, J., et al. (2017). Update on Fertility Preservation from the Barcelona International Society for Fertility Preservation-ESHRE-ASRM 2015 Expert Meeting: Indications, Results and Future Perspectives. Fertil. Sterility 108, 407–415. doi:10.1016/j.fertnstert.2017.05.024
Matzuk, M. M., Burns, K. H., Viveiros, M. M., and Eppig, J. J. (2002). Intercellular Communication in the Mammalian Ovary: Oocytes Carry the Conversation. Sci. (80- 296, 2178–2180. doi:10.1126/science.1071965
McLaughlin, M., Albertini, D. F., Wallace, W. H. B., Anderson, R. A., and Telfer, E. E. (2018). Metaphase II Oocytes from Human Unilaminar Follicles Grown in a Multi-step Culture System. Mol. Hum. Reprod. 24, 135–142. doi:10.1093/molehr/gay002
Morita, Y., Wada-Hiraike, O., Yano, T., Shirane, A., Hirano, M., Hiraike, H., et al. (2012). Resveratrol Promotes Expression of SIRT1 and StAR in Rat Ovarian Granulosa Cells: An Implicative Role of SIRT1 in the Ovary. Reprod. Biol. Endocrinol. 10, 14. doi:10.1186/1477-7827-10-14
O'Brien, M. J., Pendola, J. K., and Eppig, J. J. (2003). A Revised Protocol for In Vitro Development of Mouse Oocytes from Primordial Follicles Dramatically Improves Their Developmental Competence1. Biol. Reprod. 68, 1682–1686. doi:10.1095/biolreprod.102.013029
Ordinelli, A., Bernabò, N., Orsini, M., Mattioli, M., and Barboni, B. (2018). Putative Human Sperm Interactome: A Networks Study. BMC Syst. Biol. 12, 1–10. doi:10.1186/s12918-018-0578-6
Orisaka, M., Orisaka, S., Jiang, J.-Y., Craig, J., Wang, Y., Kotsuji, F., et al. (2006). Growth Differentiation Factor 9 Is Antiapoptotic during Follicular Development from Preantral to Early Antral Stage. Mol. Endocrinol. 20, 2456–2468. doi:10.1210/me.2005-0357
Ortega, I., Wong, D. H., Villanueva, J. A., Cress, A. B., Sokalska, A., Stanley, S. D., et al. (2012). Effects of Resveratrol on Growth and Function of Rat Ovarian Granulosa Cells. Fertil. Sterility 98, 1563–1573. doi:10.1016/j.fertnstert.2012.08.004
Palaniappan, M., and Menon, K. M. J. (2012). Luteinizing Hormone/human Chorionic Gonadotropin-Mediated Activation of mTORC1 Signaling Is Required for Androgen Synthesis by Theca-Interstitial Cells. Mol. Endocrinol. 26, 1732–1742. doi:10.1210/me.2012-1106
Paul, S., Andrew, M., Owen, O., and Nitin, S. (1971). Cytoscape: A Software Environment for Integrated Models. Genome Res. 13, 426. doi:10.1101/gr.1239303.metabolite
Qasem, R. J. (2020). The Estrogenic Activity of Resveratrol: a Comprehensive Review of In Vitro and In Vivo Evidence and the Potential for Endocrine Disruption. Crit. Rev. Toxicol. 50, 439–462. doi:10.1080/10408444.2020.1762538
Ramalingam, P. (2016). Morphologic, Immunophenotypic, and Molecular Features of Epithelial Ovarian CancerOncology (Williston Park).
Richards, J. S. (2005). Ovulation: New Factors that Prepare the Oocyte for Fertilization. Mol. Cell Endocrinol. 234, 75–79. doi:10.1016/j.mce.2005.01.004
Robinson, R. S., Woad, K. J., Hammond, A. J., Laird, M., Hunter, M. G., and Mann, G. E. (2009). Angiogenesis and Vascular Function in the Ovary. Reproduction 138, 869–881. doi:10.1530/REP-09-0283
Rosairo, D., Kuyznierewicz, I., Findlay, J., and Drummond, A. (2008). Transforming Growth Factor-β: its Role in Ovarian Follicle Development. Reproduction 136, 799–809. doi:10.1530/REP-08-0310
Rossetto, R., Saraiva, M. V. A., Bernuci, M. P., Silva, G. M., Brito, I. R., Alves, A. M. C. V., et al. (2016). Impact of Insulin Concentration and Mode of FSH Addition on the In Vitro Survival and Development of Isolated Bovine Preantral Follicles. Theriogenology 86, 1137–1145. doi:10.1016/j.theriogenology.2016.04.003
Sanfins, A., Rodrigues, P., and Albertini, D. F. (2018). GDF-9 and BMP-15 Direct the Follicle Symphony. J. Assist. Reprod. Genet. 35, 1741–1750. doi:10.1007/s10815-018-1268-4
Schube, U., Nowicki, M., Jogschies, P., Blumenauer, V., Bechmann, I., and Serke, H. (2014). Resveratrol and Desferoxamine Protect Human OxLDL-Treated Granulosa Cell Subtypes from Degeneration. J. Clin. Endocrinol. Metab. 99, 229–239. doi:10.1210/jc.2013-2692
Shah, J. S., Sabouni, R., Cayton Vaught, K. C., Owen, C. M., Albertini, D. F., and Segars, J. H. (2018). Biomechanics and Mechanical Signaling in the Ovary: a Systematic Review. J. Assist. Reprod. Genet. 35, 1135–1148. doi:10.1007/s10815-018-1180-y
Shen, M., Liu, Z., Li, B., Teng, Y., Zhang, J., Tang, Y., et al. (2014). Involvement of FoxO1 in the Effects of Follicle-Stimulating Hormone on Inhibition of Apoptosis in Mouse Granulosa Cells. Cell Death Dis 5, e1475. doi:10.1038/cddis.2014.400
Stouffer, R. L., Martínez-Chequer, J. C., Molskness, T. A., Xu, F., and Hazzard, T. M. (2001). Regulation and Action of Angiogenic Factors in the Primate Ovary. Arch. Med. Res.doi:10.1016/S0188-4409(01)00323-X
Szklarczyk, D., Franceschini, A., Wyder, S., Forslund, K., Heller, D., Huerta-Cepas, J., et al. (2015). STRING V10: Protein-Protein Interaction Networks, Integrated over the Tree of Life. Nucleic Acids Res. 43, D447–D452. doi:10.1093/nar/gku1003
Taraschi, A., Cimini, C., Capacchietti, G., Ramal-Sanchez, M., Valbonetti, L., Machado-Simoes, J., et al. (2020). Two-Player Game in a Complex Landscape: 26S Proteasome, PKA, and Intracellular Calcium Concentration Modulate Mammalian Sperm Capacitation by Creating an Integrated Dialogue-A Computational Analysis. Int. J. Mol. Sci. 21, 6256. doi:10.3390/ijms21176256
Tata, B. K., Chung, W. C. J., Brooks, L. R., Kavanaugh, S. I., and Tsai, P.-S. (2012). Fibroblast Growth Factor Signaling Deficiencies Impact Female Reproduction and Kisspeptin Neurons in Mice1. Biol. Reprod. 86, 1–7. doi:10.1095/biolreprod.111.095992
Tatone, C., Di Emidio, G., Vitti, M., Di Carlo, M., Santini, S., D’Alessandro, A. M., et al. (2015). Sirtuin Functions in Female Fertility: Possible Role in Oxidative Stress and Aging. Oxidative Med. Cell Longevity 2015, 1–11. doi:10.1155/2015/659687
Telfer, E. E., and Zelinski, M. B. (2013). Ovarian Follicle Culture: Advances and Challenges for Human and Nonhuman Primates. Fertil. Sterility 99, 1523–1533. doi:10.1016/j.fertnstert.2013.03.043
Wayne, C. M., Fan, H.-Y., Cheng, X., and Richards, J. S. (2007). Follicle-stimulating Hormone Induces Multiple Signaling Cascades: Evidence that Activation of Rous Sarcoma Oncogene, RAS, and the Epidermal Growth Factor Receptor Are Critical for Granulosa Cell Differentiation. Mol. Endocrinol. 21, 1940–1957. doi:10.1210/me.2007-0020
Wong, D. H., Villanueva, J. A., Cress, A. B., and Duleba, A. J. (2010). Effects of Resveratrol on Proliferation and Apoptosis in Rat Ovarian Theca-Interstitial Cells. Mol. Hum. Reprod. 16, 251–259. doi:10.1093/molehr/gaq002
Wright, C. S., Hovatta, O., Margara, R., Trew, G., Winston, R. M. L., Franks, S., et al. (1999). Effects of Follicle-Stimulating Hormone and Serum Substitution on the Iin-Vvitro Growth of Human Ovarian Follicles. Hum. Reprod. 14, 1555–1562. doi:10.1093/humrep/14.6.1555
Wu, J., Emery, B. R., and Carrell, D. T. (2001). In Vitro Growth, Maturation, Fertilization, and Embryonic Development of Oocytes from Porcine Preantral Follicles. Biol. Reprod. 64, 375–381. doi:10.1095/biolreprod64.1.375
Wulff, C., Wilson, H., Wiegand, S. J., Rudge, J. S., and Fraser, H. M. (2002). Prevention of Thecal Angiogenesis, Antral Follicular Growth, and Ovulation in the Primate by Treatment with Vascular Endothelial Growth Factor Trap R1R2. Endocrinology 143, 2797–2807. doi:10.1210/endo.143.7.8886
Xiao, S., Coppeta, J. R., Rogers, H. B., Isenberg, B. C., Zhu, J., Olalekan, S. A., et al. (2017). A Microfluidic Culture Model of the Human Reproductive Tract and 28-day Menstrual Cycle. Nat. Commun. 8, 1–13. doi:10.1038/ncomms14584
Xiao, S., Zhang, J., Romero, M. M., Smith, K. N., Shea, L. D., and Woodruff, T. K. (2015). In Vitro follicle Growth Supports Human Oocyte Meiotic Maturation. Sci. Rep. 5, 1–5. doi:10.1038/srep17323
Xu, M., Fazleabas, A. T., Shikanov, A., Jackson, E., Barrett, S. L., Hirshfeld-Cytron, J., et al. (2011). In Vitro Oocyte Maturation and Preantral Follicle Culture from the Luteal-phase Baboon Ovary Produce Mature Oocytes. Biol. Reprod. 84, 689–697. doi:10.1095/biolreprod.110.088674
Xu, M., West, E., Shea, L. D., and Woodruff, T. K. (2006). Identification of a Stage-specific Permissive In Vitro Culture Environment for Follicle Growth and Oocyte Development1. Biol. Reprod. 75, 916–923. doi:10.1095/biolreprod.106.054833
Xu, M., West-Farrell, E. R., Stouffer, R. L., Shea, L. D., Woodruff, T. K., and Zelinski, M. B. (2009). Encapsulated Three-Dimensional Culture Supports Development of Nonhuman Primate Secondary Follicles1. Biol. Reprod. 81, 587–594. doi:10.1095/biolreprod.108.074732
Zhang, H., Risal, S., Gorre, N., Busayavalasa, K., Li, X., Shen, Y., et al. (2014). Somatic Cells Initiate Primordial Follicle Activation and Govern the Development of Dormant Oocytes in Mice. Curr. Biol. 24, 2501–2508. doi:10.1016/j.cub.2014.09.023
Zhang, T., Du, X., Zhao, L., He, M., Lin, L., Guo, C., et al. (2019). SIRT1 Facilitates Primordial Follicle Recruitment Independent of Deacetylase Activity through Directly Modulating Akt1 and mTOR Transcription. FASEB J. 33, 14703–14716. doi:10.1096/fj.201900782R
Zhao, Y., Zhang, Y., Li, J., Zheng, N., Xu, X., Yang, J., et al. (2018). MAPK3/1 Participates in the Activation of Primordial Follicles through mTORC1‐KITL Signaling. J. Cel. Physiol. 233, 226–237. doi:10.1002/jcp.25868
Zhou, X.-L., Xu, J.-J., Ni, Y.-H., Chen, X.-C., Zhang, H.-X., Zhang, X.-M., et al. (2014). SIRT1 Activator (SRT1720) Improves the Follicle reserve and Prolongs the Ovarian Lifespan of Diet-Induced Obesity in Female Mice via Activating SIRT1 and Suppressing mTOR Signaling. J. Ovarian Res. 7, 1–12. doi:10.1186/s13048-014-0097-z
Keywords: ovarian folliculogenesis, computational biology, In vitro folliculogenesis network, hub molecules, bottleneck molecules
Citation: Bernabò N, Di Berardino C, Capacchietti G, Peserico A, Buoncuore G, Tosi U, Crociati M, Monaci M and Barboni B (2021) In Vitro Folliculogenesis in Mammalian Models: A Computational Biology Study. Front. Mol. Biosci. 8:737912. doi: 10.3389/fmolb.2021.737912
Received: 07 July 2021; Accepted: 04 October 2021;
Published: 09 November 2021.
Edited by:
Mauricio Alejandro Latorre, University of Chile, ChileReviewed by:
Pablo Smircich, Instituto de Investigaciones Biológicas Clemente Estable (IIBCE), UruguayVictor Aliaga-Tobar, University of Chile, Chile
Copyright © 2021 Bernabò, Di Berardino, Capacchietti, Peserico, Buoncuore, Tosi, Crociati, Monaci and Barboni. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alessia Peserico, YXBlc2VyaWNvQHVuaXRlLml0
†These authors have contributed equally to this work and share first authorship