# Beta-nerve growth factor stimulates spontaneous electrical activity of in vitro embryonic mouse GnRH neurons through a P75 mediated-mechanism

Caroline Pinet-Charvet, Renaud Fleurot, Flavie Derouin Tochon, Simon de Graaf, Xavier Druart, Guillaume Tsikis, Catherine Taragnat, Ana-Paula Teixeira-Gomes, Valérie Labas, Thierry Moreau, et al.

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Caroline Pinet-Charvet, Renaud Fleurot, Flavie Derouin Tochon, Simon de Graaf, Xavier Druart, et al. Beta-nerve growth factor stimulates spontaneous electrical activity of in vitro embryonic mouse GnRH neurons through a P75 mediated-mechanism. Scientific Reports, 2020, 10 (1), pp.10654. 10.1038/s41598-020-67665-4. hal-03094626

HAL Id: hal-03094626

https://hal.science/hal-03094626v1

Submitted on 2 Sep 2024

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

OPEN

# Beta-nerve growth factor stimulates spontaneous electrical activity of in vitro embryonic mouse GnRH neurons through a P75 mediated-mechanism

Caroline Pinet-Charvet<fnref n="1" />,²,⁶, Renaud Fleurot<fnref n="1" />,⁶, Flavie Derouin-Tochon<fnref n="1" />, Simon de Graaf³, Xavier Druart<fnref n="1" />, Guillaume Tsikis<fnref n="1" />, Catherine Taragnat<fnref n="1" />, Ana-Paula Teixeira-Gomes⁴, Valérie Labas<fnref n="1" />, Thierry Moreau⁵, Xavier Cayla<fnref n="1" /> & Anne H Duittoz<fnref n="1" />✉

The control of ovulation helps guarantee the success of reproduction and as such, contributes to the fitness of a species. In mammals, two types of ovulation are observed: induced and spontaneous ovulation. Recent work on camelids, that are induced ovulators, highlighted the role of a factor present in seminal plasma, beta Nerve Growth Factor (β-NGF), as the factor that triggers ovulation in a GnRH dependent manner. In the present work, we characterized alpaca β-NGF (aβ-NGF) and its 3D structure and compared it with human recombinant β-NGF (hβ-NGF). We showed that the β-NGF enriched fraction of alpaca semen and the human recombinant protein, both stimulated spontaneous electrical activity of primary GnRH neurons derived from mouse embryonic olfactory placodes. This effect was dose-dependent and mediated by p75 receptor signaling. P75 receptors were found expressed in vitro by olfactory ensheathing cells (OEC) in close association with GnRH neurons and in vivo by tanycytes in close vicinity to GnRH fibers in adult mouse. Altogether, these results suggested that β-NGF induced ovulation through an increase in GnRH secretion provoked by a glial dependent P75 mediated mechanism.

Two modes of ovulation have been described in mammals: spontaneous ovulation and induced ovulation. Spontaneous ovulation occurs in many mammalian species (mouse, rat, bovine, ovine, equine, primates). In these species, ovulation occurs spontaneously whether the female copulates or not. Induced ovulation has been described in a great variety of species belonging to various mammalian Orders. Induced ovulation has been particularly well studied in camelids<fnref n="1" />. In these species, ovulation occurs after mating in more than 95% of females. These observations drove to the hypothesis of a physical stimulation of the female genitalia and the existence of an ovulation reflex mediated through the spine medulla. Spontaneous ovulators have been intensively studied since most laboratory species belong to this category. In these species, it is acknowledged that ovulation is triggered by an increase in luteinizing hormone (LH) levels (preovulatory surge) as the consequence of an increase in gonadotropin releasing hormone (GnRH) release. The timing of ovulation is controlled by kisspeptin neurons

<footnote n="1">¹Physiologie de la Reproduction et des Comportements (PRC) UMR7247 INRA, CNRS, Centre INRA Val de Loire, Université de Tours, IFCE, 37380 Nouzilly, France. ²Physiologie de la Reproduction et des Comportements (PRC), ComUE Centre-Val de Loire, Centre INRA Val de Loire, Université de Poitiers, 37380 Nouzilly, France. ³Faculty of Science, School of Life and Environmental Sciences, The University of Sydney, Sydney, NSW 2006, Australia. ⁴Infectiologie et Santé Publique (ISP) UMR1282, INRA, Centre INRA Val de Loire, Université de Tours, 37380 Nouzilly, France. ⁵Biologie des Oiseaux et Aviculture (BOA) UMR Centre INRA Val de Loire, 37380 Nouzilly, France. ⁶These authors contributed equally: Caroline Pinet-Charvet and Renaud Fleurot. ✉email: anne.duittoz@univ-tours.fr</footnote>
located in the preoptic area (POA) (ewes) or in the antero-ventral periventricular region (AVPV) (rodents)². Kisspeptin stimulates GnRH secretion through an action on KISS1-R, also named GPR54 receptor, expressed by GnRH neurons³. In induced ovulators, ovulation is the result of an increase in LH due presumably to an increase in GnRH release, the timing of this increase is given by mating. In alpacas, intra-muscular (IM) administration of seminal plasma during the pre-ovulatory phase was able to induce ovulation in more than 90% of female alpaca⁴,⁵, this action being blocked by administration of a GnRH antagonist. Thus, seminal plasma contains one or several chemicals that can induce ovulation through a GnRH dependent mechanism, this factor was named Ovulation Inducing Factor (OIF)⁶,⁷. Early work in the bactrian camel identified two fractions of seminal plasma that were able to elicit LH secretion from in vitro culture of rat pituitary cells, suggesting the presence of a GnRH-like factor⁸. Seminal plasma was also found to be active in spontaneous ovulators such as gilts, boar 's seminal plasma placed in gilt uterine horn advance ovulation⁹. The authors identified the activity in the 1–10 kDa fraction separated by ultrafiltration and showed that it lost its activity after treatment with pronase⁹. In 2012, two groups identified the neurotrophin beta-nerve growth factor (β-NGF) as the main component of OIF¹⁰,¹¹.

β-NGF is also present in seminal plasma from spontaneously ovulating species such as bull, ram, porcine and stallion¹². Moreover OIF activity was conserved in these species since equine and porcine seminal plasma induced ovulation in prepubertal mice¹³. Although there is strong evidence for an hypothalamic effect of β-NGF, other components may act at different levels such as at the pituitary level¹⁴,¹⁵. The purpose of our study was to determine the neuroendocrine mechanism of action of β-NGF using in vitro primary cultures of GnRH neurons derived from mouse embryonic nasal placode. We used β-NGF enriched fraction of alpaca seminal plasma (aβ-NGF) and recombinant human β-NGF (hβ-NGF). We characterized the structure of purified alpaca β-NGF and compared it to hβ-NGF structure. Our study showed that recombinant hβ-NGF and enriched fraction of alpaca seminal plasma aβ-NGF increased GnRH neurons spontaneous electrical activity in a dose-dependant manner. This effect was mediated by the P75 neurotrophin receptor but not by TrkA. p75 receptors were expressed by olfactory ensheathing cells in vitro and by tanycytes of the median eminence in adult mice, at the vicinity of GnRH nerve terminals suggesting an indirect effect on GnRH neurons.

## Results

**Characterization of alpaca β-NGF.** Alpaca seminal plasma fraction enriched in β-NGF (aβ-NGF) was obtained after fractionating seminal plasma using liquid chromatography. To further determine the sequence of alpaca β-NGF, we performed bottom up and top down proteomic studies. The complete sequence of mature aβ-NGF consists of 118 amino acids displaying 97% sequence identity with llama OIF and 92% identity with human β-NGF (Fig. 1, see Supplementary Information Table S1).

**β-NGF stimulates GnRH neurons spontaneous electrical activity through a P75 mediated action.** To explore the neuroendocrine mechanisms involved in the ovulation inducing properties of OIF we tested the effect of β-NGF enriched fraction (2+3 Fig. 1) of alpaca seminal plasma (aβ-NGF) and recombinant hβ-NGF on spontaneous electrical activity of GnRH neurons cultures derived from mouse embryonic nasal placodes. In control condition, electrical activity occurred at a frequency of 42.38 ± 9.04 action potentials (AP) per minute (AP min⁻¹) (Fig. 2). This frequency was not affected by hβ-NGF 37.5 ng/mL (n = 5, 30.3 ± 3.9 AP min⁻¹) but was increased by 100% in the presence of 75 ng/mL (n = 3, 83.91 ± 9.94 AP min⁻¹, Kruskal-Wallis chi-squared = 14.588 p = 0.0022, and non parametric multiple comparisons p = 0.0164) and by 300% in the presence of hβ-NGF 150 ng/mL (n = 5, 187.09 ± 28.11 AP min⁻¹, Kruskal-Wallis chi-squared = 14.588 p = 0.0022, and non parametric multiple comparisons p = 0.0009). This effect was reversible after washout. This increase in activity was characterized by an important increase in the proportion of short intervals between action potentials. The frequency histogram of the time intervals between AP illustrates this shortening (Fig. 2). Enriched seminal plasma fraction in aβ-NGF 450 ng/mL (protein content) produced an 78% increase in AP frequency (n = 6, 67.02 ± 30.18 to 98.75 ± 28.49 AP min⁻¹, paired t-test t = 3.3048 p = 0.0213, data not shown).

Burst analysis showed that the number of APs within a burst and burst duration increased following the application of aβ-NGF 450 ng/mL. The number of APs within a burst and burst duration increased not immediately after aβ-NGF 450 ng/mL application but the effect took approximately 20 min to develop. The mean intraburst interval between APs was not affected. 3. To determine which receptor was involved in this action, we first used a specific TrkA antagonist, GW441756, at 3.2 µM (IC50 = 2nM). The application of GW441756 prior or during the application of aβ-NGF (450 ng/mL) did not prevent the increase in AP frequency (control period: 57.37 ± 42.23 AP min⁻¹, GW441756 period: 48.53 ± 23.76 AP min⁻¹, GW441756 + aβ-NGF period: 124.22 ± 38.60 AP min⁻¹, paired t-test t = -3.899, p = 0.0299) (Fig. 2).

Since no p75 specific antagonist was available, we then tested a p75 receptor agonist, LM11A31 at 10 nM (EC50 = 100 nM¹⁶) to see whether it could mimic the effect of aβ-NGF on GnRH neurons 'electrical activity. LM11A31 increased AP frequency from 53.36 ± 7.87 AP min⁻¹ in control condition up to 87.07 ± 6.22 AP min⁻¹ (n = 6, paired t-test t = -3.3908, p-value = 0.01944, data not shown). To assess the interaction of aβ-NGF with neurotrophin receptor p75, aβ-NGF was superimposed onto the human β-NGF moiety of human β-NGF-p75 complex (PDB 1SG1 entry code). Sixteen residues of site I and site II of aβ-NGF are involved in interactions with 25 residues of p75, as similarly observed in p75-hβ-NGF complex (18 aβ-NGF residues interacting with 21 p75 residues )¹⁷. The buried area upon complex formation with p75 was slightly higher in the case of aβ-NGF-p75 complex (1318 Å²) compared to hβ-NGF-p75 complex (1248 Å²). In addition, the cluster of basic residues at interaction site I is totally conserved, providing charge complementarity between this positively charged region of aβ-NGF and the corresponding negative region of p75 in this zone of interaction (Fig. 4).
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Enrichment and identification of alpaca β-NGF. (a) Representative chromatogram of seminal plasma after a Superdex 75 gel-filtration (GF) chromatography. Seven fractions were collected between 22 and 39 min of retention time. (b) Representative 12% SDS-PAGE profile of the 1–7 eluted GF fractions containing proteins of 13–14 kDa. Lane MM represents marker proteins with their molecular masses (kDa) indicated on the left. (b') Immunodetection of the 1–7 eluted fractions using an antibody raised against human β-NGF. Western blotting analysis confirmed the presence of β-NGF and the cropped blot revealed a strong enrichment in fraction 2 and 3, at the expected molecular weight. (c) On the left, MALDI-TOF mass spectrum, in the mass range 4,000–30,000 m/z, of intact β-NGF enriched by GF chromatography. The three major peaks correspond to monocharged molecular specie [M+H]+ at 13,216.34 m/z, to bicharged molecular specie [M+2H]2+ at 6,609.10 m/z and to dimer 2 × [M+H]+ at 26435.55 m/z. On the right, multicharged ions MS spectra for intact β-NGF obtained with nanoESI-LTQ Velos Orbitrap. The average molecular mass of the mature β-NGF was observed at m/z 1,202.40; 1,322.54; 1,469.38; 1,652.93 and 1,888.92 respectively for the (+ 11), (+ 10), (9+), (8+) and (+7) charge states. The monocharged molecular specie [M+H]+ was calculated (mean) at 13,216.42 ± 0.02 Daltons. Below MS/MS fragmentation and product ion map from HCD spectra from multicharged intact β-NGF at m/z 1,469 and 1,652 with or without salt adducts. Fragmentation map shows the coverage of b and y fragment ions confirming N- and C-ter sequences.</figcaption></figure>

**hθ-NGF increased intracellular calcium peak frequency.** GnRH pulsatile secretion not only depends on electrical activity but is directly related to periodic intracellular increases in [Ca²⁺]ᵢ. In mouse embryonic nasal explant cultures GnRH neurons present spontaneous [Ca²⁺]ᵢ elevations lasting a few seconds (3–10) and occurring at a rate of 3–6 per minute. We named these events: [Ca²⁺]ᵢ peaks. hβ-NGF 75 ng/mL induced a decrease in the mean duration of [Ca²⁺]ᵢ interpeaks intervals (IPI) (Fig. 5). During the control period the mean IPI was 10.42 ± 0.58 s (n = 79). After hβ-NGF application the mean IPI decreased significantly to 7.47 ± 0.17 s (Wilcoxon rank sum test, W = 2,165, p = 0.00089). This effect on IPI was either due to an increase in [Ca²⁺]ᵢ peak frequency and/or to an increase in [Ca²⁺]ᵢ peak duration. We observed both mechanisms within the same experiment.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Effect of β-NGF on GnRH neurones spontaneous electrical activity. (a) Spontaneous electrical activity recorded using the loose patch mode from a GnRH neuron (vertical bar = 40 pA, horizontal bar = 2 s). (b) Electrical activity recorded on the same neuron in (a) after the application of 150 ng/mL hβ-NGF. (c) Boxplot representing AP frequency as a function of hβ-NGF concentration. Different letters indicate a statistically significant difference (Kruskal-Wallis test, chi-squared = 14.588, p = 0.0022, and non parametric multiple comparisons p = 0.0098). (d) Control recording using the loose patch configuration of a GnRH neuron. (e) Trace of electrical activity recorded from the same neuron as in (d) after the application of 3.2 µM GW441756, TrkA antagonist. No effect on electrical activity was detected. (f) Electrical activity of the same neuron in (d,e). after application of 450 ng/mL aβ-NGF + 3.2 µM GW441756, (g) boxplot representing mean AP frequency. TrkA antagonist failed to block the increase in electrical activity induced by aβ-NGF (n = 4 experiments, paired t-test, t = -3.8992, p = 0.0299). a',b',d',e',f' corresponding frequency histogram of interspike intervals.</figcaption></figure>

**p75 but not TrkA receptors are expressed in mouse embryonic nasal placode cultures.** p75 receptors but not TrkA receptors were detected using immunohistochemistry in GnRH neuron cultures derived from mouse embryonic nasal placodes. p75 labeling was found absent on GnRH neurons, whereas bipolar cells closely associated along to GnRH neurites and known as olfactory ensheathing cells (OEC) displayed a strong membrane labeling for P75 (Fig. 6a). Higher magnification revealed that GnRH neurites are tightly associated with p75 membrane labelled OEC, some GnRH neurites wrap around p75 membrane labelled OEC, forming some ring-like structures around GnRH-labelled structures (Fig. 6a').

**p75 receptors are expressed in mouse median eminence.** Immunodetection of p75 receptor was performed on three hypothalamic sections from three adult GFAP-GFP female mice. Figure 6b shows that neither GFAP-GFP positive cells or GnRH-labeled fibers located in the ME were labeled by p75 antibody. However, we observed a strong membrane labeling of β2-tanycytes in the vicinity of GnRH axons terminals (Fig. 6b).

**hβ-NGF induced neuro-glial plasticity.** Thus, the stimulatory effect of β-NGF on GnRH neurons electrical activity should be the result of p75 signaling in glial cells forming the microenvironment of GnRH neurons. To detect a possible neuroglial plasticity, we assessed polysialic acid N-Cell Adhesion Molecule (PSA-NCAM) immunoreactivity, in control or hβ-NGF treated cultures derived from mouse embryonic nasal placodes. hβ-NGF was applied for one hour at 5 ng/mL, a concentration without effect on APs frequency and 50 ng/mL, a concentration that was effective on electrical activity. In all conditions, we found a strong membrane labeling of GnRH neurons with PSA-NCAM antibody (stars, Fig. 7) even though some GnRH somas and fibers were not labelled by PSA-NCAM antibody (arrow's head, Fig. 7). In control and hβ-NGF 5 ng/mL conditions, we observed a classical patchy PSA-NCAM labeling on GnRH neuron membrane. However, after exposure to hβ-NGF 50 ng/mL, PSA-NCAM labeling defined larger membrane patches. We assessed quantitatively these observations by measuring the volume of PSA-NCAM labelled membrane in 3D reconstructed images. In control and 5 ng/mL hβ-NGF conditions, 100% and 91.6 ± 7.4% of the total PSA-NCAM immunoreactivity volume were represented by patches smaller than 3,000 µm³. In contrast, patches with volumes greater than 3,000 µm³ formed
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Burst analysis. (a) Trace recording representing a 30 min long record (horizontal bar = 5 min, vertical bar = 40 pA) and below enlarged view of typical APs or bursts of APs (horizontal bar = 12 s, vertical bar = 40 pA). The recording was divided in 5 min blocks (Control = first 5 min control condition, NGF.05 = 5 to 10 min of recording which represent 0–5 min in the presence of aβ-NGF, NGF.10 = 5–10 min after application of aβ-NGF, NGF.15 = 10–15min, NGF.20 = 15–20 min, NGF.25 = 20–25 min) with different color codes. (b) Stripcharts representing the mean number of APs within a burst, burst duration and mean intraburst interval. Following bath application of 450 ng/mL aβ-NGF, the frequency of APs progressively increased, the number of events within a burst and the duration of burst increased. The number of bursts was 3, 2, 2, 5, 5, 5 during Control, NGF.05, NGF.10, NGF.15, NGF.20, NGF.25 periods respectively. The number of events per burst were 6 ± 1.15, 4.0 ± 0.0, 3.0 ± 0.0, 8.6 ± 3.9, 9.8 ± 3.2, 32.0 ± 16.4 during control, NGF.05, NGF.10, NGF.15, NGF.20 and NGF.25 periods respectively. The mean burst duration (s) was 7.9 ± 2.4, 3.1 ± 0.6, 1.7 ± 0.2, 8.4 ± 2.0, 15.3 ± 7.3, 41.9 ± 17.5 for control, NGF.05, NGF.10, NGF.15, NGF.20 and NGF.25 periods respectively. The mean intraburst interval between APs (s) was: 1.5 ± 0.1, 1.0 ± 0.2, 0.8 ± 0.1, 1.5 ± 0.4, 1.5 ± 0.3, 1.7 ± 0.2 during control, NGF.05, NGF.10, NGF.15, NGF.20 and NGF.25 periods respectively.</figcaption></figure>

76.8 ± 4.3% of the total PSA-NCAM immunoreactive volume in 50 ng/mL hβ-NGF treated cultures (Kruskal–Wallis rank sum test on percentage of total PSA volume arcsin transformation of class 3,000–6,000 (two classes < 3,000 and 3,000–6,000) chi-squared = 8.2927, p-value = 0.01582, Fig. 7).

Discussion

β-NGF has been identified as the OIF in llama's seminal plasma. The purified llama β-NGF was shown to induce
ovulation in female llamas after intra-muscular (IM) or intra-vaginal administration⁴. Most in vivo studies used
seminal plasma or β-NGF enriched seminal plasma for administration to female mammals (alpaca⁴,¹¹, mouse¹³,
cow¹⁸). Only one study demonstrated that purified llama β-NGF administrated IM to female llama induced
ovulation. Since other factors with a potential effect on pituitary cells have been previously detected in seminal
plasma, it is important to compare enriched β-NGF fraction of seminal plasma effects to recombinant β-NGF
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Structure of aβ-NGF complexed with p75. (a) Cartoon representation of the aβ-NGF-p75 complex where p75 (magenta) binds alongside of the alpaca homodimer model (monomer A and monomer B colored blue and cyan respectively). The solvent-accessible surface of p75 and aβ-NGF homodimer is shown in (b,c). respectively and colored according to values of electrostatic potentials (blue:positive charges; red: negative charges). (d) aβ-NGF has been rotated 90° around Y-axis to show the conserved cluster of positive charges in aβ-NGF interaction site I. Atomic coordinates of aβ-NGF were obtained by comparative modeling based on lama β-NGF using Swiss-Model server as described in Supplementary Information section. The figure was prepared with PYMOL software.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Effect of hβ-NGF on [Ca²⁺]ᵢ in GnRH neurones. (a,b) Traces of [Ca²⁺]ᵢ fluctuations in individual GnRH neurons during control condition. Boxes correspond to enlarged views of traces. a' and b' traces recorded from same neurons after application of 75 ng/mL hβ-NGF showing either an increase in [Ca²⁺]ᵢ frequency (a) and an increase in [Ca²⁺]ᵢ duration (b). (c) Boxplot of IPI values in control and hβ-NGF conditions. Different letters indicate a statistically significant difference (Wilcoxon rank sum test W = 4,076, p-value = 0.0008944, n = 79)</figcaption></figure>

effects. In the present study we aimed to investigate whether β-NGF, either purified from seminal plasma or recombinant form, could act directly on GnRH neurons using an in vitro primary culture of GnRH neurons derived from mouse E11.5 embryo nasal placodes.

We first characterized the mature structure of alpaca β-NGF using top down and bottom up mass spectrometry and found 97% sequence identity with lama OIF and 92% sequence identity with human β-NGF.

Both alpaca seminal plasma enriched in β-NGF and recombinant hβ-NGF induced an increase in spontaneous electrical activity of mouse GnRH neurones in vitro. Spontaneous electrical activity is driven by action potentials (APs). Neuronal APs involve a depolarizing step mostly due to inward sodium current provoked by the opening of voltage-gated sodium channels, this step last less than a 1 ms, followed by a repolarizing step due to potassium outward current provoked by the opening of delayed voltage potassium channels and the closure of sodium channels that are becoming refractory to depolarization. Neurons can express a variety of high voltage activated calcium channels (Cav1 (L type), Cav2 (N, P/Q and R types) which will open during the depolarizing
<figure><img src="image_9.png" /><figcaption>**Figure 6.** p75 immunodetection in vitro and in vivo. (a) In vitro immunohistochemistry. DAPI Nuclear staining. GnRH-GFP staining shows GnRH neurones (cell bodies and neurites. p75 staining shows membrane labelled cells soma and processes, GnRH neurones were not labelled. No TrkA staining was detected (bar = 10 µm). (a') In vitro immunohistochemistry of GnRH neuron was revealed with a anti-GFP (Alexa 488) here in green staining stroma and dendrites/axons. p75 was labelled with an anti P75 (Alexa 633) here in red false color. Orthographic projection shows the vicinity of GnRH neuron fibers with membranes of p75 labelled cells (bar = 20 µm). (b) In vivo immunohistochemistry: microphotographs of single confocal plane of hypothalamic section of a hGFAP-GFP female mouse at the level of bregma 2.06 according to the Franklin and Paxinos Mouse Brain Stereotaxic Coordinates atlas. DAPI staining. GFAP-GFP staining labelled dispersed glial cells close to the third ventricle and within the ME. GnRH staining labelled GnRH fibers in the lateral and medial external ME. P75 staining labelled tanycytes in the ME (bar = 20 µm).</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 7.** Effect of hβ-NGF on PSA-NCAM immunoreactivity. (a) Control condition. (b) After 1 h treatment with 5 ng/mL hβ-NGF. (c) After 1 h treatment with 50 ng/mL hβ-NGF. (a1,b1,c1) PSA-NCAM immunoreactivity. (a2,b2,c2) GnRH immunoreactivity. (a3,b3,c3) are merge images. (a4,b4,c4) are 3D reconstructed images using Imaris software. False color indicate the average volume of labelled objects (see color scale from 10 to 3,000 µm).</figcaption></figure>
step induced by sodium channels. The contribution of calcium current to the depolarizing phase of AP is small since calcium channels are expressed at low density compared to sodium channels. A third group of calcium channels (Cav3 or T type) are found to be expressed in neurons presenting regular phasic activity like GnRH neurons¹⁹, these channels are activated by low voltage and can contribute significantly to the first depolarizing step of the AP. Altogether, these channels give rise to calcium entry during the first phase of AP whose duration is less than 1 ms. Once Ca²⁺ ions reach the submembrane region of the cytoplasm, they are buffered by calcium binding proteins (calbindin, calmodulin...²⁰) or by surrounding negatively charged proteins. Calcium ions can also be transported outside by the Na+/Ca²⁺ exchanger. These mechanisms drive Ca²⁺ ions out in a few microseconds and preclude calcium ions to diffuse outside a several 100 nm radius microdomain located around the channel inner mouth²¹. At the synapse, calcium channels are clustered and form microdomains where neurotransmitter vesicles are docked and ready for release. For neuroendocrine cells, calcium channels and large dense core vesicles (LDCV) are not located in specialized regions, such as the synapse, they are spaced at larger intervals. In neuroendocrine cells, a single AP will allow only a small number of LDCVs to be released, only the one that are located within the nanodomain of a calcium channel inner mouth will be able to fuse their membrane and release their content²². In contrast, a burst of AP will maintain calcium entry and [Ca²⁺]ᵢ will build up, diffuse and recruit more LDVC.

Spontaneous electrical activity of GnRH neurons is carried by bursts of action potentials whose frequency determines the rhythmicity of GnRH pulsatile secretion²³. In contrast to other neuroendocrine neurons, GnRH neuron bursts contain only two or three APs²³,²⁴. There are also differences according to the biological model used, bursting is observed in more than 75% of mouse hypothalamic GnRH cell line GT1-7 cells and in 75% of neurones recorded from anterior hypothalamic area (AHA) slices²⁵ although this proportion varied depending on slice's thickness and orientation²⁶. Using a ventral pharyngeal approach in anesthetized GnRH-GFP mice, Constantin et al.²⁷ found a majority of neurons to have an irregular activity and only 15% were bursting. According to Lee et al.²⁸, who defined a burst as "two or more spikes clustered in a 4 s interspikes intervals", GnRH neurons derived from embryonic nasal placodes displayed bursts among irregular spiking, as described in AHA slices²³. In the presence of hβ-NGF, the median value of interspike interval was decreased in a dose-dependent manner and revealed an increase in the proportion of bursts. Bursts are important for neuroendocrine secretion since they induce a sufficient long lasting rise in intracellular calcium to allow hormone-containing LDCVs to be released (see²³ for review).

Bursting electrical activity gives rise to [Ca²⁺]ᵢ peaks lasting several seconds that allow GnRH local release. Partial synchronization of [Ca²⁺]ᵢ transients between more than 30% of GnRH neurons was shown to be correlated with in vitro episodic GnRH release²⁹. Here we showed that the inter [Ca²⁺]ᵢ peak interval (IPI) was significantly decreased after application of hβ-NGF. This decrease was either due to an increase in [Ca²⁺]ᵢ peak frequency or/and to an increase in [Ca²⁺]ᵢ peak duration. Both effects could be the result of an increase in the frequency and/or duration of bursts. Therefore, hβ-NGF induced an increase in electrical activity by increasing APs frequency and by changing APs temporal distribution in favor of a bursting pattern. These changes were associated to an increase in [Ca²⁺]ᵢ peaks frequency and/or duration and an increase in synchronization that could led to an increase in the probability of GnRH LDCVs release.

β-NGF plays multiple roles in neuronal development, function, survival, and neurite growth through an action on two types of membrane receptors: TrkA and p75³⁰. Taken independently, p75 and TrkA neurotrophin receptors binds β-NGF with low affinity, but when they are both expressed in the same cell, p75 and TrkA forms a high affinity complex³¹. It is now well established that p75 signaling in the absence of TrkA can lead to apoptosis, cell survival, myelination or neurite outgrowth according to the partner proteins recruited in different cellular contexts³²,³³. Here we showed that the specific TrkA antagonist, GW441756, used at a concentration ten times greater than its IC50, failed to block the stimulatory effect of 450 ng/mL aβ-NGF on GnRH neurones electrical activity. Moreover, the specific p75 agonist, LM11A31¹⁶, mimicked the effect of aβ-NGF on GnRH neuron activity even at a concentration 10 times lower than its EC50. Thus, in our in vitro model, β-NGF seemed to increase GnRH neuron electrical activity through a p75-mediated mechanism. In embryonic nasal placode explant cultures, we did not observed p75 labelling on GnRH neurons in contrast to³⁴. However, olfactory ensheathing cells (OECs), that form the cellular microenvironment of GnRH neurones in vitro, strongly expressed p75³⁴⁻³⁷. OECs are glial cells originating from neural crest, they ensheath GnRH neurons at their exit from the medial olfactory pit. OECs are highly plastic, they can display various glial markers and they present myelinating-like properties. OECs play a role in the nasal migration of GnRH neurons during embryogenesis³⁸. OEC are tightly associated with GnRH neurons as illustrated in Fig. 7a'. Thus, at low magnification, p75 OECs's membrane labeling might appear as belonging to GnRH neuron membrane. Similarly, the use of secondary antibodies coupled to probes whose fluorescence spectra are close may induce a fluorescence leak and an apparent membrane labelling of GnRH neurons. To prevent spectra overlap we used alexa 488- and alexa 633-coupled secondary antibodies to label anti-GnRH and anti-p75 primary antibodies, their respective spectra being sufficiently separated to avoid any overlap. Nevertheless, Raucci et al.³⁴ confirmed the expression of p75 using single cell RT-PCR on isolated GnRH neurons in vitro. Thus we cannot totally exclude the fact that p75 receptor could be expressed at low level in GnRH neurons. However the strong expression exhibited by surrounding OECs, suggests that most of effects produced by β-NGF should be mediated through OECs p75 receptors. OECs are known to be confined in the olfactory region at adulthood, however Geller et al.³⁷ demonstrated that brain lipid binding protein (BLBP)-labelled OEC progenitors accompanied GnRH neurons along peripherin fibers into the diencephalon of E15.5 mouse embryos.

Other mechanisms can be involved, for example in llama, Carrasco et al.³⁹, identified neurons expressing p75 and TrkA receptors in female llama's hypothalamus. Most hypothalamic areas contained TrkA immunoreactive cells with a greater labelling in the diagonal band of Broca. P75 labelled neurons were found in the diagonal band of Broca, the lateral preoptic area and caudal parts of the hypothalamus. P75 was also detected in ependymal cells
boarding the lateral ventricle and p75 fibers were observed in the *organum vasculorum of the lamina terminalis* (OVLT), and the ventral part of arcuate nucleus (ArcN) and the median eminence (ME). GnRH neurons were not labelled by p75 antibody and only 2.5% were labelled by TrkA antibody suggesting that they are not a direct target for β-NGF. TrkA or p75 expressing neurons could act as a relay for β-NGF action on GnRH neurons³⁹. The neuroanatomical location of TrkA and p75 neurons close to the third ventricle suggests that they could be sensitive to β-NGF in the cerebrospinal fluid. Although there is some evidence that β-NGF can cross the blood brain barrier⁴⁰, the presence of TrkA and p75 neurons in the preoptic area and the hypothalamus is not a proof of interaction between TrkA/p75 expressing neurons and GnRH neurons. However, the presence of p75 immunoreactive fibers in OVLT and in the ArcN and ME, regions located outside the brain blood barrier, is of particular interest since circulating β-NGF could reach the dense network of GnRH fibers in these areas. We particularly focused on the ME area since it is the structure where GnRH secretion take place and can be regulated by tanycytes (see Sharif et al.⁴¹ for a review) and by kisspeptin terminals found in the internal layer⁴².

We found a specific p75 expression in adult mouse tanycytes of the ME, in close vicinity to GnRH fibers as previously reported in rat⁴³. Tanycytes β1 and β2 form the tanycyte population localized within the median eminence. They are considered as stem cells and are involved in numerous neuroendocrine regulations and play a key role in the control of metabolism and seasonal functions⁴⁴,⁴⁵. Tanycytes are known to regulate GnRH neuronal activity through a mechanism involving the release of prostaglandin E2 (PGE2)⁴¹,⁴⁶. They also display a cellular plasticity across the estrus cycles in female mammals⁴⁷. To assess this cellular plasticity, we quantified PSA-NCAM labelling on GnRH neurons. The large negatively charged PSA chain of NCAM prevents homophilic interaction between cells. We used an antibody specific to the PSA group⁴⁸, highlighting membrane regions on GnRH neurons free of adhesion with surrounding cells. hβ-NGF treatment induced the apparition of enlarged regions labeled by PSA, suggesting a remodeling of cell-to-cell interactions. We have previously shown that PSA-NCAM is present on both OECs and GnRH neurons in vitro⁴⁹. A previous study showed that antagonizing p75 signaling by incubation with a p75 antibody induced the defasciculation of olfactory fibers in vitro, but also decreased the surface of GnRH neurons' soma⁴⁴. Rearranging the glial microenvironment of GnRH neurones may lead to changes in spontaneous electrical activity as it has been described with the tripartite synapse⁵⁰⁻⁵³. In astrocytes, several mechanisms are involved in the regulation of neuronal synaptic transmission. Among them are glutamate transporters (GLT1), potassium channels, hemichannels that can regulate the synaptic strength and shape the synaptic response by controlling local glutamate and K⁺ concentrations resulting from neuronal activity⁵³. Astrocytes can also release gliotransmitters such as GABA, glutamate, D-serine, adenosine/ATP that will act at the post-synaptic level (see⁵³ for a review). Whether these mechanisms are functional in other types of glial cells such as OEC and hypothalamic tanycytes remains to be ascertained.

Kisseptin terminals are found in the median eminence but in the internal layer, in contrast with GnRH terminals that are located in the external layer⁴². Attempts to detect synapses between GnRH and kisspeptin terminals in the ME have failed, thus direct regulation of GnRH secretion by kisspeptin terminals in the ME is not supported by neuroanatomy studies. However, tanycytes β spread their cytoplasm from the third ventricle floor lining to the capillaries in ME external layer⁵⁴ and could be the relay between Kp and GnRH terminals.

If this mechanism is really functional in a spontaneous ovulator species is not known. Previous work by¹³ showed that injection of alpaca seminal plasma to prepubertal mice induced ovulation. However alpaca seminal plasma contains other components that could trigger ovulation either by targeting the pituitary level or the ovary level. Gene networks governing these two modes of ovulation may have been selected in various mammalian taxa under different selection pressure. Interestingly, some mammalian species displayed the two modes of ovulation. For example the domestic cat (*Felix cati*)⁵⁵,⁵⁶ and wild felid species such as the lion (*Panthera leo*)⁵⁶,⁵⁷ displayed mixed oestrus cycles. Whether the β-NGF and the kisspeptin pathways interact is not known and should be investigated in the future, not only on rodents but also in species displaying both modes of ovulation. In conclusion, using an in vitro approach, we showed that β-NGF stimulated the electrical activity of GnRH neurons favoring a bursting behavior. By increasing the frequency and/or duration of [Ca²⁺]ᵢ peaks, this bursting behavior should increase the probability of GnRH release that could lead to a global increase in GnRH secretion driving LH release. All these effects seem to involve a neuro-glia! plasticity mediated by p75 receptor expressed by glial cells forming the microenvironment of GnRH neurons.

Methods

**Characterization of alpaca β-NGF from seminal plasma.** Technical details can be found in Supplementary Information (Methods). In order to enrich β-NGF, seminal plasma (SP) from alpaca (*Vicugna pacos*) was fractionated by liquid chromatography (LC). SP were separated from spermatozoa by centrifugation (10,000×g for 10 min at room temperature). The supernatant was centrifuged again (10,000×g for 10 min at room temperature) and stored at – 80 °C.

*Gel-filtration chromatography.* Protein concentration of SP was determined using the Uptima BC Assay kit (Interchim, Montluçon, France) according to manufacturer's instructions.

*Gel electrophoresis.* SDS-PAGE was carried out according to Laemmli's method⁵⁸ using a Mini-Protean II electrophoresis cell (BioRad, Marnes-la-Coquette, France).

*Western blotting.* Proteins were separated by SDS-PAGE and transferred onto nitrocellulose membrane (GE Healthcare Life Sciences WHATMAN). Semidry transfer of proteins was performed over 1.15 h at 0.8 mA/cm². The membrane was then incubated with rabbit polyclonal antibody directed against human β-NGF (1/5000, v/v, sc11358, Santa Cruz). After being washed with TBS with 0.5% Tween 20, the membrane was incubated 1
h at 30 °C with goat anti-rabbit conjugated with peroxidase (dilution 1:5000). Blot was developed using a mixture of two chemiluminescence substrates developing kit (GE Healthcare AmershamTH ECL SelectTH Western blotting detection Reagent RPN2235 and Supersignal West Pico #34087 Chemiluminescent Substrate Thermo Scientific). Images were digitized with a cooled CCD camera (ImageMaster VDS-CL, Amersham Biosciences, GE Healthcare Life Sciences).

**Bottom-up proteomic.** In order to evaluate the purity level of β-NGF, GF fractions, bottom-up proteomic approach was carried out using nano-liquid chromatography tandem high-resolution mass spectrometry (nanoLC-MS/MS). Proteomic experiment using bottom-up approach was performed on a dual linear ion trap Fourier Transform Mass Spectrometer LTQ Orbitrap Velos (Thermo Fisher Scientific, Bremen, Germany) coupled to an Ultimate®3000 RSLC Ultra High Pressure Liquid Chromatographer (Thermo Fisher Scientific, Bremen, Germany) controlled by Chromeleon version 6.8 SR11 software. Data were acquired using Xcalibur version 2.1 software (Thermo Fisher Scientific, San Jose, CA), in a positive data-dependent mode in the 300–1,800 m/z mass range. Resolution in the Orbitrap was set at R = 60.000. The 20 most intense peptide ions with charge states > 2 were sequentially isolated (isolation width 2 m/z, 1 microscan) and fragmented in the high-pressure linear ion trap using CID (collision induced dissociation) mode (collision energy 35%, activation time 10 ms, Qz 0.25). Dynamic exclusion was activated during 30 seconds with a repeat count of 1. The lock mass was enabled for accurate mass measurements. Proteins were identified using SEQUEST algorithm in Proteome Discoverer software (version 2.1, Thermo Fisher Scientific) against the NCBInr database with mammalia taxonomy (download January 2018). Scaffold version 4.8.4 software (Proteome Software, Portland, USA) was used to validate protein identifications. Peptide and proteins identifications were accepted if they could be established at greater than 95.0% and 99% probability as specified by the Peptide Prophet algorithm⁵⁹ and the Protein Prophet algorithm⁶⁰, respectively. Protein identifications were accepted if they contained at least two identified peptides. The abundance of identified proteins was estimated by calculating the emPAI (Exponentially Modified Protein Abundance Index)⁶¹ using Scaffold Q+ software (version 4.8.4, Proteome Software, Portland, USA).

**Top-down proteomic.** Firstly, the fraction containing β-NGF (RT 27–31 min) was analyzed by Matrix Assisted Laser Desorption/Ionization Time-Of-Flight Mass Spectrometry (MALDI-TOF MS) in order to measure whole mass of intact β-NGF and to evaluate heterogeneity between samples.

Samples were desalted spectra were acquired using a Bruker UltrafleXtreme MALDI-TOF instrument (Bruker Daltonics, Bremen, Germany) equipped with a Smartbeam laser at 2 kHz laser repetition rate following an automated method controlled by FlexControl 3.0 software (Bruker Daltonics, Bremen, Germany). Spectra were obtained in positive linear ion mode in the 1,000–20,000 m/z range. Secondly, 3 µL of β-NGF desalted by SPE and diluted with 40 µL of 50% (v/v) methanol/1% (v/v) FA was loaded into a metalized nanoelectrospray needle (PicoTip emitters, New Objective, USA). Top-Down MS and MS analyses were performed on a LTQ Orbitrap Velos instrument (Thermo Fisher Scientific, Bremen, Germany) operating in positive mode. MS and MS/MS data were acquired using Xcalibur version 2.1 software (Thermo Fisher Scientific, San Jose, CA).

**3D modeling.** The alpaca β-NGF sequence was submitted to the automated protein structure homology-modelling server Swiss-Model (http://swissmodel.expasy.org). A 3D model of alpaca β-NGF was generated based on the llama OIF X-ray structure (PDB 4EFV entry code) as a template. Because alpaca β-NGF exhibits 97% sequence identity with llama OIF, the modeling procedure was straightforward and the corresponding root mean square (rms) between both structures was very low (0.096 Å). To assess the interaction of alpaca NGF with neurotrophin receptor p75, alpaca β-NGF was superimposed onto the human β-NGF moiety of human β-NGF-p75 complex (PDB 1SG1 entry code). As expected from sequence identity of 92% between both β-NGF, their 3D structure were perfectly superimposed (not shown). Most, if not all interactions, between p75 and each β-NGF were conserved, as shown by analysis of the interface of p75 with each β-NGF using COCOMAPS tools (https://www.molnac.unisa.it/BioTools/cocomaps/).

**Animals.** hGFAP-GFP mice (Charles River), Swiss mice (Janvier Labs, Saint Genest L'isle, France), GnRH-GFP mice (donated by Daniel Spergel⁶²) were housed in the rodent facility at UEPAO, INRA centre Val de Loire, Nouzilly France. Each cage housed four mice at – 20 °C, which were used to being housed together.

**Primary culture of embryonic nasal placode explants.** Pregnant GnRH-GFP mice at embryonic day 11.5 were killed by cervical dislocation. Briefly, embryos were placed in ice cold Gey's balanced salt solution (Eurobio) enriched with 5% glucose (Sigma-Aldrich). Nasal pits were dissected out and adhered onto glass coverslips (Ø14 mm Marienfeld) coated with chicken plasma/thrombin (Sigma-Aldrich). Nasal explants were maintained in a defined serum-free medium (SFM)²⁹,⁶³.

**Electrophysiology.** Nine cultures for experiments with αβ-NGF and twenty-one cultures for experiments with hβ-NGF, at 7–14 DIV were transferred to a recording chamber (RC-25 reference 64-0232 from Warner Instruments) under a fixed-stage microscope (Leica DM IRB 275064) in a bath solution containing (in mM) NaCl 124, KCl 4.1, CaCl₂ (2H₂O) 1, MgCl₂ (6H₂O) 0.3, MgSO₄ (7H₂O) 0.4, NaHCO₃ 20.1, NaH₂PO₄ (2H₂O) 1, and glucose 35; pH (NaOH) 7.4 (see⁶³ for details). Loose patch recordings were performed at (37 °C) (Automatic Temperature controller TC-324B Warner Instruments Corp) and acquired using Clampex 10 (25 kHz sampling and low pass 5 kHz filtering) (Molecular Devices) connected to an Axopatch 200A patch-clamp amplifier and digitized with a Digidata 1550B A/D converter⁶³. Patch electrodes (5–10 MΩ)⁶³ were filled with a solution con-
taining (in mM): NaCl 120, KCl 4.7, CaCl₂ (2H₂O) 2.6, MgCl₂ (6H₂O) 2, MgSO₄ (7H₂O) 0.7, HEPES 10 and glucose 10; pH = 7.4 with NaOH⁶³. Peak detection was performed with Clampfit. Burst analysis used the Poisson surprise method in Clampfit software.

**Calcium Imaging.** Calcium imaging recordings were performed between 7 and 14 DIV (n = 79 neurons, N = 3 cultures). Explant cultures were loaded with the calcium-fluorescent dye Green-1 AM (10 µM, Thermofisher, Villebon sur Yvette, France) (see⁶³ for details). Cultures were mounted in a perfusion chamber (Warner Instruments, Hamden, CT, USA) and loaded with calcium Green-1 AM solution containing 2.5mM probenecid (blocker of organic anion transporter) and 10 mM Hepes in SFM⁶³. Fluorescence was visualized using an inverted microscope (DM-IRB; Leica Microsystems GmbH) and acquired through a × 20 objective using a cooled intensified charge-coupled device camera (CoolSNAP HQ2; Ropper Instruments, Photometrics) using Metafluor software (Molecular Devices). Excitation wavelengths were provided by a xenon lamp through a YFP-2427b filter block (BrightLine Single band filter set Semrock, excitation bandpass filter 488–512 nm, dichroic edge wavelength at 520 nm, emission bandpass filter 528.5–555.5 nm, 93% transmission)⁶³. Fluctuations in [Ca²⁺]ᵢ were analyzed with Metafluor software. Each cell was individually identified as a circular region of interest (ROI). Excitation duration was set to 100 ms and time lapse acquisition was set to 1 Hz. Calcium Green-1 mean fluorescence intensity for each region of interest was calculated and analyzed with Excel (Microsoft Corp). A calcium peak was defined as a value greater than the average of the five previous points + 1.5 × standard deviation⁶³.

**Tissue processing.** hGFAP-GFP mice were anaesthetised with pentobarbital (150 mg/kg IP) and transcarially perfused with 15 mL of 4% paraformaldehyde in PBS. Brains were dissected out of the skull and post-fixed in 4% PAF for 24 h. The brains were then soaked in 20% sucrose solution (PBS with 0.1% sodium azide) at 4 °C for 24 h. Brains were then embedded in Tissu-Tek (Sakura, Finetek Europe BV, Zoeterwoude, The Netherlands) and frozen in isopentane solution at – 50 °C and stored at – 20 °C until sectioning. Coronal 20 µm sections were made using a cryostat (CM3050S Leica Microsystemes from the beginning of the hypothalamus to the mammillary bodies. Tissue sections were stored at – 20 °C until processing for immunohistochemistry.

**Immunocytochemistry.** *Cell culture p75, TrkA, GnRH-GFP immunohistochemistry.* After 10 days in vitro (DIV), cultures were fixed in 4% formaldehyde (30–45 min). Nonspecific binding sites were blocked in PBS supplemented with 5% normal goat or horse serum or 2% BSA, 0.3% Triton X-100 and 0.1% sodium azide and incubated overnight at 4 °C with polyclonal goat anti-GFP (1:2000, Abcam ab5450), polyclonal rabbit anti-p75 (1:1000, promega ref G3231), monoclonal mouse anti-TrkA (1:1000, Merk Millipore ref 06-574). The next day, cultures were rinsed in PBS and incubated in secondary antibodies (anti-goat Alexa Fluor 488, 1:500 and anti-rabbit Alexa 633, 1:100, and anti-mouse Alexa 594, 1:500, for 2 h at room temperature), then rinsed and counterstained using the nuclear dye DAPI (1 mM), 1:5000, Sigma Aldrich, Illkirch, France) and mounted with Fluomount G (Sigma Aldrich, Illkirch, France).

*Cell culture PSA-NCAM, GnRH and GFAP immunohistochemistry (IHC).* Six cultures aged between 13 DIV and 16 DIV were exposed either to control SFM or to hβ-NGF 5 or 50 ng/mL for 1 h and fixed in 4% formaldehyde (30–45 min). Nonspecific binding sites were blocked as described above, and incubated overnight at 4 °C with monoclonal mouse IgM anti-polysialylated neural cell adhesion molecule (PSA-NCAM) (1:5000, Abcam, ref: abc0019). The following day after PBS rinses, cell cultures were incubated in goat anti mouse IgM coupled to Rhodamine (1:500, Sigma Aldrich, ref: AP128R) for two hours at room temperature. After PBS rinses, cell cultures were then incubated with monoclonal mouse IgG anti-GnRH SMI-41 (Abcam ref: ab24563) and polyclonal IgG rabbit anti-GFAP (1:500, DAKO ref: Z0334) overnight at + 4 °C. The following day after PBS rinses, cell cultures were incubated with biotinylated goat anti mouse IgG (1:200) and goat anti rabbit IgG coupled to Alexa 642 (1:500) for 2 h at room temperature. After PBS rinses, cell cultures were incubated in streptavidine coupled to AMCA (1:500) for 30 min, and mounted with Fluomount G (Sigma Aldrich, Illkirch, France).

*In vivo p75, GnRH, GFAP-GFP immunohistochemistry (IHC).* Nonspecific binding sites were blocked as previously described and incubated overnight at + 4 °C with polyclonal chicken anti-GFP (1:10000, Aves Lab), polyclonal rabbit anti-p75 (1:1000, promega ref G3231), polyclonal sheep anti-GnRH (1:5000, gift from Dr Alain Caraty)³⁷,⁶³. Sections were then rinsed in PBS and incubated for two hours at room temperature in secondary antibodies Alexa 488 donkey anti-chicken IgY (1:500), Alexa 633 donkey anti rabbit IgG (1:500) and Alexa 546 donkey anti sheep IgG (1:500). All secondary antibodies are from Thermofisher, Villebon sur Yvette, France). Sections were conterstained with DAPI, and mounted in Fluomount G (Sigma Aldrich, Illkirch, France).

**Fluorescence microscopy and image analysis.** Stack image acquisition was performed using a laser confocal microscope (LSM700 Zeiss, Germany) and the associate Zen software. Images 1024 × 1024 pixels were acquired sequentially. At the objective Plan apochromat 40×/1.3 oil DIC M27, the x, y resolutions were 0.156 µm, 0.156 µm and the z resolution was 0.5 µm. Images were analysed using FIJI software. Volumes of PSA-NCAM labelled structures from IHC image stacks were determined with Imaris software (release 9.2, Bitplane AG, Zurich, Switzerland). In image stacks 8 bits encoded, a threshold of intensity of 10 was used to eliminate background and pixel with intensity above 247 were eliminated. The segmented volumes quantified where above 2 µm³.
**Statistical analysis.** All statistical analyses were performed using R Studio⁶⁴. Two-samples comparisons were performed using Student test if the variable distribution followed a Normal law (Shapiro-Wilk test) and variances were homoscedatic (Fisher test). If these hypothesis were not verified we used non-parametric Wilcoxon ranking test). Multiple samples comparisons were performed using non parametric Kruskal-Wallis followed by post hoc non-parametric multiple comparisons (mctp function, nparcomp package). Chi-square test was used to compare the proportion of prepubertal mice having ovulated according to the treatment group. The number of COCs obtained according the treatment was analysed using a permutation multivariate analysis of variance. The comparison of average percentage of volumes for PSA-NCAM labelling between control and hβ-NGF treatments was made after arcsin transformation by Kruskal-Wallis non parametric test.

**Ethics statement.** All animals 'procedures were carried out with permission from the Val de Loire Ethics Committee of the French Ministry of Agriculture, in accordance with the European legislation for animal experimentation (Animal housing Authorization F371752, protocol authorisation no. 10378-2017052213185465 from the French Ministry of Agriculture). All experiments described were conducted in accord with accepted standards of humane animal care and all efforts were made to minimise suffering.

**Data availability**

The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

Received: 30 July 2019; Accepted: 11 June 2020

Published online: 30 June 2020

**References**

1. Adams, G. P., Ratto, M. H., Silva, M. E. & Carrasco, R. A. Ovulation-inducing factor (oif/ngf) in seminal plasma: A review and update. *Reprod. Domest. Anim.* **51** (Suppl 2), 4–17. https://doi.org/10.1111/rda.12795 (2016).

2. Caraty, A. et al. Kisspeptin synchronizes preovulatory surges in cyclical ewes and causes ovulation in seasonally acyclic ewes. *Endocrinology* **148**, 5258–5267. https://doi.org/10.1210/en.2007-0554 (2007).

3. Messenger, S. et al. Kisspeptin directly stimulates gonadotropin-releasing hormone release via g protein-coupled receptor 54. *Proc. Natl. Acad. Sci. U.S.A.* **102**, 1761–1766. https://doi.org/10.1073/pnas.0409330102 (2005).

4. Ratto, M. H., Huanca, W., Singh, J. & Adams, G. P. Local versus systemic effect of ovulation-inducing factor in the seminal plasma of alpacas. *Reprod. Biol. Endocrinol.* **3**, 29. https://doi.org/10.1186/1477-7827-3-29 (2005).

5. Silva, M. E. et al. Cetrorelix suppresses the preovulatory lh surge and ovulation induced by ovulation-inducing factor (oif) present in llama seminal plasma. *Reprod. Biol. Endocrinol.* **9**, 74. https://doi.org/10.1186/1477-7827-9-74 (2011).

6. Zhao, X. X., Li, X. L. & Chen, B. X. Isolation of ovulation-inducing factors in the seminal plasma of bactrian camel (*Camelus bactrianus*) by deae-cellulose chromatography. *Reprod. Domest. Anim.* **36**, 177–181 (2001).

7. Tanco, V. M., Ratto, M. H., Lazzarotto, M. & Adams, G. P. Dose-response of female llamas to ovulation-inducing factor from seminal plasma. *Biol. Reprod.* **85**, 452–456. https://doi.org/10.1095/biolreprod.111.091876 (2011).

8. Li, X. & Zhao, X. Separation and purification of ovulation-inducing factors in the seminal plasma of the bactrian camel (*Camelus bactrianus*). *Vet. Res. Commun.* **28**, 235–245 (2004).

9. Waberski, D. et al. Advanced ovulation in gilts by the intrauterine application of a low molecular mass pronase-sensitive fraction of boar seminal plasma. *J. Reprod. Fertil.* **105**, 247–252 (1995).

10. Ratto, M. H. et al. The nerve of ovulation-inducing factor in semen. *Proc. Natl. Acad. Sci. U.S.A.* **109**, 15042–15047. https://doi.org/10.1073/pnas.1206273109 (2012).

11. Kershaw-Young, C. M., Druart, X., Vaughan, J. & Maxwell, W. M. C. beta-nerve growth factor is a major component of alpaca seminal plasma and induces ovulation in female alpacas. *Reprod. Fertil. Dev.* **24**, 1093–1097. https://doi.org/10.1071/RD12039 (2012).

12. Druart, X. et al. Proteomic characterization and cross species comparison of mammalian seminal plasma. *J. Proteomics* **91**, 13–22. https://doi.org/10.1016/j.jprot.2013.05.029 (2013).

13. Bogle, O. A., Ratto, M. H. & Adams, G. P. Evidence for the conservation of biological activity of ovulation-inducing factor in seminal plasma. *Reproduction* **142**, 277–283. https://doi.org/10.1530/REP-11-0042 (2011).

14. Paolicchi, F., Urquieta, B., Del Valle, L. & Bustos-Obregón, E. Biological activity of the seminal plasma of alpacas: Stimulus for the production of lh by pituitary cells. *Anim. Reprod. Sci.* **54**, 203–210 (1999).

15. Bogle, O. A., Ratto, M. H. & Adams, G. P. Ovulation-inducing factor (oif) induces lh secretion from pituitary cells. *Anim. Reprod. Sci.* **133**, 117–122. https://doi.org/10.1016/j.anireprosci.2012.06.006 (2012).

16. Massa, S. M. et al. Small, nonpeptide p75ntr ligands induce survival signaling and inhibit prongf-induced death. *J. Neurosci.* **26**, 5288–5300. https://doi.org/10.1523/JNEUROSCI.3547-05.2006 (2006).

17. He, X.-L. & Garcia, K. C. Structure of nerve growth factor complexed with the shared neurotrophin receptor p75. *Science* **304**, 870–875. https://doi.org/10.1126/science.1095190 (2004).

18. Tribulo, P., Bogle, O., Mapletoft, R. J. & Adams, G. P. Bioactivity of ovulation inducing factor (or nerve growth factor) in bovine seminal plasma and its effects on ovarian function in cattle. *Theriogenology* **83**, 1394–1401. https://doi.org/10.1016/j.theriogenology.2014.12.014 (2015).

19. Nunemaker, C. S., DeFazio, R. A. & Moenter, S. M. Calcium current subtypes in gnrh neurons. *Biol. Reprod.* **69**, 1914–1922. https://doi.org/10.1095/biolreprod.103.019265 (2003).

20. Schwaller, B. Cytosolic ca2+ buffers. *Cold Spring Harb. Perspect. Biol.* **2**, a004051. https://doi.org/10.1101/cshperspect.a004051 (2010).

21. Naraghi, M., Müller, T. H. & Neher, E. Two-dimensional determination of the cellular ca2+ binding in bovine chromaffin cells. *Biophys. J.* **75**, 1635–1647. https://doi.org/10.1016/S0006-3495(98)77606-4 (1998).

22. Oheim, M., Kirchhoff, F. & Stühmer, W. Calcium microdomains in regulated exocytosis. *Cell Calcium* **40**, 423–439. https://doi.org/10.1016/j.ceca.2006.08.007 (2006).

23. Moenter, S. M. Identified gnrh neuron electrophysiology: A decade of study. *Brain Res.* **1364**, 10–24. https://doi.org/10.1016/j.brainres.2010.09.066 (2010).

24. Rossoni, E. et al. Emergent synchronous bursting of oxytocin neuronal network. *PLoS Comput. Biol.* **4**, e1000123. https://doi.org/10.1371/journal.pcbi.1000123 (2008).
25. Herbison, A. E. & Moenter, S. M. Depolarising and hyperpolarising actions of gaba(a) receptor activation on gonadotrophin-releasing hormone neurones: Towards an emerging consensus. *J. Neuroendocrinol.* **23**, 557–569. https://doi.org/10.1111/j.1365-2826.2011.02145.x (2011).

26. Constantin, S. et al. Gnrh neuron firing and response to gaba in vitro depend on acute brain slice thickness and orientation. *Endocrinology* **153**, 3758–3769. https://doi.org/10.1210/en.2012-1126 (2012).

27. Constantin, S., Iremonger, K. J. & Herbison, A. E. In vivo recordings of gnrh neuron firing reveal heterogeneity and dependence upon gabaa receptor signaling. *J. Neurosci.* **33**, 9394–9401. https://doi.org/10.1523/JNEUROSCI.0533-13.2013 (2013).

28. Lee, K., Liu, X. & Herbison, A. E. Burst firing in gonadotrophin-releasing hormone neurones does not require ionotrophic gaba or glutamate receptor activation. *J. Neuroendocrinol.* **24**, 1476–1483. https://doi.org/10.1111/j.1365-2826.2012.02360.x (2012).

29. Constantin, S., Caraty, A., Wray, S. & Duittoz, A. H. Development of gonadotropin-releasing hormone-1 secretion in mouse nasal explants. *Endocrinology* **150**, 3221–3227. https://doi.org/10.1210/en.2008-1711 (2009).

30. Chao, M. V. Neurotrophins and their receptors: A convergence point for many signalling pathways. *Nat. Rev. Neurosci.* **4**, 299–309. https://doi.org/10.1038/nrn1078 (2003).

31. Chao, C.C.-K. et al. Involvement of gas7 in nerve growth factor-independent and dependent cell processes in pc12 cells. *J. Neurosci.* **Res. 74**, 248–254. https://doi.org/10.1002/jnr.10763 (2003).

32. Roux, P. P. & Barker, P. A. Neurotrophin signaling through the p75 neurotrophin receptor. *Prog. Neurobiol.* **67**, 203–233 (2002).

33. Cragnolini, A. B. & Friedman, W. J. The function of p75ntr in glia. *Trends Neurosci.* **31**, 99–104. https://doi.org/10.1016/j.tins.2007.11.005 (2008).

34. Raucci, F., Tiong, J. D. & Wray, S. P75 nerve growth factor receptors modulate development of gnrh neurons and olfactory ensheating cells. *Front. Neurosci.* **7**, 262. https://doi.org/10.3389/fnins.2013.00262 (2013).

35. Ramón-Cueto, A., Pérez, J. & Nieto-Sampedro, M. In vitro enfolding of olfactory neurites by p75 ngf receptor positive ensheathing cells from adult rat olfactory bulb. *Eur. J. Neurosci.* **5**, 1172–1180 (1993).

36. Geller, S., Kolasa, E., Tillet, Y., Duittoz, A. & Vaudin, P. Olfactory ensheathing cells form the microenvironment of migrating gnrh-1 neurons during mouse development. *Glia* **61**, 550–566. https://doi.org/10.1002/glia.22455 (2013).

37. Geller, S. et al. Rostro-caudal maturation of glial cells in the accessory olfactory system during development: Involvement in outgrowth of gnrh neurites. *Eur. J. Neurosci.* **46**, 2596–2607. https://doi.org/10.1111/ejn.13732 (2017).

38. Barraud, P., St John, J. A., Stolt, C. C., Wegner, M. & Baker, C. V. H. Olfactory ensheathing glia are required for embryonic olfactory axon targeting and the migration of gonadotropin-releasing hormone neurons. *Biol. Open* **2**, 750–759. https://doi.org/10.1242/bio.20135249 (2013).

39. Carrasco, R. A., Singh, J. & Adams, G. P. The relationship between gonadotropin releasing hormone and ovulation inducing factor/ nerve growth factor receptors in the hypothalamus of the llama. *Reprod. Biol. Endocrinol.* **16**, 83. https://doi.org/10.1186/s1295 8-018-0402-6 (2018).

40. Pan, W., Banks, W. A. & Kastin, A. J. Permeability of the blood-brain barrier to neurotrophins. *Brain Res.* **788**, 87–94. https://doi. org/10.1016/s0006-8993(97)01525-4 (1998).

41. Sharif, A., Baroncini, M. & Prevot, V. Role of glia in the regulation of gonadotropin-releasing hormone neuronal activity and secretion. *Neuroendocrinology* **98**, 1–15. https://doi.org/10.1159/000351867 (2013).

42. Mikkelsen, J. D. & Simonneaux, V. The neuroanatomy of the kisspeptin system in the mammalian brain. *Peptides* **30**, 26–33. https://doi.org/10.1016/j.peptides.2008.09.004 (2009).

43. Gudiño-Cabrera, G. & Nieto-Sampedro, M. Schwann-like macroglia in adult rat brain. *Glia* **30**, 49–63 (2000).

44. Ebling, F. J. P. & Lewis, J. E. Tanocytes and hypothalamic control of energy metabolism. *Glia* **66**, 1176–1184. https://doi.org/10.1002/ glia.23303 (2018).

45. García-Cáceres, C. et al. Role of astrocytes, microglia, and tanocytes in brain control of systemic metabolism. *Nat. Neurosci.* **22**, 7–14. https://doi.org/10.1038/s41593-018-0286-y (2019).

46. Clasadonte, J., Sharif, A., Baroncini, M. & Prevot, V. Gliotransmission by prostaglandin e(2): A prerequisite for gnrh neuronal function?. *Front. Endocrinol. (Lausanne)* **2**, 91. https://doi.org/10.3389/fendo.2011.00091 (2011).

47. Prevot, V. et al. Neuronal-glial-endothelial interactions and cell plasticity in the postnatal hypothalamus: Implications for the neuroendocrine control of reproduction. *Psychoneuroendocrinology* **32** (Suppl 1), S46–S51. https://doi.org/10.1016/j.psyne uen.2007.03.018 (2007).

48. Rougon, G., Dubois, C., Buckley, N., Magnani, J. L. & Zollinger, W. A monoclonal antibody against meningococcus group b polysaccharides distinguishes embryonic from adult n-cam. *J. Cell Biol.* **103**, 2429–2437 (1986).

49. Franceschini, I., Desroziers, E., Caraty, A. & Duittoz, A. The intimate relationship of gonadotropin-releasing hormone neurons with the polysialylated neural cell adhesion molecule revisited across development and adult plasticity. *Eur. J. Neurosci.* **32**, 2031–2041. https://doi.org/10.1111/j.1460-9568.2010.07517.x (2010).

50. Haydon, P. G. & Carmignoto, G. Astrocyte control of synaptic transmission and neurovascular coupling. *Physiol. Rev.* **86**, 1009– 1031. https://doi.org/10.1152/physrev.00049.2005 (2006).

51. Halassa, M. M., Fellin, T. & Haydon, P. G. The tripartite synapse: roles for gliotransmission in health and disease. *Trends Mol. Med.* **13**, 54–63. https://doi.org/10.1016/j.molmed.2006.12.005 (2007).

52. Halassa, M. M. & Haydon, P. G. Integrated brain circuits: astrocytic networks modulate neuronal activity and behavior. *Annu. Rev. Phys.* **72**, 335–355. https://doi.org/10.1146/annurev-physiol-021909-135843 (2010).

53. Papouin, T., Dunphy, J., Tolman, M., Foley, J. C. & Haydon, P. G. Astrocytic control of synaptic function. *Philos. Trans. R. Soc. Lond. B Biol. Sci.* **372**. https://doi.org/10.1098/rstb.2016.0154 (2017).

54. Rodríguez, E. M. et al. Hypothalamic tanocytes: A key component of brain-endocrine interaction. *Int. Rev. Cytol.* **247**, 89–164. https://doi.org/10.1016/S0074-7696(05)47003-5 (2005).

55. Gudermuth, D. F., Newton, L., Daels, P. & Concannon, P. Incidence of spontaneous ovulation in young, group-housed cats based on serum and faecal concentrations of progesterone. *J. Reprod. Fertil.* **51**, 177–184 (1997).

56. Pelican, K. M., Wildt, D. E., Ottinger, M. A. & Howard, J. Priming with progestin, but not gnrh antagonist, induces a consistent endocrine response to exogenous gonadotropins in induced and spontaneously ovulating cats. *Domest. Anim. Endocrinol.* **34**, 160–175. https://doi.org/10.1016/j.domaniend.2007.01.002 (2008).

57. Schramm, R. D., Briggs, M. B. & Reeves, J. J. Spontaneous and induced ovulation in the lion (*Panthera leo*). *Zoo Biol.* **13**, 301–307 (1994).

58. Laemmli, U. K. Cleavage of structural proteins during the assembly of the head of bacteriophage t4. *Nature* **227**, 680–685 (1970).

59. Keller, A., Nesvizhskii, A. I., Kolker, E. & Aebersold, R. Empirical statistical model to estimate the accuracy of peptide identifications made by ms/ms and database search. *Anal. Chem.* **74**, 5383–5392 (2002).

60. Nesvizhskii, A. I., Keller, A., Kolker, E. & Aebersold, R. A statistical model for identifying proteins by tandem mass spectrometry. *Anal. Chem.* **75**, 4646–4658 (2003).

61. Ishima, R. & Torchia, D. A. Error estimation and global fitting in transverse-relaxation dispersion experiments to determine chemical-exchange parameters. *J. Biomol. NMR* **32**, 41–54. https://doi.org/10.1007/s10858-005-3593-z (2005).

62. Spergel, D. J., Krüth, U., Hanley, D. F., Sprengel, R. & Seeburg, P. H. Gaba- and glutamate-activated channels in green fluorescent protein-tagged gonadotropin-releasing hormone neurons in transgenic mice. *J. Neurosci.* **19**, 2037–2050 (1999).

63. Pinet-Charvet, C. et al. Gnrh episodic secretion is altered by pharmacological blockade of gap junctions: Possible involvement of glial cells. *Endocrinology* **157**, 304–322. https://doi.org/10.1210/en.2015-1437 (2016).
64. R Core Team. R: A language and environment for statistical computing (2017).

## Acknowledgements

This research was supported by INRA PHASE department through the IMPAIC project coordinated by C. Taragnat and IFCE NOVO project. We thank UEPAO for the animal care. This work has benefited from the facilities and expertise of the "Plateforme d'Imagerie Cellulaire" (PIC) of the UMR-PRC.

## Author contributions

C.P.-C. performed electrophysiology recordings and analysis, wrote the corresponding methods and result sections. R.F. performed cell cultures, calcium imaging, immunohistochemistry, confocal analysis. F.D.-T. participated to the writing. S.D.G. collected alpaca seminal plasma and participated to the writing. X.D. contributed to alpaca seminal plasma purification and participated to the writing. G.T. contributed to alpaca plasma seminal purification. C.T. contributed to characterize NGF enriched fraction binding to p75. A.-P.T.-G. contributed to alpaca plasma seminal purification. V.L. did the proteomics, analysed the results and wrote the corresponding material and methods and results sections. T.M. did the 3D alpaca NGF-P75 complex modelling, wrote the corresponding material and methods. X.C. did the 3D analysis of PSA-NCAM labelling and contributed to the writing. A.D. designed the experiments and managed the project, wrote the paper, analysed calcium imaging data. All authors reviewed the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41598-020-67665-4.

Correspondence and requests for materials should be addressed to A.H.D.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020