Selective YAP activation in Procr cells is essential for ovarian stem/progenitor expansion and epithelium repair
Abstract
Ovarian surface epithelium (OSE) undergoes recurring ovulatory rupture and OSE stem cells rapidly generate new cells for the repair. How the stem cell activation is triggered by the rupture and promptly turns on proliferation is unclear. Our previous study has identified that Protein C Receptor (Procr) marks OSE progenitors. In this study, we observed decreased adherent junction and selective activation of YAP signaling in Procr progenitors at OSE rupture site. OSE repair is impeded upon deletion of Yap1 in these progenitors. Interestingly, Procr+ progenitors show lower expression of Vgll4, an antagonist of YAP signaling. Overexpression of Vgll4 in Procr+ cells hampers OSE repair and progenitor proliferation, indicating that selective low Vgll4 expression in Procr+ progenitors is critical for OSE repair. In addition, YAP activation promotes transcription of the OSE stemness gene Procr. The combination of increased cell division and Procr expression leads to expansion of Procr+ progenitors surrounding the rupture site. These results illustrate a YAP-dependent mechanism by which the stem/progenitor cells recognize the murine ovulatory rupture, and rapidly multiply their numbers, highlighting a YAP-induced stem cell expansion strategy.
Editor's evaluation
This is a well conducted and interesting study that increases our knowledge of mechanisms governing healing after ovarian rupture. This study specifically demonstrates the importance of Procr+ progenitors that are positively regulated by YAP signalling.
https://doi.org/10.7554/eLife.75449.sa0Introduction
During the adult reproductive cycles, the ovarian surface epithelium (OSE) undergoes recurring ovulatory rupture and repair (Gaytán et al., 2005; Auersperg et al., 2001). After ovulation, to maintain the physiological function and morphology of the ovary, the wound is completely closed within 12 hr to 3 days following rupture (Burdette et al., 2006; Ng and Barker, 2015; Tan and Fleming, 2004). Cells surrounding the damaged sites are required to respond to the wound by turning on cell proliferation to supply sufficient cells as building blocks for regeneration (Wang et al., 2017). Our previous study has identified that Procr+ OSE stem/progenitor cells are the major contributor for ovulatory rupture repair. Targeted ablation of these cells hampers the repair (Wang et al., 2019). Interestingly, we observed that Procr+ cells expand instantly upon ovulation, reminiscent of a result of symmetric division (Wang et al., 2019). It remains unknown how the stem cell is triggered by the ovulation event, and what is the signal that links the rupture to the instant stem cell expansion.
The cue for this stem/progenitor cell amplification likely comes from a particular extracellular signal occurring upon ovulation. One possibility is that the follicular fluid expelled during ovulation consists of Wnts and other potential niche signals (Ahmed et al., 2006; Boyer et al., 2010; Nilsson et al., 2001; Parrott et al., 2000), which may regulate Procr+ stem/progenitor cell expansion. Another possibility is the involvement of mechanical force-induced signals, either through the change of adherent junctions or via a directed mechanism during ovulation, resulting in Procr+ stem/progenitor cell expansion.
YAP (Yes-associated protein, also known as YAP1) signaling is an evolutionarily conserved pathway and a master regulator of organ size and tissue growth during animal development (Moya and Halder, 2019). As a downstream effector, YAP is critical for regeneration in different organs, through triggering cell proliferation, cell survival, or expansion of stem and progenitor cell compartments (Johnson and Halder, 2014; Lin et al., 2017b; Moya and Halder, 2016; Panciera et al., 2016; Patel et al., 2017; Xiao et al., 2016; Zhang and Del Re, 2017). YAP is a transcriptional coactivator protein that shuttles between the cytoplasm and nucleus, and regulates the expression of target genes, such as Ccn1 and Birc5, through binding with TEAD transcription factors (Chen et al., 2001; Lin et al., 2017a; Piccolo et al., 2014; Totaro et al., 2018; Zhao et al., 2008). Vgll4, a member of Vestigial-like proteins, serves as a transcriptional repressor of YAP through direct interactions with TEADs (Chen et al., 2004). Previous studies from us and others have demonstrated the important roles of Vgll4 plays during development and regeneration in various tissues (Feng et al., 2019; Suo et al., 2020; Lin et al., 2016; Yu et al., 2019). Cell–cell junctions links cells to each other in epithelial tissues, and is an upstream negative regulator of YAP (Ramos and Camargo, 2012; Yang et al., 2015). Mechanical forces regulate cell–cell adhesion stability, and cell–cell adhesion junctions may be intrinsically weak at high forces (Pinheiro and Bellaϊche, 2018). It has been shown that disruption of adherent junctions turns on YAP nuclear activities in lung stem/progenitor cells (Zhou et al., 2018). However, whether YAP signaling is implicated in ovulatory rupture repair is unknown.
In this study, we investigated how OSE stem/progenitor cells are triggered by the rupture postovulation and divide subsequently. We found that, in the proximity of rupture site, decreased adherent junction is associated with increased incidence of Yap1 nuclear localization in OSE cells. Interestingly, only Procr+ OSE cells displayed a low level of Vgll4, allowing YAP signaling activation, and conditional deletion of Yap1 in Procr+ cells hampers OSE repair. We generated a new TetO-Vgll4 mouse. Ectopic expression of Vgll4 in the stem/progenitor cells using Procr-rtTA;TetO-Vgll4 mice blocked OSE ovulatory repair. Moreover, we found that YAP signaling activation resulted in Procr+ cells expansion at the rupture site, through the combination of inducing cell division, and directly activating Protein C Receptor (Procr) transcription. The activation of Procr is essential, as when Procr was deleted, stemness property was lost and OSE repair was hindered.
Results
Decreased E-cadherin expression at the rupture site and selective activation of YAP signaling in Procr+ cells
To investigate what could be the potential extracellular stimuli at the rupture site, we performed immunostaining of various adherent or tight junction components on ovarian sections. To increase rupture incidences, superovulation was induced by injection of pregnant mare serum gonadotropin (PMSG) and human chorionic gonadotropin (HCG), and the ovaries were harvested at 0.5 days after HCG injection when ovulation just occurred (Figure 1—figure supplement 1A). Interestingly, we found that E-cadherin, α-catenin, and ZO-1 staining is markedly decreased at the proximal region of rupture (defined as within 20 cells on one side of the rupture in section) compared to other regions, that is, rupture distal region (Figure 1A and Figure 1—figure supplement 1C) and nonrupture region (Figure 1—figure supplement 1B,D). As adherent junction has been implicated as a modulator of YAP signaling (Kim et al., 2011; Schlegelmilch et al., 2011; Varelas et al., 2010; Yang et al., 2015), we examined YAP activities at the rupture area by immunostaining. We observed an increased incidence of nuclear Yap1 at the proximal region of rupture compared to other regions (Figure 1B–D). These results suggest that compromised adherent junctions resulting from ovulatory rupture are associated with Yap1 nuclear localization in OSE cells surrounding the wound.
Our previous study has established that Procr+ progenitor cells surrounding the wound instantly proliferate upon rupture and are responsible for OSE repair (Wang et al., 2019). We therefore investigated whether Procr+ cells close to the rupture site are associated with YAP signaling activities. We performed Yap1 immunostaining using Procr-rtTA;tetO-H2B-GFP reporter, in which H2B-GFP signal marks Procr-expressing cells. Superovulation was induced in these animals by PMSG and HCG injections, and ovaries were harvested 0.5 days after HCG injection. We found that Procr+ (H2B-GFP+) cells at the rupture proximal region (referred to as rupture site from here on) have significantly higher nuclear Yap1 staining (75.9% ± 1.7%) compared to Procr− cells (39.6% ± 1.0%) (Figure 1E, G), or compared to Procr+ cells at the nonrupture region (Figure 1F, G). This was further validated by RNA double in situ hybridization with Procr and a YAP target gene Ccn1. We found that, at the rupture site, Ccn1 is preferentially activated in Procr+ OSE cells, with 44.8% ± 2.9% of Procr+ cells being Ccn1+, which is markedly higher than that of Procr− cells (8.9% ± 0.8%) (Figure 1H, J). At the nonrupture site, both Procr+ and Procr− cells had a rather low expression of Ccn1 expression (Figure 1I, J). Furthermore, we FACS-isolated Procr+ and Procr− OSE cells from the rupture sites (Figure 1—figure supplement 2A), validated the epithelial identity by Krt8 (K8) staining (Figure 1—figure supplement 2B). qPCR analysis indicated that Procr+ cells at rupture sites have a higher level of Birc5 expression compared to Procr− cells (Figure 1—figure supplement 2C). Together, these results suggest that YAP signaling was specifically activated in Procr+ cells at the rupture site. Considering the role of YAP signaling in promoting cell proliferation, these results are in line with our previous observations that only Procr+, but not Procr−, cells at the rupture site displayed increased proliferation (Wang et al., 2019).
Deletion of Yap1 in Procr+ cells hinders rupture repair and progenitor proliferation
To investigate whether YAP signaling is important for OSE repair, we deleted Yap1 specifically in Procr+ cells using Procr-CreER;Yap1fl/fl mice (Yap1-cKO). Yap1fl/fl mice was used as control (Ctrl). Tamoxifen (TAM) was administered in 4-week-old mice, following by superovulation at 2 days after TAM injection (Figure 2A). The impact on OSE repair by Yap1 deletion was analyzed by ovary whole-mount imaging. At 4.5 dpi (ovulation), the two groups had similar ruptures (Figure 2B, E). At 6 dpi, Ctrl ovaries underwent rapid repairing (Figure 2C, E), and the OSE was completely recovered by 7.5 dpi (Figure 2D, E). In contrast, the OSE repair in Yap1-cKO ovaries was significantly delayed at both 6 and 7.5 dpi (Figure 2C–E). The efficacy of Yap1 deletion and the reduced expression of the target gene Ccn1 in OSE cells were validated by qPCR analyses (Figure 1—figure supplement 2D).
To analyze the proliferative capacity of Procr+ OSE cells, mice were subjected to 12 hr of 5-ethynyl-29-deoxyuridine (EdU) incorporation before harvesting the ovaries (Figure 2A). When analyzed at 4.5 dpi (ovulation), the number of proliferating OSE cells at rupture site (defined as 20 cells on one side from the opening) was significantly decreased from 5.5 ± 0.3 EdU+ in Ctrl to 1.4 ± 0.2 EdU+ in Yap1-cKO (Figure 2F–H). The impact to cell proliferation was further analyzed in vitro. Our previous study has established that Procr+, but not Procr−, OSE cells can form colonies in vitro (Wang et al., 2019). At 4.5 dpi, total OSE cells were isolated from both Ctrl and Yap1-cKO mice (Figure 2I), and placed in culture as previously described (Wang et al., 2019). Deletion of Yap1 in Procr+ cells drastically inhibited OSE colony formation (Figure 2J, K).
To visualize the contribution of Procr+ progenitors toward the repair in the presence or absence of Yap1, we performed in vivo lineage tracing. TAM was administered to 4-week-old mice to simultaneously delete Yap1 and initiate lineage tracing in Procr+ cells (Figure 2L). At 4.5 dpi, control (Procr-CreER;R26-mTmG) ovary displayed a zone of mGFP+ cells that are the progeny of Procr+ progenitors surrounding the rupture sites (Figure 2M). In contrast, Yap1-cKO (Procr-CreER;Yap1fl/fl;R26-mTmG) ovaries have markedly fewer mGFP+ cells around the wound (Figure 2M, O), supporting the notion that the activity of Procr+ progenitors was hampered at the beginning of the repairing process. At 7 dpi, control ovaries had generated patches of mGFP+ cells covering the newly formed corpus luteum (Figure 2N). However, Yap1-cKO ovaries still had obvious openings with few mGFP+ cells (Figure 2N, O). Together, these results suggest that YAP signaling activation is crucial for the proliferation of Procr+ progenitor cells and the timely repair of OSE after rupture.
An intrinsic lower level of Vgll4 in Procr+ cells is essential for their progenitor property and OSE rupture repair
Next, we investigated what could be the reason that YAP signaling is specifically activated in Procr+ cells. Vgll4 is a negative regulator of YAP by inhibiting the binding of YAP and TEAD4 (Feng et al., 2019). We FACS-isolated Procr+ cells and Procr− cells from the rupture sites (Figure 3A). qPCR analysis indicated that Procr+ cells have lower level of Vgll4 compared to Procr− cells (Figure 3B). This was further validated by Vgll4 immunostaining using Procr-rtTA;tetO-H2B-GFP reporter, in which H2B-GFP signals mark Procr-expressing cells. Consistent with the qPCR results, Procr+ cells also exhibited lower Vgll4 protein expression compared to Procr− cells (Figure 3C–E).
To examine whether the reduced level of Vgll4 is significant for the selective YAP signaling activation in Procr+ cells and rupture repair, we set to overexpress Vgll4 specifically in Procr+ cells. A new TetO-Vgll4 mouse line was generated, by inserting a tetO-Vgll4-Flag-wpre-polyA cassette behind the 3′UTR of the Col1a1 gene (Figure 3F and Figure 3—figure supplement 1A–C). Subsequently, Procr-rtTA;TetO-Vgll4 (Vgll4-OE) mice were generated by genetic crosses with TetO-Vgll4 as control (Ctrl). The efficacy of overexpression was validated by western blotting and qPCR, showing increased expression of Vgll4 and decreased expression of Ccn1 in Vgll4-OE cells (Figure 3—figure supplement 1D, E). Furthermore, immunostaining confirmed the increased number of Vgll4 high cells in the OSE layer of Vgll4-OE mice (Figure 3G). For this experiment, superovulation was performed to 4-week-old mice and doxycycline hyclate (DOX) was fed throughout the process (Figure 3H). The impact of Vgll4 overexpression was analyzed throughout the repairing process, at 4.5 dpi (ovulation), 6 dpi (OSE repair ongoing), and 7.5 dpi (OSE repair completed) by ovary whole-mount imaging. We found that the rupture in Ctrl and Vgll4-OE ovaries are comparable at 4.5 dpi (Figure 3I, L). At 6 dpi, while Ctrl ovaries had sights of wound closure, Vgll4-OE ovaries still showed larger areas of rupture (Figure 3J, L). At 7.5 dpi, Ctrl ovaries displayed complete OSE, whereas the repair in Vgll4-OE ovaries was obviously delayed (Figure 3K, L).
Next, we examined whether overexpression of Vgll4 affects progenitor proliferation. At 4.5 dpi (ovulation), the number of proliferating OSE cells at rupture site was significantly decreased from 3.7 ± 0.3 in Ctrl to 1.0 ± 0.2 in Vgll4-OE (Figure 3M–O). At 4.5 dpi, total OSE cells were isolated and cultured in vitro for 7 days (Figure 3P). Consistently, overexpression of Vgll4 inhibits cell proliferation and colony formation (Figure 3Q, R).Together, these results suggest that overexpression of Vgll4 in Procr+ cell impaired Procr+ cell proliferation and ovulatory rupture repair.
YAP signaling promotes Procr+ cells expansion at rupture site
We have previously found that Procr+ progenitor cells expand instantly at the periphery of the rupture site upon ovulation (Wang et al., 2019). To investigate whether YAP signaling activation is linked to the expansion of Procr+ progenitor cells, TAM was administered to 4-week-old Procr-CreER;Yap1fl/fl (Yap1-cKO) and Yap1fl/fl (Ctrl) mice for two times, followed by superovulation. At 4.5 dpi (ovulation), FACS analysis showed a dramatic decrease of Procr+ progenitor population when Yap1 was deleted (Figure 4A–C).
To better visualize the change ofProcr+ progenitor cells under the influence of YAP signaling, we generated Procr-rtTA;TetO-H2B-GFP+/−;TetO-Vgll4+/− mice (Vgll4-OE). Superovulation was performed to 4-week-old mice and DOX was fed throughout the experiments to maintain the expression of H2B-GFP in Procr+ cells (Figure 4D). When analyzed at 4.5 dpi (ovulation), at the wound edge (defined as 20 cells on one side from the opening) of control ovary (Procr-rtTA;TetO-H2B-GFP+/−), there were about 7.4 ± 0.3 H2B-GFP+ cells expressing the peak level of GFP (Figure 4E, G). In contrast, in Vgll4-OE ovary (Procr-rtTA;TetO-Vgll4;TetO-H2B-GFP+/−), only 2.4 ± 0.2 H2B-GFP+ cells were observed at the wound edge (Figure 4F, G). FACS analysis also showed that the percentage of Procr+ progenitor population decreased significantly from 18.5% ± 1.0% in Ctrl to 7.4% ± 1.2% in Vgll4-OE at ovulation stage (Figure 4H, I).
The proliferative activity of Procr+ cells was further evaluated in vitro. We isolated OSE cells from control (Procr-rtTA;TetO-H2B-GFP+/−) and Vgll4-OE (Procr-rtTA;TetO-Vgll4;TetO-H2B-GFP+/−) mice and placed in culture, followed by live imaging to document the division of H2B-GFP+ (Procr+) cells (Figure 4—figure supplement 1A, B). In control cells, we observed frequent division of Procr+ cells, and in most cases, it was one Procr+ cell dividing into two Procr+ cells (Figure 4—figure supplement 1A, C). But in Vgll4-OE, we could hardly observe cell division (Figure 4—figure supplement 1B, C). Together, these results suggest that inhibition of YAP signaling, by either Yap1-deletion or Vgll4-OE, impairs the expansion of Procr+ progenitors upon rupture.
YAP signaling enhances Procr expression
It is unclear how YAP maintains Procr expression during or after cell division. Thus, we investigated the association of YAP activation and Procr expression. OSE cells were isolated from Procr-rtTA;TetO-H2B-GFP+/− mice, and cultured on glass (YAP activation) or soft condition (0.48 kPa, YAP inactivation) (Figure 4—figure supplement 2A, B). DOX was added 2 days before harvest. Consistent with the notion, we found that, in soft condition, Yap1 was mostly cytoplasmic and most OSE cells are H2B-GFP− (Figure 4—figure supplement 2A). In contrast, most OSE cells are H2B-GFP+ in stiff condition and Yap1 was found in the nucleus (Figure 4—figure supplement 2B). These observations suggest that YAP activation might induce Procr expression. We verified by qPCR that Procr expression is upregulated in stiff conditions (Figure 4—figure supplement 2C). Our results support the notion that YAP activation induces Procr expression.
To further investigate whether Yap1 regulates Procr expression, we knocked down Yap1 by shRNA in OSE culture and found that this inhibits Procr expression (Figure 4—figure supplement 2D, E). Furthermore, blocking YAP activation by Verteporfin (VP) or Vgll4 overexpression also resulted in lower Procr expression (Figure 4—figure supplement 2F, G). These results suggest that inhibiting YAP signaling suppresses Procr expression.
To investigate whether YAP/TEAD4 directly regulate Procr expression, we analyzed the promoter of Procr. A Tead4-binding motif (5′-CATTCC-3′) was found at the proximal promoter of Procr (−1486 to −1481 bp) (Figure 4J). Chromatin immunoprecipitation-qPCR (ChIP-qPCR) showed that Tead4 could directly bind to the Procr promoter (Figure 4K). Therefore, we examined whether this Tead4-binding motif is responsible for induction of Procr expression by Yap1. While Yap1 induced the wild-type promoter luciferase in a dose-dependent manner (Figure 4L and Figure 4—figure supplement 2H, I), it could not activate the mutant reporter with the deletion of the Tead4-binding motif (Figure 4M and Figure 4—figure supplement 2H, I). These results suggest that Yap1 directly promotes Procr expression. Together, our data support a model that YAP signaling promotes expansion of Procr+ cells at rupture site through a combination of increased cell division and Procr expression (Figure 4N).
Procr is essential for the progenitor property
The upregulation of Procr expression coupled with YAP-induced cell division implies that the expression of Procr may be important for keeping the stem cell property in OSE. To assess the significance of Procr, we utilized a Procr-flox allele (Liu and Zeng, unpublished) and specifically deleted Procr in the progenitor using ProcrCreER/fl (Procr-cKO) mice. TAM was administered in 4-week-old mice for two times, followed by superovulation at 2 days after TAM injection (Figure 5A), and the phenotype was analyzed by ovary whole-mount imaging. Ctrl and Procr-cKO ovaries formed comparable ruptures at 4.5 dpi (ovulation) (Figure 5B, E). At 6 dpi (OSE repair ongoing), control ovaries showed smaller openings compared to Procr-cKO (Figure 5C, E). At 7.5 dpi (OSE repair completed), Ctrl ovaries were covered by complete OSE, whereas Procr-cKO ovaries still had regions with unrepaired OSE (Figure 5D, E). Furthermore, at 4.5 dpi (ovulation), the ovaries were harvested after 12 hr of EdU incorporation. The number of proliferated OSE cells at rupture site decreased from 4.7 ± 0.4 cells in control (Procrfl/+) to 1.6 ± 0.2 cells in Procr-cKO (Figure 5F–H). After deletion of Procr in vivo, total OSE cells were isolated and cultured in vitro for 7 days (Figure 5I). We found that deletion of Procr inhibits the proliferation of progenitor cells, resulting in reduced colony sizes (Figure 5J–L). Overall, these data suggest that Procr is essential for progenitor property upon rupture.
Discussion
In this study, we addressed the molecular mechanism which links the ovulatory rupture to OSE stem/progenitor cells activation, thus promptly turning on proliferation and repairing the wound. Our findings support the following model. Procr+ OSE progenitors have intrinsically lower levels of Vgll4. Upon ovulatory rupture, the decreased adherent junction at the proximity of the rupture site promotes Yap1 nuclear localization. These intrinsic and extrinsic factors together lead to YAP signaling activation in Procr+ progenitors around the wound, which sequentially stimulates the proliferation of the progenitors. Importantly, YAP activation directly upregulates Procr expression in the dividing cells, resulting in the expansion of Procr+ progenitors around the wound (Figure 5M). Blocking YAP signaling in the progenitors by Yap1-cKO or Vgll4-OE impairs the progenitors’ activities and hinders OSE repair. Furthermore, Procr function is essential for these progenitors. When Procr was deleted, stem cell property was lost hindering OSE repair.
While we uncovered the significance of selective activation of YAP in OSE progenitors, it is still unclear how ovarian rupture is sensed and how YAP signaling is induced by injury. We observed reduced cell-adhesion protein E-cadherin is lowered at the site of the ovarian rupture, and it is previously known that E-cadherin signals through the Hippo pathway to block YAP (Yang et al., 2015). The current missing link is why and how E-cadherin would be reduced at the ovarian rupture site. We speculate that during the late stage of follicle development, the pre-ovulatory follicle forms a protrusion toward OSE. Subsequently, ovulation generates a rupture on OSE. These contiguous events likely induce the thinning of OSE surrounding the pre-ovulatory follicles and at the proximity of the rupture site, resulting in the reduced adherent junction proteins. Yet, there are other potentially more direct possibilities, i.e. mechanical stretching induces YAP (Halder et al., 2012). First, the pre-ovulatory follicle protrusion or the release of oocytes induces a mechanical force on the OSE surrounding the wound, activating YAP signaling. Second, at the rupture site, epithelial cells no longer become packed together because the epithelium has been denuded, therefore potentially cells flanking the rupture site would become ‘stretched’, consequently actives YAP.
YAP signaling promotes Procr+ cell expansion at the rupture site through a combination of increased cell division and Procr expression. In the current study, YAP is particularly activated in Procr+ progenitor cells at the rupture site. We observed that at the rupture sites, Vgll4 is highly expressed in Procr− cells, preventing YAP pathway activation in those cells around the rupture sites. Our findings demonstrate that the reduced levels of Vgll4 in Procr+ progenitors likely contribute to the selective activation of YAP signaling in these cells. Further study should investigate what mechanism determines the lower expression of Vgll4 in Procr+ progenitor cells.
In the current study, we generated a new TetO-Vgll4 mouse model that enables the overexpression of Vgll4 in a specific cell type. The overexpression of Vgll4 in the progenitor of OSE has been validated using Procr-rtTA. The advantages brought by our TetO-Vgll4 reporter will be of broad value in studies of Hippo-Yap signaling across all tissues.
Procr expression is initially found on the surface of vascular cells exerting an anticoagulation role, by binding and activating protein C (PC) in the extracellular compartment (Fukudome et al., 1998). More recently, studies from us and others have identified Procr as a stem cell surface marker in multiple tissues (Wang et al., 2015; Wang et al., 2020; Wang et al., 2019; Yu et al., 2016), but less is known regarding the function of Procr in stem/progenitor cells. In the current study, we demonstrate that, Procr is essential for the proliferation of Procr+ progenitor cells and OSE repair upon rupture. Our previous report indicated that PROCR concomitantly activates multiple pathways including ERK, PI3K–Akt–mTOR, and RhoA–Rock–P38 signaling in breast cancer cells (Wang et al., 2018). We speculate that similar intracellular pathways might be involved in the Procr+ OSE cells. Procr is regarded as a Wnt target gene from an in vitro screen in mammary stem cell culture (Wang et al., 2015). In this study, we identify YAP as a novel upstream regulator of Procr. ChIP-qPCR and promoter luciferase experiments demonstrate that Procr transcription can be directly upregulated by YAP activation.
The phenomena of YAP promoting stem/progenitor cell expansion have been reported in various tissues (Beverdam et al., 2013; Camargo et al., 2007; Cao et al., 2008; Ramos and Camargo, 2012; Schlegelmilch et al., 2011; Zhou et al., 2018). Yet, in this process, less is known about how YAP maintains stem cell properties. To the best of our knowledge, this is the first report illustrating a mechanism through which YAP promotes cell proliferation, and simultaneously upregulates the expression of an essential stemness gene to maintain cell fate, leading to a rapid expansion of stem cell numbers around the wound. In summary, our study provides new evidence and molecular insights into how ovulatory rupture triggers the activation of OSE stem cells, resulting inpromptly expanding their numbers for repair. This may have a broad implication to understand the action of tissue stem cells during would healing in other tissue.
Materials and methods
Lead contact and materials availability
Request a detailed protocolFurther information and requests for reagents should be directed to and will be fulfilled by the Lead Contact, Yi Arial Zeng (yzeng@sibcb.ac.cn). All unique/stable reagents generated in this study are available from the Lead Contact with a completed Materials Transfer Agreement.
Experiment animals
Request a detailed protocolTetO-H2B-GFP+/− (Stock: 005104), R26-mTmG+/− (Stock: 007576) from Jackson Laboratories, ProcrCreER (Wang et al., 2015), ProcrrtTA (Wang et al., 2019), YAPfl/+ (Feng et al., 2019), Procrfl/+(Liu and Zeng, unpublished), TetO-Vgll4 were used in this study. The TetO-Vgll4 mouse line was generated by knocking in a cassette of TetO-Vgll4-Flag-wpre-polyA behind 3′UTR of Col1a1 gene (Figure 3—figure supplement 1). All mice were housed in the SIBCB animal facility under IVC standard with a 12 hr light/dark cycle at room temperature. Both ovaries were used per mice and the number of mice per experiment was shown in figure legends. For targeted knockout in vivo, 4–5 weeks mice were administered with TAM diluted in sunflower oilby intraperitoneal (IP) injection at a concentration of 2 mg per 25 g body weight for two or three times (on every second day). For superovulation experiments, 4- to 5-week-old mice were injected with 10 IU of PMSG by IP, followed by IP injection of 10 IU of HCG about 48 hr later. For DOX feeding, DOX was dissolved in drink water at a concentration of 1 mg/ml. Experimental procedures were approved by the Animal Care and Use Committee of Shanghai Institute of Biochemistry and Cell Biology, Chinese Academy of Sciences, with a project license number of IBCB0065.
OSE cells isolation and flow cytometry
Request a detailed protocolOvaries from superovulated or 4- to 12-week-old female mice were isolated, and the oviduct and bursa were carefully cleared out under dissect microscope. The ovaries were minced into pieces as small as possible, and then placed in 10 ml digest buffer (RPMI 1640 with 5% fetal bovine serum (FBS), 1% penicillin–streptomycin, 25 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES), and 300 U/ml collagenase IV). After digestion at 37°C, 100 rpm for about 1 hr, ovarian cells were obtained after centrifugation at 1000 rpm for 5 min. The red blood cells were lysed with buffer at room temperature for 5 min, and then single cells were obtained with 0.25% trypsin treatment at 37°C for 5 min, followed by 0.1 mg/ml DNaseI incubation at 37°C for 5 min with gently pipetting before filtering through 70 μm cell strainers. The single cells were incubated on the ice and in dark with the following antibodies at a dilution of 1:200: FITC conjugated, PE conjugated, or biotinylated CD31, CD45, EpCAM-APC, Procr-PE, Procr-Biotin, Streptavidin-APC-Cy7, and Streptavidin-V450. All analysis and sorting were performed using a FACSJazz (Becton Dickinson). The purity of sorted population was routinely checked and ensured to be >95%.
OSE cells 3D culture assay
Request a detailed protocolFACS sorted OSE cells were resuspended with 60 μl 100% growth factor-reduced Matrigel and placed around the rim of a well of a 24-well plate, and allowed to solidify for at least 15 min at 37°C in a 5% CO2 incubator before adding 0.5–1 ml culture medium. Colonies were grown for 7–9 days and the medium was changed every other days. The culture medium was prepared by adding 5% FBS, 4 mM L-glutamine, 1 mM sodium pyruvate, 10 ng/ml epidermal growth factor, 500 ng/ml hydrocortisone, 5 mg/ml insulin, 5 mg/ml transferrin, 5 ng/ml sodium selenite, 0.1 mM Minimum Essential Medium (MEM) nonessential amino acids, 10−4 M 2-mercaptoethanol into Dulbecco’s modified essential medium (DMEM)/F12. The organoid images were captured by Zeiss inverted microscope at days 7–9.
Immunohistochemistry
Request a detailed protocolFor section staining, ovarian tissues were fixed in 4% PFA at room temperature for 15 min, following by washed with phosphate-buffered saline (PBS) for three times, dehydrated in 30% sucrose at 4°C overnight and embedded with Optimum Cutting Temperature. 16–18 μm tissue sections were incubated in 0.1% or 0.5% Triton X-100 diluted with PBS (PBST) for 20 min and then 1 hr blocking using 10% FBS in PBST. Then sections were incubated with primary antibodies diluted in blocking buffer at 4°C overnight, followed by washes for three times (20 min per time). After wash, the sections were further incubated with secondary antibodies and 4′,6-diamidino-2-phenylindole (DAPI) diluted in blocking buffer for 2 hr at room temperaturein dark, followed by washes for three times (20 min per time) and mounted with mounting medium.
For staining of cultured colonies, colonies were released from Matrigel by incubating with dispase for 20–30 min. Then the colonies were fixed in 4% PFA on ice for 10 min, followed by cytospin (Thermo Fisher) into slides and staining protocol described above.
For whole mouse ovary immunohistochemistry, mouse ovaries that cleared without bursa and oviduct were fixed with fresh 4% PFA at room temperature for 15 min in 4 ml Eppendorf tubes, followed by washing with PBST for three times (20 min per time). The staining of whole ovaries was then transferred into the 2 ml Falcon tubes using a dropper carefully. Ovaries were blocked for 1 hr using 10% FBS in PBST. Then, the ovaries were incubated with primary antibodies diluted in blocking buffer at 4°C for 48 hr on a transference shaker with 10rpm, followed by washing for three times (20min per time) at room temperature. After washing, the ovaries were incubated with secondary antibodies diluted in blocking buffer for 24 hr at 4°C in dark, and counterstained with DAPI on a transference shaker with 10rpm, followed by washing for three times (20min per time) at room temperature. The ovaries could be stored in PBST at 4°C for at least 2 weeks.
For Yap1 staining in vivo, tyramide signal amplification assay (TSA staining) with Yap1 antibody from CST (Cat# 14074) was used. Briefly, paraffin sections were rehydrated in histoclear and gradual ethanol (100%, 100%, 95%, 85%, 75%, 50%, and 30%) and the TSA staining was performed using the Opal 4-Color Automation IHC Kit (PerkinElmer) following the manufacturer’s instructions. After TSA staining for Yap1, staining for GFP and Krt8 was performed following protocol described above.
Tissue sections and organoids fluorescent images were captured using Leica DM6000 TCS/SP8 laser confocal scanning microscope with a ×20/0.75 or ×40/0.75 or ×63/0.75 IMM objective with 1–3 μm z-step. Confocal images were processed with maximum intensity projections.
Whole mouse ovarian fluorescent images were captured with inverted Leica TCS SP8 WLL at a ×10/0.75 objective, z-stack was ~50–80 layers with 6–7 μm per layer, and the area was about 1.5 mm × 1.5 mm, which was about 1/6–1/4 of the ovary surface.
Western blotting
Request a detailed protocolDigested cells were lysed in sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS–PAGE) loading buffer and boiled for 10 min. Proteins were separated by SDS–PAGE and transferred to nitrocellulose membrane (GE Company). Bolts were blocked with 3% Bovine Serum Albumin (BSA) in Tris-Buffered Saline with 0.5% Tween 20 (TBST) (50 mM Tris–HCl, 150 mM NaCl, 0.05% Tween-20, pH 7.5) for 1 hr and incubated with primary antibodies at 4°C overnight, followed by incubated with secondary IgG-HRP antibodies for 2 hr at room temperature. Protein bands were visualized with chemiluminescent reagent and exposed to Mini Chemiluminescent Imager.
RNA in situ
Request a detailed protocolIn situ hybridization was performed using the RNA scope kit (Advanced Cell Diagnostics) following the manufacturer’s instructions. Procr probes (REF#410321) and Ccn1 probes (REF#429001) were ordered from Advanced Cell Diagnostics. After in situ hybridization, TSA method was used for Krt8 staining following the manufacturer’s instructions using the Opal 4-Color Automation IHC Kit (PerkinElmer). The images were captured using Leica DM6000 TCS/SP8 laser confocal scanning microscope with a ×63/0.75 IMM objective.
EdU labeling assays
Request a detailed protocolThe proliferation of OSE cells in vivo was measured by EdU uptake. Briefly, mice were injected with 100 μl EdU (2.5 mg/ml in dimethyl sulfoxide) for 12 hr. Then ovaries were harvested for section, following by EdU color staining using Click-iTEdU Alexa Fluor Imaging Kit (prepared according to the manufacturer’s instructions). After washed with PBS for three times (10 min per time), EdU color development was performed following the manufacturer’s protocol. After EdU signal developing, sections were blocked in blocking buffer for 1 hr at room temperature followed by antibody staining and mounted with mounting medium for imaging and quantification.
Living image of cultured OSE cells
Request a detailed protocolOSE cells were isolated from the mice and cultured on glass for 3–4 days. DOX was added into the medium 1 day and Hoechst 33342 was added 30 min before image. Live-cell imaging was performed at 37°C on a Zeiss Cell discoverer seven with perfect focus system. Cells were imaged at 1 time per 5 min for 24 hr with 70% laser power.
Chromatin immunoprecipitation-qPCR
Request a detailed protocolCultured primary OSE cells were crosslinked in a final concentration of 1% formaldehyde (Sigma) PBS buffer for 15 min at 37°C, then added glycine to stop crosslinking. Chromatin from nuclei was sheared to 200–600 bp fragments using ultrasonic apparatus, then immunoprecipitated with antibody of TEAD4 (ab58310, Abcam) or normal mouse IgG (sc-2025, Santa Cruz) overnight. Antibody/antigen complexes were recovered with Protein A/G PLUS-Agarose (sc-2003, Santa Cruz Biotechnology) for 2 hr at 4°C. After washing, the chromatin was eluted, decrosslinked and digested. The immunoprecipitated DNA was collected with QIAQIUCK PCR Purification Kit (QIAGEN). Purified DNA was performed with ChIP-qPCR. Assessing the enrichment of the proteins of interest on the targeting region by calculating the value of ‘fold over IgG’.ChIP-qPCR primers used were as follows.
Negative Ctrl CHIP-R, TATCCCCACTGCCCAGAAGA.
Negative Ctrl CHIP-F, GATCAACGCAGGGGAGAGAG.
Procr CHIP-R, GTGAATGCACACACACACCC.
Procr CHIP-F, ATATCCGAGCTACACACGGC.
Ctfg CHIP-R, GAACTGAATGGAGTCCTACACA.
Ctfg CHIP-F, TGTGCCAGCTTTTTCAGACG.
Preparation of Procr promoter luciferase reporter and luciferase assay
Request a detailed protocolThe DNA sequence of Procr promoter containing TEAD4-binding sites (about 2 kb before the initiation codon) were amplified by PCR, separated by agarose gel, purified by Gel Extraction Kit, and then cloned into pGL3-promoter vector. Luciferase assays were performed in 293T cells with the pGL3-Procr promoter luciferase reporter described above 0.2 mg reporter plasmid were transfected together with CMV-Renilla (0.005 mg) to normalize for transfection efficiency. For luciferase assays in overexpression plasmid-transfected cells, cells were transfected with the indicated plasmids and reporter plasmid together, and then the luciferase activity was measured 36 hr later using Dual-Luciferase Reporter Assay System Technical Manual kit following the manufacturer’s protocol.
Cell culture, viral production, and infection
Request a detailed protocolHEK 293T was obtained from American Type Culture Collection (ATCC) and cultured in DMEM supplemented with 10% FBS plus 1% penicillin and streptomycin antibiotics at 37°C in 5% CO2 (vol/vol). For cells cultured on different modulus of elasticity, hydrogel substrates with tunable mechanical properties were prepared following the previous protocol (Tse and Engler, 2010), and the glass was as solid control. HEK 293T cells were used to produce lentivirus. When cells were up to 80%–90%, indicated constructs and packaging plasmids transfection was performed in Opti-MEM, and the media were replaced 12 hr later. Viral supernatants were collected 48–72 hr after medium change and filtered through a 0.45-μm filter, followed by concentration. For primary OSE cells infection, concentrated virus was diluted in the culture medium along with 1:100 polybrene.
RNA isolation and quantitative real-time PCR
Request a detailed protocolTotal RNA was isolated from fresh OSE cells or cultured cells lysed with Trizol according to the manufacturer’s instructions. The cDNA was generated from equal amounts of RNA using the SuperScriptIII kit. qPCR was performed on a StepOne Plus (Applied Biosystems) with Power SYBR Green PCR Master Mix. RNA level was normalized to Gapdh. The cycling condition was as 10 min at 95°C for initial denaturing, 40 cycles of 15 s at 95°C for denaturing, 1 min at 60°C for annealing and extension, following by melt curve test.
Quantification and statistical analysis
Request a detailed protocolFor quantification of nuclear Yap1+, Vgll4+, and EdU+ cells, 40 OSE cells at the both edges of ruptured sites (20 OSE cells at one side of rupture site) was identified as rupture regions, while other regions as nonrupture regions. At least 30 rupture regions and 30 nonrupture regions were counted. For quantification of the diameter of rupture, the longest diameter was counted, and at least 20 rupture sites were counted. For quantification of mG+ clone sizes, about 0.3 mm2 circle centered on ruptured sites was identified as rupture regions. At least 30 rupture regions were counted. For quantification of colonies size, diameters of the colonies were measured using Zeiss software.
Statistical analyses were calculated in GraphPad Prism (Student’s t-test or one-way analysis of variance). For all experiments with error bars, the standard error of measurement was calculated to indicate the variation within each experiment. All the p values were calculated using GraphPad PRISM six with the following significance: n.s. p > 0.05; *p < 0.05; **p < 0.01; ***p < 0.001. Statistical details for each experiment can be found in the figures and the legends.
Appendix 1
Data availability
All data generated or analysed during this study are included in the manuscript and supporting file; Source Data files have been provided for Figure 1-source data 1, Figure 2-source data 1, Figure 3-source data 1, Figure 4-source data 1, Figure 5-source data 1, Figure 1-figure supplement 2-source data 1, Figure 2-figure supplement 1-source data 1,Figure 2-figure supplement 1-source data 2,Figure 4-figure supplement 1-source data 1;Figure 4-figure supplement 2-source data 1;Figure 4-figure supplement 2-source data 2.
References
-
Molecular pathways regulating EGF-induced epithelio-mesenchymal transition in human ovarian surface epitheliumAmerican Journal of Physiology. Cell Physiology 290:C1532–C1542.https://doi.org/10.1152/ajpcell.00478.2005
-
Ovarian surface epithelium: biology, endocrinology, and pathologyEndocrine Reviews 22:255–288.https://doi.org/10.1210/edrv.22.2.0422
-
Yap controls stem/progenitor cell proliferation in the mouse postnatal epidermisThe Journal of Investigative Dermatology 133:1497–1505.https://doi.org/10.1038/jid.2012.430
-
WNT signaling in ovarian follicle biology and tumorigenesisTrends in Endocrinology and Metabolism 21:25–32.https://doi.org/10.1016/j.tem.2009.08.005
-
YAP1 increases organ size and expands undifferentiated progenitor cellsCurrent Biology 17:2054–2060.https://doi.org/10.1016/j.cub.2007.10.039
-
YAP regulates neural progenitor cell number via the TEA domain transcription factorGenes & Development 22:3320–3334.https://doi.org/10.1101/gad.1726608
-
The angiogenic factor Cyr61 activates a genetic program for wound healing in human skin fibroblastsThe Journal of Biological Chemistry 276:47329–47337.https://doi.org/10.1074/jbc.M107666200
-
Vgl-4, a novel member of the vestigial-like family of transcription cofactors, regulates alpha1-adrenergic activation of gene expression in cardiac myocytesThe Journal of Biological Chemistry 279:30800–30806.https://doi.org/10.1074/jbc.M400154200
-
Dual function of VGLL4 in muscle regenerationThe EMBO Journal 38:e101051.https://doi.org/10.15252/embj.2018101051
-
Activation mechanism of anticoagulant protein C in large blood vessels involving the endothelial cell protein C receptorThe Journal of Experimental Medicine 187:1029–1035.https://doi.org/10.1084/jem.187.7.1029
-
Cyclic changes of the ovarian surface epithelium in the ratReproduction (Cambridge, England) 129:311–321.https://doi.org/10.1530/rep.1.00401
-
Transduction of mechanical and cytoskeletal cues by YAP and TAZNature Reviews. Molecular Cell Biology 13:591–600.https://doi.org/10.1038/nrm3416
-
The two faces of Hippo: targeting the Hippo pathway for regenerative medicine and cancer treatmentNature Reviews. Drug Discovery 13:63–79.https://doi.org/10.1038/nrd4161
-
Regulation of the Hippo Pathway Transcription Factor TEADTrends in Biochemical Sciences 42:862–872.https://doi.org/10.1016/j.tibs.2017.09.003
-
The Hippo pathway in cellular reprogramming and regeneration of different organsCurrent Opinion in Cell Biology 43:62–68.https://doi.org/10.1016/j.ceb.2016.08.004
-
Hippo-YAP/TAZ signalling in organ regeneration and regenerative medicineNature Reviews. Molecular Cell Biology 20:211–226.https://doi.org/10.1038/s41580-018-0086-y
-
Ovary and fimbrial stem cells: biology, niche and cancer originsNature Reviews. Molecular Cell Biology 16:625–638.https://doi.org/10.1038/nrm4056
-
The biology of YAP/TAZ: hippo signaling and beyondPhysiological Reviews 94:1287–1312.https://doi.org/10.1152/physrev.00005.2014
-
The Hippo signaling pathway and stem cell biologyTrends in Cell Biology 22:339–346.https://doi.org/10.1016/j.tcb.2012.04.006
-
YAP/TAZ upstream signals and downstream responsesNature Cell Biology 20:888–899.https://doi.org/10.1038/s41556-018-0142-z
-
Preparation of hydrogel substrates with tunable mechanical propertiesCurrent Protocols in Cell Biology Chapter 10:Unit.https://doi.org/10.1002/0471143030.cb1016s47
-
Protein C receptor stimulates multiple signaling pathways in breast cancer cellsThe Journal of Biological Chemistry 293:1413–1424.https://doi.org/10.1074/jbc.M117.814046
-
Hippo/Yap Signaling in Cardiac Development and RegenerationCurrent Treatment Options in Cardiovascular Medicine 18:38.https://doi.org/10.1007/s11936-016-0461-y
-
A growing role for the Hippo signaling pathway in the heartJournal of Molecular Medicine (Berlin, Germany) 95:465–472.https://doi.org/10.1007/s00109-017-1525-5
-
TEAD mediates YAP-dependent gene induction and growth controlGenes & Development 22:1962–1971.https://doi.org/10.1101/gad.1664408
-
Claudin-18-mediated YAP activity regulates lung stem and progenitor cell homeostasis and tumorigenesisThe Journal of Clinical Investigation 128:970–984.https://doi.org/10.1172/JCI90429
Article and author information
Author details
Funding
National Key Research and Development Program of China (2019YFA0802000)
- Yi Arial Zeng
National Key Research and Development Program of China (2020YFA0509002)
- Yi Arial Zeng
National Natural Science Foundation of China (31625020)
- Yi Arial Zeng
Chinese Academy of Sciences (XDA16020200)
- Yi Arial Zeng
National Natural Science Foundation of China (32030025)
- Lei Zhang
National Natural Science Foundation of China (31625017)
- Lei Zhang
Fundamental Research Funds for the Central Universities (2020XZZX002-22)
- Junfen Fu
National Key Research and Development Program of China (2021YFC2701901)
- Junfen Fu
China Postdoctoral Science Foundation (2020TQ0260)
- Jingqiang Wang
Zhejiang Provincial Preferential Postdoctoral Foundation (ZJ2020150)
- Jingqiang Wang
Postdoctoral Exchange Fellowship Program 2021 (No. 160 Document of OCPC 2021)
- Jingqiang Wang
The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.
Acknowledgements
Thank to Dr. Chi Chung Hui of University of Toronto for helpful comments on the manuscript. This research was supported by grants from National Key Research and Development Program of China (2019YFA0802000 and 2020YFA0509002 to YAZ); National Natural Science Foundation of China (31625020 to YAZ; 32030025 and 31625017 to LZ); the Chinese Academy of Sciences (XDA16020200 to YAZ); Fundamental Research Funds for the Central Universities (2020XZZX002-22 to JF); the National Key Research and Development Programme of China (2021YFC2701901 to JF); China Postdoctoral Science Foundation (2020TQ0260 to JW); Zhejiang Provincial Preferential Postdoctoral Foundation (ZJ2020150 to JW); Postdoctoral Exchange Fellowship Program 2021 (No.160 Document of OCPC, 2021) to JW; Shanghai Leading Talents Program to LZ.
Ethics
All mice were housed in the SIBCB animal facility under IVC standard with a 12-hr light/dark cycle at room temperature. Experimental procedures were approved by the Animal Care and Use Committee of Shanghai Institute of Biochemistry and Cell Biology, Chinese Academy of Sciences, with a project license number of IBCB0065.
Copyright
© 2022, Wang et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 802
- views
-
- 194
- downloads
-
- 0
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
Premature infants with bronchopulmonary dysplasia (BPD) have impaired alveolar gas exchange due to alveolar simplification and dysmorphic pulmonary vasculature. Advances in clinical care have improved survival for infants with BPD, but the overall incidence of BPD remains unchanged because we lack specific therapies to prevent this disease. Recent work has suggested a role for increased transforming growth factor-beta (TGFβ) signaling and myofibroblast populations in BPD pathogenesis, but the functional significance of each remains unclear. Here, we utilize multiple murine models of alveolar simplification and comparative single-cell RNA sequencing to identify shared mechanisms that could contribute to BPD pathogenesis. Single-cell RNA sequencing reveals a profound loss of myofibroblasts in two models of BPD and identifies gene expression signatures of increased TGFβ signaling, cell cycle arrest, and impaired proliferation in myofibroblasts. Using pharmacologic and genetic approaches, we find no evidence that increased TGFβ signaling in the lung mesenchyme contributes to alveolar simplification. In contrast, this is likely a failed compensatory response, since none of our approaches to inhibit TGFβ signaling protect mice from alveolar simplification due to hyperoxia while several make simplification worse. In contrast, we find that impaired myofibroblast proliferation is a central feature in several murine models of BPD, and we show that inhibiting myofibroblast proliferation is sufficient to cause pathologic alveolar simplification. Our results underscore the importance of impaired myofibroblast proliferation as a central feature of alveolar simplification and suggest that efforts to reverse this process could have therapeutic value in BPD.
-
- Developmental Biology
We have generated single cell transcriptomic atlases of vomeronasal organs (VNO) from juvenile and adult mice. Combined with spatial molecular imaging, we uncover a distinct, previously unidentified class of cells that express the vomeronasal receptors (VRs) and a population of canonical olfactory sensory neurons in the VNO. High-resolution trajectory and cluster analyses reveal the lineage relationship, spatial distribution of cell types, and a putative cascade of molecular events that specify the V1r, V2r, and OR lineages from a common stem cell population. The expression of vomeronasal and olfactory receptors follow power law distributions, but there is high variability in average expression levels between individual receptor and cell types. Substantial co-expression is found between receptors across clades, from different classes, and between olfactory and VRs, with nearly half from pairs located on the same chromosome. Interestingly, the expression of V2r, but not V1r, genes is associated with various transcription factors, suggesting distinct mechanisms of receptor choice associated with the two cell types. We identify association between transcription factors, surface axon guidance molecules, and individual VRs, thereby uncovering a molecular code that guides the specification of the vomeronasal circuitry. Our study provides a wealth of data on the development and organization of the accessory olfactory system at both cellular and molecular levels to enable a deeper understanding of vomeronasal system function.