SEMA6A drives GnRH neuron-dependent puberty onset by tuning median eminence vascular permeability ================================================================================================= * A. Lettieri * R. Oleari * M.H. van den Munkhof * E.Y. van Battum * M.G. Verhagen * C. Tacconi * M. Spreafico * A.J.J. Paganoni * R. Azzarelli * V. Andre’ * F. Amoruso * L. Palazzolo * I. Eberini * L. Dunkel * SR. Howard * A. Fantin * R.J. Pasterkamp * A. Cariboni ## ABSTRACT Innervation of the hypothalamic median eminence by Gonadotropin-Releasing Hormone (GnRH) neurons is vital to ensure puberty onset and successful reproduction. However, the molecular and cellular mechanisms underlying median eminence development and pubertal timing are incompletely understood. Here we show that Semaphorin-6A is strongly expressed by median eminence-resident oligodendrocytes positioned adjacent to GnRH neuron projections and fenestrated capillaries, and that Semaphorin-6A is required for GnRH neuron innervation and puberty onset. *In vitro* and *in vivo* experiments reveal an unexpected function for Semaphorin-6A, via its receptor Plexin-A2, in the control of median eminence vascular permeability to maintain neuroendocrine homeostasis. To support the significance of these findings in humans, we identify patients with delayed puberty carrying a novel pathogenic variant of *SEMA6A*. In all, our data reveal an undescribed role for Semaphorin-6A in regulating GnRH neuron patterning by tuning the median eminence vascular barrier and thereby controlling puberty onset. ## INTRODUCTION Puberty onset and reproductive capacity depend on the projection of hypothalamic Gonadotropin-Releasing Hormone (GnRH) neurons to the median eminence (ME) and the release of GnRH into the hypothalamic-pituitary-gonadal (HPG) axis1,2. During development, GnRH neurons are born in the nasal placode and migrate along the terminal nerve (TN) and in association with olfactory ensheathing cells (OECs) to the medial preoptic area (MPOA) of the hypothalamus3,4. Here, GnRH neurons cease to migrate and extend their terminal projections to contact the ME, where complex neuro-glial-endothelial (NGE) interactions occur to regulate GnRH homeostasis5,6. Defective GnRH neuron migration and/or innervation of the ME leads to reduced GnRH secretion, and consequently delayed or absent puberty. These are the hallmarks of reproductive disorders such as Congenital Hypogonadotropic Hypogonadism (HH), and also contribute to the etiology of familial Delayed Puberty (DP)7,8. The genetic determinants underlying these reproductive disorders remain unknown in many congenital cases, suggesting the existence of unexplored mechanisms controlling the development and function of the GnRH neuronal system relevant to these conditions9,10. Recent studies have highlighted key roles for semaphorins and their receptors in the GnRH neuronal system11–18. These molecules have been implicated in GnRH neuron development and adult GnRH neuronal plasticity, which are both required to ensure GnRH secretion19. Consistently, pathogenic variants in various semaphorin signaling genes have been identified in patients with HH and DP12,20–24. Semaphorin-6A (SEMA6A) is a transmembrane semaphorin mainly known for its roles during axon guidance25,26, neuronal positioning27,28 and angiogenesis29. To exert its biological functions, SEMA6A signals as a ligand through Plexin-A2 and Plexin-A4 receptors28,30,31. Consistent with a role in the GnRH neuronal system, SEMA6A is expressed in the developing mouse olfactory system and hypothalamus32,33. However, its function in the GnRH system or in mechanisms underlying puberty onset are unknown. Here, by combining mouse and human genetics with tailored *in silico*, *in vitro* and *in vivo* models, we unveil a previously unexplored role for SEMA6A in GnRH neuron innervation and ME vascular permeability, and identify SEMA6A as an additional player in the regulation of puberty timing. ## RESULTS ### SEMA6A is expressed along GnRH neuron migration route during brain development To investigate the potential role of SEMA6A in GnRH neuron development and function, we first studied SEMA6A expression in *wild-type* (WT) mouse tissues at stages relevant for GnRH neuron development and function. Double immunofluorescence staining in the nasal compartment at embryonic day (E) 12.5 and E14.5, corresponding to the start and peak of GnRH neuron migration, respectively, showed co-localization of SEMA6A and Peripherin, a TN marker (Fig. 1a). In contrast, GnRH neurons in the nasal compartment, identified by anti-GnRH staining, did not express SEMA6A at E12.5 and E14.5 (Fig. 1b). Similarly, no SEMA6A expression was found in OECs, marked with an anti-Brain lipid-binding protein (BLBP) antibody (Fig. 1c).  [Figure 1.](http://medrxiv.org/content/early/2023/10/26/2023.10.24.23297449/F1) Figure 1. SEMA6A is expressed in territories relevant for embryonic GnRH neuron development. **(a)** Coronal sections of E12.5 (left) and E14.5 (right) mouse heads at the nose level were immunolabelled for SEMA6A (green) and Peripherin (red). Solid arrowheads indicate examples of Peripherin-positive axons that express SEMA6A at the VNO level and in the nasal parenchyma. **(b)** Coronal sections of E12.5 (left) and E14.5 (right) mouse heads at the nose level were immunolabelled for SEMA6A (green) and GnRH (red). Empty arrowheads indicate examples of GnRH neurons that lack SEMA6A expression at the VNO level and in the nasal parenchyma. **(c)** Coronal sections of E12.5 (left) and E14.5 (right) mouse heads at the nose level were immunolabelled for SEMA6A (green) and BLBP (red). Empty arrowheads indicate examples of OECs that lack SEMA6A expression at the VNO level and in the nasal parenchyma. **(d)** Coronal sections of E18.5 mouse brains at the MPOA (left) and ME (right) levels were immunolabelled for SEMA6A (green) and GnRH (red). Empty arrowheads indicate examples of GnRH neurons that lack SEMA6A expression in the cell body (left) and neurite (right). **(e)** Sagittal view of a CS19 human embryo in a schematic drawing (left) and sagittal sections of schematics squared area immunolabeled for SEMA6A (right). Black arrows indicate examples of SEMA6A-expressing axons in the nasal parenchyma. All sections were counterstained with DAPI. Dotted boxes indicate areas shown at higher magnification on the right of the corresponding panel. Abbreviations: VNO, vomeronasal organ; NS, nasal septum; FB, forebrain; OE, olfactory epithelium; OB, olfactory bulb; MPOA, medial preoptic area; ME, median eminence; 3v, third ventricle; CS, Carnegie Stage. Scale bars: 500 μm (**e**, low magnification), 250 μm (**e**, second panel), 200 μm (**a-d**, low magnifications), 125 μm (**e**, high magnifications), 100 μm (**d**, low magnification of right panel), 50 μm (**a-d**, high magnifications). At E14.5, when the majority of GnRH neurons are scattered in the MPOA of the hypothalamus along TN fibers14,34, SEMA6A was strongly expressed in this region but it did not co-localize with GnRH+ neurons or Peripherin+ TN fibers (Supplementary Fig. 1a, b). Next, expression of SEMA6A was examined at E18.5, when GnRH neurons have attained their final position in the hypothalamus and start to project to the ME35. Although GnRH+ cell bodies in the MPOA and GnRH+ neurites in the ME did not express SEMA6A (Fig. 1d), strong and punctate SEMA6A expression was detected in the ME in close proximity to GnRH+ axonal terminals (Fig. 1d, right panels). Our observations in mouse embryos were confirmed in human embryonic tissue at Carnegie Stage (CS) 19, corresponding to E12.5 in mouse, showing SEMA6A co-localization with axons in the nasal parenchyma (Fig. 1e). Together, these data show that during GnRH neuron migration from the nasal compartment to the MPOA SEMA6A initially localizes in TN fibers, while during late gestation SEMA6A is expressed in close proximity to GnRH+ axons in the ME. ### Loss of SEMA6A does not affect GnRH neuron migration but reduces GnRH innervation of the ME Our expression analysis suggests that SEMA6A may regulate GnRH neuron migration along nasal axons either directly or indirectly, by controlling nasal axon development, consistent with its well-characterized role as an axon guidance cue25,26. To explore this model, nasal axon patterning and GnRH neuron migration were analyzed in E14.5 *Sema6a*+/+ and *Sema6a*-/- mouse embryos using immunohistochemistry, as described previously14. Unexpectedly, no obvious defects in the patterning of Peripherin+ axons were observed in *Sema6a*-/- embryos as compared to *Sema6a*+/+ littermates (Fig. 2a). Also, no significant differences in the number and positioning of GnRH neurons were found (Fig. 2b, c and Supplementary Table 1).  [Figure 2.](http://medrxiv.org/content/early/2023/10/26/2023.10.24.23297449/F2) Figure 2. SEMA6A loss does not impair nasal axon patterning and GnRH neuron migration, but reduces GnRH innervation of the ME. **(a)** Coronal sections of E14.5 mouse heads of the indicated genotypes at the level of the nose (left) and MPOA (right) were immunolabeled for Peripherin (green) to mark nasal axons. White arrows indicate normal nasal axon development in the nose and normal TN fibers in the MPOA. Sections were counterstained with DAPI. **(b)** Coronal sections of E14.5 mouse heads with the indicated genotypes were immunolabelled for GnRH. Black arrows indicate examples of normal GnRH neuron distribution in the nasal parenchyma (left) at the CP (middle) and in the forebrain (right). **(c)** Quantification of the distribution and total number of GnRH neuron in E14.5 heads (n = 4 per group; nose *p* = 0.7974, CP *p* = 0.6360, forebrain *p* = 0.8596, TOT *p* = 0.5670). **(d)** Quantification of GnRH neuron number in E18.5 brains (from OB to MPOA) (*Sema6a+/+*: n = 3; *Sema6a-/-*: n = 5; *p* = 0.7819). **(e)** Coronal sections of E18.5 brains with the indicated genotypes at the level of the ME were immunolabelled for GnRH. Black arrows indicate normal presence of GnRH neuron axon terminals in *Sema6a+/+* embryos, whereas the decreased presence of GnRH neuron axon terminals is indicated with Δ. **(f)** Quantification of GnRH intensity at ME in E18.5 brains (*Sema6a+/+*: n = 3; *Sema6a-/-*: n = 5; *p* = 0.0189). * *p* < 0.05 and NS (not significant) after two-tailed unpaired Student’s t test. Abbreviations: VNO, vomeronasal organ; OE, olfactory epithelium; OB, olfactory bulb; CP, cribriform plate; FB, forebrain; MPOA, medial preoptic area; 3v, third ventricle. Scale bars: 200 μm (**a**), 125 μm (**b, e**). Data are presented as mean ± SD. Source data are provided as a Source Data file. Since the GnRH system appeared intact at E14.5, we next analyzed 1) the final positioning of GnRH neurons in the MPOA, and 2) ME innervation in E18.5 embryos, when GnRH axons can be found at or near the ME35,36. Although *Sema6a*-/- and *Sema6a*+/+ mice displayed a similar number of GnRH neurons in the forebrain (FB) (Fig. 2d and Supplementary Table 1), GnRH neuron innervation of the ME was significantly reduced in *Sema6a*-/- mice by 2-fold (Fig. 2e, f; *p* = 0.0189; two-tailed unpaired Student’s t test). However, the canonical SEMA6A receptors Plexin-A2 and Plexin-A428,30,31 (Supplementary Fig. 2a) were not detected in MPOA-resident GnRH neurons by immunolabelling (Supplementary Fig. 2b, c). This suggests that SEMA6A does not modulate ME innervation by acting on GnRH neurons but rather through functions in other cell types. Together these findings indicate that SEMA6A is dispensable for GnRH neuron migration and positioning, but is required for innervation of the ME by GnRH+ axons. ### Female and male mice lacking SEMA6A exhibit delayed puberty onset and sexual maturation defects Proper innervation of the ME by GnRH axon terminals during development is crucial to ensure GnRH secretion and thereby puberty onset and reproductive function at postnatal stages1,35. To analyze whether SEMA6A expression in the ME is maintained postnatally, we assessed SEMA6A localization relative to GnRH axon terminals in the ME of pre- (P24) and post- (P60) pubertal mice. Interestingly, strong punctate SEMA6A expression was observed in close proximity to GnRH axon terminals in the ME at both stages (Fig. 3a, b). This suggests that SEMA6A might also play a role in maintaining the GnRH system.  [Figure 3.](http://medrxiv.org/content/early/2023/10/26/2023.10.24.23297449/F3) Figure 3. SEMA6A is expressed in the postnatal ME in close proximity to GnRH neuron terminals. **(a,b)** Coronal sections of P24 (**a**) and P60 (**b**) ME were immunolabelled for SEMA6A (green) and GnRH (red) to reveal GnRH axon terminals. White boxes in the left panels indicate areas shown at higher magnification on the right together with single channels for SEMA6A and GnRH. Empty arrowheads indicate examples of GnRH neuron axon terminals that lack SEMA6A expression. White arrows indicate the expression of SEMA6A in the basal part of the ME. All sections were counterstained with DAPI. Abbreviations: 3v, third ventricle. Scale bars: 100 μm (low magnifications), 50 μm (high magnifications). Since SEMA6A was strongly expressed in the ME just before puberty (P24), we next investigated the role of SEMA6A in pubertal timing. The day of vaginal opening (VO), a proxy measurement for mouse pubertal activation of the HPG axis35,37,38, was assessed in *Sema6a+/+*, *Sema6a+/-* and *Sema6a**-/-*** female mice. *Sema6a-/-* females displayed significantly delayed VO compared to WT mice, with 3 out of 9 *Sema6a-/-* females displaying a closed vaginal phenotype at the moment of termination at P45, 54 and 57 (Fig. 4a, red dots and Supplementary Fig. 3a). VO of *Sema6a**-/-*** mice was delayed by 12.1 ± 2.6 days (mean difference ± SEM) compared to *Sema6a+/+*mice, and by 5.9 ± 0.9 days when the 3 female mice still displaying a closed vagina at the time of sacrifice were excluded (Fig. 4a; *Sema6a+/- p* = 0.9376, *Sema6a-/- p* < 0.0001 and delayed *Sema6a-/-*only *p* < 0.0001 vs. *Sema6a+/+*; One-way ANOVA followed by Dunnett’s post-hoc test). *Sema6a-/-* mice had significantly increased weight at VO as they were on average older than their WT littermates at puberty onset (Supplementary Fig. 3b). In subsequent studies, we assessed the appearance of the first estrous (FE) after VO, which corresponds to the first ovulation and indicates the reaching of sexual competence35. We confirmed delayed VO in this new group of mice (Fig. 4b; *Sema6a+/- p* = 0.5186 and *Sema6a-/- p* = 0.0002 vs. *Sema6a+/+*; Two-way ANOVA followed by Dunnett’s post-hoc test) and detected a significant delay in FE (12.9 ± 2.1 days compared to *Sema6a+/+*). One female out of 9 had not reached FE at the time of termination (P58, 24 days after VO) (Fig. 4b, red triangle; *Sema6a+/- p* = 0.5550 and *Sema6a-/- p* < 0.0001 vs. *Sema6a+/+*; Two-way ANOVA followed by Dunnett’s post-hoc test, Supplementary Fig. 3c).  [Figure 4.](http://medrxiv.org/content/early/2023/10/26/2023.10.24.23297449/F4) Figure 4. SEMA6A loss delays puberty and gonadal maturation in female and male mice. **(a)** Age at VO (*Sema6a+/+* n = 15; *Sema6a+/-* n = 17; *Sema6a-/-* n = 9). Red dots indicate *Sema6a-/-* females without VO at sacrifice. **(b)** Ages at VO and FE (*Sema6a+/+* n = 9; *Sema6a+/-* n = 18; *Sema6a-/-* n= 9). Red triangle indicates *Sema6a-/-* female without FE at sacrifice. **(c)** Age at BPS (*Sema6a+/+* n = 11; *Sema6a+/-* n = 22; *Sema6a-/-* n= 14). Red dots indicate *Sema6a-/-*males without BPS at sacrifice. **(d)** Quantification of GnRH intensity at ME (as in e,f) in female and male brains at VO/BPS (females: n = 3 per group; males: n = 5 per group). Red dots indicate *Sema6a-/-* females without VO. **(e, f)** Coronal sections of female (e) and male (f) ME immunolabelled for GnRH. Presence of GnRH neuron axon terminals in *Sema6a+/+* brains (black arrows) is reduced in *Sema6a-/-* (Δ). **(g)** Quantification of ovarian weight (*Sema6a+/+* n = 20, *Sema6a-/-* n= 12). Red dots indicate *Sema6a-/-*female without FE. **(h,i)** H&E histological analysis of mouse ovary folliculogenesis (h) quantified in (i) (n = 6 per group). **(j)** Quantification of testis weight (*Sema6a+/+* n = 10, *Sema6a-/-* n= 12). Red dots indicate *Sema6a-/-* males without BPS. (**k,l**) Immunological analysis of Leydig cells marked by CYP17A1 expression in mouse testes (k), quantified in (l) (n = 6 per group). Red dots indicate *Sema6a-/-* males without BPS. **(m)** RT-qPCR analysis for *Lhb* transcript in male pituitaries (*Sema6a+/+* n = 7, *Sema6a-/-* n= 9) calculated relative to controls using *Gapdh*-normalized Cq threshold values. **p* < 0.05, ***p* < 0.01, \***|*p* < 0.001, \**\*|\**p* < 0.0001 and NS (not significant) after One- (a, c) or Two-way (b) ANOVA followed by Dunnett’s post hoc test, or Two- tailed unpaired Student’s t test (d, g, i, j, l, m) Abbreviations: VO, vagina opening; FE, first estrous; BPS, balanopreputial separation; 3v, third ventricle; del., delayed. Scale bars: 500 μm (**h**), 250 μm (**e, f, k**). Data are presented as mean ± SD. Source data are provided as a Source Data file. In addition to female sexual development, male puberty onset was assessed by examining balanopreputial separation (BPS)37,38. Similar to mutant females, *Sema6a-/-* male mice exhibited a delayed pubertal onset (of 3.3 ± 0.9 days), with 2 out of 14 males not exhibiting BPS at the day of sacrifice (P30) (Fig. 4c, red dots; *Sema6a+/- p* = 0.0615 and *Sema6a-/- p* < 0.0001 vs. *Sema6a+/+*; One-way ANOVA followed by Dunnett’s post-hoc test). *Sema6a-/-* and *Sema6a+/+* males had a similar weight at BPS (Supplementary Fig. 3d). To assess whether the reduced ME innervation observed at E18.5 persisted at later stages, we next analyzed GnRH axon terminals in the ME of adult *Sema6a-/-* mice. Consistent with the observed delayed puberty phenotype, all *Sema6a-/-*females displayed a 2.6-fold decrease in ME innervation compared to age-matched *Sema6a+/+* females (Fig. 4d, e; *p* = 0.0311; two-tailed unpaired Student’s t test). Notably, two *Sema6a-/-* female mice that failed to achieve VO (Fig. 4a) displayed a stronger reduction in ME innervation as compared to *Sema6a-/-* mice exhibiting delayed VO. Similar to our observations at E18.5, the number of GnRH neurons in the MPOA was unaffected in adult *Sema6a-/-* females (Supplementary Table 1; *p* = 0.3258). Similarly, *Sema6a-/-* males had a significant 2.2-fold decrease in ME innervation compared to age- matched *Sema6a+/+* mice (Fig. 4d, f; *p* = 0.0389; two-tailed unpaired Student’s t test). We further found that adult *Sema6a-/-* females exhibited smaller ovaries but normal folliculogenesis (Fig. 4g-i; ovary weight: *p* = 0.0049; follicle subtype number/total number: primary *p* = 0.7666, secondary *p* = 0.6837, pre-antral *p* = 0.3500, early antral *p* = 0.8205; two-tailed unpaired Student’s t). *Sema6a-/-* males displayed decreased gonadal size and a reduced number of Leydig cells, as detected by Cytochrome P450 17A1 (CYP17A1) immunostaining (Fig. 4j-l; testis weight: *p* = 0.0182; CYP17A1+ area: *p* = 0.0002; two-tailed unpaired Student’s t). Consistent with this observation, *Lhb* mRNA levels were significantly reduced in *Sema6a-/-* pituitary glands (Fig. 4m; *p* = 0.0161). Finally, to assess general fertility, breeding records of *Sema6a* mutant mice were compared to those of WT *C57Bl/6J* mice for a period of 9 months (Supplementary Table 2). In this period, *Sema6a+/-* males bred with *Sema6a+/-* females produced fewer litters compared to the breeding of WT mice (68.6% successful in *Sema6a+/-* vs. 80.3% successful in WT). Crosses of *Sema6a-/-* males with either Sema6a+/- or WT females were only successful in 35% of the cases. Out of 9 *Sema6a-/-* males in our colony, 4 were able to produce a litter, which occurred in only 50% of their matings. Overall litter size was not affected in the different crosses. These data support the other changes found in the reproductive system of *Sema6a* mutants described above. Together, these results show a requirement for SEMA6A in the innervation of the ME by GnRH neurons, and in the regulation of puberty onset and sexual maturation. ### SEMA6A is expressed by oligodendrocytes in the proximity of blood capillary loops To establish the mechanisms through which SEMA6A controls GnRH neuron innervation and puberty onset, we first defined the cellular source of SEMA6A at the ME. Interestingly, SEMA6A signal was prominent in the basal part of the ME, which is rich in fenestrated blood capillary loops39. These mediate the delivery of GnRH to the pituitary, but also the entry of cues from peripheral tissues that fine tune ME homeostasis and consequently modulate GnRH structural plasticity40–42. Hence, we evaluated SEMA6A localization relative to blood vessels by immunolabelling SEMA6A in combination with endothelial cell (EC)12 marker Isolectin B4 (IB4). Interestingly, SEMA6A was not detected in ECs but was concentrated and aligned to ME fenestrated blood capillary loops (Fig. 5a). To identify the source of SEMA6A in the ME, additional staining with markers for different cell types (secretory neurons, astrocytes, β-tanycytes and oligodendrocytes (OLs)) involved in ME homeostasis were performed43,44. SEMA6A was not detected in the Synaptophysin+ terminals of neurosecretory neurons (Fig. 5b), but surrounded these structures, as shown for GnRH- enriched terminals. Further, SEMA6A did not localize in Glial fibrillary acidic protein (GFAP)+ astrocytes (Fig. 5c) or in Vimentin+ β-tanycytes (Fig. 5d). Instead, and in agreement with previous reports showing that SEMA6A is required for OL differentiation45, SEMA6A was detected in Oligodendrocyte transcription factor 2 (Oligo2)+ cells (Fig. 5e). SEMA6A selectively localized to Glutathione S-transferase P1 (Gst P1)+ maturing OLs but not oligodendrocyte precursor cells (OPCs) that expressed Platelet-derived growth factor receptor alpha (PDGFR-α) (Supplementary Fig. 4a,b).  [Figure 5.](http://medrxiv.org/content/early/2023/10/26/2023.10.24.23297449/F5) Figure 5. SEMA6A is detected in OLs in the ME and not by other cell types composing the NGE unit. **(a-e)** Coronal sections of adult ME were immunolabelled for SEMA6A (green) together with IB4, Synaptophysin, GFAP, Vimentin and Oligo2 (red) to label blood vessels, neuroendocrine neurons, astrocytes, tanycytes and oligodendrocytes, respectively. White boxes indicate areas shown at higher magnification on the right of the corresponding panel together with single channels for SEMA6A and different cell markers. Empty arrowheads indicate examples of blood vessels penetrating in the ME that lack SEMA6A expression **(a)** and the lack of SEMA6A expression on neuroendocrine axon terminals (**b**), astrocytes (**c**), tanycytes (**d**). Solid arrowheads indicate the expression of SEMA6A on oligodendrocytes (**e**). All sections were counterstained with DAPI. Abbreviations: 3v, third ventricle. Scale bars: 100 μm (**a,b** low magnifications), 50 μm (**a,b**, high magnifications; **c-e**, low magnifications), 25 μm (**c,d**, high magnifications), 12 μm (**e**, high magnification). Protein localization patterns were confirmed by analysis of the Tabula Muris single-cell RNA sequencing (scRNA-seq) database46. Dimensionality reduction combined with Uniform Manifold Approximation and Projection (UMAP) of the FACS/SMARTseq2-based dataset from adult mouse brain identified seven main clusters of cell populations: OPCs (enriched for *Olig2*, *Pdgfra*, *Cspg4* and *Sox10*), maturing OLs (enriched for *Olig2*, *Gstp1*, *Plp1, Cnp, Mog, Mag* and *Sox10*), OLs (enriched for *Mbp*), tanycytes/pericytes (enriched for *Vim*, *Cspg4* and *Pdgfrb*), astrocytes (enriched for *Gfap*), neurons (enriched for *Syp*, *Eno2* and *Rbpfox3*) and ECs (enriched for *Emcn*) (Supplementary Fig. 4c). *Sema6a* transcripts were highly enriched in premyelinating maturing OLs, with fewer transcripts in OPCs and ECs. Expression of *Sema6a* in mature OLs, tanycytes, astrocytes and neurons was very rare (Supplementary Fig. 4c). These findings therefore suggest a possible role for SEMA6A produced by maturing OLs in the regulation of ME vascular permeability, which is crucial for proper GnRH homeostasis5,6,39. ### Abluminal SEMA6A induces vascular permeability To explore the hypothesis that SEMA6A regulates ME vascular permeability, we tested the ability of SEMA6A to modulate the barrier function of EC monolayers. Trans-endothelial electrical resistance (TEER) was measured in response to SEMA6A stimulation from either the apical (luminal) or basal (parenchymal) side47 (Fig. 6a). Although SEMA6A is produced as a transmembrane protein, its extracellular domain can be released to generate a secreted ectodomain48. Therefore, we transfected COS-7 cells with a vector encoding mouse SEMA6A ectodomain (SEMA6AECTO) and used conditioned medium to stimulate ECs.  [Figure 6.](http://medrxiv.org/content/early/2023/10/26/2023.10.24.23297449/F6) Figure 6. Abluminal SEMA6A regulates vascular permeability via Plexin- A2. **(a)** Schematic of vascular permeability measurement by TEER. Apical (red) and basal (blue) compartments represent luminal and parenchymal sides of a blood vessel, respectively. Created with BioRender.com. **(b)** TEER quantification in mBECs apically (red) or basally (blue) treated with conditioned media from SEMA6AECTO or mock-transfected COS-7 cells. TEER values are expressed as percentage of TEER at T0, ranging 46-80 Ωcm2. Graph shows one out of n = 2 independent experiments (n = 3 per group) **(c,d)** scRNA-seq analysis of adult mBECs from the EC atlas dataset. UMAP plots show distinct cell types (**c,** left panel) and *Plxna2* transcript levels (**c**, right panel). Violin plots compare *Plxna2* transcript levels in different EC subsets (**d**, left panel). *Plxna1*, *Plxna3* and *Plxna4* transcripts could not be detected (ND) in any EC subpopulation (**d**, right panel). **(e)** X,y maximum intensity projections (MIP) of mBECs immunolabelled for Plexin-A2 (green) and luminal EC marker ICAM-2 (red) counterstained with DAPI, including control staining for the secondary antibody used to detect Plexin-A2 (anti-goat, left). Yellow line indicates the single optical y,z cross section displayed on the side. Bottom panel shows 3D surface rendering of the z-stack used to generate the x,y MIP. **(f)** TEER quantification in HUVECs treated basally with conditioned media from SEMA6AECTO or mock-transfected COS-7 cells, following knockdown of *PLXNA2* (*sh*RNA B, green; C, orange) or in controls (*sh*Scr, blue). TEER values are expressed as percentage of TEER at T, ranging 36-44 Ωcm2. Graph shows one out of n = 2 independent experiments (n = 3 per group). * refers to SEMA6AECTO vs mock, # refers to SEMA6AECTO *shScr* vs SEMA6AECTO *shPLXNA2*. * *p* < 0.05, ** or ## *p* < 0.01, \***| or ### *p* < 0.001 and \**\*|\* or #### *p* < 0.0001 after Two-way ANOVA followed by Tukey post-hoc test. Scale bars: 8 μm (**e**). Abbreviations: mBECs, mouse brain endothelial cells; aECs, arterial ECs; cECs, capillary ECs; vECs, venous ECs, fECs, fenestrated capillary ECs; UMAP, uniform manifold approximation and projection; ND, not detected. Data are presented as mean ± SD. Source data are provided as a Source Data file. The ability of SEMA6AECTO to affect permeability was initially examined in human umbilical vein endothelial cells (HUVECs), a well-established human model of ECs, whose sensitivity to hyperpermeability-inducing agents was confirmed by the significant TEER reduction induced by apical exposure to vascular endothelial growth factor A (VEGFA), a well-known potent inducer of vascular hyperpermeability49–51 (Supplementary Fig. 5a). Exposure of the apical surface of the HUVEC monolayer to SEMA6AECTO did not affect TEER (Supplementary Fig. 5b; SEMA6AECTO apical vs mock apical, 60 min *p* = 0.0011, 120 min *p* = 0.9998, 180 min *p* = 0.9810; Two-way ANOVA followed by Tukey’s post-hoc test). In contrast, TEER was significantly reduced, and hence permeability was increased, following treatment of the basal compartment with SEMA6AECTO, as compared to HUVECs treated with control medium (Supplementary Fig. 5b; SEMA6AECTO basal vs mock basal: 60 min *p* = 0.0365, 120 min *p* = 0.0006, 180 min *p* < 0.0001; Two-way ANOVA followed by Tukey’s post-hoc test). This indicates that SEMA6A can selectively induce vascular permeability in an EC monolayer. Next, to better recapitulate ME vascular barrier *in vitro*, mouse brain endothelial cells (mBECs) freshly isolated from adult WT mice were analyzed50,52. mBECs were first validated for purity by double immunostaining for the pan-endothelial markers Platelet endothelial cell adhesion molecule (PECAM1) and Cadherin-5 (Supplementary Fig. 5c), and for responsiveness to VEGFA (Supplementary Fig. 5d). In agreement with results obtained in HUVECs, application of SEMA6AECTO to the apical compartment did not affect mBEC TEER (Fig. 6b; SEMA6AECTO apical vs mock apical: 30 min *p* = 0.3798, 60 min *p* = 0.8985, 120 min *p* = 0.7426, 180 min *p* = 0.8563, 240 min *p* = 0.3279; Two-way ANOVA followed by Tukey’s post-hoc test). Conversely, TEER was significantly decreased when SEMA6AECTO was added to the basal compartment (Fig. 6b; SEMA6AECTO basal vs mock basal: 30 min *p* = 0.9982, 60 min *p* = 0.0653, 120 min *p* = 0.0314, 180 min *p* = 0.0078, 240 min *p* = 0.0508; Two-way ANOVA followed by Tukey’s post-hoc test). This indicates that parenchymal, but not luminal, SEMA6A acts as a potent inducer of endothelial permeability *in vitro*. ### SEMA6A regulates vascular permeability via Plexin-A2 As a first step to determine the signaling pathway by which SEMA6A controls EC permeability, the expression of putative SEMA6A receptors in ECs was examined in two publicly available scRNA-seq databases. As Plexin-A2 and Plexin-A4 are canonical SEMA6A receptors28,30,31 (Supplementary Fig. 2a), our analysis focused on class A plexins. For analysis of the Tabula Muris dataset, the FACS/SMARTseq2-based dataset was used and the cluster containing *Emcn*+ ECs was selected (Supplementary Fig. 4c). *Plxna2*, but not *Plxna1*, *Plxna3* or *Plxna4*, mRNA was present at high levels and frequency (Supplementary Fig. 5e,f). Analysis of the EC atlas53, obtained by droplet-based sequencing of FACS-isolated PECAM1+ ECs, allowed separation of two discrete EC clusters thanks to the higher number of EC transcriptomes obtained despite a lower sequencing depth compared to Tabula Muris. One cluster was composed of arterial, capillary and venous ECs, whereas a smaller subset comprised *Plvap*+ ECs, likely representing fenestrated capillary ECs (fECs) such as those in the ME (Fig. 6c). *Plxna1*, *Plxna3* and *Plxna4* transcripts were not detected in any EC subpopulation (Fig. 6d). In contrast, *Plxna2* was detected in the different EC subtypes and was enriched in fECs (Fig. 6c,d). We validated these results by RT-qPCR in HUVECs and confirmed *Plxna2* expression (Supplementary Fig. 5g). Together, these observations identify Plexin-A2 as the likely binding partner of SEMA6A on ECs. In line with the observation that SEMA6A induced vascular permeability only when exposed to EC monolayers from the abluminal side (Fig. 6b), Plexin-A2 did not colocalize with the luminal/apical membrane marker Intercellular adhesion molecule 2 (ICAM-2)54, but rather localized to the basal side of the monolayer (Fig. 6e). Finally, to establish a functional requirement for Plexin-A2 in SEMA6A-induced vascular permeability, *PLXNA2* expression was knocked down in HUVECs using *sh*RNA approaches (Supplementary Fig. 5h). As expected, basal SEMA6AECTO exposure significantly and steadily reduced TEER compared to mock treatment in HUVEC monolayers infected with scrambled control short hairpins (*shScr*) (Fig. 6f; *shScr*, blue symbols: SEMA6AECTO basal vs. mock basal, 30 min *p* = 0.0001, 60-240 min *p* < 0.0001; Two-way ANOVA followed by Tukey’s post-hoc test). Instead, treatment with SEMA6AECTO from the basal compartment in HUVECs with decreased levels of Plexin-A2 no longer affected TEER (Fig. 6f; *shPLXNA2* (B), green symbols: SEMA6AECTO basal vs. mock basal, 30 min *p* = 0.0517, 60 min *p* = 0.9994, 120 min *p* = 0.3132, 180 min *p* = 0.8817, 240 min *p* = 0.1738; Two-way ANOVA followed by Tukey’s post-hoc test) or caused only a slight and delayed TEER reduction (Fig. 6f; *shPLXNA2* (C), orange symbols: SEMA6AECTO basal vs. mock basal, 30 min *p* > 0.9999, 60 min *p* = 0.5015, 120 min *p* = 0.7981, 180 min *p* = 0.0149, 240 min *p* = 0.0089; Two-way ANOVA followed by Tukey’s post-hoc test). Still, SEMA6AECTO efficacy in reducing TEER was significantly reduced in both *shPLXNA2*-infected HUVECs compared to *shScr* HUVECs at all time points (Fig. 6f, hashtags; SEMA6AECTO *shScr* vs SEMA6AECTO *shPLXNA2* (B): 30 min *p* = 0.2205, 60-180 min *p* < 0.0001, 240 min *p* = 0.0051; SEMA6AECTO *shScr* vs SEMA6AECTO *shPLXNA2* (C): 30 min *p* = 0.0006, 60-180 min *p* < 0.0001, 240 min *p* = 0.0049; Two-way ANOVA followed by Tukey’s post-hoc test). Together, these experiments show that SEMA6A can regulate barrier opening in EC monolayers via Plexin-A2. ### SEMA6A regulates ME capillary fenestration and vascular permeability in vivo Having found that parenchymal SEMA6A was able to bind PLXNA2 on the abluminal side of EC monolayers *in vitro* to promote endothelial permeability suggests that *in vivo* SEMA6A, which is expressed in close proximity of blood capillary loops (Fig. 5a), may play a role in regulating the ME vascular barrier. In agreement with this hypothesis, PECAM1+ capillary loops in the ME expressed Plexin-A2 (Fig. 7a). Further, transcardial perfusion of Evans Blue (EB) dye, to mimic peripheral signals diffusing through the ME barrier *in vivo*, led to EB extravasation in the ME of control mice. In contrast, EB leakage was significantly less in *Sema6a-/-*MEs (Fig. 7b, c; *p* = 0.0212; two-tailed unpaired Student’s t test). To assess whether this reduction in ME barrier permeability was associated with structural changes, the number of Plasmalemma Vesicle Associated Protein (PLVAP)+ fenestrated capillary loops was quantified, as previously reported55–57. Notably, *Sema6a-/-* mice had significantly fewer PLVAP+ loops invading the ME parenchyma, as compared to WT mice (Fig. 7d, e; *p* = 0.023; two-tailed unpaired Student’s t test). Double staining for PLVAP and PECAM1 showed reduced PLVAP immunostaining and a decrease in PECAM1+ capillary loops in the ME of *Sema6a-/-* mice (Supplementary Fig. 6a). This indicates that reduced ME barrier permeability caused by SEMA6A loss is accompanied by defective structural remodeling of ME capillaries. The morphology of Vimentin+ β-tanycytes, known to modulate NGE unit plasticity and consequently GnRH secretion17, was unchanged in *Sema6a-/-* mice (Supplementary Fig. 6b).  [Figure 7.](http://medrxiv.org/content/early/2023/10/26/2023.10.24.23297449/F7) Figure 7. SEMA6A loss prevents ME barrier opening and induces structural rearrangements in the ME capillary bed. **(a)** Coronal sections of adult ME were immunolabelled for Plexin-A2 (green) together with PECAM1 (red) to label blood vessels. The white box indicates the area shown at higher magnification on the right of the corresponding panel together with single channels for Plexin-A2 and PECAM1. Solid arrowheads indicate the expression of Plexin-A2 on capillaries of the basal portion of the ME. **(b)** Coronal sections of brains from EB-injected females of the indicated genotypes. White boxes indicate areas shown at higher magnification on the right of the corresponding panel. Δ indicates areas of no EB leakage in *Sema6a-/-* ME parenchyma. **(c)** Quantification of EB diffusion in the ME parenchyma of female brains (n = 5 per group; *p* = 0.0212). **(d)** Coronal sections of female brains with the indicated genotypes were immunolabelled for PLVAP (green) to reveal capillary loops invading the ME parenchyma. White boxes indicate areas shown at higher magnification on the right of the corresponding panel together with a single channel for PLVAP. White arrows or Δ indicate the normal or decreased presence of capillary loops in *Sema6a+/+* and *Sema6a-/-*ME parenchyma, respectively. **(e)** Quantification of capillary loops in the ME of female brains (n = 3 per group; *p* = 0.0230). All sections were counterstained with DAPI. * *p* < 0.05 after two-tailed unpaired Student’s t test. Abbreviations: 3v, third ventricle; EB, Evans Blue. Scale bars: 100 μm (**b,d** low magnifications), 50 μm (**a**, low magnification; **b** and **d**, high magnifications), 25 μm (**a**, high magnification). Data are presented as mean ± SD. Source data are provided as a Source Data file. Together, these findings show that SEMA6A is required *in vivo* to maintain ME barrier permeability and the formation of fenestrated capillary loops in the ME. ### Exome sequencing of patients with DP identifies a loss-of-function SEMA6A variant in six individuals In view of the observation that *Sema6a-/-* mice present with delayed pubertal timing, we interrogated exome sequencing data from a large cohort of patients with self-limited DP and identified a potentially pathogenic, variant in *SEMA6A* (ENSG00000092421, NCBI gene ID: 57556). This variant (NM_020796; c.1268T>C; p.I423T) was carried by six affected individuals from our cohort, including four members of one family, and all individuals were heterozygous for this variant (Fig. 8a). The presence of this *SEMA6A* variant in probands and family members was confirmed by Sanger sequencing (Supplementary Fig. 7). The affected subjects did not carry any other predicted pathogenic variants in known HH/DP causing genes nor other pituitary hormone deficiencies were excluded.  [Figure 8.](http://medrxiv.org/content/early/2023/10/26/2023.10.24.23297449/F8) Figure 8. Exome sequencing identifies a *SEMA6A* variant with effects on protein stability in DP patients. **(a)** Genealogy of family presenting DP. Symbols: DP phenotype (black), unknown phenotype (grey); unaffected individual (white). Horizontal black lines indicate the presence of heterozygous p.I423T variant identified by WES and verified by Sanger sequencing. US signifies that DNA was unavailable for sequencing. **(b)** Chromosome (Chr) position, nucleotide substitution (nt sub), amino acid substitution (aa sub), and bioinformatic predictions of the p.I423T *SEMA6A* variant according to Poly-Phen, SIFT, REVEL and CADD software. **(c)** Partial protein sequence alignment of vertebrate SEMA6A orthologs shows that I423 conservation in mouse and zebrafish (left). Genomic evolutionary rate profiling of sequence constraint for SEMA6A variant using GERP++ analysis (RS score = 6.0, right). **(d)** Side (left) and top (right) 3D views of modeled human SEMA6A monomer. Red asterisks indicate the I423 residue within the central cavity of SEMA domain 6th beta-propeller blade (d). **(e)** High magnification images of I423 (left) and T423 (right) residues and their interactions. **(f)** Lysates from untransfected (UT) COS-7 cells and cells transfected with an empty control expression vector (CTRL) or vectors encoding c-myc tagged human SEMA6AWT or mutant SEMA6AI423T were immunoblotted for c-myc to recognize human (SEMA6A: 120 kDa). GAPDH (37 kDa) was used as a loading control. **(g)** Quantification of relative abundance of SEMA6AI423T protein compared to SEMA6AWT after normalization with GAPDH (n = 3 independent experiments; *p* < 0.0001). **(h)** COS-7 cells transfected with vectors encoding c-myc tagged SEMA6AWT (left) or SEMA6AI423T (right) were immunolabeled for c-myc (green) and F-actin (red). Solid arrowheads indicate cell surface SEMA6AWT expression (left), whereas empty arrowheads indicate lack of surface SEMA6AI423T but intracellular and perinuclear expression (right). **(i)** COS-7 cells co-transfected with vectors encoding SEMA6AI423T and ER-emerald, visualizing ER (green), were immunolabeled for c-myc (red) and counterstained with DAPI. Solid arrowheads indicate examples of SEMA6AI423T ER retention. \**\*|\* *p* < 0.0001 after two-tailed unpaired Student’s t test. Abbreviations: ER, endoplasmic reticulum. Scale bars: 50 μm (h, i). Data are presented as mean ± SD. Source data are provided as a Source Data file. All the individuals (both males and females) carrying the *SEMA6A* variant showed classical clinical and biochemical features of self-limited DP, with delayed onset of Tanner stage 2 and delayed peak height velocity. One boy from this family (proband A) presented between 11 and 15 years with testicular volumes of 3.7 mL, bone age delayed by 3.6 years and short stature (height SD score −2.75) (see Table 1 for detailed auxological data). He had normal olfaction and no evidence of chronic illness. His short stature was related to pubertal delay, as height SDS at both 7 years and at adult height was −1.3, in keeping with his mid-parental target height. His delayed puberty was inherited from his affected relative, who had menarche aged between 11 and 15 years; three of his relatives also had delayed onset of puberty and growth spurt. Both the proband and his relatives had spontaneous pubertal development without testosterone therapy on follow up, thus excluding HH. The second male proband (B) presented between 11 and 15 years with delayed pubertal onset (testicular volumes 3.5 mL and lack of secondary sexual characteristics). He had a family history of delayed puberty with affected relatives. Proband C, also male, was first seen between age 15 and 17 years with testicular volumes of 4.5 mL and bone age delay of 2.8 years. Both probands B and C had progressed through spontaneous puberty by the age of 18-19 years, confirming a self-limited DP phenotype. DNA was not available from family members of these two probands. View this table: [Table 1](http://medrxiv.org/content/early/2023/10/26/2023.10.24.23297449/T1) Table 1 Clinical data of auxological measurements and biochemical assessments obtained from childhood to adult age for one of the probands from the index family. ### The SEMA6AI423T variant is predicted to affect protein stability The identified *SEMA6A* variant was not found in the gnomAD browser (accessed 28.06.2023). The substitution of the non-polar amino acid isoleucine (Ile, I) with a neutral polar amino acid, threonine (Thr, T), is predicted to be damaging to protein function or disease-causing by 3/3 bioinformatic tools (Poly-Phen, SIFT and REVEL) with a combined annotation dependent depletion (CADD) score of 31 (Fig. 8b). The affected I423 residue lies in the SEMA domain and is highly conserved between species (Fig. 8c). To further predict the possible damaging impact of the I423T variant on protein stability, an *in silico* analysis was performed starting from a CryoEM-resolved structure of SEMA6A in complex with the *C. sordellii* lethal toxin, TcsL (PDB ID: 6WTS)58. Within the SEMA domain, I423 lies close to the central cavity of the beta-sheet structure of the 6th beta-propeller blade (Fig. 8d), which has no direct contact to the SEMA6A dimerization or SEMA6A::Plexin-A2 interaction surfaces58. Since I423 forms hydrophobic interactions with some of the surrounding residues (e.g., A298, T300, L312, A313, T314), an Ile to Thr substitution would reduce the size of the side chain and its hydrophobicity, thereby reducing the extent of these interactions (hydrophobicity from 4.5 for Ile to −0.7 for Thr, according to previous reports59 (Fig. 8e). Additionally, substitution with a smaller residue would make this local niche part of the large solvent accessible cavity of the SEMA domain, leading to a possible decrease in protein stability. Therefore, the impact of Ile to Thr substitution on protein stability of both SEMA6A chains was evaluated *in silico* with the BioLuminate *Residue Scan* tool. This analysis showed that the stability of the mutated SEMA6A homodimer was substantially decreased compared to the WT SEMA6A homodimer (ΔStability: 16.22 ± 0.14 kcal/mol for each single SEMA6A protomer), supporting an effect of the identified variant on SEMA6A stability. ### The SEMA6AI423T variant affects protein levels and localization in vitro To functionally validate the *in silico* modelling predictions, we generated myc-tagged vectors encoding WT and mutated forms of human *SEMA6A* (*SEMA6AWT* and *SEMA6AI423T* respectively) followed by expression in COS-7 cells. Then, protein synthesis and cellular localization were analyzed as previously described for other semaphorins at 48 hours post-transfection12,18. Western blotting analysis revealed that both SEMA6AWT and SEMA6AI423T proteins were successfully synthetized, although expression of the mutant protein was significantly reduced as compared to SEMA6AWT (Fig. 8f, g; *p* < 0.0001; two-tailed unpaired Student’s t test). This suggests that the mutation affects protein expression and/or stability. To examine whether the variant could also affect protein localization, we performed immunofluorescence using an anti-c-myc antibody to visualize SEMA6AWT or SEMA6AI423T in transfected cells (Fig. 8h). SEMA6AWT protein was detected at the cell membrane, as previously reported60, while SEMA6AI423T protein mainly localized in the perinuclear region. To understand the subcellular localization of mutant SEMA6A, COS-7 cells were co-transfected with SEMA6AI423T and a fluorescent endoplasmic reticulum (ER)-marker (mEmerald-ER-3). High level of signal colocalization suggests that SEMA6AI423T was retained in the ER (Fig. 8i). Together, the *in vitro* results strongly support our *in silico* predictions and indicate that the I423T variant might affect SEMA6A stability and impair intracellular protein trafficking. The consequently defective membrane localization would hamper any paracellular interaction with SEMA6A receptors, suggesting a loss-of-function effect of the SEMA6AI423T variant in patients. ## DISCUSSION Our study provides so far undescribed insights into the mechanisms that regulate puberty onset by implicating SEMA6A in the control of ME homeostasis and GnRH neuron plasticity and by identifying human patients with delayed onset of puberty carrying a pathogenic *SEMA6A* variant. Similar to other members of the semaphorin family, SEMA6A regulates axon guidance during brain development25,61. In line with this, we detected strong SEMA6A expression in nasal axons of developing mouse and human embryos (Fig. 1). However, SEMA6A loss did not affect GnRH axon pathfinding or GnRH neuron migration (Fig. 2). These results indicate a dispensable role for SEMA6A in these processes, as already observed for other semaphorins (e.g. SEMA3F)11. In addition to their key roles in axon guidance and neuronal migration at embryonic stages, semaphorins are expressed in the postnatal hypothalamus where they control GnRH neuron structural plasticity15–17. We found that SEMA6A is stably expressed in the ME parenchyma of perinatal (E18.5) and postnatal mice (P24 and P60) (Figs. 1,3), a region that undergoes several cellular rearrangements to regulate neurohormonal homeostasis39. Intriguingly, neither SEMA6A transcript nor protein were detected in neuroendocrine neurons, β-tanycytes or astrocytes (Fig. 5 and Supplementary Fig. 4). These cells are known to play key roles in GnRH neuron plasticity43,44. Instead, we identified maturing OLs as the source of SEMA6A in the ME (Fig. 5 and Supplementary Fig. 4). Thus, in addition to its previously described role in OL differentiation45, SEMA6A derived from OLs in the ME serves as a signal for regulating ME vascular permeability and GnRH neuron plasticity. *Sema6a-/-* mice displayed decreased GnRH axon innervation of the ME starting from E18.5 (Fig. 2), when the first GnRH+ fibers project to the ME35, until adult stages (P60) (Fig. 4), when GnRH neuron plasticity must be tightly regulated to guarantee proper GnRH secretion1. In line with these observations, both male and female *Sema6a*-/- mice showed delayed puberty onset and sexual maturation defects (Fig. 4). The decrease in GnRH neuron innervation of the ME was not accompanied by a reduction in the number of GnRH neurons reaching the MPOA (Fig. 2). This excludes a decrease in neuronal number as the cause for the reduced density of GnRH-positive axons in the ME. Moreover, our data also argue against a direct effect of SEMA6A on GnRH neuron axons, as no expression of the SEMA6A receptors Plexin-A2 and Plexin-A4 was detected in MPOA-resident GnRH neurons (Supplementary Fig. 2). The expression of SEMA6A in the external portion of the ME adjacent to fenestrated blood capillary loops invading the brain parenchyma, suggested a possible role in the NGE unit. This is a fundamental structure formed by a dynamic network of capillaries, β-tanycytes and neuroendocrine projections5,39. Interestingly, previous studies have described roles for semaphorins, such as SEMA3A and SEMA7A, in the control of GnRH neuron plasticity by modulating GnRH neurite extension15 or tanycyte endfeet sprouting17, respectively. However, evidence for direct control of ME vascular permeability by semaphorins had not been reported. Here, by combining *in vitro* vascular permeability experiments and *in vivo* functional studies, we unveiled a thus far undiscovered role for SEMA6A in ME vascular barrier homeostasis. Specifically, we found that a secreted form of SEMA6A increases permeability of cultured EC monolayers (Fig. 6 and Supplementary Fig. 5), and that *Sema6a-/-* mice display a less permissive ME vascular barrier. This was associated with fewer fenestrated loop capillaries, but normal β-tanycyte morphology (Fig. 7 and Supplementary Fig. 6). These findings strongly support a direct role for SEMA6A in the regulation of ME vascular permeability. The effect of SEMA6A on ECs is likely to be mediated by its canonical receptor Plexin-A2. Plexin-A2 is expressed by fenestrated capillaries of the ME *in vivo* (Fig. 7) and *PLXNA2* knockdown attenuated SEMA6A-induced vascular permeability *in vitro* (Fig. 6). These data show that SEMA6A can act via Plexin-A2 to modulate vascular permeability. Still, we cannot exclude that SEMA6A also interferes with Vascular endothelial growth factor receptor 2 (VEGFR-2) signaling, as previously reported in HUVECs29,62, possibly through the binding with other SEMA6A interactors such as members of the Ena/VASP protein family33,63. VEGFR-2 is expressed by fenestrated ECs and has previously been shown to sustain ME vascular permeability through VEGF-A signaling56. Our findings provide support for a key role of SEMA6A in the regulation of NGE unit. The NGE unit regulates the outflux of neurohormones, including GnRH, to target tissues, as well as the influx of molecules from peripheral tissues. Therefore the ability of SEMA6A to regulate ME vascular barrier could modulate trafficking of molecules leading to changes in factor secretion or in uptake of circulating peptides39,64. In particular, decreased ME vascular permeability in the absence of SEMA6A (Fig. 7) may reduce the uptake of circulating factors that regulate GnRH neurite extension and/or GnRH secretion. For example, blood-borne Insulin-like growth factor 1 (IGF-1) and Fibroblast growth factor 21 (FGF-21) can induce the release of GnRH65 and stimulate GnRH neurite outgrowth with subsequent peptide release40, respectively (Figure 9). Our observation that loss of SEMA6A reduces the diffusion of EB dye from the periphery through the ME barrier supports this model. Of note, mice lacking IGF-1 or FGF-21 receptors also exhibit DP phenotypes like *Sema6a-/-*mice 40,66.  [Figure 9.](http://medrxiv.org/content/early/2023/10/26/2023.10.24.23297449/F9) Figure 9. SEMA6A modulates GnRH neuron homeostasis by tuning ME vascular permeability. Schematic drawing representing the proposed cellular mechanisms through which SEMA6A is expressed by OLs to modulate GnRH homeostasis. In the ME of WT mice, SEMA6A (green dots) is produced by OLs (green cells) and distributed close to fenestrated ECs (red cells). In our model, SEMA6A acts on fenestrated ECs expressing Plexin-A2 (blue receptors) to increase the number of capillary loops in the brain parenchyma and maintain a permissive vascular barrier (1). Such permissive vascular barrier may ensure optimal GnRH peptide (gray dots) release in ME blood vessels either directly (2) or indirectly, via entry of molecules from the periphery (e.g. IGF-1, FGF-21; black triangles) in the ME (3) where they can act on GnRH neurons (in gray), by modulating the extension of axon terminals towards ME and ECs (4). This mechanism will ensure a normal release of GnRH peptide (5) (gray dots) and therefore a physiological puberty onset. Conversely, in *Sema6a*-/-mice, the lack of SEMA6A production by OLs (grey cells) in the basal part of the ME leads to fewer capillary loops and reduced vascular permeability. Consequently, direct GnRH release to the pituitary or entry of peripheral signals that normally sustain axon terminal elongation and GnRH secretion do not occur, resulting in a less innervated ME and ultimately in delayed puberty onset. Created with BioRender.com. The identification of a predicted pathogenic variant in *SEMA6A* in 6 individuals with a DP phenotype serves to translate our experimental observations to the clinic and provides a mechanistic link between SEMA6A deficiency and GnRH functionality in humans. DP is a highly heritable condition affecting up to 2% of the population8,67. Yet, despite recent progress, the genetic determinants of DP remain largely unknown10. The analysis of exome sequencing from 100 probands from a self-limited DP cohort led to the identification of this novel missense variant in *SEMA6A*, segregating with the DP trait in 4 individuals from the same family and two further probands from unrelated families. Interestingly, our patients do not exhibit other pituitary hormone deficiencies, strongly supporting a selective effect of SEMA6A on the regulation of the HPG axis. Our work provides strong support for a deleterious effect of this *SEMA6AI423T* variant on protein stability and localization. The DP patients in our cohort are heterozygous carriers of the *SEMA6AI423T* variant, whilst in the mouse model *Sema6a-/-* but not *Sema6a+/-* mice displayed DP. However, previous work on several related genes has demonstrated that even when patients carry heterozygous pathogenic mutants, only full knock-out, but not heterozygous, mice targeting the corresponding genes exhibit a pathological reproductive phenotype (e.g. *Sema3a*, *Sema7a*, *Plxna1*)11,14,68. The same phenomenon has been demonstrated for other developmental genes69. This would suggest that patient-derived mutations can exhibit a dominant effect on protein function, such as the altered stability and localization of the *SEMA6AI423T* variant. In contrast, gene targeting strategies in mice do not usually affect the stability of the protein encoded by the *WT* allele in heterozygous animals. In conclusion, our results uncover an additional role for the axon guidance cue SEMA6A in the control of puberty onset and sexual maturation and help to understand the pathogenic mechanisms underlying a novel *SEMA6A* variant in patients affected by DP. Specifically, these data unveil an unprecedented mechanism through which OL-produced SEMA6A acts on Plexin-A2+ ECs to regulate ME homeostasis and, as a consequence, GnRH neuron innervation. Our findings provide a framework for further understanding the genetic and molecular mechanisms that underlie puberty onset as well as for guiding diagnostic and therapeutic approaches in pubertal delay. Further, the outlined mechanism of disrupted ME vascular permeability could potentially have broader implications on various neuroendocrine functions beyond pubertal maturation. ## METHODS ### Ethics Human embryonic material was provided by The MRC-Wellcome Trust Human Developmental Biology Resource (HDBR, [http://hdbr.org](http://hdbr.org)) that has ethical approval to collect and distribute embryonic and fetal material to individual researchers, upon project approval (n° 200343). The human study protocol was approved by the Ethics Committee for Pediatrics, Adolescent Medicine and Psychiatry, Hospital District of Helsinki and Uusimaa (570/E7/2003). Ethical approval for human studies was granted by the UK London-Chelsea NRES committee (13/LO/0257). All participants provided written informed consent prior to study participation. The study was conducted in accordance with the guidelines of The Declaration of Helsinki. Participants did not receive compensation. ### Animals and tissue preparation *HA-Sema6afl/fl* mice, which harbor a 3xHA tag in exon 2 of *Sema6a* and lox sites surrounding exon 3, were a kind gift of Alex Kolodkin (Johns Hopkins University School of Medicine, Baltimore, USA). Germline deletion of *Sema6a* exon 3 was performed by crossing these mice with *EIIa-Cre* mice (JAX stock #003724)70 to create *Sema6a*-null mice (here referred to as *Sema6a-/-*) that were maintained on a *C57Bl/6J* background (JAX strain, strain code 632, Charles River). Genotyping was performed by PCR using DNA from mouse tail or ear tissue as a template and primers flanking the HA tag sequence (fw ‘5-CCTTGTCACACATGCAGTTG-3’ and rev ‘5- TGAAGTAGAGAAGCA-3’), and by primers controlling for germline deletion of exon 3 (fw ‘5-GCTCTTTGGCAGGGTGAAGTA-3’ and rev ‘5-GTCCAAGTCCACGTGACCC-3’). For all experiments, mice were housed socially with chow diet (CRM (E) SQC 811437, Special Diet Services) and water ad libitum, on a 12 hour light/dark cycle, humidity of 45-65% and a temperature between 20-24 degrees. The morning of vaginal plug observation was defined as embryonic day 0.5 (E0.5). The day of birth was defined as postnatal day 0 (P0). Timed pregnant females were sacrificed by cervical dislocation to collect embryos at E12.5, E14.5 and E18.5. Embryo heads were fixed for 3-6 h in 4% PFA at 4°C. Adult mice were euthanized by overdosing Euthanimal (200 mg/mL, Alfasan) intra-peritoneal (i.p.), and transcardially perfused with ice-cold 4% PFA. Subsequently, brains were isolated and post-fixed in 4% PFA overnight at 4°C. All samples for cryosectioning were cryoprotected overnight in 30% sucrose prior to OCT-embedding, whereas samples for paraffin embedding were dehydrated by immersion through crescent alcohol series and xylene. All *in vivo* experiments were conducted by certified personnel in compliance with the Animal Ethics Committee of Utrecht University (Dierexperimenten Ethische Commissie) (CCD licence: AVD115002016532 and AVD11500202114777) and the Animal Welfare Body of the University of Milan, and in agreement with Dutch laws (Wet op de Dierproeven, 1996; revised 2014), the Italian Minister of Health and European regulations (Guideline 86/609/EEC; Directive 2010/63/EU). Animal welfare was checked daily. ### Immunofluorescence labelling 20 μm cryostat sections of PFA-fixed embryos were incubated with serum-free protein block (DAKO) after permeabilization of sections with 0.1% TritonX-100. We used as primary antibodies mouse anti-nTUBB3 (1:500, clone Tuj1; Covance, cat. MMS-435P), rabbit anti-Peripherin (1:100; Millipore, cat. AB1530), rabbit anti-GnRH, (1:400; Immunostar, cat. 20075)14, rabbit anti-BLBP (1:200; Millipore, cat. ABN14)71, goat anti-SEMA6A (1:200; R&D Systems, cat. AF1615)29, rabbit anti-GFAP (1:500; Dako, cat. Z0334)72, rabbit anti-Synaptophysin (1:200; Synaptic System, cat. 101002), mouse anti-Vimentin (1:200; Santa Cruz Biotechnology, cat. sc-373717)73, rabbit anti-Oligo2 (1:300; Millipore, cat. AB9610)74, rabbit anti-Gst P1 (1:200; MBL, cat. 312)75, rat anti-PDGFR-α (1:100, clone APA5; Invitrogen, cat. 14-1401-82)76, goat anti-Plexin-A2 (R&D Systems, cat. AF5486)77, rabbit anti-Plexin-A4 (Novus Biologicals, cat. NBP1-85128)78, mouse anti-Plexin-A2 (1:100, clone A-2; Santa Cruz Biotechnology, cat. sc- 393939), rabbit anti-PECAM1 (1:200; Abcam, cat. ab28364)79 and rat anti- PECAM1 (1:200; BioLegend, cat. 102502). To visualize GnRH neurons in the MPOA when another antibody raised in rabbit was used, we employed goat anti-Plexin-D1 antibody (1:200; R&D Systems, AF4160) as previously done80. For blood capillary loop analysis, 16 μm cryostat sections of PFA-fixed brains at ME level were first boiled in an antigen retrieval solution (0.1M citric acid, 1M sodium citrate) for 20 min and then incubated in a blocking solution (PBS 0.1% TritonX-100, 0.05% Tween-20, 2.5% donkey serum, 5% BSA) and primary rat anti-PLVAP antibody (1:100, clone Meca32; BD Biosciences, cat. 550563)57. Secondary antibodies used were 488-, Cy3- or 647-conjugated donkey anti- goat, anti-mouse, anti-rabbit and anti-rat Fab fragments (1:200; Jackson Immunoresearch). In some experiments, blood vessels were labelled using biotinylated IB4 (1:100, Vector laboratories) followed by Cy3-conjugated streptavidin. Nuclei were counterstained with DAPI (1:10000; Sigma Aldrich). ### Immunohistochemistry 20 µm cryostat sections of PFA-fixed samples were incubated with hydrogen peroxide to quench endogenous peroxidase activity, and sequentially incubated with 10% normal goat serum in PBS or serum-free blocking solution (Agilent, cat. X0909) and then immunostained with rabbit anti- GnRH (1:1000; Immunostar), previously validated to recognize both the pre- hormone and the processed hormone81. 8 µm microtome sections of PFA-fixed paraffin-embedded testes were similarly processed and incubated with rabbit anti-CYP17A1 antibody (1:200; Proteintech, cat. 14447-1-AP)18. Sections were then incubated with anti-rabbit biotinylated antibody (1:400; Vector Laboratories) and developed with the ABC kit (Vector Laboratories) and 3,3- diaminobenzidine (DAB; Sigma). To determine the total number of GnRH neurons at E14.5, E18.5 and postnatal stages, 20 µm coronal sections through each entire head (E14.5), FB (E18.5) or MPOA (postnatal) were immunolabelled for GnRH and all GnRH-positive cells quantified as previously reported14. To compare the abundance of GnRH-positive neurites in the ME, we measured the pixel intensity of GnRH staining as previously described82. To compare the abundance of Leydig cells in testes, we quantified for each section the percentage of CYP17A1+ area using ImageJ software (v1.52a, NIH). ### Hematoxylin & eosin staining and follicle quantification 8 µm microtome sections of PFA-fixed paraffin-embedded ovaries were rehydrated with xylene/decrescent alcohol series, stained with Mayer’s hematoxylin (Bio- Optica) for 2 min and Eosin Y alcoholic solution (Bio-Optica) for 3 min and finally dehydrated with crescent alcohol series/xylene. Number and subtype of ovarian follicles were quantified as previously described83. ### Determination of male balanopreputial separation, female vaginal opening and first estrous *Sema6a-/-*, *Sema6A*+/-, and *Sema6a+/+*male and female mice were weighed and checked daily for BPS and VO37,38, respectively, after weaning from P22 onwards. When female mice showed VO, vaginal smears were obtained every 24 h to assess the start of the estrous cycle by gently pipetting and resuspending 30 µL saline solution (0.9% NaCl) in the VO using a P200 pipette. The collected cells were transferred to a glass microscope slide and air-dried. Slides were then treated with 70% EtOH for 3 min, followed by 0.1% Cresyl Violet stain in 50% EtOH for 1 min, and briefly washed with 100% EtOH. Subsequently, smears were evaluated using light microscopy and the relative amount of nucleated, anucleated epithelial cells and leukocytes present was used to determine the phase of the estrous cycle (Supplementary fig. 3c)37. Gross absence of leukocytes and presence of multitudes of anucleated epithelial cells was considered as estrous. Tissues for analysis were collected from male and female mice at the day of BPS and FE, respectively. Mice were euthanized by cervical dislocation, brains and gonads were isolated and fixed overnight in 4% PFA at 4°C, and pituitaries were snap-frozen and stored at −80°C. Males were collected from 13 litters, females were collected from 14 litters. ### Evans Blue (EB) *in vivo* permeability assay At P30*, Sema6a-/-*, *Sema6A*+/-, and *Sema6a+/+* female mice were euthanized by overdose of Euthanimal i.p. (200 mg/mL, Alfasan). Subsequently, animals were transcardially perfused with 20 mL sterile-filtered 68 mg/mL EB dye (Sigma-Aldrich) dissolved in 4% BSA in PBS, followed by perfusion with 25 mL PBS. Perfusion speed was kept at 3 mL/min. Following perfusion, brains were isolated and cryoprotected by incubation in ascending concentrations of sucrose (10%-20%-30%) in PBS for 30 min each at 4°C. Brains were then embedded in OCT and snap-frozen in dry-ice. *Sema6A*+/- and *Sema6a+/+*females who did not yet show VO at P30 were excluded from analysis. *WT* and *Sema6a* mutant animals were collected from 4 litters. For analysis, 30 µm cryosections were directly imaged under laser scanning confocal microscope to visualize EB staining excited at 640 nm. EB diffusion in the ME parenchyma was assessed by quantifying the background-normalized mean gray values of three equally sized fields per section, randomly placed in the ME region, using ImageJ software (v1.52a, NIH). ### RNA extraction and RT-qPCR Total RNA was obtained from HUVEC cell pellets and mouse pituitaries using Trizol/NucleoZOL (Macherey-Nagel). 500 ng of RNA were then retrotranscribed using All-in-one 5x RT MasterMix kit (Abm, cat. G592) and random primers. RT-qPCR was performed on CFXOplus96 thermocycler (Bio-Rad) using Master Mix with SYBR Green (Takara Bio) and specific primers (mouse *Lhb*: fw 5’- CAAGAATGGAGAGGCTCCAG-3’, rev 5’-ACTGGGCAGAACTCATTCTCTG-3’; mouse *Gapdh*: fw 5’-CATCCCAGAGCTGAACG-3’, rev 5’- CTGGTCCTCAGTGTAGCC-3’; human *PLXNA1*: fw 5’- GTGGAGAGGTACTATGCAGAC-3’, rev 5’-CATGCTGTTGAA CTGGCTC-3’; human *PLXNA2*: fw 5’-AAGTATAGTGAGGAGCTCATCG-3’, rev 5’- CTGCTCCACCTTATAAGCCA-3’; human *PLXNA3*: fw 5’-GTGCCAGTTTGTAAGGAGAG-3’, rev 5’-CAATGGCAAGTGTGATGGG-3’; human *PLXNA4*: fw 5’-AATGACCGCATTAAGGAGC-3’, rev 5’-AATGGTTAAGAGCGCACTG-3’; human *RPLP0*: fw 5’-CAGATTGGCTACCCAACTGTT-3’, rev 5’-GGGAAGGTGTAATCCGTCTCC-3’). The quantification cycle (ΔCq) value and the ΔΔCq were calculated relative to control samples using Cq threshold values normalized to the housekeeping gene *RPLP0*/*Gapdh*. ### scRNA-seq analysis Raw data for adult brain scRNA-seq datasets from Tabula Muris and EC atlas were obtained from NCBI gene Expression Omnibus (GSE109774; [https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE109774](https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE109774)) and EBI ArrayExpress (E-MTAB-8077; [https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-8077?query=E-MTAB-8077](https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-8077?query=E-MTAB-8077)), respectively, and analyzed with RStudio v.1.3.1056. The raw gene expression matrices (UMI counts per gene per cell) were filtered, normalized and clustered using the R package Seurat v.3.2.0. Cells containing less than 200 feature counts and genes detected in less than 3 cells were removed. Downstream analysis included data normalization (“LogNormalize” method and scale factor of 10,000) and variable gene detection (“vst” selection method, returning 2,000 features per dataset). The principal components analysis (PCA) was performed on variable genes, and the optimal number of principal components, PCs, for each sample was chosen using the elbow plot. The selected PCs were used for Louvain graph-based clustering at a resolution of 0.3 followed by manual curation based on marker genes to obtain the final seven and four clusters from Tabula Muris and EC atlas brain datasets, respectively. Uniform manifold approximation and projection (UMAP) was chosen as a non-linear dimensionality reduction method, and each relevant gene was then examined using the *FeaturePlot* and *VlnPlot* functions. Cluster cell identity was assigned by manual annotation based on known marker genes. Clustering was then manually curated based on marker genes to obtain the final seven and four clusters from Tabula Muris and EC atlas brain datasets, respectively, followed by subset selection in the Tabula Muris dataset of the cluster containing *Emcn*+ ECs. ### Cell culture, transfection and lentivirus infection COS-7 cells (American Type Culture Collection, cat. CRL-1651) were grown as a monolayer at 37°C in a humidified CO2 incubator in complete DMEM (Euroclone) and supplemented with 10% fetal bovine serum (FBS; Gibco). Subconfluent cells were harvested by trypsinization and cultured in 57 cm2 dishes. For transfection, COS-7 cells (at 80% confluence) were grown in culture plates in complete culture medium for 24 h and incubated for 48 h with the selected expression vector (1 μg/mL) in the presence of Lipofectamine 3000 (Invitrogen) according to the manufacturer’s instructions. To visualize Endoplasmic Reticulum (ER), the mEmerald-ER-3 vector (gift from Michael Davidson; Addgene plasmid #54082) was co-transfected with the *SEMA6A*I423T vector (ratio 1:4). The conditioned media of COS-7 cells either mock transfected or expressing SEMA6AECTO were collected after 48 h and stored at −80°C. HUVECs from pooled donors (Promocell, cat. FB60C12203) were cultured on fibronectin (Promocell) or rat collagen I coated vessels (Corning, cat. DLW354236) in EGM-2 complete medium (Promocell). To generate stable *PLXNA2* knockdown cells, HUVECs were separately infected with a set of lentiviral particles expressing four different *sh*RNAs (A-D) against human *PLXNA2* (Origene, cat. TL302405V) with a calculated MOI of 10; cells infected with lentiviral particles expressing scramble *sh*RNA were used as controls. At 48 hours post infection, cells were selected with 1 μg/mL puromycin. Knockdown efficiency was evaluated by RT-qPCR. ### Plasmids The *pAPtag-5* vector (6.6 kb) (GenHunter, cat. Q202) was chosen for the cloning of the shorter isoform of human *SEMA6A* gene (NM_020796.5), encoding for a protein of 1030 amino acids (expected molecular weight 120 kDa). Two different vectors have been prepared with human *SEMA6A*WT and *SEMA6A*I423T genes inserted between the NheI and the XhoI restriction enzyme sites, in frame with c-myc tag, cutting off the AP protein. The empty vector has been used as a control vector. The *pCAGG-AP-SEMA6A-ecto-Fc-His* vector was employed to produce mouse SEMA6A soluble extracellular domain ranging from amino acid 18 to 688 (SEMA6AECTO)31. ### Isolation and culture of ECs from adult mouse brain mBECs were isolated from 3-5 brains of adult WT C57/Bl6 mice as previously reported52. Briefly, brains were cut along the longitudinal fissure to separate the two hemispheres and the cerebellum was removed from the rest of the brain. The meninges were removed by rolling the individual parts on 3M paper. Brain hemispheres were cut and minced using a scalpel, resuspended in dissection buffer (Optimem, 2% FCS, 2 ng/mL heparin and pen/strep), and passed through a 10 mL serological pipet until no evident clumps were present. The homogenate was centrifuged (300 g, 5 min) and the pellet was digested in dissection buffer containing 1 mg/ml of collagenase/dispase (Sigma Aldrich) and 50 µg/mL of DNase (Sigma Aldrich) on a benchtop shaker (900 RPM) for 45 min at 37°C. To remove the myelin, the pellet was resuspended in 22% (w/v) of bovine serum albumin (BSA) in PBS and centrifuged (1620 g for 15 min at 4°C) for 3 times. The combined pellets were then resuspended in EGM complete medium (Promocell) supplemented with 5µg/mL of puromycin to kill non-ECs and seeded on collagen coated (rat tail collagen I 50 µg/mL in PBS; Corning, cat DLW354236) 0.4 μm pore-size Transwell Permeable Supports (Merck, cat. MHCT24H48). Culture medium was replaced after 48 hours with complete medium and changed every other day until confluence. mBECs purity was evaluated by morphological inspection and immunofluorescence staining for the EC markers PECAM1 and Cadherin-5 as described below. ### Immunofluorescence Cells were fixed with 4% paraformaldehyde (PFA) fixative solution for 15 min at room temperature and washed with PBS. COS-7 cells were blocked and permeabilized with 10% goat serum 0.1% TritonX-100 (all from Sigma Aldrich) solution and incubated with mouse anti c-myc (1:500, clone 9E10; ThermoFisher, cat. 13-2500)18 antibody overnight at 4°C. mBECs were blocked and permeabilized with 5% donkey serum, 5% BSA, 0.1% Triton X-100 (all from Sigma Aldrich) solution for 30 min at room temperature. After blocking, mBECs were incubated with goat anti-Cadherin-5 (1:200; Santa Cruz Biotechnology, cat. sc-6458), rat anti-PECAM1 (1:200; BioLegend, cat. 102502)84, rat anti-ICAM-2 (1:200; clone 3C4; BD Biosciences, cat. 553326)54 and goat anti-Plexin-A2 (1:200; R&D Systems, cat. AF5486) antibodies overnight at 4°C. Secondary antibodies used were 488-, Cy3 or 647-conjugated donkey anti-goat, anti-mouse and anti-rat Fab fragments (1:200; Jackson Immunoresearch). F-actin staining was performed as previously described85. Nuclei were counterstained with DAPI (1:10000; Sigma Aldrich). ### ECs permeability assay HUVECs or mBECs were seeded on fibronectin or collagen I coated 0.4 μm pore-size Transwell Permeable Supports (Merck, cat. MHCT24H48), cultured in complete culture EGM-2 or EGM medium (Promocell), respectively, and let to form a monolayer. After the establishment of a stable monolayer, cells were starved overnight with DMEM without phenol red supplemented with 2% FBS and antibiotics. The next day, HUVECs or mBECs were stimulated with conditioned medium containing murine SEMA6AECTO or a mock control added in the apical or basal side of the transwell. DMEM without phenol red supplemented with 2% FBS with and without VEGF-A (50 ng/mL) was used as negative and positive control, respectively. TEER was assessed using a Millicall-ERS Voltohmeter (World Precision Instruments)86,87. TEER measurements were performed using the same conditions to allow comparison of different groups and calculated as area-integrated resistance (Ωcm2). The percentage of TEER change was calculated with respect to time 0 for each indicated time point. ### Human genetic analysis Patients from a large cohort with self-limited DP were analyzed for this study. The cohort, has been described in previous reports from our group88. Inclusion criteria were the onset of Tanner genital stage II (testicular volume > 3 mL) >13.5 years in boys or Tanner breast stage II > 13.0 years in girls (i.e., two SD later than average pubertal development) and their unaffected relatives. Normal values of thyroid function tests (TFTs) and IGF-1 for bone age were also used as inclusion criteria. Chronic illness as a cause for functional HH was excluded by medical history, clinical examination, and biochemical investigations. Congenital or acquired HH, if suspected, was excluded by spontaneous pubertal development by 18 years of age at follow-up. Sex was considered in the study design, and sex of participants was determined based on self-report. ### Whole exome sequencing analysis Whole exome sequencing was performed on DNA extracted from peripheral blood leukocytes of 100 probands, using an Agilent V5 platform and Illumina HiSeq 2000 sequencing. The exome sequences were aligned to the UCSC hg19 reference genome using the Burrows-Wheeler Aligner software (BWA-MEM [bwa-0.7.12]). Picard tools (picard-tools-1.119) was used to sort alignments and mark PCR duplicates. The genome analysis toolkit (GATK-3.4-46) was used to realign around indels and recalibrate quality scores using dbSNP, Mills and 1000 genomes as reference resources. Variant calling and joint genotyping using pedigree information was performed using HaplotypeCaller in GVCF mode from the genome analysis toolkit. The resulting variants were filtered using the variant quality score recalibration (VQSR) function from GATK. An analysis of the called variants was performed using Ingenuity Variant Analysis (QIAGEN Redwood City, [www.qiagen.com/ingenuity](http://www.qiagen.com/ingenuity)). Filtering for potential causal variants was carried out using filters for quality control (read depth and Phred strand bias), minor allele frequency (MAF <0.5% in the ExAC and Genome Aggregation Database (gnomAD) databases v2.0.2), predicted functional annotation (Poly-Phen, SIFT, REVEL, CADD score) and conservation score (GERP). Potential causal variants were confirmed by Sanger sequencing. ### *In silico* mutagenesis The CryoEM structure of the SEMA6A homodimer in complex with *C. sordellii* lethal toxin TcsL (PDB ID: 6WTS; [https://doi.org/10.2210/pdb6WTS/pdb](https://doi.org/10.2210/pdb6WTS/pdb))58 was retrieved from the RCSB Protein Data Bank. After a structure preparation procedure *via* Schrödinger BioLuminate Protein Preparation Tool, including an energy minimization using the OPLS3e force field89, the TcsL toxin was removed. The evaluation of the impact of I423T mutation on the protein stability has been carried out with the Schrödinger BioLuminate Residue Scanning Tool. Stability is computed from a thermodynamic cycle which considers the relative stability of each entity: unfolded SEMA6A, folded SEMA6A, unfolded mutated SEMA6A and folded mutated SEMA6A90. ### Immunoblotting Transfected COS-7 cells were lysed in 150 mM NaCl, 50 mM Tris-HCl (pH 7.4), and 1% TritonX-100, supplemented with protease and phosphatase inhibitors (Roche). Lysates were centrifuged at 10,000 g for 10 min at 4°C and protein concentration determined with the Bradford assay (Bio-Rad). Protein lysates (20 μg) were used for SDS-PAGE. Proteins were transferred to nitrocellulose membranes (Bio-Rad) and immunoblotted with mouse anti c-myc (1:500, clone 9E10; ThermoFisher, cat. 13-2500) or rabbit anti-GAPDH (1:1000; Cell Signaling, cat. 5174S) followed by goat anti-mouse (1:4000; Agilent, cat. P044701-2) or anti-rabbit (1:2000; Sigma-Aldrich, cat. A4914) HRP-conjugated antibodies. Images were acquired by using Image Lab software (v5.2.1, Biorad). For quantification, 3 independent experiments were performed for each condition, the OD of the signal measured with ImageJ software (v1.52a, NIH) and the mean pixel intensity calculated. ### Image acquisition Brightfield images were acquired with IS Capture software (v3.6.7, TiEsseLab). Cells and mouse tissues were examined with a Nikon A1R/AX laser scanning confocal microscope equipped with a Nikon A1/AX plus camera and the following objectives (Nikon): Plan Fluor 10X DIC L N1 (NA 0.3), Plan Fluor 20X DIC N2 (NA 0.5), 20X Plan Apo λD OFN25 DIC N2 (NA 0.80), 40X Apo LWD WI λS DIC N2 (NA 1.15), 60X Plan Apo λD OFN25 DIC N2 (NA 1.42). DAPI, Alexa488, Cy3 and Alexa647 were excited at 405, 488, 561 nm and 640 nm and observed at 425-475, 500-550, 570-620 and 663-738 nm, respectively. Plexin-A2 polarization and EB *in vivo* permeability assays were examined with a Zeiss LSM900 Airyscan laser scanning confocal microscope equipped with a Zeiss Axiocam 305 camera and the following objectives (Zeiss): EC Plan-Neofluar 20X M-27 (NA 0.5), 63X Plan-Apochromat Oil DIC M27 (NA 1.40). DAPI, Alexa488, Cy3 and Alexa647 were excited at 405, 488, 561 nm and 640 nm and observed at 400-605, 410-545, 535-617 nm and 645-700 nm. 1024x1024 pixels images were captured in a stepwise fashion over a defined z-focus range corresponding to all visible fluorescence within the sample. Maximum projections of the z-stack were performed post-acquisition by using NIS Elements AR (v5.21.03, Nikon) or ZEN 3.0 Suite (Zeiss) software. Adobe Photoshop CS6 (Adobe) and Fiji (v2.3.0/1.53f, NIH) software were used to obtain the presented images. ### Statistics and reproducibility Statistical tests employed and number of samples analyzed are outlined in figure legends and were conducted when the experiment had been performed a minimum of 3 times on a minimum of 3 distinct samples, using Prism software (v8.2.1, GraphPad). Expression studies has been performed a minimum of 2 independent times. Data are presented as mean ± SD; the differences between means are reported as mean difference ± SEM. Results are considered significant with a *p* value less than 0.05. ## DATA AVAILABILITY All data supporting the findings described in this manuscript are available in the article and its Supplementary Information files, and from the corresponding authors, upon request. Source data are provided with this paper. All unique materials are readily available from the authors upon request. The exome sequencing data are available under restricted access due to ethics restrictions on storing and sharing our pediatric patient exome data. These data will have controlled access and will be limited to individuals who enter a research agreement; use of these genomic data will be restricted to those named on the agreement, and exome data will be patient de-identified. Access can be obtained by contacting the corresponding authors, ensuring a response within 1 month. The SEMA6A variant has been deposited in the ClinVar database under accession code VCV002571587.2 ([https://www.ncbi.nlm.nih.gov/clinvar/variation/2571587/](https://www.ncbi.nlm.nih.gov/clinvar/variation/2571587/)). The single-cell transcriptomic data Tabula Muris ([https://tabula-muris.ds.czbiohub.org/](https://tabula-muris.ds.czbiohub.org/)) and EC atlas ([https://endotheliomics.shinyapps.io/ec_atlas/](https://endotheliomics.shinyapps.io/ec_atlas/)) used in this study are available in the GEO and EBI ArrayExpress database under accession code GSE109774 [https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE109774](https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE109774)) E-MTAB-8077 ([https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-8077?query=E-MTAB-8077](https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-8077?query=E-MTAB-8077)). In addition, the following databases were used: gnomAD ([https://gnomad.broadinstitute.org/](https://gnomad.broadinstitute.org/)), VarSome ([https://varsome.com](https://varsome.com)) and PDB ([https://www.rcsb.org](https://www.rcsb.org)). ## Supporting information Supplementary information [[supplements/297449_file03.pdf]](pending:yes) ## Data Availability All data supporting the findings described in this manuscript are available in the article and its Supplementary Information files, and from the corresponding authors, upon request. Source data are provided with this paper. All unique materials are readily available from the authors upon request. The exome sequencing data are available under restricted access due to ethics restrictions on storing and sharing our pediatric patient exome data. These data will have controlled access and will be limited to individuals who enter a research agreement; use of these genomic data will be restricted to those named on the agreement, and exome data will be patient de-identified. Access can be obtained by contacting the corresponding authors, ensuring a response within 1 month. The SEMA6A variant has been deposited in the ClinVar database under accession code VCV002571587.2 (https://www.ncbi.nlm.nih.gov/clinvar/variation/2571587/). The single-cell transcriptomic data Tabula Muris (https://tabula-muris.ds.czbiohub.org/) and EC atlas (https://endotheliomics.shinyapps.io/ec_atlas/) used in this study are available in the GEO and EBI ArrayExpress database under accession code GSE109774 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE109774) E-MTAB-8077 (https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-8077?query=E-MTAB-8077). In addition, the following databases were used: gnomAD (https://gnomad.broadinstitute.org/), VarSome (https://varsome.com) and PDB (https://www.rcsb.org). ## AUTHOR CONTRIBUTION A.L. R.O., M.H.v.d.M., E.Y.v.B., M.G.V., S.R.H, A.F., R.J.P. and A.C. conceived and designed experiments, collected and analyzed data, and wrote the manuscript; A.L. and R.O., performed immunohistochemical experiments; M.H.v.d.M., E.Y.v.B., M.G.V. performed *in vivo* experiments; A.L., A.J.J.P., R.A., V.A. and F.A. performed *in vitro* experiments; S.R.H and L.D. obtained and analyzed patient data; I.E and L.P. performed *in silico* modeling; C.T and M.S. performed *in vitro* EC studies; R.J.P. and A.C. performed supervision and provided final approval of the manuscript. ## COMPETING INTERESTS Authors declare no competing interests. ## ACKNOWLEDGEMENTS This work received support from: the Italian Ministry of Health (GR-2016-02362389) and the Italian Fondazione Telethon (GMR22T1054) to A.C.; the Netherlands Organization for Scientific Research (ENW-VICI) to R.J.P; the National Institute for Health Research (CL-2017-19-002), Wellcome Trust (222049/Z/20/Z), Barts Charity [MGU0552] and the Rosetrees Trust [M222-F1] to S.R.H.; Fondazione Cariplo (2018-0298) and the Fondazione Italiana per la Ricerca sul Cancro (AIRC, 22905) to A.F.. A.L. was partially supported by Fondazione Veronesi and R.O by Fondazione Collegio Ghislieri and Fondazione Telethon. We would like to thank Chiara Parravicini, Silvia Manfro and Marta Frisa for technical help, Diego De Stefani for sharing the ER-vector, Filippo Casoni, Giulia Magni and Davide Marangon for sharing antibodies, Alex Kolodkin for sharing *HA-Sema6afl/fl* mice, Chiara Colletto for assistance with scRNA-seq analysis and Nicky van Kronenburg for help with mouse breeding. Part of this work was carried out at NOLIMITS, an advanced imaging facility established by the Università degli Studi di Milano. The human embryonic and fetal material was provided by the Joint MRC/Wellcome Trust (grant # 099175/Z/12/Z) Human Developmental Biology Resource ([www.hdbr.org](http://www.hdbr.org)). ## Footnotes * $ These authors jointly supervised this work * Received October 24, 2023. * Revision received October 24, 2023. * Accepted October 26, 2023. * © 2023, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## REFERENCES 1. 1.Herbison, A. E. Control of puberty onset and fertility by gonadotropin-releasing hormone neurons. Nat. Rev. Endocrinol. 12, 452–466 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrendo.2016.70&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27199290&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) 2. 2.Abreu, A. P. & Kaiser, U. B. Pubertal development and regulation. lancet. Diabetes Endocrinol. 4, 254–264 (2016). 3. 3.Wierman, M. E., Kiseljak-Vassiliades, K. & Tobet, S. Gonadotropin-releasing hormone (GnRH) neuron migration: initiation, maintenance and cessation as critical steps to ensure normal reproductive function. Front. Neuroendocrinol. 32, 43–52 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.yfrne.2010.07.005&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20650288&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000286644600004&link_type=ISI) 4. 4.Cho, H.-J., Shan, Y., Whittington, N. C. & Wray, S. Nasal Placode Development, GnRH Neuronal Migration and Kallmann Syndrome. Front. Cell Dev. Biol. 7, 1–27 (2019). 5. 5.Prevot, V. et al. Definitive evidence for the existence of morphological plasticity in the external zone of the median eminence during the rat estrous cycle: implication of neuro-glio-endothelial interactions in gonadotropin-releasing hormone release. Neuroscience 94, 809–819 (1999). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0306-4522(99)00383-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10579572&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000083250900014&link_type=ISI) 6. 6.Prevot, V. et al. GnRH nerve terminals, tanycytes and neurohaemal junction remodeling in the adult median eminence: functional consequences for reproduction and dynamic role of vascular endothelial cells. J. Neuroendocrinol. 22, no-no (2010). 7. 7.Boehm, U. et al. Expert consensus document: European Consensus Statement on congenital hypogonadotropic hypogonadism--pathogenesis, diagnosis and treatment. Nat. Rev. Endocrinol. 11, 547– 64 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrendo.2015.112&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26194704&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) 8. 8.Howard, S. R. & Dunkel, L. Delayed Puberty - Phenotypic Diversity, Molecular Genetic Mechanisms and Recent Discoveries. Endocr. Rev. 1285–1317 (2019). doi:10.1210/er.2018-00248 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/er.2018-00248&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31220230&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) 9. 9.Oleari, R., Massa, V., Cariboni, A. & Lettieri, A. The Differential Roles for Neurodevelopmental and Neuroendocrine Genes in Shaping GnRH Neuron Physiology and Deficiency. Int. J. Mol. Sci. 22, 9425 (2021). 10. 10.Saengkaew, T. & Howard, S. R. Genetics of pubertal delay. Clin. Endocrinol. (Oxf*).* (2021). doi:10.1111/cen.14606 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/cen.14606&link_type=DOI) 11. 11.Cariboni, A. et al. Defective gonadotropin-releasing hormone neuron migration in mice lacking SEMA3A signalling through NRP1 and NRP2: implications for the aetiology of hypogonadotropic hypogonadism. Hum. Mol. Genet. 20, 336–44 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddq468&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21059704&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000285626900012&link_type=ISI) 12. 12.Cariboni, A. et al. Dysfunctional SEMA3E signaling underlies gonadotropin-releasing hormone neuron deficiency in Kallmann syndrome. J. Clin. Invest. 125, 2413–28 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1172/JCI78448&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25985275&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) 13. 13.Marcos, S. et al. Defective signaling through plexin-A1 compromises the development of the peripheral olfactory system and neuroendocrine reproductive axis in mice. Hum. Mol. Genet. 26, 2006–2017 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddx080&link_type=DOI) 14. 14.Oleari, R. et al. PLXNA1 and PLXNA3 cooperate to pattern the nasal axons that guide gonadotropin-releasing hormone neurons. Development 146, (2019). 15. 15.Giacobini, P. et al. Brain endothelial cells control fertility through ovarian-steroid-dependent release of semaphorin 3A. PLoS Biol. 12, e1001808 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pbio.1001808&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24618750&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) 16. 16.Parkash, J. et al. Suppression of β1-integrin in gonadotropin-releasing hormone cells disrupts migration and axonal extension resulting in severe reproductive alterations. J. Neurosci. 32, 16992–7002 (2012). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjExOiIzMi80Ny8xNjk5MiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzEwLzI2LzIwMjMuMTAuMjQuMjMyOTc0NDkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 17. 17.Parkash, J. et al. Semaphorin7A regulates neuroglial plasticity in the adult hypothalamic median eminence. Nat. Commun. 6, 6385 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms7385&link_type=DOI) 18. 18.Oleari, R. et al. A Novel SEMA3G Mutation in Two Siblings Affected by Syndromic GnRH Deficiency. Neuroendocrinology 111, 421–441 (2021). 19. 19.Oleari, R., Lettieri, A., Paganoni, A., Zanieri, L. & Cariboni, A. Semaphorin Signaling in GnRH Neurons: From Development to Disease. Neuroendocrinology 109, 193–199 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000495916&link_type=DOI) 20. 20.Hanchate, N. K. et al. SEMA3A, a gene involved in axonal pathfinding, is mutated in patients with Kallmann syndrome. PLoS Genet. 8, e1002896 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1002896&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22927827&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) 21. 21.Young, J. et al. SEMA3A deletion in a family with Kallmann syndrome validates the role of semaphorin 3A in human puberty and olfactory system development. Hum. Reprod. 27, 1460–5 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/humrep/des022&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22416012&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) 22. 22.Känsäkoski, J. et al. Mutation screening of SEMA3A and SEMA7A in patients with congenital hypogonadotropic hypogonadism. Pediatr. Res. 75, 641–644 (2014). 23. 23.Kotan, L. D. et al. Prevalence and associated phenotypes of PLXNA1 variants in normosmic and anosmic idiopathic hypogonadotropic hypogonadism. Clin. Genet. 95, 320–324 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/cge.13482&link_type=DOI) 24. 24.Kotan, L. D. et al. Loss-of-function variants in SEMA3F and PLXNA3 encoding semaphorin-3F and its receptor plexin-A3 respectively cause idiopathic hypogonadotropic hypogonadism. Genet. Med. (2021). doi:10.1038/s41436-020-01087-5 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41436-020-01087-5&link_type=DOI) 25. 25.Mitsogiannis, M. D., Little, G. E. & Mitchell, K. J. Semaphorin-Plexin signaling influences early ventral telencephalic development and thalamocortical axon guidance. Neural Dev. 12, 6 (2017). 26. 26.Okada, T., Keino-Masu, K., Suto, F., Mitchell, K. J. & Masu, M. Remarkable complexity and variability of corticospinal tract defects in adult Semaphorin 6A knockout mice. Brain Res. 1710, 209–219 (2019). 27. 27.Renaud, J. & Chédotal, A. Time-lapse analysis of tangential migration in Sema6A and PlexinA2 knockouts. Mol. Cell. Neurosci. 63, 49–59 (2014). 28. 28.Hatanaka, Y. et al. Semaphorin 6A-Plexin A2/A4 Interactions with Radial Glia Regulate Migration Termination of Superficial Layer Cortical Neurons. iScience 21, 359–374 (2019). 29. 29.Segarra, M. et al. Semaphorin 6A regulates angiogenesis by modulating VEGF signaling. Blood 120, 4104–15 (2012). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImJsb29kam91cm5hbCI7czo1OiJyZXNpZCI7czoxMToiMTIwLzE5LzQxMDQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8xMC8yNi8yMDIzLjEwLjI0LjIzMjk3NDQ5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 30. 30.Janssen, B. J. C. et al. Structural basis of semaphorin-plexin signalling. Nature 467, 1118–22 (2010). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature09468&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20877282&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000283548600050&link_type=ISI) 31. 31.Suto, F. et al. Interactions between plexin-A2, plexin-A4, and semaphorin 6A control lamina-restricted projection of hippocampal mossy fibers. Neuron 53, 535–47 (2007). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuron.2007.01.028&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17296555&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000245126900012&link_type=ISI) 32. 32.Zhou, L. et al. Cloning and expression of a novel murine semaphorin with structural similarity to insect semaphorin I. Mol. Cell. Neurosci. 9, 26–41 (1997). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/mcne.1997.0607&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9204478&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997XB90800003&link_type=ISI) 33. 33.Klostermann, A., Lutz, B., Gertler, F. & Behl, C. The orthologous human and murine semaphorin 6A-1 proteins (SEMA6A-1/Sema6A-1) bind to the enabled/vasodilator-stimulated phosphoprotein-like protein (EVL) via a novel carboxyl-terminal zyxin-like domain. J. Biol. Chem. 275, 39647–53 (2000). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyNzUvNTAvMzk2NDciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8xMC8yNi8yMDIzLjEwLjI0LjIzMjk3NDQ5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 34. 34.Taroc, E. Z. M., Prasad, A., Lin, J. M. & Forni, P. E. The terminal nerve plays a prominent role in GnRH-1 neuronal migration independent from proper olfactory and vomeronasal connections to the olfactory bulbs. Biol. Open 6, 1552–1568 (2017). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiYmlvbG9wZW4iO3M6NToicmVzaWQiO3M6OToiNi8xMC8xNTUyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTAvMjYvMjAyMy4xMC4yNC4yMzI5NzQ0OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 35. 35.Prevot, V. Puberty in Mice and Rats. in Knobil and Neill’s Physiology of Reproduction 2, 1395–1439 (Elsevier, 2015). 36. 36.Gill, J. C. et al. The Gonadotropin-Releasing Hormone (GnRH) Neuronal Population Is Normal in Size and Distribution in GnRH-Deficient and GnRH Receptor-Mutant Hypogonadal Mice. Endocrinology 149, 4596–4604 (2008). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/en.2008-0403&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18499748&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000258660800043&link_type=ISI) 37. 37.Caligioni, C. S. Assessing Reproductive Status/Stages in Mice. Curr. Protoc. Neurosci. 48, Appendix 4I (2009). 38. 38.Hoffmann, H. M. Determination of Reproductive Competence by Confirming Pubertal Onset and Performing a Fertility Assay in Mice and Rats. J. Vis. Exp. 1–10 (2018). doi:10.3791/58352 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3791/58352&link_type=DOI) 39. 39.Ciofi, P. et al. Brain-Endocrine Interactions: A Microvascular Route in the Mediobasal Hypothalamus. Endocrinology 150, 5509–5519 (2009). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/en.2009-0584&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19837874&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000272042300035&link_type=ISI) 40. 40.Xu, C., et al. KLB, encoding β-Klotho, is mutated in patients with congenital hypogonadotropic hypogonadism. EMBO Mol. Med. 9, 1379– 1397 (2017). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZW1ib21tIjtzOjU6InJlc2lkIjtzOjk6IjkvMTAvMTM3OSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzEwLzI2LzIwMjMuMTAuMjQuMjMyOTc0NDkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 41. 41.Daftary, S. S. & Gore, A. C. IGF-1 in the brain as a regulator of reproductive neuroendocrine function. Exp. Biol. Med. (Maywood*).* 230, 292–306 (2005). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/153537020523000503&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15855296&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000228786800003&link_type=ISI) 42. 42.Quennell, J. H. et al. Leptin indirectly regulates gonadotropin-releasing hormone neuronal function. Endocrinology 150, 2805–12 (2009). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/en.2008-1693&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19179437&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000266256700043&link_type=ISI) 43. 43.Smedlund, K. B. & Hill, J. W. The role of non-neuronal cells in hypogonadotropic hypogonadism. Mol. Cell. Endocrinol. 518, 110996 (2020). 44. 44.Desroziers, E. Unusual suspects: Glial cells in fertility regulation and their suspected role in polycystic ovary syndrome. J. Neuroendocrinol. 34, 1–12 (2022). 45. 45.Bernard, F. et al. Role of transmembrane semaphorin Sema6A in oligodendrocyte differentiation and myelination. Glia 60, 1590–1604 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/glia.22378&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22777942&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000307568900012&link_type=ISI) 46. 46.Rashighi, M. & Harris, J. E. Single-cell transcriptomics of 20 mouse organs creates a Tabula Muris. Nature 562, 367–372 (2018). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-018-0590-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30283141&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) 47. 47.Srinivasan, B. et al. TEER Measurement Techniques for In Vitro Barrier Model Systems. J. Lab. Autom. 20, 107–126 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/2211068214561025&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25586998&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) 48. 48.St. Clair, R. M. et al. <lt;scp>PKC</scp> induces release of a functional ectodomain of the guidance cue semaphorin6A. FEBS Lett. 593, 3015– 3028 (2019). 49. 49.Esser, S., Lampugnani, M. G., Corada, M., Dejana, E. & Risau, W. Vascular endothelial growth factor induces VE-cadherin tyrosine phosphorylation in endothelial cells. J. Cell Sci. 111 **(** **Pt 1**, 1853–65 (1998). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NToiam9jZXMiO3M6NToicmVzaWQiO3M6MTE6IjExMS8xMy8xODUzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTAvMjYvMjAyMy4xMC4yNC4yMzI5NzQ0OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 50. 50.Fantin, A. et al. VEGF165-induced vascular permeability requires NRP1 for ABL-mediated SRC family kinase activation. J. Exp. Med. 214, 1049–1064 (2017). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamVtIjtzOjU6InJlc2lkIjtzOjEwOiIyMTQvNC8xMDQ5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTAvMjYvMjAyMy4xMC4yNC4yMzI5NzQ0OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 51. 51.Brash, J. T., Ruhrberg, C. & Fantin, A. Evaluating Vascular Hyperpermeability-inducing Agents in the Skin with the Miles Assay. J. Vis. Exp. 2018, 2–7 (2018). 52. 52.Bernard-Patrzynski, F. et al. Isolation of endothelial cells, pericytes and astrocytes from mouse brain. PLoS One 14, e0226302 (2019). 53. 53.Kalucka, J. et al. Single-Cell Transcriptome Atlas of Murine Endothelial Cells. Cell 180, 764–779.e20 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2020.01.015&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32059779&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) 54. 54.Stenzel, D. et al. Endothelial basement membrane limits tip cell formation by inducing Dll4/Notch signalling in vivo. EMBO Rep. 12, 1135–1143 (2011). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NToiZW1ib3IiO3M6NToicmVzaWQiO3M6MTA6IjEyLzExLzExMzUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8xMC8yNi8yMDIzLjEwLjI0LjIzMjk3NDQ5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 55. 55.Mullier, A., Bouret, S. G., Prevot, V. & Dehouck, B. Differential distribution of tight junction proteins suggests a role for tanycytes in blood-hypothalamus barrier regulation in the adult mouse brain. J. Comp. Neurol. 518, 943–62 (2010). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/cne.22273&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20127760&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000274834200001&link_type=ISI) 56. 56.Langlet, F. et al. Tanycytic VEGF-A Boosts Blood-Hypothalamus Barrier Plasticity and Access of Metabolic Signals to the Arcuate Nucleus in Response to Fasting. Cell Metab. 17, 607–617 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cmet.2013.03.004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23562080&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000326265600015&link_type=ISI) 57. 57.Jiang, H. et al. MCH Neurons Regulate Permeability of the Median Eminence Barrier. Neuron 107, 306–319.e9 (2020). 58. 58.Lee, H. et al. Recognition of Semaphorin Proteins by P. sordellii Lethal Toxin Reveals Principles of Receptor Specificity in Clostridial Toxins. Cell 182, 345–356.e16 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2020.06.005&link_type=DOI) 59. 59.Kyte, J. & Doolittle, R. F. A simple method for displaying the hydropathic character of a protein. J. Mol. Biol. 157, 105–132 (1982). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0022-2836(82)90515-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7108955&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1982NQ28000006&link_type=ISI) 60. 60.Perez-Branguli, F. et al. Reverse Signaling by Semaphorin-6A Regulates Cellular Aggregation and Neuronal Morphology. PLoS One 11, e0158686 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0158686&link_type=DOI) 61. 61.Jongbloets, B. C. & Pasterkamp, R. J. Semaphorin signalling during development. Development 141, 3292–7 (2014). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZGV2ZWxvcCI7czo1OiJyZXNpZCI7czoxMToiMTQxLzE3LzMyOTIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8xMC8yNi8yMDIzLjEwLjI0LjIzMjk3NDQ5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 62. 62.Dhanabal, M. et al. Recombinant semaphorin 6A-1 ectodomain inhibits in vivo growth factor and tumor cell line-induced angiogenesis. Cancer Biol. Ther. 4, 659–668 (2005). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1158/1535-7163.MCT-04-0290&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15917651&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000236042600024&link_type=ISI) 63. 63.Zink, J. et al. EVL regulates VEGF receptor-2 internalization and signaling in developmental angiogenesis. EMBO Rep. 22, 1–19 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.15252/embr.202051851&link_type=DOI) 64. 64.Miyata, S. Advances in Understanding of Structural Reorganization in the Hypothalamic Neurosecretory System. Front. Endocrinol. (Lausanne*).* 8, 1–18 (2017). 65. 65.Ojeda, S. R., Lomniczi, A. & Sandau, U. S. Glial-Gonadotrophin Hormone (GnRH) Neurone Interactions in the Median Eminence and the Control of GnRH Secretion. J. Neuroendocrinol. 20, 732–742 (2008). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-2826.2008.01712.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18601696&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000256582900017&link_type=ISI) 66. 66.DiVall, S. A. et al. Divergent roles of growth factors in the GnRH regulation of puberty in mice. J. Clin. Invest. 120, 2900–2909 (2010). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1172/JCI41069&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20628204&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000280492100025&link_type=ISI) 67. 67.Sedlmeyer, I. L. & Palmert, M. R. Delayed puberty: analysis of a large case series from an academic center. J. Clin. Endocrinol. Metab. 87, 1613–20 (2002). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/jc.87.4.1613&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11932291&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000174963100030&link_type=ISI) 68. 68.Messina, A. et al. Dysregulation of Semaphorin7A/β1-integrin signaling leads to defective GnRH-1 cell migration, abnormal gonadal development and altered fertility. Hum. Mol. Genet. 20, 4759–74 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddr403&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21903667&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000297242100001&link_type=ISI) 69. 69.Chen, H. et al. Limb and kidney defects in Lmx1b mutant mice suggest an involvement of LMX1B in human nail patella syndrome. Nat. Genet. 19, 51–55 (1998). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng0598-51&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9590288&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000073346700016&link_type=ISI) 70. 70.Lakso, M. et al. Efficient in vivo manipulation of mouse genomic sequences at the zygote stage. Proc. Natl. Acad. Sci. U. S. A. 93, 5860–5 (1996). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMDoiOTMvMTIvNTg2MCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzEwLzI2LzIwMjMuMTAuMjQuMjMyOTc0NDkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 71. 71.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–66 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/glia.22455&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23404564&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) 72. 72.Chi, S. et al. Astrocytic Piezo1-mediated mechanotransduction determines adult neurogenesis and cognitive functions. Neuron 110, 2984–2999.e8 (2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuron.2022.07.010&link_type=DOI) 73. 73.Ramalho, A. F. et al. Dietary fats promote functional and structural changes in the median eminence blood/spinal fluid interface—the protective role for BDNF. J. Neuroinflammation 15, 10 (2018). 74. 74.Bonfanti, E. et al. Abnormal Upregulation of GPR17 Receptor Contributes to Oligodendrocyte Dysfunction in SOD1 G93A Mice. Int. J. Mol. Sci. 21, 2395 (2020). 75. 75.Marangon, D. et al. In vivo silencing of miR-125a-3p promotes myelin repair in models of white matter demyelination. Glia 68, 2001–2014 (2020). 76. 76.Han, S. et al. Proneural genes define ground-state rules to regulate neurogenic patterning and cortical folding. Neuron 109, 2847–2863.e11 (2021). 77. 77.Plein, A. et al. Neural crest-derived SEMA3C activates endothelial NRP1 for cardiac outflow tract septation. J. Clin. Invest. 125, 2661– 2676 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1172/JCI79668&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26053665&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) 78. 78.Poltavski, D. M. et al. Venous endothelin modulates responsiveness of cardiac sympathetic axons to arterial semaphorin. Elife 8, 1–23 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/elife.39324&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) 79. 79.Méndez-Gómez, H. R. et al. Transcytosis in the blood–cerebrospinal fluid barrier of the mouse brain with an engineered receptor/ligand system. Mol. Ther. - Methods Clin. Dev. 2, 15037 (2015). 80. 80.Oleari, R. et al. Autism-linked NLGN3 is a key regulator of gonadotropin-releasing hormone deficiency. Dis. Model. Mech. 16, (2023). 81. 81.Taroc, E. Z. M., Lin, J. M., Tulloch, A. J., Jaworski, A. & Forni, P. E. GnRH-1 Neural Migration From the Nose to the Brain Is Independent From Slit2, Robo3 and NELL2 Signaling. Front. Cell. Neurosci. 13, 70 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fncel.2019.00070&link_type=DOI) 82. 82.Macchi, C. et al. Iron overload induces hypogonadism in male mice via extrahypothalamic mechanisms. Mol. Cell. Endocrinol. 454, 135–145 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.mce.2017.06.019&link_type=DOI) 83. 83.Lee, E. B., Chakravarthi, V. P., Wolfe, M. W. & Rumi, M. A. K. ERβ Regulation of Gonadotropin Responses during Folliculogenesis. Int. J. Mol. Sci. 22, 10348 (2021). 84. 84.Domingues, A. & Fantin, A. Neuropilin 1 Regulation of Vascular Permeability Signaling. Biomolecules 11, 666 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/biom11050666&link_type=DOI) 85. 85.Cannarella, R. et al. Anti-Müllerian Hormone, Growth Hormone, and Insulin-Like Growth Factor 1 Modulate the Migratory and Secretory Patterns of GnRH Neurons. Int. J. Mol. Sci. 22, 2445 (2021). 86. 86.Tacconi, C. et al. Vascular Endothelial Growth Factor C Disrupts the Endothelial Lymphatic Barrier to Promote Colorectal Cancer Invasion. Gastroenterology 148, 1438–1451.e8 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1053/j.gastro.2015.03.005&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25754161&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) 87. 87.Ungaro, F. et al. MFSD2A Promotes Endothelial Generation of Inflammation-Resolving Lipid Mediators and Reduces Colitis in Mice. Gastroenterology 153, 1363–1377.e6 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1053/j.gastro.2017.07.048&link_type=DOI) 88. 88.Mancini, A. et al. LGR4 deficiency results in delayed puberty through impaired Wnt/β-catenin signaling. JCI Insight 5, (2020). 89. 89.Harder, E. et al. OPLS3: A Force Field Providing Broad Coverage of Drug-like Small Molecules and Proteins. J. Chem. Theory Comput. 12, 281–296 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/acs.jctc.5b00864&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26584231&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F10%2F26%2F2023.10.24.23297449.atom) 90. 90.Beard, H., Cholleti, A., Pearlman, D., Sherman, W. & Loving, K. A. Applying Physics-Based Scoring to Calculate Free Energies of Binding for Single Amino Acid Mutations in Protein-Protein Complexes. PLoS One 8, e82849 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0082849&link_type=DOI)