Next Article in Journal
Ethylene and Nitric Oxide Involvement in the Regulation of Fe and P Deficiency Responses in Dicotyledonous Plants
Next Article in Special Issue
SUL-151 Decreases Airway Neutrophilia as a Prophylactic and Therapeutic Treatment in Mice after Cigarette Smoke Exposure
Previous Article in Journal
Hypergravity Load Modulates Acetaminophen Nephrotoxicity via Endoplasmic Reticulum Stress in Association with Hepatic microRNA-122 Expression
Previous Article in Special Issue
Role of Endocrine-Disrupting Chemicals in the Pathogenesis of Non-Alcoholic Fatty Liver Disease: A Comprehensive Review
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 9
10.3390/ijms22094902
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fibroinflammatory Signatures Increase with Age in the Human Ovary and Follicular Fluid

by
Jordan H. Machlin
1,
Seth J. Barishansky
1,
John Kelsh
2,
Megan J. Larmore
3,
Brian W. Johnson
3,
Michele T. Pritchard
4,
Mary Ellen Pavone
1 and
Francesca E. Duncan
1,*
1
Department of Obstetrics and Gynecology, Feinberg School of Medicine, Northwestern University, Chicago, IL 60611, USA
2
Department of Anatomy and Cell Biology, University of Kansas Medical Center, Kansas City, KS 66160, USA
3
Department of Comparative Medicine, Histology and Imaging core, University of Washington, Seattle, WA 98195, USA
4
Department of Pharmacology, Toxicology, & Therapeutics, University of Kansas Medical Center, Kansas City, KS 66160, USA
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(9), 4902; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22094902
Submission received: 18 March 2021 / Revised: 29 April 2021 / Accepted: 30 April 2021 / Published: 5 May 2021
(This article belongs to the Special Issue Inflammation in Health and Disease: New Insights and Therapeutic Avenues)
Browse Figures
Versions Notes

Abstract

:
The female reproductive system ages before any other organ system in the body. This phenomenon can have tangible clinical implications leading to infertility, miscarriages, birth defects and systemic deterioration due to estrogen loss. “Fibroinflammation” is a hallmark of aging tissues; there is an increase in inflammatory cytokines and fibrotic tissue in the aging ovarian stroma. We systematically evaluated immunomodulatory factors in human follicular fluid, which, like the stroma, is a critical ovarian microenvironment directly influencing the oocyte. Using a cytokine antibody array, we identified a unique fibroinflammatory cytokine signature in follicular fluid across an aging series of women (27.7–44.8 years). This signature (IL-3, IL-7, IL-15, TGFβ1, TGFβ3 and MIP-1) increased with chronologic age, was inversely correlated to anti-Müllerian hormone (AMH) levels, and was independent of body mass index (BMI). We focused on one specific protein, TGFβ3, for further validation. By investigating this cytokine in human cumulus cells and ovarian tissue, we found that the age-dependent increase in TGFβ3 expression was unique to the ovarian stroma but not other ovarian sub-compartments. This study broadens our understanding of inflammaging in the female reproductive system and provides a defined fibroinflammatory aging signature in follicular fluid and molecular targets in the ovary with potential clinical utility.

1. Introduction

Aging is associated with a general decline of tissue and cellular function [1]. The reproductive system is unique in that it is the first organ system in the female body to exhibit overt signs of aging. Ovarian function ceases at menopause, which occurs at ~51 years of age [2]. However, even before this, fertility decreases when women are in their mid-30 s due to a decline in egg quantity and quality [2]. Female reproductive aging has important clinical implications because it can lead to infertility, miscarriages, birth defects and systemic deterioration due to reduced estrogen [3]. Evidence from assisted reproductive technology (ART) cycles strongly suggests that the age-dependent decrease in fertility is largely due to defects at the level of the egg. For example, if women use their own eggs to conceive, there is a strong maternal age effect with the likelihood of women in their mid-forties having a live birth being close to zero percent [4]. However, if women use donor eggs from young, healthy women to conceive, the maternal age effect is essentially abrogated [5]. In fact, there are significant age-related changes in the female mammalian gamete, including defects in chromosome structure and cohesion, microtubule and kinetochore function, mitochondrial function and bioenergetics, nucleolar architecture and ribosome numbers, inorganic signaling, gene expression, and epigenetics [6,7,8,9]. However, whether these changes are inherent to the gamete or due to indirect effects of the aging microenvironment are not known.
Fibroinflammation is a hallmark of aging tissues [10,11]. Inflammation is a natural immune response to defend against pathogens or sterile injury, is typically short-lived, and ultimately resolves [12]. However, if acute inflammation persists and becomes chronic, tissue damage and disease pathogenesis can ensue [13]. Chronic inflammation causes fibrosis in which normal wound healing processes are skewed towards overproduction of extracellular matrix and accumulation of scar tissue [14]. We and others have recently demonstrated that the aging ovarian microenvironment or stroma becomes inflammatory and fibrotic with age [15,16]. Reproductive aging in the ovary is associated with increased production and secretion of proinflammatory and pro-fibrotic cytokines and growth factors and a shift in macrophage populations to multinucleated macrophage giant cells [15,17]. In some mouse strains, this inflammation is accompanied by fibrosis characterized by increased collagen I and III levels in mice of advanced reproductive age [15]. Fibrosis and additional changes in the extracellular matrix are also conserved in the aging human ovary, demonstrating the translational relevance of this aging phenomenon [16,18,19].
Although these prior studies focused on the inflammatory and fibrotic nature of the aging ovarian stroma, another critical ovarian microenvironment is the follicular fluid (FF). FF directly surrounds the oocyte as it grows and matures within the cumulus oocyte complex of the antral follicle. The FF is derived from blood plasma macromolecules flowing through the thecal capillaries that cross the blood-follicle barrier and from granulosa and thecal cell secretions themselves [20,21,22]. FF consists of large polysaccharides, hormones, inflammatory factors, reactive oxygen species, and growth factors [21,23]. FF plays a critical role in the development and the eventual quality of the oocyte, and its composition likely reflects the larger circulating biochemical environment [22]. Additionally, the cytokine composition of FF has been shown to vary with age [24].
Given the age-dependent increase in inflammation and fibrosis in the ovarian stroma, we hypothesized that the follicular fluid milieu also becomes fibroinflammatory with age. To test this hypothesis, we used cytokine antibody arrays to semiquantitatively profile the expression of 80 fibroinflammatory cytokines in human follicular fluid from women across an aging series (27.7–44.8 years). We identified a cytokine signature (IL-3, IL-7, IL-15, TGFβ1, TGFβ3 and MIP-1) that increases with chronologic age is inversely correlated to anti-Müllerian hormone (AMH) levels and is independent of body mass index (BMI). Thus, this represents a unique and specific fibroinflammatory reproductive aging cytokine profile in human follicular fluid. We then performed a comprehensive analysis of one fibroinflammatory factor, TGFβ3, in human cumulus cells and ovarian tissue. We found that it increased specifically in the ovarian stroma with age, which may suggest a source for its corresponding increased levels in follicular fluid. These findings advance our understanding of how the human ovarian microenvironment—both the follicular fluid and stromal milieu—change with age and may impact the gametes and follicles that grow within it.

2. Results

2.1. Chronologic and Reproductive Age Are Associated with an Increase in Specific Subsets of Fibroinflammatory Cytokines

To examine whether fibroinflammatory cytokines increase in human follicular fluid with age, we obtained samples from 30 participants undergoing oocyte aspiration for ART. We interrogated aging from two perspectives: chronologic (i.e., biologic age) as well as reproductive age. These women ranged in age from 27.7 years to 44.8 years (average: 35.7 ± 0.9 y) (Table 1). Anti-Müllerian hormone (AMH) is a hormone produced by growing follicles and is often used as an indirect marker of the ovarian reserve and thus a marker of reproductive age [25,26,27]. In our study, population, serum AMH levels ranged from 0.1 ng/mL to 6.2 ng/mL (average: 2.6 ± 0.33 ng/mL) (Table 1). As expected, because follicle number decreases with age, there was an inverse correlation between age and AMH (Figure 1A) (r2 = 0.16, p = 0.03). Given that we excluded participants with a BMI ≥ 30 kg/m2, the BMI in our participants ranged from 17.5 kg/m2 to 29.6 kg/m2 (average 23.3 ± 0.57 kg/m2) (Table 1). There was no correlation between age and BMI in this population (Figure 1B) (r2 = 0.001, p = 0.86).
Using the follicular fluid samples from these participants, we performed an antibody array that allowed the semiquantitative detection of 80 cytokines involved in fibroinflammatory pathways (Figure 1C, Supplemental Figure S2). Of the cytokines interrogated, 61 had levels above a background intensity threshold. To determine whether there was a relationship between these cytokine levels and age, we plotted the intensity values versus chronologic age (years) and reproductive age (AMH). Six cytokines, IL-3 (p = 0.0022), IL-7 (p = 0.0170), IL-15 (p = 0.0084), TGFβ1 (p = 0.0108), TGFβ3 (p = 0.0193), and MIP-1β (p = 0.0315), all showed a positive correlation with chronologic age (Figure 2A). Conversely, the same cytokines were negatively correlated with AMH: IL-3 (p = 0.0457), IL-7 (p = 0.0154), IL-15 (p = 0.0127), TGF-β1 (p = 0.0296), TGF-β3 (p = 0.0461), and MIP-1β (p = 0.0413) (Figure 2B). However, these cytokines did not show any association with BMI (Supplemental Figure S3). Thus, these results demonstrate a unique subset of cytokines that are tightly associated with advanced age defined in two different ways (Table 2). An additional subset of 16 cytokines exhibited a positive correlation with age but not AMH (Table 2, Supplemental Figure S4, and data not shown). The remaining 39 cytokines were not correlated with age (Table 2 and Supplemental Figure S5). Of note, there were no cytokines that showed a significant correlation with AMH and not age (data not shown).

2.2. Cytokines That Increase with BMI Are Unique from Those That Change with Age

Although we selected an upper BMI limit of ≤ 30 kg/m2 to avoid confounding effects of obesity on the fibroinflammatory profile in human follicular fluid, there were two cytokines, leptin (p = 0.0048) and IL-8 (p = 0.0173), that nevertheless were positively correlated with BMI in our sample population (Figure 3A). These cytokines did not correlate with age or AMH (Figure 3B,C). These results validate the efficacy of our cytokine antibody array approach because leptin was previously shown to increase follicular fluid with BMI [28,29]. Moreover, these results demonstrate that the cytokine signature that changes with age is independent of that related to BMI.

2.3. TGFβ3 Content in Human Cumulus Cells Does Not Increase with Age

To validate and extend our findings in human follicular fluid, we prioritized a more in depth analysis of TGFβ3, since knowledge of this specific TGFβ isoform is lacking in the human ovary, especially in the context of aging. Transforming growth factor-beta (TGFβ) is a family of proteins with unique immunoregulatory and fibrosis-fostering properties, and in the human, there are three isoforms: TGFβ-1, 2, and 3 [30]. Normally the actions of TGFβ cytokines contribute to tissue repair, and TGFβ action is shut down through regulatory mechanisms that feedback to limit excess TGFβ signaling. However, TGFβ signaling resulting in excessive ECM accumulation can cause fibrosis and scarring [31,32]. Several pathologies stem from TGFβ dysregulation in various organs, including kidney disease, lung, liver and cardiac fibrosis and dermal scarring [31,33]. TGFβ dysregulation has been implicated in female reproductive system pathologies, including PCOS and the pathogenesis of endometriosis and uterine fibroids [32,34].
Components of human follicular fluid are derived in part from granulosa cell secretions [22]. Therefore, to determine if the increase in TGFβ3 observed in follicular fluid was paralleled in the somatic compartment of the follicle, we performed immunoblot analysis on intact cumulus cell masses (Figure 4A). As a positive control to validate the antibody, we ran protein extracts from MCF7 cells, a hormone-responsive breast cancer cell line known to express TGFβ3, alongside protein extracts from human cumulus cells [35]. Like the control, TGFβ3 was detected as a single immunoreactive band at 50 kDa in whole cumulus cell masses (Figure 4B). To examine whether there were age-dependent changes in TGFβ3 expression in the cumulus cells, we performed immunoblot analysis on protein extracts from cumulus cell masses from reproductively young (<35 y) and old (≥35 y) participants (Figure 4C–E). To minimize confounding variables, we analyzed cumulus cells from reproductively young and old participants that surrounded oocytes that produced the same embryo outcomes on day 5 or day 6, indicating similar quality (Table 3). Densitometry was performed to determine band intensities of the young and old age cohorts, and we normalized TGFβ3 expression to GAPDH. Although the ratio of TGFβ3:GAPDH tended to increase with age, this was not significant (Figure 4D). Moreover, the fold change of TGFβ3 expression between reproductively young and old samples was not consistent across experimental replicates, suggesting that there was no age-dependent change in cumulus cell expression of TGFβ3 (Figure 4E).

2.4. TGFβ3 Is Expressed in Follicles, Vasculature and Stroma in the Human Ovary

Although we did not see age-associated changes in TGFβ3 expression in cumulus cells, it is possible that its expression increased in the ovarian microenvironment especially given that TGFβs are anchored in the extracellular matrix (ECM) in a latent form. Thus, we next examined TGFβ3 localization in the human ovary via immunohistochemistry on tissue sections. TGFβ3 expression was observed throughout the ovarian stroma with more dense staining at the cortex than the medulla (Figure 5A). This staining, except for that at the tunica albuginea, was attenuated when the tissue was stained with the same concentration of nonimmune IgG, confirming specificity (Figure 5B). A further examination of sub-compartments within the ovary demonstrated that TGFβ3 is specifically expressed in follicles, surrounding the vasculature, and throughout the stroma (Figure 5C–E).

2.5. Stromal and Vascular TGFβ3 Content Increases with Age in the Human Ovary

To investigate age-associated differences in ovarian TGFβ3 expression, we performed immunohistochemistry using a TMA containing cortical tissue from 60 participants in two age cohorts: 1–20 y and 39–58 y [16]. (Supplemental Figure S1). TGFβ3 staining qualitatively appeared more intense in the samples from the older cohort, so we performed further quantification (Supplemental Figure S1B). First, we created two additional cohorts and thus quantified expression in four age cohorts across an aging continuum: 0–10, 11–20, 39–50, 51–60 (years) using an average of both the core 1 and core 2 data per participant (Figure 6A). There was a trend towards increased TGFβ3 expression across these four age cohorts (Figure 6A). This increased expression was significant when comparing TGFβ3 content between the 0–20 y and 39–60 y (p = 0.012) and 0–10 y and 51–60 y (p = 0.006) cohorts (Figure 6B).
Although each of the cores from the TMA was obtained from the pathologist-confirmed ovarian cortex, there was heterogeneity in the tissue cores. Some cores contained follicles, others had blood vessels, and there were varying amounts of stroma throughout. We, therefore, analyzed TGFβ3 content taking into account the presence of dominant structures within the tissue. Based on H&E staining, cores were classified into those containing follicles, blood vessels, or stroma. We plotted the average TGFβ3 per total area across the same four age cohorts, stratified by the three ovarian structures (Figure 7). Cores containing follicles showed a trend towards a decrease in TGFβ3 levels with age, but this analysis is limited by the age-dependent decrease in the number of cores containing follicles (Figure 7A). In contrast, cores containing blood vessels and stroma demonstrated increased TGFβ3 content with age when comparing the 0–10 y vs. 51–60 y cohorts (Figure 7B,C, p = 0.042 blood vessels and p = 0.004 stroma). Therefore, the age-associated increase in TGFβ3 expression in the ovary occurs in tissue regions containing stroma and vasculature.

3. Discussion

In this study we interrogated the fibroinflammatory signature of human follicular fluid. Both maternal age and elevated BMI can negatively impact fertility and affect the circulating levels of cytokines, but the precise mechanisms by which these factors impact oocyte quality are not fully elucidated [36,37,38,39,40,41,42,43]. We identified six cytokines, IL-3, IL-7, IL-15, TGFβ1, TGFβ3, and MIP-1, which were positively correlated with chronologic age and inversely correlated with AMH, an indirect marker of the ovarian reserve or reproductive age [10]. These specific cytokines, however, were not correlated with BMI. Obesity alters the levels of certain adipokines, and although there was an upper BMI limit of ≤30 kg/m2 in our study, we identified leptin and IL-8 as being positively correlated with BMI. These findings are consistent with others, demonstrating that leptin is a principal cytokine that increases with BMI in human follicular fluid [28,29]. Leptin is produced by granulosa and cumulus cells of ovarian follicles [44] and can impair follicle development, ovulation, and oocyte maturation [45,46]. Like leptin, IL-8 is also part of the immunologic profile associated with excessive body weight [47]. The observation that IL-8 and leptin levels correlated with BMI, but not age and that IL-3, IL-7, IL-15, TGFβ1, TGFβ3, and MIP-1 levels correlated with age, but not BMI suggest distinct fibroinflammatory responses in obesity and aging, and thereby unique mechanisms by which these factors negatively impact oocyte quality.
There are several potential origins of the unique fibroinflammatory signature seen with aging. These cytokines could be due to the known age-associated shift of macrophage plasticity and ontogeny in the ovary [48]. For example, ovarian aging is associated with increased multinucleated macrophage giant cells, which may be responsible for a highly inflammatory cytokine milieu [15,17]. An alternative possibility is that senescent ovarian cells, which may or may not be immune cells, may produce an inflammatory milieu. “Inflammaging,” a leading theory on aging, refers to the low-grade, chronic inflammation that underlies aging tissues and pathologies in the absence of overt infection [49,50]. Senescent cells, which increase with age, are a key source of inflammaging [51]. Cellular senescence is a process in which cells cease dividing but remain metabolically active and undergo distinctive phenotypic alterations, fueling aging [52]. These changes include secretion of cytokines, chemokines, growth factors, and proteases, which are collectively referred to as the senescence-associated secretory phenotype (SASP) [53].
Our data are consistent with the well-documented SASP in other tissues by way of composition. Specifically, two of the cytokines that we found to increase with age in follicular fluid, IL-7 and IL-15, are well-characterized senescent cell markers [54]. IL-7 is implicated in aging due to its role in the generation and maintenance of naive and memory CD8 and CD4 T cells [55,56]. Increased concentrations of IL-15 in follicular fluid from women undergoing in vitro fertilization are negatively correlated with achieving a clinical pregnancy [57]. MIP-1β, another fibroinflammatory chemokine that increased with age in our samples, is also a potential biomarker for inflammaging [58]. On the other hand, we did not identify an age-dependent increase in other conserved SASP components, including IL-6, IL-8, GROα, or GROβ. This could be because SASP composition is diverse and highly tissue and cell-type-dependent, suggesting that the ovary may have a unique signature [54,59,60]. However, future studies are necessary and ongoing to determine whether senescent cells exist and drive inflammaging in the ovary.
Our results also provide insight into other chronic inflammatory regulators that can induce specific fibrotic pathways, such as those regulated by the transforming growth factor-beta (TGFβ) family. The three TGFβ family isoforms have unique immunoregulatory properties that normally contribute to tissue repair [30]. However, excess action of TGFβ involving surplus ECM deposition can cause fibrosis and scarring [31,32]. This type of pathogenic dysregulation is more widely documented in the TGFβ1 and 2 isoforms [31,32,61]. Interestingly, our results demonstrated that TGFβ3 increased specifically with age, was inversely correlated with AMH, and was independent of BMI. Although TGFβ3 is considered a proliferative growth factor and a scar-preventive protein, there is evidence that excessive TGFβ3 signaling can cause tissue fibrosis and disease progression [62,63,64]. For example, TGFβ3 levels can be elevated 3–5 times more in uterine fibroids than in normal myometrium [62].
In the human ovary, TGFβ3 localized throughout the stroma, blood vessels, and surrounding follicles. These findings are consistent with previous observations in the bovine ovary, where TGFβ isoforms are expressed in follicles and likely act as important paracrine and autocrine signaling molecules regulating ovarian follicle growth and proliferation [65]. Interestingly, TGFβ3 expression decreased with age in histological sections of the ovarian cortex that contained follicles. Since this factor is important for growth and proliferation, a reduced expression with age may be correlated with changes in follicle quality. Finally, we observed increased expression of TGFβ3 with age in tissue sections containing stroma and vasculature but not in the cumulus granulosa cells. The TGFβ family has angiogenic and angiostatic properties depending on expression levels and the tissue context [66]. TGFβ3 may propagate fibroinflammation by increasing vasculature giving rise to additional sources where cytokines and inflammatory molecules may accumulate in the ovary with advancing age. Interestingly, we also saw increased vascular endothelial growth factor (VEGF) expression with age in human follicular fluid. This is consistent with data that women of advanced reproductive age undergoing IVF have higher FF levels of VEGF, which could reflect a hypoxic follicular environment [67]. Thus, increased vasculature and pro-angiogenic factors may act as a compensatory mechanism for promoting blood supply and oxygen through the dense fibrotic tissue.
Overall, our findings advance our knowledge of age-related inflammation in the ovary by defining a unique subset of fibroinflammatory cytokines associated with aging. Whether these cytokines are a cause or consequence of female reproductive aging is not known. The TGFβ pathway is currently being investigated as a therapeutic target for pathologies in the reproductive system, and our findings suggest that these strategies may also have relevance in ovarian aging. However, future studies are needed to decipher the relationship between the TGFβ3 that increases with age in ovarian tissue and in the follicular fluid. Translational studies are also ongoing to determine whether these particular cytokines influence gamete quality and whether these effects are direct on the oocyte or indirect via the supporting somatic cells. The fibroinflammatory cytokine/chemokine signature in the follicular fluid may also be an important biomarker of the quality and developmental potential of the oocyte, which may have predictive value for assisted reproductive technologies outcomes. Moreover, the aging follicular fluid signature may reflect systemic aging, which would broaden our understanding of fertility as a marker of overall health.

4. Materials and Methods

4.1. Human Follicular Fluid, Cumulus Cells, and Ovarian Tissue

Follicular fluid was obtained from women undergoing oocyte retrieval for a noncancerous diagnosis at Fertility and Reproductive Medicine (FRM) (Northwestern University, Chicago, IL, USA) through the Northwestern University Reproductive Tissue Library (NU-RTL), an Institutional Review Board-approved tissue repository protocol, following written informed consent. Inclusion criteria were women undergoing oocyte retrieval ≥ 18 y, ≤30 kg/m2, and English speaking. Participants with multiple controlled ovarian stimulation cycles, a history of reproductive pathologies or other conditions associated with inflammation were excluded. To avoid blood contamination, follicular fluid aspirated only from the first follicle from each ovary was saved from each participant. Upon collection, a Halt protease inhibitor cocktail (Thermo Fisher Scientific, Waltham, MA, USA) was added to the follicular fluid samples to a final concentration of 1×. The follicular fluid was then centrifuged at 300× g for 20 min at 4 °C to pellet cellular debris. The supernatant was then transferred to sterile tubes, snap-frozen in liquid nitrogen and stored at −80 °C until utilized in the cytokine antibody arrays. In addition to follicular fluid, participant-specific health information was obtained, including age, body mass index (BMI), and serum levels of anti-Müllerian hormone (AMH) (Table 1).
We also obtained human cumulus cells that would have otherwise been discarded from women undergoing oocyte retrieval at the FRM clinic through the NU-RTL. The inclusion and exclusion criteria were the same as above, but for cumulus cells, participants had to be undergoing intracytoplasmic sperm injection and preimplantation genetic testing. The oocytes retrieved were used clinically, whereas a portion of the cumulus cells from the larger cumulus–oocyte complexes (COC) was designated for research purposes. Each cumulus mass was assigned a number corresponding to the oocyte to track the oocyte and its respective embryo outcome in a one-to-one manner. All cumulus cells were transferred from the clinic in 4-well dishes containing HEPES buffered human tubal fluid (HTF) media (Irvine, CA, USA, Catalog #90,126) and placed in an incubator (37 °C, 5% CO2) before processing. Each cumulus cell mass was imaged on an EVOS FL auto microscope (Life Technologies Corporation, Grand Island, NY, USA). Next, cumulus cell masses were rinsed through six large drops of Leibovitz’s L-15 Medium (Thermo Fisher Scientific, Waltham, MA, USA) supplemented with 3 mg/mL polyvinylpyrrolidone (PVP) and 0.5% penicillin–streptomycin (PS) to remove any exogenous protein from the HTF media. Cumulus cells were then snap-frozen in a minimal volume of L-15 medium and stored at −80 °C. Deidentified information for each participant was obtained, including age, diagnosis, cumulus cell number corresponding to the oocyte retrieved, the oocyte’s maturity status and embryo outcome (Supplementary Table S1). We received an average of 7 cumulus cell masses per participant and utilized cumulus cells from a total of 12 participants samples for this study.
De-identified, formalin-fixed, paraffin-embedded (FFPE) human ovarian tissue samples used for this study were obtained through two established tissue archives at Northwestern University: the Oncofertility Consortium’s National Physicians Cooperative (OC-NPC) and the NU-RTL. All procedures for tissue collection were performed under IRB-approved protocols after obtaining informed consent from each patient. Tissues obtained through the OC-NPC were from patients undergoing ovarian tissue cryopreservation for fertility preservation purposes. Under these protocols, 80% of the tissue was cryopreserved for future patient use to restore endocrine function and/or fertility, and 20% was donated to research. A portion (10%) of the research tissue was fixed, paraffin-embedded and stored in a fixed tissue archive, accessed for this study. Tissue blocks were obtained and either used to generate an ovarian tissue microarray (TMA) as described previously or were sectioned at 5 µm thickness for immunohistochemistry as also described below [16].

4.2. Cytokine Antibody Arrays

Cytokine antibody arrays were performed on human follicular fluid using the Ray Bio C-series human cytokine antibody array C5 (Ray Biotech Inc, Norcross, GA, USA) according to the manufacturer’s instructions. Human follicular fluid samples from each participant were run in duplicate, and 500 μL of follicular fluid was used on each array. Although media contamination in the follicular fluid samples was anticipated to be low because only the first follicle from each ovary was used, we performed cytokine antibody arrays on media alone to determine the background and, thus, establish signal thresholds. The cytokine antibody arrays were visualized using chemiluminescence according to the manufacturer’s instructions. The cytokine antibody arrays were imaged, and the Protein Array Analyzer Plugin for ImageJ was used to quantify signal intensity using densitometry [68]. The region of interest for each cytokine measured was kept constant across all experiments. Intensity values were exported into Ray Biotech’s analysis tool, and intensity values were normalized to the oldest participant. Relative intensity values for each cytokine were plotted versus age, BMI, and AMH values.

4.3. Antibodies

To detect TGFβ3, we used a rabbit polyclonal antibody (Abcam, AB15537, Cambridge, MA, USA) that recognizes the pro-form of the TGFβ3 molecule (50 kDa) for immunohistochemistry at a 1:100 dilution (0.01 µg/µL) and immunoblot analysis at a 1:500 dilution (0.002 µg/µL). A nonimmune rabbit IgG antibody (Vector Laboratories, Inc, Burlingame, CA, USA) was used as a negative control at the same concentrations as the TGFβ3 antibody. We used a rabbit monoclonal antibody against GAPDH (36 kDa) (Catalog 5174S, Cell Signaling Technology, Danvers, MA, USA) at a 1:2000 dilution as a protein loading control to normalize TGFβ3 immunoreactivity in immunoblots. A horseradish peroxidase-conjugated donkey anti-rabbit secondary antibody was purchased from GE Healthcare Life Sciences (Pittsburgh, PA, USA) and was used to detect both primary antibodies in the immunohistochemistry and immunoblot analyses.

4.4. Protein Extraction and Immunoblot Analysis (Non-Reducing Conditions)

Cumulus cells were thawed in a 1× cell lysis buffer (Cell Signaling Technology, Danvers, MA, USA) supplemented with Halt™ Protease and Phosphatase inhibitor Cocktail (100x) (Thermo Fischer Scientific, Waltham, MA, USA). Cells were homogenized using a Pellet Pestles cordless motor (Sigma Aldrich, St. Louis, MO, USA) for 30 s intervals and spun down for 20 min at 10,000 g at 4 °C. The supernatant was divided into 5 µL aliquots, which were either used immediately for immunoblot analysis or stored at −80 °C. To run immunoblots, cumulus cell protein extract (run individually with one extract from one cumulus cell per lane) and molecular weight standards (Bio-Rad, Hercules, CA, USA) were loaded onto precast 10% MiniPROTEAN TGX gels (Bio-Rad, Hercules, CA, USA) and transferred onto 0.45 µm PVDF membranes (GE Healthcare Life Sciences). All steps were done at room temperature unless otherwise specified. Membranes were blocked with 3% ECL Blocking Agent (GE Healthcare Life Sciences) in 1× Tris-buffered saline (TBS; 20 mM Tris-HCL and 150 mM NaCl (pH 7.6)) containing 0.1% Tween-20 (TBS-T) for 3–4 h. Blots were incubated overnight in primary antibodies (TGFβ3 and GAPDH) diluted in a 3% ECL block at 4 °C. Blots were rinsed three times with 1× TBS-T for 3 × 30 min and incubated with secondary antibody diluted in 3% ECL block for 1 h. After rinsing three times again with 1× TBS-T for 30 min each, blots were developed using GE Healthcare Amersham ECL detection reagent. Blots were then reprobed with an antibody against GAPDH, which served as a loading control. To do this, blots were rinsed in 1× TBS-T and stripped in 5 mL of Restore stripping buffer (ThermoFisher, Waltham, MA, USA) for 15 min, rinsed 4 × 5 min in 1× TBS-T, and processed as described above for the TGFβ3 antibody. Films were scanned, and densitometry was performed using ImageJ software (National Institutes of Health, Bethesda, MD, USA). Individual band intensities were measured and normalized to the GAPDH band intensities. All graphs were created using GraphPad Prism 7 software.

4.5. Immunohistochemistry (IHC)

To detect TGFβ3 in human ovarian tissue, we performed IHC on whole ovarian tissue sections and on TMAs. Slides containing human FFPE ovarian tissue sections were deparaffinized in Citrisolv (Thermo Fisher Scientific, Waltham, MA, USA) and rehydrated through a graded ethanol series (100% EtOH 2 × 3 min, 95% EtOH 1 × 3 min, 85% EtOH 1 × 3 min, 70% EtOH 1 × 3 min, 50% EtOH 1 × 3 min ddH2O 1 × 3 min). Slides were rinsed in TBS-T and then incubated in 3% H2O2 for 15 min to quench endogenous peroxidases. A Pap pen (Vector Laboratories, Inc., Burlingame, CA, USA) was used to experimental and control sections on each slide. Slides were rinsed in 1× TBS, and endogenous biotin and avidin binding was blocked with an avidin/biotin blocking kit (Vector Laboratories, Inc., Burlingame, CA, USA). Nonspecific binding was attenuated by incubating slides in blocking solution according to the VECTASTAIN Elite ABC System (Vector Laboratories, Inc., Burlingame, CA, USA) with 1× TBS as the buffer. Experimental slide sections were incubated with the primary antibody, and control sections were incubated in nonimmune rabbit IgG overnight at 4 °C. Slides were rinsed again in TBS-T and incubated in biotinylated antibody according to manufacturer’s instructions (VECTASTAIN Elite ABC System, Vector Laboratories, Inc., Burlingame, CA, USA) for two hours at room temperature. Slides were washed in TBS-T and incubated in an ABC reagent (Vector Laboratories, Inc., Burlingame, CA, USA) for 30 min at room temperature. Slides were then washed again in TBS-T, and staining was visualized using the DAB peroxidase (HRP) substrate kit (Vector Laboratories, Inc., Burlingame, CA, USA), which resulted in a brown precipitate. The DAB reaction was stopped by placing slides in dH2O for 5 min as soon as the signal was first visible in the negative control section. Slides were then counterstained in hematoxylin (Harris Hematoxylin, Mercedes Medical, Sarasota, FL, USA) and then coverslipped with Cytoseal XYL (Thermo Fisher Scientific, Waltham, MA, USA). Stained slides were imaged on an EVOS FL auto microscope (Life Technologies Corporation).

4.6. Human Ovarian Tissue Microarray (TMA) TGFβ3 Analysis

We used a previously generated human ovarian TMA that included FFPE tissue blocks from 60 participants ranging from 0–20 and ≥39 [16]. The cores containing all 60 samples were arranged in duplicate on the TMA block (Supplemental Figure S1). Samples from each tissue archive were randomly distributed concerning age in the TMA block to avoid staining bias. Cores from the human liver (fibrotic and non-fibrotic) and ovarian cancer tissues were included as controls. IHC for TGFβ3 was performed on a TMA section as described above. Due to the precious nature of the human ovarian tissue, we did not want to apply a nonimmune IgG on one of the TMA slides as a negative control. Instead, as a positive control and an indicator for stopping the DAB reaction, we ran a slide containing human ovarian tissue in parallel with the TMA section. Additionally, a sequential TMA section was stained with hematoxylin and eosin (H&E) according to the standard procedure to correlate TGFβ3 expression with specific histological structures in the ovarian cores.
To quantify TGFβ3 expression in each core, whole slides were scanned in bright-field using a 20X objective on a NanoZoomer digital pathology system (Hamamatsu City, Japan) at the University of Washington Histology and Imaging Core (Seattle, Washington, DC, USA). The digital images were then imported into Visiopharm software (Hørsholm, Denmark) for analysis. Using the Visiopharm image analysis module, regions of interest (ROI) were manually drawn around each TMA core. The digital images were converted into grayscale values using two-color space bands, FastRed_DAB—DAB and RGB—B, with the mean filter of 3 × 3 pixels. The software training to determine positive TGFβ3 staining and background tissue counterstain (hematoxylin) was performed on three ROIs from two serial sections of control tissue staining in conjunction with the TMA slides. From the control slide, one section was stained with TGFβ3 and the other with IgG isotype-matched control. To be conservative in our thresholding, we averaged three different ROIs from the IgG isotype-matched control. These represented areas of low, medium and high background intensity (Supplemental Figure S1D). Individual thresholds of pixel values were created for each of the three ROIs, averaged to account for tissue heterogeneity and then the average threshold was applied to the TMA for analysis. The TMA cores were processed in batch mode using this configuration to generate the desired per area outputs and analyzed at 100%. Based on this strategy, the green label visible in the images corresponds to the thresholding for tissue that is not considered positively stained, whereas the blue label is any tissue considered positively stained. Once the analysis was complete, we received a spreadsheet indicating the ratio of positive TGFβ3 staining per total core area for each core run in duplicate. These values were used for all quantification and graphical analysis on the TMA.
To stratify the patterns of TGFβ3 localization based on the presence of specific cellular features within each tissue core, we analyzed the sequential section of the TMA that was H&E-stained. We classified all 120 cores as being enriched in either follicles, blood vessels, or stroma. If any follicle was present, the core was classified as the follicle. The vasculature category was assigned if 3 or more medium (>75 µm) to large blood vessels were present. If most of the tissue (~50% or greater) was stroma and it lacked a follicle, then the core was classified as stroma. We then stratified the ratio of TGFβ3 expression per core area by specific structure-based categories or “landmark features.”

4.7. Statistical Analysis

To examine the relationship between maternal age, AMH, BMI and cytokine relative intensities in follicular fluid, simple linear regression analysis (Pearson’s coefficient) was used. Goodness-of-fit was quantified using R-squared and sum-of-squares. Statistically significant differences between two independent groups were analyzed using unpaired t-test, and comparisons of more than two independent groups were analyzed using one-way analysis of variance (ANOVA) followed by multiple comparisons test (Tukey post hoc test). p < 0.05 was considered statistically significant. Statistical analysis was performed using GraphPad Prism Version 8.1.2 software (GraphPad Software Inc., La Jolla, CA, USA).

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijms22094902/s1, Figure S1: Human ovarian tissue microarray (TMA) layout with sample age breakdown and threshold intensity method, Figure S2: RayBio C-series human cytokine antibody array C5 decoder of cytokines in Figure 1C, Figure S3: Cytokine profiles that were significantly different with age or AMH, but did not change with BMI, Figure S4: Follicular fluid cytokine profiles that increased with chronologic age, but not with reproductive age (AMH), Figure S5: Cytokine profiles in the follicular fluid that showed no change with chronologic age, Table S1: Participant cumulus cell information.

Author Contributions

The author contributions for this study are as follows: Conceptualization, F.E.D., M.E.P., and J.H.M.; methodology, F.E.D., M.E.P., J.K. and J.H.M.; software, J.K., M.J.L., B.W.J.; validation, J.H.M.; formal analysis, J.K., J.H.M., M.J.L.; investigation, all authors.; resources, F.E.D., M.E.P., M.T.P.; data curation, J.K. and J.H.M.; writing—original draft preparation, J.H.M., S.J.B.; writing—review and editing, F.E.D. and M.T.P.; visualization, J.H.M., F.E.D.; supervision, F.E.D., M.E.P., M.T.P.; project administration, F.E.D.; funding acquisition, F.E.D., M.E.P., M.T.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from Friends of Prentice (F.E.D. and M.E.P.), the Eunice Kennedy Shriver National Institutes of Child Health and Development (R01 HD093726, F.E.D. and M.T.P.), and the Master of Science in Reproductive Science and Medicine program at Northwestern University.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board of Northwestern University (STU00205662, 1 August 2017 and STU00200166, 13 Janurary 2015).

Informed Consent Statement

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

Data Availability Statement

Data are available from the corresponding author upon reasonable request.

Acknowledgments

We thank Shawn Briley for his assistance with performing the cytokine antibody arrays, John Zhang, Saurabh Malpani, and Rafael Confino for assistance with human follicular fluid and cumulus cell collection; Sharrón Manuel and Jian Jun Wei for generating the ovarian tissue microarray, and Farners Amargant i Riera for her critical input and discussions.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Navarro-Costa, P. A (micro)environmental perspective on the evolution of female reproductive aging. J. Assist. Reprod. Genet. 2018, 35, 2129–2131. [Google Scholar] [CrossRef]
  2. Hall, J.E. Endocrinology of the Menopause. Endocrinol. Metab. Clin. N. Am. 2015, 44, 485–496. [Google Scholar] [CrossRef] [PubMed]
  3. Igarashi, H.; Takahashi, T.; Nagase, S. Oocyte aging underlies female reproductive aging: Biological mechanisms and therapeutic strategies. Reprod. Med. Biol. 2015, 14, 159–169. [Google Scholar] [CrossRef] [PubMed]
  4. Cimadomo, D.; Fabozzi, G.; Vaiarelli, A.; Ubaldi, N.; Ubaldi, F.M.; Rienzi, L. Impact of Maternal Age on Oocyte and Embryo Competence. Front. Endocrinol. 2018, 9, 327. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Ameratunga, D.; Weston, G.; Osianlis, T.; Catt, J.; Vollenhoven, B. In vitro fertilisation (IVF) with donor eggs in post-menopausal women: Are there differences in pregnancy outcomes in women with premature ovarian failure (POF) compared with women with physiological age-related menopause? J. Assist. Reprod. Genet. 2009, 26, 511–514. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Liang, X.; Ma, J.; Schatten, H.; Sun, Q. Epigenetic changes associated with oocyte aging. Sci. China Life Sci. 2012, 55, 670–676. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Eichenlaub-Ritter, U.; Wieczorek, M.; Lüke, S.; Seidel, T. Age related changes in mitochondrial function and new approaches to study redox regulation in mammalian oocytes in response to age or maturation conditions. Mitochondrion 2011, 11, 783–796. [Google Scholar] [CrossRef]
  8. Hornick, J.E.; Duncan, F.E.; Sun, M.; Kawamura, R.; Marko, J.F.; Woodruff, T.K. Age-associated alterations in the micromechanical properties of chromosomes in the mammalian egg. J. Assist. Reprod. Genet. 2015, 32, 765–769. [Google Scholar] [CrossRef]
  9. Duncan, F.E.; Jasti, S.; Paulson, A.; Kelsh, J.M.; Fegley, B.; Gerton, J.L. Age-associated dysregulation of protein metabolism in the mammalian oocyte. Aging Cell 2017, 16, 1381–1393. [Google Scholar] [CrossRef] [PubMed]
  10. López-Otín, C.; Blasco, M.A.; Partridge, L.; Serrano, M.; Kroemer, G. The Hallmarks of Aging. Cell 2013, 153, 1194–1217. [Google Scholar] [CrossRef] [Green Version]
  11. Thannickal, V.J.; Murthy, M.; Balch, W.E.; Chandel, N.S.; Meiners, S.; Eickelberg, O.; Selman, M.; Pardo, A.; White, E.S.; Levy, B.D.; et al. Blue Journal Conference. Aging and Susceptibility to Lung Disease. Am. J. Respir. Crit. Care Med. 2015, 191, 261–269. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Weiss, U. Inflammation. Nature 2008, 454, 427. [Google Scholar] [CrossRef] [Green Version]
  13. Nathan, C.; Ding, A. Nonresolving Inflammation. Cell 2010, 140, 871–882. [Google Scholar] [CrossRef] [Green Version]
  14. Sahin, H.; Wasmuth, H.E. Chemokines in tissue fibrosis. Biochim. Biophys. Acta (BBA)-Mol. Basis Dis. 2013, 1832, 1041–1048. [Google Scholar] [CrossRef] [Green Version]
  15. Briley, S.M.; Jasti, S.; McCracken, J.M.; Hornick, J.E.; Fegley, B.; Pritchard, M.T.; Duncan, F.E. Reproductive age-associated fibrosis in the stroma of the mammalian ovary. Reproduction 2016, 152, 245–260. [Google Scholar] [CrossRef] [Green Version]
  16. Amargant, F.; Manuel, S.L.; Tu, Q.; Parkes, W.S.; Rivas, F.; Zhou, L.T.; Rowley, J.E.; Villanueva, C.E.; Hornick, J.E.; Shekhawat, G.S.; et al. Ovarian stiffness increases with age in the mammalian ovary and depends on collagen and hyaluronan matrices. Aging Cell 2020, 19, e13259. [Google Scholar] [CrossRef] [PubMed]
  17. Foley, K.G.; Pritchard, M.T.; Duncan, F.E. Macrophage-derived multinucleated giant cells: Hallmarks of the aging ovary. Reproduction 2021, 161, V5–V9. [Google Scholar] [CrossRef]
  18. Huang, Y.; Hu, C.; Ye, H.; Luo, R.; Fu, X.; Li, X.; Huang, J.; Chen, W.; Zheng, Y. Inflamm-Aging: A New Mechanism Affecting Premature Ovarian Insufficiency. J. Immunol. Res. 2019, 2019, 8069898. [Google Scholar] [CrossRef] [Green Version]
  19. Umehara, T.; Richards, J.S.; Shimada, M. The stromal fibrosis in aging ovary. Aging 2018, 10, 9–10. [Google Scholar] [CrossRef] [Green Version]
  20. Edwards, R. Follicular Fluid. Reproduction 1974, 37, 189–219. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Rodgers, R.J.; Irving-Rodgers, H.F. Formation of the Ovarian Follicular Antrum and Follicular Fluid1. Biol. Reprod. 2010, 82, 1021–1029. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Revelli, A.; Piane, L.D.; Casano, S.; Molinari, E.; Massobrio, M.; Rinaudo, P. Follicular fluid content and oocyte quality: From single biochemical markers to metabolomics. Reprod. Biol. Endocrinol. 2009, 7, 40. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Emori, M.M.; Drapkin, R. The hormonal composition of follicular fluid and its implications for ovarian cancer pathogenesis. Reprod. Biol. Endocrinol. 2014, 12, 60. [Google Scholar] [CrossRef] [Green Version]
  24. Bouet, P.-E.; Boueilh, T.; de la Barca, J.M.C.; Boucret, L.; Blanchard, S.; Ferré-L’Hotellier, V.; Jeannin, P.; Descamps, P.; Procaccio, V.; Reynier, P.; et al. The cytokine profile of follicular fluid changes during ovarian ageing. J. Gynecol. Obstet. Hum. Reprod. 2020, 49, 101704. [Google Scholar] [CrossRef] [Green Version]
  25. de Vet, A.; Laven, J.S.; de Jong, F.H.; Themmen, A.P.; Fauser, B.C. Antimüllerian hormone serum levels: A putative marker for ovarian aging. Fertil. Steril. 2002, 77, 357–362. [Google Scholar] [CrossRef]
  26. Bentzen, J.G.; Forman, J.L.; Johannsen, T.H.; Pinborg, A.; Larsen, E.C.; Andersen, A.N. Ovarian Antral Follicle Subclasses and Anti-Müllerian Hormone During Normal Reproductive Aging. J. Clin. Endocrinol. Metab. 2013, 98, 1602–1611. [Google Scholar] [CrossRef] [PubMed]
  27. Van Rooij, I.; Broekmans, F.; Velde, E.T.; Fauser, B.; Bancsi, L.; Jong, F.; Themmen, A. Serum anti-Mullerian hormone levels: A novel measure of ovarian reserve. Hum. Reprod. 2002, 17, 3065–3071. [Google Scholar] [CrossRef] [Green Version]
  28. Welt, C.K.; Schneyer, A.L.; Heist, K.; Mantzoros, C.S. Leptin and Soluble Leptin Receptor in Follicular Fluid. J. Assist. Reprod. Genet. 2003, 20, 495–501. [Google Scholar] [CrossRef] [Green Version]
  29. Llaneza-Suarez, D.; Llaneza, P.; González, C.; De-La-Fuente, P.; García-Ochoa, C.; Garrido, P.; Castañón, V.; Pérez-López, F.R. Assessment of follicular fluid leptin levels and insulin resistance as outcome predictors in women undergoing in vitro fertilization–intracytoplasmic sperm injection. Fertil. Steril. 2014, 102, 1619–1625. [Google Scholar] [CrossRef]
  30. Letterio, J.J.; Roberts, A.B. Regulation of Immune Responses by TGF-β. Annu. Rev. Immunol. 1998, 16, 137–161. [Google Scholar] [CrossRef] [Green Version]
  31. Border, W.A.; Ruoslahti, E. Transforming growth factor-beta in disease: The dark side of tissue repair. J. Clin. Investig. 1992, 90, 1–7. [Google Scholar] [CrossRef] [Green Version]
  32. Raja-Khan, N.; Urbanek, M.; Rodgers, R.J.; Legro, R.S. The role of TGF-beta in polycystic ovary syndrome. Reprod. Sci. 2014, 21, 20–31. [Google Scholar] [CrossRef] [Green Version]
  33. Dewidar, B.; Meyer, C.; Dooley, S.; Meindl-Beinker, A.N. TGF-β in Hepatic Stellate Cell Activation and Liver Fibrogenesis—Updated 2019. Cells 2019, 8, 1419. [Google Scholar] [CrossRef] [Green Version]
  34. Monsivais, D.; Matzuk, M.M.; Pangas, S.A. The TGF-beta Family in the Reproductive Tract. Cold Spring Harb. Perspect Biol. 2017, 9, 10. [Google Scholar] [CrossRef] [Green Version]
  35. Hachim, M.Y.; Hachim, I.Y.; Dai, M.; Ali, S.; Lebrun, J.J. Differential expression of TGFbeta isoforms in breast cancer highlights different roles during breast cancer progression. Tumour Biol. 2018, 40, 1010428317748254. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Janny, L.; Menezo, Y.J.R. Maternal age effect on early human embryonic development and blastocyst formation. Mol. Reprod. Dev. 1996, 45, 31–37. [Google Scholar] [CrossRef]
  37. Steuerwald, N.M.; Bermúdez, M.G.; Wells, D.; Munné, S.; Cohen, J. Maternal age-related differential global expression profiles observed in human oocytes. Reprod. Biomed. Online 2007, 14, 700–708. [Google Scholar] [CrossRef]
  38. Munné, S.; Sandalinas, M.; Escudero, T.; Márquez, C.; Cohen, J. Chromosome mosaicism in cleavage-stage human embryos: Evidence of a maternal age effect. Reprod. Biomed. Online 2002, 4, 223–232. [Google Scholar] [CrossRef]
  39. Hou, Y.-J.; Zhu, C.-C.; Duan, X.; Liu, H.-L.; Wang, Q.; Sun, S.-C. Both diet and gene mutation induced obesity affect oocyte quality in mice. Sci. Rep. 2016, 6, 18858. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Sneed, M.L.; Uhler, M.L.; Grotjan, H.E.; Rapisarda, J.J.; Lederer, K.J.; Beltsos, A.N. Body mass index: Impact on IVF success appears age-related. Hum. Reprod. 2008, 23, 1835–1839. [Google Scholar] [CrossRef] [Green Version]
  41. Roubenoff, R.; Harris, T.B.; Abad, L.W.; Wilson, P.W.F.; Dallal, G.E.; Dinarello, C.A. Monocyte cytokine production in an elderly population: Effect of age and inflammation. J. Gerontol. Ser. A Boil. Sci. Med Sci. 1998, 53, M20–M26. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Kiecolt-Glaser, J.K.; Preacher, K.J.; Maccallum, R.C.; Atkinson, C.; Malarkey, W.B.; Glaser, R. Chronic stress and age-related increases in the proinflammatory cytokine IL-6. Proc. Natl. Acad. Sci. USA 2003, 100, 9090–9095. [Google Scholar] [CrossRef] [Green Version]
  43. Khaodhiar, L.; Ling, P.-R.; Blackburn, G.L.; Bistrian, B.R. Serum Levels of Interleukin-6 and C-Reactive Protein Correlate With Body Mass Index Across the Broad Range of Obesity. J. Parenter. Enter. Nutr. 2004, 28, 410–415. [Google Scholar] [CrossRef]
  44. Cioffi, J.A.; Van Blerkom, J.; Antczak, M.; Shafer, A.; Wittmer, S.; Snodgrass, H.R. The expression of leptin and its receptors in pre-ovulatory human follicles. Mol. Hum. Reprod. 1997, 3, 467–472. [Google Scholar] [CrossRef]
  45. Duggal, P.S.; Van der Hoek, K.H.; Milner, C.R.; Ryan, N.K.; Armstrong, D.T.; Magoffin, D.A.; Norman, R.J. The in vivo and in vitro effects of exogenous leptin on ovulation in the rat. Endocrinology 2000, 141, 1971–1976. [Google Scholar] [CrossRef]
  46. Swain, J.E.; Dunn, R.L.; McConnell, D.; Gonzalez-Martinez, J.; Smith, G.D. Direct Effects of Leptin on Mouse Reproductive Function: Regulation of Follicular, Oocyte, and Embryo Development. Biol. Reprod. 2004, 71, 1446–1452. [Google Scholar] [CrossRef] [Green Version]
  47. Sharabiani, M.T.A.; Vermeulen, R.; Scoccianti, C.; Hosnijeh, F.S.; Minelli, L.; Sacerdote, C.; Palli, D.; Krogh, V.; Tumino, R.; Chiodini, P.; et al. Immunologic profile of excessive body weight. Biomarkers 2011, 16, 243–251. [Google Scholar] [CrossRef]
  48. Zhang, Z.; Schlamp, F.; Huang, L.; Clark, H.; Brayboy, L.M. Inflammaging is associated with shifted macrophage ontogeny and polarization in the aging mouse ovary. Reproduction 2020, 159, 325–337. [Google Scholar] [CrossRef] [PubMed]
  49. Franceschi, C.; Campisi, J. Chronic inflammation (inflammaging) and its potential contribution to age-associated diseases. J. Gerontol. Ser. A Biol. Sci. Med. Sci. 2014, 69 (Suppl. 1), S4–S9. [Google Scholar] [CrossRef] [PubMed]
  50. Xia, S.; Zhang, X.; Zheng, S.; Khanabdali, R.; Kalionis, B.; Wu, J.; Wan, W.; Tai, X. An Update on Inflamm-Aging: Mechanisms, Prevention, and Treatment. J. Immunol. Res. 2016, 2016, 1–12. [Google Scholar] [CrossRef] [PubMed]
  51. Sanada, F.; Taniyama, Y.; Muratsu, J.; Otsu, R.; Shimizu, H.; Rakugi, H.; Morishita, R. Source of Chronic Inflammation in Aging. Front. Cardiovasc. Med. 2018, 5, 8426874. [Google Scholar] [CrossRef] [Green Version]
  52. van Deursen, J.M. The role of senescent cells in ageing. Nature 2014, 509, 439–446. [Google Scholar] [CrossRef] [Green Version]
  53. Chung, H.Y.; Ki, W.C.; Lee, E.K.; Chung, K.W.; Chung, S.; Lee, B.; Seo, A.Y.; Chung, J.H.; Jung, Y.S.; Im, E.; et al. Redefining Chronic Inflammation in Aging and Age-Related Diseases: Proposal of the Senoinflammation Concept. Aging Dis. 2019, 10, 367–382. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Coppé, J.-P.; Desprez, P.-Y.; Krtolica, A.; Campisi, J. The Senescence-Associated Secretory Phenotype: The Dark Side of Tumor Suppression. Annu. Rev. Pathol. Mech. Dis. 2010, 5, 99–118. [Google Scholar] [CrossRef] [Green Version]
  55. Schluns, K.S.; Kieper, W.C.; Jameson, S.C.; Lefrançois, L. Interleukin-7 mediates the homeostasis of naïve and memory CD8 T cells in vivo. Nat. Immunol. 2000, 1, 426–432. [Google Scholar] [CrossRef]
  56. Seddon, B.; Tomlinson, P.; Zamoyska, R. Interleukin 7 and T cell receptor signals regulate homeostasis of CD4 memory cells. Nat. Immunol. 2003, 4, 680–686. [Google Scholar] [CrossRef]
  57. Vujisi, S.; Lepej, S.I.; Emedi, I.; Bauman, R.; Remenar, A.; Tiljak, M.K.; Vujisić, S. Ovarian follicular concentration of IL-12, IL-15, IL-18 and p40 subunit of IL-12 and IL-23. Hum. Reprod. 2006, 21, 2650–2655. [Google Scholar] [CrossRef] [Green Version]
  58. Morrisette-Thomas, V.; Cohen, A.A.; Fülöp, T.; Riesco, É.; Legault, V.; Li, Q.; Milot, E.; Dusseault-Bélanger, F.; Ferrucci, L. Inflamm-aging does not simply reflect increases in pro-inflammatory markers. Mech. Ageing Dev. 2014, 139, 49–57. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Freund, A.; Orjalo, A.V.; Desprez, P.-Y.; Campisi, J. Inflammatory networks during cellular senescence: Causes and consequences. Trends Mol. Med. 2010, 16, 238–246. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Faget, D.V.; Ren, Q.; Stewart, S.A. Unmasking senescence: Context-dependent effects of SASP in cancer. Nat. Rev. Cancer 2019, 19, 439–453. [Google Scholar] [CrossRef] [PubMed]
  61. Ingman, W.V.; Robker, R.L.; Woittiez, K.; Robertson, S.A. Null mutation in transforming growth factor beta1 disrupts ovarian function and causes oocyte incompetence and early embryo arrest. Endocrinology 2006, 147, 835–845. [Google Scholar] [CrossRef] [Green Version]
  62. Halder, S.; Al-Hendy, A. Hypovitaminosis D and high serum transforming growth factor beta-3: Important biomarkers for uterine fibroids risk. Fertil. Steril. 2016, 106, 1648–1649. [Google Scholar] [CrossRef] [Green Version]
  63. Klar, R.M.; Duarte, R.; Dix-Peek, T.; Ripamonti, U. The induction of bone formation by the recombinant human transforming growth factor-β3. Biomaterials 2014, 35, 2773–2788. [Google Scholar] [CrossRef]
  64. Ohno, S.; Hirano, S.; Kanemaru, S.-I.; Tateya, I.; Kitani, Y.; Kojima, T.; Nakamura, T.; Ito, J. Prevention of buccal mucosa scarring with transforming growth factor β3. Laryngoscope 2011, 121, 1404–1409. [Google Scholar] [CrossRef] [PubMed]
  65. Nilsson, E.E.; Doraiswamy, V.; Skinner, M.K. Transforming growth factor-beta isoform expression during bovine ovarian antral follicle development. Mol. Reprod. Dev. 2003, 66, 237–246. [Google Scholar] [CrossRef] [PubMed]
  66. Guerrero, P.A.; McCarty, J.H. TGF-β Activation and Signaling in Angiogenesis. In Physiologic and Pathologic Angiogenesis; Simionescu, D., Simionescu, A., Eds.; IntechOpen: Rijeka, Croatia, 2017. [Google Scholar]
  67. Friedman, C.I.; Danforth, D.R.; Herbosa-Encarnacion, C.; Arbogast, L.; Alak, B.M.; Seifer, D.B. Follicular fluid vascular endothelial growth factor concentrations are elevated in women of advanced reproductive age undergoing ovulation induction. Fertil. Steril. 1997, 68, 607–612. [Google Scholar] [CrossRef]
  68. Schneider, C.A.; Rasband, W.S.; Eliceiri, K.W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 2012, 9, 671–675. [Google Scholar] [CrossRef] [PubMed]
Ijms 22 04902 g001 550
Figure 1. Relationship between age and anti-Müllerian hormone (AMH) and body mass index (BMI) and sample cytokine antibody array (A) Graph demonstrating that with increasing age, there was a statistically significant inverse correlation with AMH in our sample population (total n = 30). (B) Graph demonstrating that with increasing age, there was no statistically significant correlation with BMI in our population. (C) Representative human cytokine array C5 performed on the follicular fluid from the first follicle aspirated from either the right (participant 04) or left ovary (participant 04′) from a single participant. One array was incubated in media alone to determine the background and thus, establish signal thresholds. The lines in (A) and (B) correspond to the fitted regression equation. * p < 0.05, N.S. = no statistical significance.
Figure 1. Relationship between age and anti-Müllerian hormone (AMH) and body mass index (BMI) and sample cytokine antibody array (A) Graph demonstrating that with increasing age, there was a statistically significant inverse correlation with AMH in our sample population (total n = 30). (B) Graph demonstrating that with increasing age, there was no statistically significant correlation with BMI in our population. (C) Representative human cytokine array C5 performed on the follicular fluid from the first follicle aspirated from either the right (participant 04) or left ovary (participant 04′) from a single participant. One array was incubated in media alone to determine the background and thus, establish signal thresholds. The lines in (A) and (B) correspond to the fitted regression equation. * p < 0.05, N.S. = no statistical significance.
Ijms 22 04902 g001
Ijms 22 04902 g002 550
Figure 2. A unique subset of fibroinflammatory cytokines tightly associated with reproductive and chronological aging. (A) Graphs showing relative cytokine intensities in the follicular fluid of the first follicle in women undergoing assisted reproductive technology (ART) vs. age. Fibroinflammatory markers IL-3, IL-7, IL-15, TGF-β1, TGF-β3 and MIP-1β were significantly correlated with increasing chronological age (years). (B) Graphs showing cytokine markers IL-3, IL-7, IL-15, TGF-β1, TGF-β3 and MIP-1β were also significantly negatively correlated with anti-Müllerian hormone (ng/mL) levels, an indicator for decreased ovarian reserve. The lines correspond to the fitted regression equation. * p < 0.05, ** p < 0.01.
Figure 2. A unique subset of fibroinflammatory cytokines tightly associated with reproductive and chronological aging. (A) Graphs showing relative cytokine intensities in the follicular fluid of the first follicle in women undergoing assisted reproductive technology (ART) vs. age. Fibroinflammatory markers IL-3, IL-7, IL-15, TGF-β1, TGF-β3 and MIP-1β were significantly correlated with increasing chronological age (years). (B) Graphs showing cytokine markers IL-3, IL-7, IL-15, TGF-β1, TGF-β3 and MIP-1β were also significantly negatively correlated with anti-Müllerian hormone (ng/mL) levels, an indicator for decreased ovarian reserve. The lines correspond to the fitted regression equation. * p < 0.05, ** p < 0.01.
Ijms 22 04902 g002
Ijms 22 04902 g003 550
Figure 3. Cytokines that increase with BMI are unique from those that change with age. (A) Graphs showing that cytokines leptin and IL-8 were significantly correlated with increasing BMI (kg/m2). (B) Graphs showing that leptin and IL-8 did not significantly change with increasing age. (C) Graphs showing that leptin and IL-8 did not significantly change with increasing AMH. The lines correspond to the fitted regression equation. * p < 0.05, ** p < 0.01, N.S. = no statistical significance.
Figure 3. Cytokines that increase with BMI are unique from those that change with age. (A) Graphs showing that cytokines leptin and IL-8 were significantly correlated with increasing BMI (kg/m2). (B) Graphs showing that leptin and IL-8 did not significantly change with increasing age. (C) Graphs showing that leptin and IL-8 did not significantly change with increasing AMH. The lines correspond to the fitted regression equation. * p < 0.05, ** p < 0.01, N.S. = no statistical significance.
Ijms 22 04902 g003
Ijms 22 04902 g004 550
Figure 4. Specific detection of TGFβ3 in human cumulus cells shows no significant difference with age. (A) Three representative images of the varying sizes of cumulus cell masses from smaller to larger are shown from participant 42 (39 y). All cumulus cell images were taken upon receipt to the lab. (B) We observed a specific immunoreactive band for TGFβ3 at 50 kDa representing the pro-form of the TGFβ3 molecule. MCF7 is the positive control. Two intact cumulus cell masses from participant 31 are shown in lanes 2 and 3. The arrow indicates TGFβ3 (50 kDa). (C) Representative blot probed with antibodies recognizing TGFβ3 and reprobed for GAPDH (36 kDa). (D) Ratio of TGFβ3 to GAPDH in young and old cumulus cells for all replicates. Young and old groups N.S. (E) Fold-change of TGFβ3 in old cumulus cells compared to young cumulus cells for all replicates. More detailed information about each participant can be found in Supplementary Table S1. Scale bar = 100 µm.
Figure 4. Specific detection of TGFβ3 in human cumulus cells shows no significant difference with age. (A) Three representative images of the varying sizes of cumulus cell masses from smaller to larger are shown from participant 42 (39 y). All cumulus cell images were taken upon receipt to the lab. (B) We observed a specific immunoreactive band for TGFβ3 at 50 kDa representing the pro-form of the TGFβ3 molecule. MCF7 is the positive control. Two intact cumulus cell masses from participant 31 are shown in lanes 2 and 3. The arrow indicates TGFβ3 (50 kDa). (C) Representative blot probed with antibodies recognizing TGFβ3 and reprobed for GAPDH (36 kDa). (D) Ratio of TGFβ3 to GAPDH in young and old cumulus cells for all replicates. Young and old groups N.S. (E) Fold-change of TGFβ3 in old cumulus cells compared to young cumulus cells for all replicates. More detailed information about each participant can be found in Supplementary Table S1. Scale bar = 100 µm.
Ijms 22 04902 g004
Ijms 22 04902 g005 550
Figure 5. TGFβ3 is found throughout the ovarian stroma, within the follicles and around blood vessels. (A) Tissue section with specific staining throughout the stroma and around key structures in the ovary. (B) IgG negative control; (C) primordial follicles (asterisks); (D) blood vessels (arrowheads); (E) dense stroma. Panels (A,B) are representative histological sections performed in n = 6 patients, and (C–E) are cores from the TMA with representative images shown. Arrowheads: blood vessels, ta: tunica albuginea, c: cortex, m: medulla, asterisk: primordial follicles. (A,B) Scale bar = 400 µm; (C,D) scale bar = 100 µm.
Figure 5. TGFβ3 is found throughout the ovarian stroma, within the follicles and around blood vessels. (A) Tissue section with specific staining throughout the stroma and around key structures in the ovary. (B) IgG negative control; (C) primordial follicles (asterisks); (D) blood vessels (arrowheads); (E) dense stroma. Panels (A,B) are representative histological sections performed in n = 6 patients, and (C–E) are cores from the TMA with representative images shown. Arrowheads: blood vessels, ta: tunica albuginea, c: cortex, m: medulla, asterisk: primordial follicles. (A,B) Scale bar = 400 µm; (C,D) scale bar = 100 µm.
Ijms 22 04902 g005
Ijms 22 04902 g006 550
Figure 6. Quantitative analysis of human ovarian tissue microarray (TMA) reveals increased TGFβ3 content with age. (A) Average ratio of TGFβ3 expression per total area for core 1 and core 2 (unique cores from the same patient) broken down into four age cohorts: 0–10, 11–20, 39–50, 51–60. Images from cores closest to the average data point in each age cohort are shown with TGFβ3 staining (brown) in the panel above and quantitative data analysis in green and blue below. The dark blue area represents a positive TGFβ3 signal above the threshold, while green is TGFβ3 negative. (B) Ratio of TGFβ3 per total area in the 0–20 vs. 39–60 age cohorts. (C) Ratio of TGFβ3 per total area in the 0–10 vs. 51–60 age cohorts. Asterisks (*) indicate statistical significance with panel (B) p = 0.01 and panel (C) p = 0.005.
Figure 6. Quantitative analysis of human ovarian tissue microarray (TMA) reveals increased TGFβ3 content with age. (A) Average ratio of TGFβ3 expression per total area for core 1 and core 2 (unique cores from the same patient) broken down into four age cohorts: 0–10, 11–20, 39–50, 51–60. Images from cores closest to the average data point in each age cohort are shown with TGFβ3 staining (brown) in the panel above and quantitative data analysis in green and blue below. The dark blue area represents a positive TGFβ3 signal above the threshold, while green is TGFβ3 negative. (B) Ratio of TGFβ3 per total area in the 0–20 vs. 39–60 age cohorts. (C) Ratio of TGFβ3 per total area in the 0–10 vs. 51–60 age cohorts. Asterisks (*) indicate statistical significance with panel (B) p = 0.01 and panel (C) p = 0.005.
Ijms 22 04902 g006
Ijms 22 04902 g007 550
Figure 7. Analysis of ovarian cores from the TMA reveal changes in TGFβ3 content with age depending on the presence of landmark tissue structures. Each core was categorized into one of three major tissue structures: follicles, blood vessels or stroma. Representative images of youngest and oldest age cohort’s cores are shown with the first panel: TGFβ3 staining in the brown, middle panel: hematoxylin and eosin stain and last panel: quantitative data measurement with positive TGFβ3 signal in dark blue and negative TGFβ3 stain in green. (A) Average core 1 and 2 ratios of TGFβ3 per total area across age cohorts with follicles. (B) Average core 1 and 2 ratios of TGFβ3 per total area across age cohorts with blood vessels. (C) Average core 1 and 2 ratios of TGFβ3 per total area across age cohorts with stroma. Asterisks (*) indicate statistical significance with panel (B) p = 0.04 and panel (C) p = 0.004.
Figure 7. Analysis of ovarian cores from the TMA reveal changes in TGFβ3 content with age depending on the presence of landmark tissue structures. Each core was categorized into one of three major tissue structures: follicles, blood vessels or stroma. Representative images of youngest and oldest age cohort’s cores are shown with the first panel: TGFβ3 staining in the brown, middle panel: hematoxylin and eosin stain and last panel: quantitative data measurement with positive TGFβ3 signal in dark blue and negative TGFβ3 stain in green. (A) Average core 1 and 2 ratios of TGFβ3 per total area across age cohorts with follicles. (B) Average core 1 and 2 ratios of TGFβ3 per total area across age cohorts with blood vessels. (C) Average core 1 and 2 ratios of TGFβ3 per total area across age cohorts with stroma. Asterisks (*) indicate statistical significance with panel (B) p = 0.04 and panel (C) p = 0.004.
Ijms 22 04902 g007
Table
Table 1. Participant-specific health information.
Table 1. Participant-specific health information.
Participant IDAge (Years)AMH (ng/mL)BMI (kg/m2)
Participant 0127.71.826.3
Participant 0227.73.225.8
Participant 0328.82.119.9
Participant 0430.55.920.5
Participant 0530.50.920.2
Participant 0631.65.219.8
Participant 0731.92.724.1
Participant 0832.44.622.7
Participant 0932.62.024.0
Participant 1032.76.222.9
Participant 1132.91.221.2
Participant 1232.91.929.6
Participant 1333.63.621.1
Participant 1433.74.317.5
Participant 1534.03.925.7
Participant 1635.54.122.9
Participant 1735.60.223.0
Participant 1836.03.729.2
Participant 1936.81.423.6
Participant 2039.00.425.5
Participant 2139.04.328.3
Participant 2239.11.119.7
Participant 2339.95.321.1
Participant 2440.02.124.0
Participant 2540.21.028.4
Participant 2642.13.521.3
Participant 2742.40.621.3
Participant 2843.30.121.0
Participant 2944.20.326.1
Participant 3044.80.822.3
Average (± SE)35.7 ± 0.902.6 ± 0.3323.3 ± 0.57
Table
Table 2. Relationship between fibroinflammatory cytokines and biologic and reproductive aging.
Table 2. Relationship between fibroinflammatory cytokines and biologic and reproductive aging.
ParametersFibroinflammatory Cytokines
Age and AMHIL-3, IL-7, IL-15, MIP-1β, TGF-β1, TGF-β3
AgeM-CSF, SDF-1, EGF, oncostatin M, VEGF, PDGF-BB, BLC, CK β 8-1, BDNF, FGF-9, GDNF, IGFBP-3, LIF, MIF, PIGF, TGF-β2
Non-significantGRO, GRO-α, IL-β, IL-8, IL-12 p40, MCP-1, MCP-2, MCP-3, MIP-1δ, RANTES, SCF, IGF-1, angiogenin, leptin, eotaxin 1, eotaxin 2, eotaxin 3, FGF-4, FGF-6, FGF-7, Flt-3 ligand, fractalkine, GCP-2, HGF, IGFBP-1, IGFBP-2, IGFBP-4, IL-16, IP-10, LIGHT, MIP-3α, NAP-2, NT-3, NT-4, osteopontin, osteoprotegerin, PARC, TIMP-1, TIMP-2
Table
Table 3. Young vs. old participant cumulus cell information.
Table 3. Young vs. old participant cumulus cell information.
GroupParticipant IDAgeCumulus Mass IDEmbryo Outcome
1Participant 3232.5aBlastocyst day 6
bBlastocyst day 6
Participant 3344.7aBlastocyst day 6
bBlastocyst day 5
2Participant 3433.9aHatching blastocyst day 5
bHatching blastocyst day 5
cHatching blastocyst day 6
Participant 3541.5aBlastocyst day 5
bBlastocyst day 6
cBlastocyst day 6
3Participant 3634.4aBlastocyst day 6
bBlastocyst day 6
Participant 3742.3a6 cell day 6
bMature (not fert.)
4Participant 3834.8aBlastocyst day 5
Participant 3940.1aBlastocyst day 5
5Participant 4034.9bEarly blastocyst day 6
cMorula day 6
Participant 4143a10 cell day 6
b6 cell day 6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Machlin, J.H.; Barishansky, S.J.; Kelsh, J.; Larmore, M.J.; Johnson, B.W.; Pritchard, M.T.; Pavone, M.E.; Duncan, F.E. Fibroinflammatory Signatures Increase with Age in the Human Ovary and Follicular Fluid. Int. J. Mol. Sci. 2021, 22, 4902. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22094902

AMA Style

Machlin JH, Barishansky SJ, Kelsh J, Larmore MJ, Johnson BW, Pritchard MT, Pavone ME, Duncan FE. Fibroinflammatory Signatures Increase with Age in the Human Ovary and Follicular Fluid. International Journal of Molecular Sciences. 2021; 22(9):4902. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22094902

Chicago/Turabian Style

Machlin, Jordan H., Seth J. Barishansky, John Kelsh, Megan J. Larmore, Brian W. Johnson, Michele T. Pritchard, Mary Ellen Pavone, and Francesca E. Duncan. 2021. "Fibroinflammatory Signatures Increase with Age in the Human Ovary and Follicular Fluid" International Journal of Molecular Sciences 22, no. 9: 4902. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22094902

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

Article Metrics

Back to TopTop