Note: This rebuttal was posted by the corresponding author to Review Commons. Content has not been altered except for formatting.
Learn more at Review Commons
Reply to the reviewers
RC-2022-01661
Response to reviewers:
Review Commons questions and Reviewers’ comments verbatim in plain text.
Authors’ responses in bold text. Line numbers refers to numbers in the marked-up manuscript. In text citations in this document – see bibliography at bottom of this document.
Reviewer #1 (Evidence, reproducibility and clarity (Required)):
Summary:
Cells within multicellular organisms are mutually dependent on each other - cells of one type or in one location provide signals that can regulate the health and differentiation of the target cells that receive those signals. Such signalling can operate bi-directionally, emphasizing the co-dependence of cells upon each other. The ovarian follicle provides an excellent model system to study intercellular signaling and its consequences, in this case between the oocyte and the somatic granulosa cells that surround it. Oocytes secrete members of the TGFbeta growth factor family that are required for normal differentiation of the granulosa cells, which in turn is necessary for normal development of the oocyte. Here the autohors show that adding TGFB-type growth factors (cumulin or BMP15) to the cuture medium during in vitro maturation increases the fraction of oocytes that can reach the blastocyst stage (improved developmental competence) and alters the pattern of protein landscape in both the (cumulus) granulosa cells and the oocyte. Changes in the mitochondria and parameters relevant to energy metabolism are also altered. They conclude that these changes underpin the acquisition of developmental competence by the oocytes.
Major issues:
The authors are world leaders in this field and therefore exceptionally well-qualified to carry out the proposed work. There are a number of issues, however, that limit the confidence with which conclusions may be drawn.
First, the experimental strategy makes drawing inferences about the role of cumulin and BMP15 challenging. Maturing oocytes express GDF9 and BMP15 (the components of cumulin). Thus, the experiments are not comparing presence vs absence of cumulin and BMP15, but rather comparing oocytes and cumulus cells exposed to supra-physiological levels of these factors to controls that are exposed to physiological levels. In other words, the experimental setup detects changes that occur in response to higher than normal levels of the factors. Ideally, one would have complementary experiments where GDF9 and BMP15 were deleted from the system, to illustrate the effects of their absence. This would be a massive additional undertaking, however. Yet, without such experiments, relying on the results of the overexpression approach to understand the functions of cumulin and BMP15 at physiological levels is risky.
RESPONSE #1 We appreciate these insightful perspectives. We apologise for not making it clear that the model used is not in fact an overexpression model. This is because, by removing the cumulus-oocyte complex from the follicle and studying it in vitro (oocyte IVM), secretion of these growth factors by the oocyte is notably compromised, so the controls are not exposed to normal physiological levels as suggested by the reviewer. This loss of normal secretion ex vivo is evidenced by: 1) in Mester B. et al _[1]_; Figure 2, we showed the mouse oocytes matured in vitro (i.e. as per the current study) are essentially devoid of the mature domain BMP15 protein, which will therefore be likewise for cumulin as cumulin contains one subunit of BMP15, and 2) mammalian cumulus-oocyte complexes explanted and cultured in vitro by IVM benefit (in terms of developmental competence) from the addition of exogenous oocyte-secreted factors such as BMP15, GDF9 and cumulin, demonstrating that they are rate-limiting under IVM conditions. We were the first to demonstrate this in 2006 _[2]_ which has been subsequently verified in many papers, including in the current paper for cumulin. The exact extent to which the controls are deficient in BMP15 and cumulin is unclear, as there are not yet reliable mouse ELISAs for these, but the model is an add-back model rather than an overexpression model. We have now added text at lines 150-152 and in the Fig 1 legend, to make this point clearer.
Re using complimentary deletion, knock-out or antagonist-type experiments: we agree this would be ideal. However, this is likely impossible as cumulin is a non-covalent heterodimer of BMP15 and GDF9 (as first named and characterised by us: Mottershead DG et al ____[3]____). Hence, to knockout cumulin one needs to knockout either or both of BMP15 and GDF9, making it impossible to discriminate the actions of the heterodimer from the homodimers. In support of this, reviewer #3 made exactly this point, and stated “Such functional analysis cannot be done using gene knockout mouse lines…… only functional work as the one presented in this manuscript can find the mechanisms of action of these hormones”. This issue is further complicated by the fact that BMP15 and GDF9 are thought to exist as homodimers, as well as monomers, including in equilibrium in heterodimeric form as cumulin (also noted by Reviewer #3). Furthermore, there is no cumulin-specific antagonist, e.g. a cumulin-specific neutralizing antibody. Small molecule signaling inhibitors (e.g. Smad2/3 or Smad1/5/8 antagonists) certainly block cumulin actions, but therefore simultaneously also block GDF9 or BMP15 actions. Collectively, these unique (with the TGFβ superfamily) structural peculiarities of cumulin make it complex to interrogate its mechanisms of action, to the extent that others have largely focused on BMP15 or GDF9 homodimer actions only, when in reality, cumulin is likely the key natural protagonist responsible for oocyte paracrine signalling. We have added a paragraph to this effect to the discussion, at lines 417-423, including acknowledging the experimental limitations of the study dictated by having to deal with a noncovalent heterodimer.
Second, the granulosa cells and oocytes interact throughout the prolonged period of growth, and this is the time when the beneficial effects of the granulosa cells on the oocyte have been most clearly documented. Yet the experiments focus on the much shorter period of meiotic maturation. This is when oocyte-granulosa cell interaction is being down-regulated, even if not entirely disrupted.
RESPONSE #2: Indeed, oocyte-granulosa interaction is absolutely essential during oocyte growth, development and meiotic maturation, for healthy oocyte function, including the orchestrated down-regulation of oocyte-granulosa interactions during the latter phase. As pioneered by John Eppig and others, including ourselves ____[4]____ (ref has 673 citations), the master conductor of this dynamic oocyte-granulosa interaction during oocyte meiotic maturation are the oocyte-secreted factors. Hence, these factors are critical at this stage, and we maintain that this is a very important phase of oocyte development to study.
Third, the data reported illustrate associations or correlations, but no experiments test the function of the changes in the proteome or of the changes in the morphology of the mitochondria or ER. Which if any of these is linked to the improved development of the oocytes after fertilization is unknown. Moreover, no experiments address how the growth factors cause the observed changes, which occur over a period of a few hours.
RESPONSE #3 This is true. The study is already very large and has many functional experiments (e.g. oocyte respiration, oocyte MS, etc), that follow-up the findings from the proteomic analysis. Hence, the study has taken a global cellular metabolism approach, e.g. we show that cumulin downregulates oxidative phosphorylation globally, c.f. pathways within OXPHOS. We found an abundance of individual proteins altered in this period (see figure 4) and to follow up on the actions and consequences of individual proteins would: 1) at best show small incremental effects, as metabolism of such a cellular syncytium is vastly complex and inter-connected, 2) further increase the size of what is already a large study, and 3) detract from the more important wholistic effects on cumulus-oocyte complex metabolism, which must act as whole, interacting entity, to support the complexities of supporting early life post-fertilization.
__Taken together, these issues unfortunately limit the potential impact of the work. But the amount of work required to address them would be substantial and not really feasible for this manuscript. The best route may be to present the work as an overexpression study that has identified associations, with a discussion that acknowledges the limitations of this approach.
__RESPONSE #4 This is not an over-expression study – see RESPONSE #1 above. We have added text in the discussion at lines 417-423, that acknowledges the limitations of the study by the impossibility to conduct a killer knockout experiment of cumulin.
Minor issues:
The text of the manuscript should be revised in a number of places.
32: We characterized the molecular mechanisms by which two model OSFs, cumulin and BMP15, regulate oocyte maturation and cumulus-oocyte cooperativity.
--Mechanistic studies were not performed.
RESPONSE #5 The scope of this work was to; (a) identify global changes to protein expression, and (b) to use this data to implement follow-up experiments on some of the lead indicators, such as metabolism (respiration, small molecule metabolic markers) and cellular morphology. This work provides the groundwork, insight and rationale for future additional studies of specific mechanisms of COC interactions. As discussed at RESPONSE# 1, these studies are as close as anyone can probably get currently to mechanistic studies of a NOVEL noncovalent heterodimer, when the noncovalent homodimers are also in play, as also noted by reviewer #3 who specifically references mechanisms: “…… only functional work as the one presented in this manuscript can find the mechanisms of action of these hormones”.
In some instances, in the interests of brevity, we made remarks based on our data, but without specifying details in the text. To redress this, we have now added specific details which illustrate and justify our statements based on the data collected (see RESPONSES #6, #7, #9 below). For greater clarity, we have also restructured our supplementary data set to cover the analysis progression from full raw proteomic data to differentially expressed proteins, to use of differentially expressed proteins in network analysis. The supplementary data set now includes the full proteomics lists for both cells and treatments (Supplementary Tables S1, S2, S3, S4), protein sequences confidently identified by both proteomic software platforms (Supplementary Tables S5, S6), differentially expressed proteomics lists for both cells and treatments (Supplementary Tables S7, S8, S9, S10), differentially expressed protein list used for the network analysis (supplementary Table S11). The Table S11 lists are intended to facilitate use by readers to perform their own analyses, if they so wish, since they can simply copy and paste the list to the on-line STRING platform. Finally, the reanalysed network analysis output, based on the differentially expressed proteins shown in supplementary Table S11, are shown in supplementary Tables S12 and S13.
__40: Collectively, these data demonstrate that OSFs remodel cumulus cell metabolism during oocyte maturation in preparation for ensuing fertilization and embryonic development.
--No mechanistic studies demonstrate this.
__RESPONSE #6 There is no mention of mechanism in this sentence at line 40 and we have provided exhaustive evidence that cumulus cell metabolism is remodelled as stated (Figures 4B and 4C). For example, of the 59 upregulated proteins in the cumulus cells of cumulin treated COC (Figure 4C and supplementary Table S11), 38 (i.e. 64%) are involved in primary metabolic processes (supplementary Table S12), including amino acid metabolism (GOT2, SHMT1, CTH, MAT2B), lipid and steroid metabolism (CERS5, DHCR7, HSD17b4), aldehydes metabolism (RDH11), nucleotides biosynthesis (RRM1, GMPR2), glycans biosynthesis and protein glycosylation (UGDH, GFPT2, GALNT2), respiratory chain (mt-ND1). The cellular macromolecule metabolic process is also a significantly enriched network, involving 26 out of the 59 upregulated proteins (i.e. 44%, Figure 4C and supplementary Table S11) and includes processes such as protein complex assembly (TM9sF4, DHX30, AP2M1), RNA metabolism and mRNA processing (DDX17, DDX5, DDX39bPRPF19, PRPF6, HNRNPF, CPSF6). To help clarify the specificity of our findings, we have added this text to the revised manuscript (lines 465-474).
__46: Oocyte-secreted factors downregulate protein catabolic processes, and upregulate DNA binding, translation, and ribosome assembly in oocytes.
--No direct evidence is provided.
__RESPONSE #7 The proteomic data provides direct evidence that these processes are involved. Sentence modified at lines 47-48 to be more specific re processes. Additional text has been included (revised manuscript lines 434-443) to provide specific details of the differentially expressed proteins involved in each of these processes.
48: Oocyte-secreted factors alter mitochondrial number...
--Need to establish that the MitoTracker is a suitable tool to measure the number of mitochondria.
RESPONSE #8____ We recognise that total mitochondrial uptake of the MitoTracker Orange dye could be a reflection of either mitochondrial function (polarity) and/or mitochondrial number, given the manufacturer’s (Thermo Fischer) statement that “MitoTracker™ Orange CMTMRos is an orange-fluorescent dye that stains mitochondria in live cells and its accumulation is dependent upon membrane potential”, as we specified in several places in the original manuscript (Lines 354-355, 366-367 and 235 of the marked up manuscript version) . However, we agree that in several places in the manuscript we also indicated that MitoTracker was being used as a measure of mitochondrial number. To avoid this ambiguity, we have made some clarifications in the text (revised manuscript lines 235, 351-352, 377, 481-482, and in Figure 5B legend). Given the extensive and diverse metabolic changes indicated by the proteomic data, our aim was to explore the potential role of mitochondria in response to cumulin and BMP15 treatment of COCs, which we did by use of EM morphology studies (figure 5A), mitochondrial respiration (figures 6B and 6C), quantification of energy metabolites, such as ATP, NAD and related compounds, by mass spectrometry (figure 6D), metabolites identified in multispectral unmixing studies (figure 7) and mitochondrial function using MitoTracker (figure 5B). Collectively this data suggested a modest downturn of energy metabolism, particularly in cumulin treated COCs. This downturn did not cause a change in net energy charge in COCs (figure 6D) despite a reduction in redox ratio in both cells (figure 7A and 7B) and respiration in COCs (Figure 6B and 6C), and could reflect adaptive changes in response to cumulin and BMP15, reflecting metabolic plasticity/Warburg effect, as explained in the discussion (revised manuscript lines 453-551).
79: ...for maintaining genomic stability and integrity of the oocyte...
83: ...minimizing secondary production of potentially DNA damaging free radicals.
--Please provide supporting references from the literature.
RESPONSE #9 References have been added (lines 82 and 85 of the revised manuscript)
373: This study provides a detailed exploration of the mechanisms by which oocyte-secreted factors...
--No mechanistic studies were performed.
RESPONSE #10 We respectfully disagree. One of many mechanisms we have studied here is OXPHOS. We have shown this is how OSFs change metabolism – that is a mechanism. As discussed at RESPONSE #1, these studies are as close as anyone can probably get currently to mechanistic studies of a noncovalent heterodimer, when the noncovalent homodimers are also in play, as also noted by reviewer #3 who specifically references mechanisms: “…… only functional work as the one presented in this manuscript can find the mechanisms of action of these hormones”. Please also refer to the comments in RESPONSE #5.
383: Collectively, these data demonstrate that oocyte paracrine signaling remodels COC metabolism in preparation for ensuing fertilization and embryonic development.
--Studies do not show that the differences observed between control and treatment groups are related to fertilizability or embryonic development.
RESPONSE #11 The data in Fig 2C, 2D show exactly that; that the difference between control and treatment (cumulin) is an increase in embryonic development. It does not show fertilizability, so we removed that at lines 41 and 415.
396: suggesting that cumulin affects meiosis in the oocyte and may increase meiotic fidelity...
--This statement is highly speculative.
RESPONSE #12 We accept this critique - reference to meiosis and meiotic fidelity removed, line 435 (revised manuscript).
409: ...lacks the machinery for amino acid uptake...
--Is the oocyte unable to take up any amino acids or only certain amino acids?
RESPONSE #13 Thank you for noting this as this sentence is too absolute. Oocytes have a very poor capacity to take up most or even all AAs, which are instead supplied to the oocyte via cumulus cells. Sentence modified at lines 455-456 to be less absolute.
In general, the manuscript is written clearly. However, in several places, technical terms or jargon will make tough going for readers who are not already familiar with the techniques being used. These should be explained using language that will be understood by journal readers who are unfamiliar with the details of the techniques. Examples include:
51: define metabolic workload using scientific terms.
RESPONSE #14____ “metabolic workload” rephrased to “metabolic processes”. Lines 52-53.
67: metabolically 'inept' requires more precision.
RESPONSE #15 “metabolically inept” rephrased to “metabolically dependent on surrounding granulosa cells” ____[5]____. Line 69
262: explain 'multispectral analysis'
RESPONSE #16 A citation has been added, which explains the technique (ref ____[6]____ at the end of this response letter, which is the same paper as citation [34] in the revised manuscript; lines 111 and 217; revised manuscript). A detailed explanation of this technique has also been added in the supplementary information, under the section “Multispectral microscopy”.
268: how is 'limited' overlap defined.
RESPONSE #17 The phrase “distinct profiles, with limited overlap between…” has been rephrased to “distinct profiles, between…” (line 279 of the revised manuscript), as the main point is that the patterns/profiles across treatments are different, and we did not quantify the extent of overlap.
318: define higher workload
RESPONSE #18 the phrase “…implying a higher workload for both organelles” has been replaced with a more specific explanation; “We suggest that such changes in morphology may be related to the remarkable increase in a diversity of metabolic processes which we observed (Figure 4C and supplementary Table S12), since ER morphology and architecture is known to be highly dynamic in response to environmental and developmental factors which affect cells” ____[7]____ (Lines 342-345).
324: provide documentation or citations to support the assertion that the intensity of MitoTracker staining is an accurate proxy for the number of mitochondria.
RESPONSE #19____ Please refer to explanation under RESPONSE #8
358: Multispectral discrimination modelling utilised cellular image features from the autofluorescent profiles of oocytes and cumulus cells.
--Please clarify this merthodology and provide support for its utility.
RESPONSE #20____ The supplementary information section (Multispectral microscopy, lines 239-258) has been expanded and clarifications provided as to the wavelengths of the channels, the features used and the unsupervised nature of algorithms.
360: intersection of union of 5-22%
RESPONSE #21____ This is a measure of the extent of overlap of data distribution for each class (treatment), i.e. of how different they are. The ellipse (Fig 3D) represents one standard deviation around the central mean value for that data set. The overlap of these ellipses is quantified by their intersection over union (IoU) value, which is the ratio of the area of the two-ellipse intersections, divided by the area of their union (the shape created by their overlap being treated as creating one continuous object). IoU values range from 0 to 100% for fully separated and fully overlapping, respectively. Hence, a 5% IoU represents a low level of overlap of data distribution between treatments. Brief explanatory text has now been added at line 387-388.
Comments on Figures.
Fig. 3A, B. The total number of proteins and the number of differentially expressed proteins among the treatment groups don't match between A and B. For example, A (Mascot-Sheffield) indicates that 17 proteins were differentially expressed between untreated and cumulin-treated oocytes. B shows (138 + 74) expressed un the untreated but not cumulin-treated and (156 + 87) expressed in the cumulin-treated but not untreated. Please account for this difference.
RESPONSE #22 The panels in Fig 3A and Fig 3B each contain different representations of the information contained within the proteomics dataset, and explain different aspects of the data. The Venn diagram panels in Figure 3B display the level of overlap of specific proteins identified in each cell, treatment and software subgroup. The degree of overlap in each cluster is high (i.e., 76 – 78% for Mascot/scaffold and 95 - 97 % for PD2.4) as would be expected within the same cell type and analysis approach, where the main variable is cell treatment. We agree that the total numbers in the Venn diagrams did not exactly match the total numbers in Figure 3A, which likely resulted from using slightly different parameters during data processing. We have now used exactly the same data set in panels A and B (the full PD2.4 and Mascot/scaffold datasets are shown in the supplementary proteomics summary Excel spreadsheet), so that total numbers are now identical, and will hopefully avoid any confusion in comparing across panels. However, the main conclusion to be drawn from Fig 3B remains unchanged, in that it shows that by far the majority of identified proteins overlap between treatments (control, BMP, cumulin), regardless of cell type or data analysis approach. However, it should be noted that Figure 3B has no information about protein fold change/differential expression, and only represents numbers of proteins confidently identified, and the level of overlap of identified proteins between treatments. Only panel 3A shows differential protein expression relative to the respective control groups.
Fig. 3D. What do the circles represent and how were their parameters (size, position) established?
RESPONSE #23 The separation of data distributions for each class is shown by an ellipse for each cluster, which encompasses one standard deviation around the central mean values. This text has now been added to the Fig 3 legend.
Reviewer #1 (Significance (Required)):
These studies identify changes in cumulus cells and oocytes that occur in response to addition of cumulin or BMP15 to the culture medium during in vitro maturation. While the data are new, the significance of the advance is limited by (i) the fact that the control group were exposed to physiological levels of GDF9 and BMP15, so this is essentially an over-epxression study and (ii) no mechanistic studies experimentally tested how the observed changes (eg, in quantity of a specific protein) affect the developmental potential of the oocytes or cumulus cells.
RESPONSE #24 We thank the reviewer for their perspectives however we respectfully disagree on all accounts. We have rebutted these 2 concerns: point (i) at RESPONSE #1, and point (ii) at RESPONSE #5 above.
Reviewer expertise: growth and meiotic maturation of the mammalian oocyte
Reviewer #2 (Evidence, reproducibility and clarity (Required)):
Summary:
The report by Richani et al, presents a research carried out in mice, in which they treated cumulus-oocyte complexes with either BMP15 and cumulin. Upon treatment they evaluated a series of biologically relevant parameters in oocytes and cumulus cells. Their findings indicate that the treatment with these molecules alter the molecular composition of oocytes and cumulus cells (proteome and metabolome), mitochondrial morphology in cumulus cells and overall oxygen consumption in COCs.
Major comments:
- Are the key conclusions convincing?
- Should the authors qualify some of their claims as preliminary or speculative, or remove them altogether?
* part of the discussion related to metabolic pathways being up regulated due to the treatments need to the revised because For instance, It is hard for me to grasp how a pathway with 2 proteins achieved FDR significance below 0.01, as I see in figure 4c
RESPONSE #25____ Network enrichment was performed using the open access software STRING ( ____https://string-db.org/ [8]____), and we have now provided additional information on how we utilised STRING in the supplementary information section, under “Gene Ontology Network Enrichment Analysis” (lines 176-217). STRING utilises information available in the Gene Ontology (GO) database ( ____http://geneontology.org/docs/ontology-documentation/____ ) to determine; (a) how many of the differentially expressed proteins identified in the proteomics experimental data fall into specific networks, (b) how much enrichment this represents relative to a random network of the same size, and (c) whether the enrichment is statistically significant based on the FDR statistic. The size of each GO network within the background set (whole genome or other) will therefore be a major determinant of whether the number of proteins identified in the proteomics experiment represents significant enrichment of a particular network. A few proteins identified within a small background network will represent greater enrichment (and lower FDR score) than the same number of proteins in a much larger network. In fact the “count in network” is often approximately the inverse of the enrichment strength (see supplementary Table S12, within the supplementary dataset Excel spreadsheet). Note that only significantly differentially expressed proteins were used for the network analysis presented in this paper, so even in the case where just 2 proteins are significantly enriched in a network (e.g., “Farnesyl diphosphate metabolic process” identified in the GO biological process section of BMP15 treated cumulus cells) they represent two upregulated proteins in a small network, so the functional/biological significance of this is likely quite high.
In revision of the manuscript we noticed that we had likely originally used the full lists of differentially expressed proteins for network analysis, rather than separating up and downregulated proteins as intended. Furthermore an updated version of STRING is now available, with improvements in the method of correction for multiple testing within the FDR output (STRING version 11.5, current since August 12, 2021). We have therefore revised the STRING network analyses, and have provided a list of the STRING input proteins (supplementary Table S11), STRINGv11.5 gene ontology (GO) functional enrichments for up and downregulated proteins in BMP and cumulin treated cumulus cells and oocytes respectively (supplementary Tables S12 and S13), and replaced the very large Figure 4C and D heatmaps (submitted version) with a summary (new Figure 4C; revised version). The updated heat maps can still be viewed in supplementary Tables S12 and S13 (the heatmaps now being the updated ones, deriving from our review response).
* In the discussion the authors use the term "oocyte secreted factors" a lot (one example lanes 490, 515, 516, 517), but they should specify BMP15 and cumulin, because these were their treatments.
*Including in the title, you did not evaluate all oocyte paracrine factors, just BMP15 and cumulin
RESPONSE #26 “Oocyte secreted factors (OSFs)” replaced with BMP15 and cumulin throughout the manuscript where we refer specifically to our treatments, results or discussion of results, except where we refer to “these OSFs” (eg line 34), and not where we refer to the principal of OSF signalling more generically. Re the latter, hence we wish to retain the title as is, as BMP15 and cumulin are prototypical oocyte secreted factors, as the title refers to the principal more generally.
- Would additional experiments be essential to support the claims of the paper? Request additional experiments only where necessary for the paper as it is, and do not ask authors to open new lines of experimentation.
NA
-
Are the suggested experiments realistic in terms of time and resources? It would help if you could add an estimated cost and time investment for substantial experiments.
NA
-
Are the data and the methods presented in such a way that they can be reproduced?
*no, in some instances, the methods are not described, see my comment below about enrichment analysis.
RESPONSE #27 Addressed next below
-
Are the experiments adequately replicated and statistical analysis adequate?
*I was not able to access enrichment analysis.
RESPONSE #28____ The method of Network Enrichment is now described in more detail in the supplementary methods section. See previous explanation under RESPONSE #25 above.
*lines 241-242: "MitoTracker staining and data from metabolite analysis by mass spectrometry were analysed by one-way ANOVA with Tukey's (parametric data) or Kruskal-Wallis (non-
parametric data) post-hoc tests. "
Specify which test was used for which data
RESPONSE #29 Post-hoc test for MitroTracker data was Tukey’s, as already stated in Figure 5 legend. Post-hoc test for metabolite analyses was Kruskal-Wallis – text now added to Figure 6 legend.
Minor comments:
- Specific experimental issues that are easily addressable.
NA
-
Are prior studies referenced appropriately?
Yes
-
Are the text and figures clear and accurate?
*lines 178-180: "expressed proteins list was further analyzed using STRING software to explore clustering and enrichment of specific molecular functions, and biological pathways. Detailed methodology and rationale for this approach is provided in the supplementary methods."
I did not read text in the supplementary materials indicating how enrichment analysis was carried out.
RESPONSE #30____ Our apologies for this oversight. We have now provided additional information on how we utilised STRING in the supplementary information, in a new section titled “Gene Ontology Network Enrichment Analysis” (lines 176-217).
* What was the concentration of treatment for the samples used for proteome and mascot/scaffold experiments?
RESPONSE #31____ The two bioinformatic analyses were conducted on common biological samples, so naturally the treatment concentrations were also the same. Text modified at line 175 to make this clearer.
* lanes 263 and 264: "Cell types and treatment conditions can be clearly distinguished based on these orthogonal global approaches."
I did not see what is the basis for this statement
RESPONSE #32____ The sentences immediately following this (i.e. lines 274-281) elaborated the basis for this statement, particularly where it is explicitly stated “____Proteomic heat maps (Fig. 3C) and multispectral analysis plots (Fig. 3D) both show distinct profiles, between controls, BMP15 and cumulin treated COCs, in both cell types.____”, at lines 277-281.
The data for the two global approaches are shown in Figure 3C (heat maps generated by PD2.4 comparing differences in protein abundance across treatments, shown separately for cumulus cells and oocytes), and Figure 3D (linear discriminant analysis comparing differences in multispectral imaging data across treatments, shown separately for cumulus cells and oocytes). Both of these global analyses show clear differences in distribution pattern between controls (untreated) and treated samples (BMP15 and cumulin), in both oocytes and cumulus cells. The approaches are (a) global, since each relates to analysis of the complete cell extracts (as opposed to targeting a specific component/analyte), and (b) orthogonal because different and unrelated measurement techniques are used i.e., proteomics (mass spectrometry) and multispectral imaging (spectroscopy).____
*I did not understand the discrepancy between the numbers observed in Figure 3A and Figure 3B.
RESPONSE #33____ Refer to RESPONSE #22 above. We have checked the data, and revised the Venn diagrams (Figure 3B) with data analysed using identical parameters, for both Figures 3A and 3B, to avoid confusion over protein numbers. We also noticed and corrected a discrepancy with regard to the number of differentially expressed oocyte proteins under the merged data column of Figure 3A.____
*I could not make sense of the shades of green or red that were used in 4C and 4D. Is the reader only supposed to make those comparisons within column?
RESPONSE #34 Note: Figures 4C and 4D are now Supplementary Tables S12 and S13. The red shades represent network enrichment analysis of upregulated proteins, while the green shades represent network enrichment analysis of downregulated proteins. The colour gradients in each case follow the numerical values for “count in network”, enrichment strength, and lower FDR, with greater colour intensity for higher numbers (and lower FDR). However, we agree that the original four panels (A, B, C and D) comprising figure 4, made for a very large and potentially overwhelming figure. To simplify the data presentation we have reprocessed the data in STRING (see details under RESPONSE #25 above) and have moved the now considerably shorter network lists (originally displayed as Figures 4C and 4D) to supplementary Tables S12 and S13, and the new Figure 4C provides a network enrichment summary instead. This is likely easier to comprehend, with the marked contrast in networks identified between oocytes and cumulus cells easier to see. The numbers of up and downregulated proteins on which the network analysis is based are also shown in Figure 4C, while the specific proteins used and networks identified are shown in supplementary tables S11, S12 and S13 (original colour coding retained, and also explained within each table).
- Do you have suggestions that would help the authors improve the presentation of their data and conclusions?
*Figure 4 is really hard to process. At least in my pdf it spanned 4 pages.
RESPONSE #35____ Indeed Figure 4 was large and has now been shortened. We made considerable effort to attempt to present in Fig 4 the vast amount of proteomic data in a summarized, hopefully comprehensible fashion. We have now moved Figs 4C and 4D to the supplementary, and replaced it with the simplified new Fig 4C (tabular format). Pease also see comments under RESPONSES#25 and #34 re this.
*I did not understand why put networks that are not significant as up-regulated or down-regulated. Besides, as mentioned above, I do not know how significance was assessed..
RESPONSE #36 Network analysis was performed using only those proteins which were significantly differentially expressed and had a consistent direction of fold change in both mascot/scaffold spectral counting and PD2.4 peak intensity proteomics quantitative approaches. Proteins with no significant expression change (i.e., the majority of proteins, which represented proteins with
__Reviewer #2 (Significance (Required)):
-
Describe the nature and significance of the advance (e.g. conceptual, technical, clinical) for the field. - Place the work in the context of the existing literature (provide references, where appropriate).
*This paper is significant because it provided a variety of measurements following the treatment of cumulus cells with BMP15 and cumulin. The authors show that these two oocyte factors can impact the molecular structure, physiology and structure of organelles in cumulus cells. The work is well contextualized with the current literature.
RESPONSE #37 We thank the reviewer for these positive remarks.
-
State what audience might be interested in and influenced by the reported findings.
*Researchers in the field of developmental biology would be most interested in this report.
-
Define your field of expertise with a few keywords to help the authors contextualize your point of view. Indicate if there are any parts of the paper that you do not have sufficient expertise to evaluate.
* I do not have expertise in hyperspectral analysis. I have been working with cumulus-oocyte complexes for over a decade, mixing technologies in cell biology, microscopy, high-throughput genome, and proteome analysis. We do all our bioinformatics work in-house.
Reviewer #3 (Evidence, reproducibility and clarity (Required)):
The work is interesting. Cumulin is a heterodimer hormone formed of GDF9 and BMP15. It is the main oocyte secreted factor. Being an heterodimer, gene knockout provides very little information about its mechanism of action. The team has a unique form of cumulin that is stable. This is why I think this work is important. However, I found two technical issues: one regarding mitochondrial count using MitoTracker and the other about comparing gene lists between the two cell types when protein input submitted to mass spectrometry differ between the two cell types. It is expected to find more with more input material. The text would need to be adjusted accordingly. Also, there is a lot of free statements and a lack of precision that is annoying. In my opinion, there are many overstatements that are not supported by the data because the work was not designed to test what is stated. The Discussion is very circular as the same statements come back on the next pages.
RESPONSE #38 See specific responses below
Detailed review:
The manuscript entitled "Oocyte and cumulus cell cooperativity and metabolic plasticity under the direction of oocyte paracrine factors" reports an in depth analysis of the exposure of cumulus oocyte complexes to either BMP15 or cumulin, the GDF9-BMP15 heterodimer. Impact assessment was done by determining developmental competence of the exposed oocytes, comparative profiling of the proteomes and spectral emissions as well as testing a potential impact at the ultrastructure level by electron microscopy imagery. Mitochondrial respiration as well as abundance of related metabolites was contrasted between the two treatments.
Overall, the work is interesting. It is very difficult to study hormonal heterodimers because they originate from two different genes and they can naturally be found in a monomeric as well as a dimeric state. Such functional analysis cannot be done using gene knockout mouse lines. Genetic disruption provided the background that GDF9 and BMP15 are key oocyte secreted factors however only functional work as the one presented in this manuscript can find the mechanisms of action of these hormones.
RESPONSE #39 We thank the reviewer for these positive comments, especially in relation to the difficulty of getting to the mechanism of actions of a non-covalent heterodimer, and hence the importance of functional experiments in providing mechanistic insights.
Comments:
I really appreciated the reference to auto-symbiosis. We often see the reference to a cellular syncytium but this one is interesting.
RESPONSE #40 Thank you.
Although I appreciated the work, two important technical issues (between cell types comparisons and mitochondrial count) have been raised and there is a bit of unnecessary overselling throughout the manuscript. Sticking to the results would keep the value of the work high and wouldn't give that impression of overstatement.
RESPONSE #41 Technical issues – see responses below, as well as responses to other reviewers. We have provided additional methodological information for greater clarity, and added specific observations from our data, to support all statements, to avoid the impression of unsubstantiated overstatements.
Technical issues:
While the gene/protein enrichment analysis can be influenced by the input material submitted to mass spectrometry, the gene network analysis is influenced by the number of gene/proteins available for the enrichment analysis. It is thus difficult to compare both cell types.
RESPONSE #42 We agree that shorter protein lists might be expected to result in fewer networks. However, it is interesting to consider the possible reasons for the shorter list:
(1) In our case the amount of protein extracted from oocytes (2-3____m____g) was much less than from cumulus cells (15-17____m____g) as explained in the “Mass Spectrometry for proteomic analysis” section in the Supplementary Information. This is because COCs have many more cumulus cells than oocytes by number as well as total mass. Consequently it was possible to load a larger ____m____g amount of total peptides from cumulus cells onto the nanoLCMSMS system, but it should be noted that on-column loading is not only determined by the total amount of material injected, but also by the limits in capacity of the C18 peptide capture cartridge upstream from the column (which is 1 – 1.5 ug estimated from the binding capacity of C18 with a bed volume of 0.35____m____L, since the trap cartridges have dimensions of 300____m____m ID and 5mm length; ____http://tools.thermofisher.com/content/sfs/manuals/Man-M5001-LC-Nano-Capillary-Micro-Columns-ManM5001-EN.pdf____ and ____https://www.optimizetech.com/opti-pak-trap-columns/____ ). Consequently, the different initial loading of oocyte vs cumulus cell proteins/peptides are likely to have made little if any contribution to proteome coverage, since 2-17____m____g all exceed the trap cartridge binding capacity, and consequently 1 - 1.5____m____g was captured and transferred to the nano-column, while the excess was transferred to waste. Based on the capacity limits of the capture cartridge, there was likely enough peptides/proteins in both oocyte and cumulus cell extracts to reach the saturation point, and therefore much more consistent on-column loadings than the initial ____m____g loadings would imply. We have added some additional information re this to the method section (see the section “Mass Spectrometry for proteomic analysis” in the Supplementary Information).
(2) The expressed proteomes of different cell types may be expected to differ not only in specific proteins expressed but also in the number of different proteins. In a recent study by Marei et al ____[9]____, equal amounts of total protein (9.5ug) from bovine oocytes and matching cumulus cells were prepared for their proteomics comparisons, and interestingly these authors also report about half as many proteins identified in oocytes as compared with cumulus cells, despite equal amounts of total protein used; “A total of 1703 and 1185 proteins were identified in CCs and oocytes, respectively, 679 of which were common.” Furthermore, a transcriptomic study of bovine oocytes and cumulus cells by Moorey et al ____[10]____, showed 69 and 128 differentially expressed genes (DEGs) in oocytes and cumulus cells respectively (comparing small vs large cells in each case), pointing to about double the differential gene expression in cumulus cells than oocytes, again implying a larger cumulus cell vs oocyte transcriptome. Our data support these observations, which collectively suggest a real difference in proteome size between oocytes and cumulus cells. If the difference in proteome size is real, then differences in network enrichment are also likely to have biological relevance, despite differences in size of the differentially expressed proteins lists.
(3) Even if initial protein loading was a contributing factor to the size of the oocyte vs cumulus cell proteomes, it is of note that we observed approximately 2 fold fewer total proteins identified in oocytes as in cumulus cells (Figure 3A, 3B and new Figure 4C), yet the difference between number of identified networks is multiple-fold (a cumulative total of 2 networks identified in BMP15 and cumulin treated oocytes vs 143 networks identified in BMP15 and cumulin treated cumulus cells – see new Figure 4C). Furthermore, there does not seem to be a strictly linear relationship between the number of proteins submitted for network analysis and the numbers of networks identified. For example, 34 upregulated proteins in cumulin treated oocytes identified a single enriched network, while a similar number of 38 upregulated proteins in BMP15 treated cumulus cells, identified a total of 42 networks (new Figure 4C), and similarly cumulin treated cumulus cells had 59 upregulated proteins and 58 downregulated proteins, which resulted in 57 and 23 enriched networks respectively.
Also, when performing GO terms analysis, the level of "branching" can explain the results. In other words, GO terms are organized in a tree like structure where general elements (e.g. nucleus) are delineated in finer elements (e.g. nuclear function) leading to finer ones (e.g. DNA binding)... to finer ones (e.g. DNA repair)... etc. The number of genes/proteins available in the initial list directly dictates to which level of precision the analysis can reach. In the present work, the number of identified network may simply reflect the number of elements available in the initial lists. With more info on the cumulus cells side, it is logical to be able to reach finer branches that contain only a few genes. I have looked in the supplemental data files but could not find more info about the background used. Was it all known proteins? Was it all identified proteins where the differentially expressed proteins are compared to the detected proteins? Using the list of detected proteins as background for the analysis could help. Proteome Discoverer generated much less differentially expressed proteins between treatments than Mascot/Scaffold (2-17 vs. 74-390). Maybe use the Mascot/Scaffold data using the same number of top genes (e.g. 87) between both cell types. Then it would be much more comparable.
RESPONSE #43 Please also refer to the explanations under RESPONSE #34 and #42 above. We have added an additional explanation of how we performed the enrichment analysis, in the supplementary information section under the heading “Gene Ontology Network Enrichment Analysis”. In the data presented here we used the whole mouse genome as our background set. The number of total proteins identified by Mascot/Scaffold and ProteomeDiscoverer were similar, but actually considerably more differentially expressed proteins were identified using ProteomeDiscoverer (Fig 3A), as expected using peak intensity vs spectral counting ____[11]____. The spectral counting approaches usually identify fewer differentially expressed proteins, but also with a lower quantitative false positive rate, while peak intensity approaches tend to identify more differentially expressed proteins, but with a higher quantitative false positive rate ____[11]____. Our reasoning was therefore to combine proteins which vary in common across both platforms, to maximise the differentially expressed proteins list while simultaneously minimising the quantitative false positive rate. We thank the reviewer for the suggestion of using our full protein list as the background set. Initially we revised our network enrichment analysis (see comments under RESPONSE #25) still using the mouse whole genome, resulting in fewer overall networks, but improved contrast between oocytes and cumulus cells (see summary in new Figure 4C, and network analysis details in supplementary Tables S12 and S13). We then repeated the network analyses using our full protein list (4450 proteins identified in both oocytes and cumulus cells; see background list in Supplementary Table S11) as the background set. With this we similarly found no enriched GO networks for BMP15 and cumulin treated oocytes, and only 6 and 1 enriched network in BMP15 and cumulin treated cumulus cells. We suggest that detecting network enrichment against a cell specific background list may not give us the same level of “contrast” as can be achieved when comparing against the whole mouse genome.
Line 226 and 324-328 and line 350: I have never seen the use of MitoTracker Orange to count mitochondria. According to the manufacturer: "MitoTracker Orange CMTMRos is an orange-fluorescent dye that stains mitochondria in live cells and its accumulation is dependent upon membrane potential. The dye is well-retained after aldehyde fixation." It is indicative of mitochondrial potential but it is not a method to count the number of mitochondria within a cell. I do not agree that more fluorescence means more mitochondria.
RESPONSE #44 We agree and in places we used ambiguous language re mitochondrial function vs mitochondrial number. We have now clarified and corrected this - please refer to detailed comments and manuscript changes under RESPONSE #8.
I understand that the MitoTracker data is counterintuitive to the oxygen consumption rate and stable levels of energetic metabolites. However, as the authors mention, mitochondria are known to be capable of switching from aerobic to anaerobic energy production. In some cases, heterogeneity in the mitochondrial population (such as the one in the oocyte) could mean that a mosaic respiratory potential exists where some mitochondria are more aerobic than others... To change the number of mitochondria, either fission or mitophagy must occur. Although mitochondrial DNA replication is done in approximatively 2 h and fission/division can occur over 1 h, and protein ubiquitination is done over 12 h-18 h during mitophagy, TEM micrographs (figure 5) do not show elongated mitochondria in the process of division. To detect active mitophagy, protein markers and association with lysosome would be needed. A shift in mitochondrial number may not be the suitable interpretation of the data.
RESPONSE #45 Please refer to comments under RESPONSE #8
For the spectral data analysis (Figure 3D), how did the three replicates perform? The figure does not show the replication variance relative to the treatment variance.
RESPONSE #46 A version of Figure 3D but with the replicates colour-coded has been added to Supplementary Material (Supplementary Figure 2) and the manuscript text has been revised with the information that data from the three replicates are shown, added to the caption to Figure 3D.
Wording/interpretation issues
Lines 114-116: "This intercellular cooperativity facilitates oocyte maturation while simultaneously protecting germ-line genomic integrity, in a manner which could not be achieved by a single cell." This is an overstatement because genomic integrity was not assessed. Why consider that the nuclear function found in the proteome contrast is necessarily associated with genomic integrity. Miosis requires in dept chromatin handling. What evidence provided from the results is associated with cellular numbers. The presence of cumulus cells is known to support meiosis but it doesn't mean that some of the cellular processes have been imparted to the surrounding somatic cells. The work done for this manuscript does not test any of this claim.
RESPONSE #47 We accept this point and agree that, especially the claim re germ-line genomic integrity, is an overstatement. This has been removed. We maintain however that there is ample evidence in our results that there is clear inter-cellular metabolic cooperativity between oocyte and cumulus cells and that this ultimately leads to an oocyte with improved developmental competence. The sentence has been modified to reflect this, line 117-118.
On numerous occasions, the statements are imprecise. For example: Line 274: "More than double the number..." Since doubling a minute value does not mean the same thing as doubling a large value, values, measurements with units and ideally with SEM should be added.
RESPONSE #48 Has been rephrased (see line 284 of the revised manuscript)
Line 287: "... and almost a third more significant networks..." Please add values.
RESPONSE #49 Section has been deleted (line 291-300 of the revised manuscript)
On the same statement, since sample input material to the mass spectrometry is vastly different between cumulus cells and oocytes, is it truly comparable? Could these differences between the two cell types be associated with the amounts of proteins in the extracted samples? Typically, more variable results are obtained with the low input. It sometimes lead to apparently more difference between treatments simply because of low count numbers. On line 292, authors mentioned that protein loading was considered. How was that done? Low input cannot be compensated or normalized. The following statement on line 293 indicate that more proteins were identified in cumulus cells. This is probably due to more input material submitted to mass spectrometry. It is not necessarily a difference in protein diversity between cumulus cells and oocytes.
RESPONSE #50 Please refer to detailed explanations under RESPONSES #42 and #43
Line 293: "... resulted in the identification of about double the number..." Please add values.
RESPONSE #51 Values added at lines 305-306, and additional detail has been added to this section of the manuscript (lines 305-317 revised manuscript).
Line 294: "However, there were 4-5 times as many differentially expressed proteins..." Please add values.
RESPONSE #52 Values added and additional detail added to this section (new lines 309-312 of the revised manuscript).
Line 298: "...difference was quite marked..." More factual info should be added.
RESPONSE #53 Values added and additional detail has been provided (lines 314-317 of the revised manuscript), as follows; “____Cumulin appeared to have a greater impact on proteomic differential expression in both cell types than BMP15 did, with 59 vs 38 and 34 vs 27 upregulated proteins in cumulin vs BMP15 treated cumulus cells and oocytes respectively, and similarly 14 vs 6 downregulated proteins in cumulin vs BMP15 treated cumulus cells and oocytes respectively (Figure 4C)”.
Line 305: Again, the whole comparison between cell types could be argued from the standpoint of input material subjected to the analysis. Given the point is to state that cumulin has a profound impact on cumulus cells, maybe it is not necessary to compare with the oocyte data. It is logical that an oocyte secreted factor targets the neighbouring cells. The point can be made without raising the question about the potential issue of input material.
RESPONSE #54 We agree with the reviewers point that it is logical that OSFs should target cumulus cells, with lesser impact on the oocyte. Nonetheless the treatments were performed on COCs, and even though the OSFs are targeting the cumulus cells, however ultimately the cumulus cells response is expected to impact oocytes. Therefore, it is relevant to look at proteomic changes to both cell types and also the related network analysis. We have however rephrased this section, to be more specific as to which data we are reporting, and have included additional citations (lines 325-334 of the revised manuscript).
__Line 317-317: "... exhibited more rounded and swollen mitochondria..." How was that determined? In the periphery of the oolemma, mitochondria aggregates in clusters which can be quite different from one another. Maybe proportions of different shapes of mitochondria could be provided if enough mitochondria are counted from the EM micrographs.
__RESPONSE #55 These are subjective observations of the typical morphological features seen in response to the different treatments. This is the typical application of TEM. Quantitative features of mitochondria are better assessed using confocal than TEM, which is the complimentary approach we took using MitroTracker in the companion figure 5B, the text for which immediately follows the TEM results. We altered the text at the sentence in question to note that these are subjective observations (line 340).
Line 169: What do you mean by "The results were merged based on consistency..."? This seems to be a trivial way to analyse the data.
RESPONSE #56 The majority of published papers reporting data dependent analysis (DDA) proteomics results utilise just a single quantitative method (i.e., either spectral counting or peak intensity). This certainly simplifies reporting, and avoids confronting uncomfortable discrepancies between different analytical approaches. However, we reasoned that robust expression change data would maintain consistency, despite the orthogonal quantitative methods. We consider it a notable strength of the approach used here that we have utilised a differentially expressed proteins list which includes only those proteins with consistent direction of fold-change in both the spectral counting and peak intensity workflows. Please also refer to comments under RESPONSE #43, re spectral counting vs peak intensity quantitative methods in data dependent analysis (DDA) proteomics.
Line 170: "A further requirement was that at least one, if not both methods..." Again, when did you decide to use one method or to use both? Why not use the common ground from both methods?
RESPONSE #57 Refer also to RESPONSE #43. In fact the main question being asked in many/most proteomics experiments is whether there is a real expression change between treatment groups. Therefore fold-change is the most pertinent common ground across disparate quantitative methodology, and indeed commonality of fold-change was the basis for merging the datasets. Since integrating peak areas is a very different approach to counting the number of spectra, then this difference in approach can make a big difference to the p-values, and is the reason why spectral counting is less sensitive to detect differential expression. For similar reasons the fold-change ratio may differ somewhat between these quantitative methodologies. However direction of fold-change is a minimum requirement for demonstrating consistent trends, hence we used this as the common ground for merging the datasets.
Line 384: Is the paracrine signaling remodeling COC metabolism or is it enhancing the rate at which it is done? I believe this switch in metabolism occurs in untreated COCs.
RESPONSE #58 We see the reviewers point in this subtle difference in wording. We agree that there is a switch in metabolism in untreated COCs during maturation – our point is that that process of changing metabolism is further remodelled by oocyte paracrine signals, to the overall betterment of the oocyte in terms of competence. We have edited this sentence to make this point clearer (line 413-415). Our data on energy charge, respiration, energy metabolite levels (Figure 6), redox potential (Figure 7) and mitotracker intensity (Figure 5) are all presented in comparison with “untreated” cells, and our conclusion that there is remodelling of metabolism is therefore relative to “untreated” COCs.
__
__The Discussion is somewhat circular. Section will need to be adjusted if the Mitotracker-based mitochondrial count and between cell types gene/protein lists comparisons are removed.
Accounts for mitochondrial counts: (lines 387-393) (lines 424-427) (line 463).
RESPONSE #59 All reference to Mitotracker in the context of mitochondrial counts only have been altered to Mitotracker being an indicator of mitochondrial function/polarity and/or counts.
Accounts for comparisons of gene lists length between cell types: Lines 389-391 and 475-477 and 496-499).
RESPONSE #60 Please see comments under RESPONSE # 53 and the new Figure 4C.
Line 395: "... a substantial number of oocyte upregulated proteins... Please provide number.
RESPONSE #61 Additional specific proteins have been listed to support our claims of effects on specific processes (see lines 435-443 of the revised manuscript). Also see comments under RESPONSE # 7.
Line 397: The data was not designed to test the potential of cumulin to preserve meiotic fidelity. This is an overstatement since DNA binding is part of the normal course of even during meiosis. Again, cumulin could accelerate the kinetic of meiosis.
RESPONSE #62 Reference to meiosis and meiotic fidelity removed, line 435.
__
__Line 402-405: the work was not designed to determine if cumulin would shift work allocation between the oocyte and the cumulus cells. Showing that cumulin drives meosis is interesting by itself.
__RESPONSE #63____ Not clear that any change is requested or needed. This sentence is interpreting the significance of the results, as required in a Discussion.
__Line 453-455: the link with the epigenome is an overstatement. RNA and DNA processing pathways are general cellular processes.
RESPONSE #64____ The link to the epigenome was a reference to some published work. However it was linked to observations in the current data, and additional information has now been added to the updated manuscript to explain this further, as follows (currently lines 509-516);
"These included significantly enriched networks of RNA binding, helicase activity, ribonucleoprotein complex biogenesis, and mRNA processing (supplementary Tables S11 and S12; upregulated proteins RNF20, SHMT1, DHX30, DDX17, DDX5, PRPF19, RPS4X, NOP58, DDX39b, HNRNPF, RPS271, NOP56, PRPF6, POLR2b, CPSF6, OOEP), as well as upregulation of key epigenetic regulators (HDAC2 and UHRF1; see supplementary Table S11), histone modifying protein MTA2, and significant network enrichment of the spliceosomal complex (supplementary Table S12; proteins DDX5, PRPF19, HNRNPF, PRPF6, POLR2B), which has been linked to epigenetic regulation ____[12]____.
Minor details
Line 36: I suggest to be more precise on the "nuclear" function that is affected in the oocyte. Given that oocytes are transcriptionaly quiescent at this stage, some might argue that it is a vague statement.
RESPONSE #65____ Information relating to specific oocyte upregulated proteins and their cellular roles has been added to the updated manuscript (currently lines 434-443).
DNA binding and ribosomal constituents (Fig. 4A, 4C),
In vitro should be in italic because it is Latin.
RESPONSE #6____6 corrected throughout
__Lines 125-126: are the batch numbers relevant to anything?
__RESPONSE #6____7 We would assume so – for the historical record. These are in-house produced proteins, cumulin is complex to produce and only a few labs worldwide have made it.
__Line 168: Mascor = Mascot
__RESPONSE #6____8 Corrected
__Line 168: a reference for the software?
__RESPONSE #6____9 URL and published references added (lines 172-175 revised manuscript)
Line 178: need a reference for the software?
RESPONSE #70 URL and published references added (line 185)
__Line 187: Need a complete source for "Procure, 812"
__RESPONSE #71 Added
Line 188: Need a complete source for "Diatome"
RESPONSE #72 Added
Line 197: Need a complete source for "Cell-Tak"
RESPONSE #73 Added
Line 232: though = through
RESPONSE #74 Corrected
Line 243: define OCR
RESPONSE #75 Added
Line 268: If I am not mistaking, it is not a multispectral analysis. The multispectral values were analysed through a principal component analysis.
RESPONSE #7____6 Data was analysed through linear discriminative analysis (LDA). This information has been added in Line 278.
Line 363: What is the "behaviour" of an oocyte and cumulus cells?
RESPONSE #77 replaced with “function”
Line 512-513: Maybe add more on the fact that most clinics use ovulated eggs and do not perform IVM. However, IVM is needed is specific contexts such as PCOS.
RESPONSE #78 Edited accordingly; lines 575-577.
Reviewer #3 (Significance (Required)):
Cumulin is the most potent oocyte secreted factor. Its mecanism of action is still unknown.
I have been working on the mammalian oocyte for the past 25 years.
References
- Mester, B., et al., Oocyte expression, secretion and somatic cell interaction of mouse bone morphogenetic protein 15 during the peri-ovulatory period. Reprod Fertil Dev, 2015. 27(5): p. 801-11.
- Hussein, T.S., J.G. Thompson, and R.B. Gilchrist, Oocyte-secreted factors enhance oocyte developmental competence. Dev Biol, 2006. 296(2): p. 514-21.
- Mottershead, D.G., et al., Cumulin, an Oocyte-secreted Heterodimer of the Transforming Growth Factor-beta Family, Is a Potent Activator of Granulosa Cells and Improves Oocyte Quality. J Biol Chem, 2015. 290(39): p. 24007-20.
- Gilchrist, R.B., M. Lane, and J.G. Thompson, Oocyte-secreted factors: regulators of cumulus cell function and oocyte quality. Hum Reprod Update, 2008. 14(2): p. 159-77.
- Sugiura, K., F.L. Pendola, and J.J. Eppig, Oocyte control of metabolic cooperativity between oocytes and companion granulosa cells: energy metabolism. Dev Biol, 2005. 279(1): p. 20-30.
- Campbell, J.M., et al., Multispectral autofluorescence characteristics of reproductive aging in old and young mouse oocytes. Biogerontology, 2022. 23(2): p. 237-249.
- Schwarz, D.S. and M.D. Blower, The endoplasmic reticulum: structure, function and response to cellular signaling. Cell Mol Life Sci, 2016. 73(1): p. 79-94.
- Szklarczyk, D., et al., STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res, 2019. 47(D1): p. D607-D613.
- Marei, W.F.A., et al., Proteomic changes in oocytes after in vitro maturation in lipotoxic conditions are different from those in cumulus cells. Sci Rep, 2019. 9(1): p. 3673.
- Moorey, S.E., et al., Differential Transcript Profiles in Cumulus-Oocyte Complexes Originating from Pre-Ovulatory Follicles of Varied Physiological Maturity in Beef Cows. Genes (Basel), 2021. 12(6).
- Ramus, C., et al., Benchmarking quantitative label-free LC-MS data processing workflows using a complex spiked proteomic standard dataset. J Proteomics, 2016. 132: p. 51-62.
- Luco, R.F., et al., Epigenetics in alternative pre-mRNA splicing. Cell, 2011. 144(1): p. 16-26.