Embryonic origins of forebrain oligodendrocytes revisited by combinatorial genetic fate mapping
eLife assessment
In this study the authors revisited the question of the embryonic origin of telencephalic oligodendrocytes using some new and powerful genetic tools. There is convincing evidence to support previous suggestions of a predominantly cortical origin of oligodendrocytes in the cerebral cortex, however the new studies suggest that LGE/CGE-derived oligodendrocytes make a modest contribution in some areas, while MGE/POA-derived oligodendrocytes make a small but enduring contribution. The findings are valuable and should be of interest to developmental and myelin biologists.
https://doi.org/10.7554/eLife.95406.3.sa0Valuable: Findings that have theoretical or practical implications for a subfield
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Convincing: Appropriate and validated methodology in line with current state-of-the-art
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Abstract
Multiple embryonic origins give rise to forebrain oligodendrocytes (OLs), yet controversies and uncertainty exist regarding their differential contributions. We established intersectional and subtractional strategies to genetically fate map OLs produced by medial ganglionic eminence/preoptic area (MGE/POA), lateral/caudal ganglionic eminences (LGE/CGE), and dorsal pallium in the mouse brain. We found that, contrary to the canonical view, LGE/CGE-derived OLs make minimum contributions to the neocortex and corpus callosum, but dominate piriform cortex and anterior commissure. Additionally, MGE/POA-derived OLs, instead of being entirely eliminated, make small but sustained contribution to cortex with a distribution pattern distinctive from those derived from the dorsal origin. Our study provides a revised and more comprehensive view of cortical and white matter OL origins, and established valuable new tools and strategies for future OL studies.
Introduction
Oligodendrocytes (OLs) are an important class of macroglia responsible for producing the myelin sheaths that insulate and protect neuronal axons. Forebrain OLs arise from multiple embryonic origins. Previous fate-mapping study using Nkx2.1Cre (Xu et al., 2008), Gsh2Cre (Kessaris et al., 2006), and Emx1Cre (Gorski et al., 2002) reported consecutive and competing waves of OLs derived from medial ganglionic eminence/preoptic area (MGE/POA), lateral/caudal ganglionic eminences (LGE/CGE), and dorsal pallium (Kessaris et al., 2006). The first wave of OLs generated by MGE/POA (MPOLs) was believed to be eliminated postnatally, while those from the second and third waves (LCOLs and dOLs) survive and populate the cortex and corpus callosum at comparable proportions. Several other studies provided both supporting and contradicting evidence to this model (Nakahira et al., 2006; Tsoa et al., 2014; Naruse et al., 2016; Orduz et al., 2019; Liu et al., 2021; Shen et al., 2021; Tripathi et al., 2011; Winkler et al., 2018). Moreover, Gsh2 was recently found to be expressed in dorsal progenitors (Zhang et al., 2020), casting doubt on the interpretation of lineage tracing data from Gsh2Cre.
In this study, we generated new genetic tools and combinatorial fate-mapping strategies which allow direct visualization and comparison among OLs derived from different origins. We found that neocortical OLs are primarily composed of dOLs, rather than similar proportions of LCOLs and dOLs. In contrast, LCOLs and dOLs made comparable contributions to piriform cortex. We also found that although MPOLs only make a small contribution, they do persist in the cortex beyond adulthood with a unique spatial pattern distinct from that of the dOLs. In the two major white matter commissure tracts, dOLs are the vast majority in corpus callosum but make little contribution to anterior commissure, while LCOLs behaved the opposite. These findings significantly revised the classical view and provided a new and more comprehensive picture of cortical and white matter OL origins.
Results and discussion
To unambiguously track OLs from different embryonic origins, we first generated a knock-in driver, OpalinP2A-Flpo-T2A-tTA2 (Figure 1), orthogonal to Cre drivers that label dorsal or ventral progenitors (ProgenitorCre). Opalin (also known as Tmem10) encodes oligodendrocytic myelin paranodal and inner loop protein that are specifically expressed in differentiated OLs (Kippert et al., 2008; Yoshikawa et al., 2008; Jiang et al., 2013; Marques et al., 2016). In OpalinP2A-Flpo-T2A-tTA2, Flpo and tTA2 were inserted before the STOP codon and linked by self-cleavage peptide P2A and T2A (Figure 1A–C), allowing co-transcription and translation with Opalin. Flp-mediated recombination by this driver (hereinafter referred to as OpalinFlp for simplicity) enables highly specific, efficient, and irreversible OL labeling, while the tTA2 component offers the flexibility for OL-specific labeling in tunable densities (Figure 1D–F).
-
Figure 1—source data 1
- https://cdn.elifesciences.org/articles/95406/elife-95406-fig1-data1-v1.zip
-
Figure 1—source data 2
- https://cdn.elifesciences.org/articles/95406/elife-95406-fig1-data2-v1.zip
-
Figure 1—source data 3
- https://cdn.elifesciences.org/articles/95406/elife-95406-fig1-data3-v1.zip
-
Figure 1—source data 4
- https://cdn.elifesciences.org/articles/95406/elife-95406-fig1-data4-v1.zip
-
Figure 1—source data 5
- https://cdn.elifesciences.org/articles/95406/elife-95406-fig1-data5-v1.xlsx
Next, we established two types of genetic combinatorial fate-mapping strategies to directly visualize OLs from different embryonic origins (Figure 2): (1) combining OpalinFlp and ProgenitorCre with intersectional reporters Ai65 to label OLs derived from Cre+ progenitor domain by RFP (Figure 2A); (2) combining OpalinFlp and ProgenitorCre with RC::FLTG (Plummer et al., 2015) to simultaneously label OLs derived from Cre+ progenitors by green fluorescent protein (GFP) and OLs derived from the complementing Cre− progenitors by RFP (Figure 2B, Figure 2—figure supplement 1A). The first approach allowed us to track dOLs and MPOLs (Figure 2C, E). The second approach empowered us to observe and compare OLs generated from dorsal and ventral origins (Figure 2—figure supplement 1B), or those from Gsh2+ and Gsh2− progenitors, in the same brain (Figure 2—figure supplement 1C). Importantly, the subtraction power enabled us to target OLs derived from LGE/CGE progenitors that express neither Emx1 nor Nkx2.1 (Figure 2D and Figure 2—figure supplement 1D). In addition, these strategies greatly facilitated the identification of OLs derived from specific origin which exist at relatively low density in certain regions.
Deploying these strategies, we assessed the differential contributions of dOLs, LCOLs, and MPOLs by analyzing RFP+ cells in the following mice: OpalinFlp::Emx1Cre::Ai65 (Figure 2C), OpalinFlp::Emx1Cre::Nkx2.1Cre::RC::FLTG (Figure 2D), and OpalinFlp::Nkx2.1Cre::Ai65 (Figure 2E). To better assess their contributions to the total OL population (Figure 2F), we co-stained RFP with the mature OL marker aspartoacylase (ASPA) (Huang et al., 2023; Figure 2G–I) and quantified the ratio of co-localization (Figure 2J). Notably, all RFP+ cells are ASPA+, reassured the specificity of our label strategies. We observed two significant differences from the traditional model in the neocortex. The first major deviation is that, instead of comparable contributions by dOLs and LCOLs, the vast majority of neocortical OLs were dOLs but not LCOLs. The densities (Figure 2G–I and Figure 2—figure supplement 2A–F) and ASPA ratios (Figure 2J) of dOLs are much higher than those of LCOLs. Considering the possibility of incomplete recombination in combinatorial reporters, and the relatively low Cre activity in the dorsal MGE of Nkx2.1Cre (Xu et al., 2008), the genuine contribution of LCOLs to the neocortex could be even lesser than our current observation. Therefore, the large quantity of neocortical OLs labeled by Gsh2Cre in previous study (Kessaris et al., 2006) or by GFP in OpalinFlp::Gsh2Cre::RC::FLTG (Figure 2—figure supplement 1C) most likely were predominantly dOLs generated by Gsh2+ dorsal progenitors (Zhang et al., 2020), rather than bona fide LCOLs.
The second major deviation is that cortical MPOLs are not completely depleted postnatally. Instead, they make a small but continued contribution with a unique spatial distribution pattern (Figure 2I, J and Figure 2—figure supplement 2D–G). MPOLs display a clear rostrocaudal density decline (Figure 2—figure supplement 2E, F), a higher density in somatosensory cortex (SS) than motor cortex (Mo) (Figure 2I), and a laminar preference toward layer 4 (L4) in SS (Figure 2—figure supplement 2G). In contrast, the distribution of dOLs and LCOLs do not vary significantly across the rostrocaudal axis (Figure 2—figure supplement 2E, F) or between Mo and SS (Figure 2G, H), but exhibits increased density toward deeper layers (Figure 2—figure supplement 2G). Importantly, we have observed cortical MPOLs in mice as old as 1 year (Figure 2—figure supplement 2H), well beyond the age analyzed in previous reports (Kessaris et al., 2006; Orduz et al., 2019; Liu et al., 2021), suggesting a persisted contribution.
We then turned our attention to the lateral three-layer archicortex, piriform cortex (Pir). Different from the neocortex, Pir contains higher proportions (Figure 2J) of LCOLs than dOLs. MPOLs make the lowest contribution (Figure 2J) at a density similar to SS and higher than Mo (Figure 2I).
These combinatorial models also grant us the opportunity to revisit the differential contributions of dOLs, LCOLs, and MPOLs to the two commissural white matter tracts, corpus callosum (cc), and anterior commissure (ac), which contain high density of OLs (Figure 2K–M). We found that, similar to the neocortex, cc is mainly populated by dOLs and supplemented by very low proportions of LCOLs and MPOLs (Figure 2K–N). Interestingly, LCOLs and MPOLs seem to show preferential distribution in the lateral and medial regions of cc (cc-l and cc-m), respectively (Figure 2L–N). Different from cc, ac is mainly populated by LCOLs and MPOLs and supplemented by very low proportion of dOLs (Figure 2K–N).
To substantiate the above results, we further breed OpalinFlp::Emx1Cre::Nkx2.1Cre::Ai65 to label dOLs together with MPOLs by RFP and co-stained them with ASPA (Figure 2—figure supplement 3). RFP−ASPA+ cells were difficult to find in Mo, SS, and cc, but were more easily observed in Pir and ac, consistent with the respective low and high LCOL contributions in these regions.
In summary, our findings significantly revised the canonical model of forebrain OL origins (Figure 3A), and provided a new and more comprehensive view (Figure 3B). We demonstrated that neocortical OLs are mainly derived from dorsal origin with small but lasting contribution from the ventral origin (Figure 2, Figure 2—figure supplements 1B and 2). Our data showed that LGE/CGE makes little contribution to neocortex and cc, but makes major contribution to piriform cortex and ac (Figure 2 and Figure 2—figure supplement 3). This finding is supported by another report in which in utero electroporation failed to label LGE-derived cortical OLs in both embryonic and early postnatal brains, and an exclusion strategy revealed very low percentage of LGE/CGE-derived cortical OLs in neonatal brains (Li et al., 2023). The lack of adult labeling in our study together with the lack of developmental labeling in the other study suggests that the lack of LCOL in neocortex is less likely caused by competitive postnatal elimination, but more likely due to limited production and/or allocation. We further discovered that MGE/POA makes a small but persistent contribution to the neocortex with a distinct distribution pattern featured by a rostral-high to caudal-low gradient and a preference toward L4 in SS (Figure 2—figure supplement 2). Whether their enduring existence and highly biased localization has functional implications awaits future exploration. In addition, we found that the cc showed a similar OL composition as the neocortex, but the Pir and the ac each exhibited distinct OL compositions in term of their embryonic origins. LCOLs are the major contributor to both regions, while dOLs and MPOLs mainly contribute to Pir and ac, respectively (Figure 2).
In addition to the new framework of forebrain OL origins (Figure 3), we also generated a new driver (Figure 1) and established multiple combinatorial genetic models (Figure 2) for efficient tracking and direct visualization of OLs from different embryonic origins without interference from other cells types sharing the same progenitor domains such as OL precursors, astrocytes, and neurons (Figures 1–2). These tools set up a firm foundation and will provide reliable experimental access for future inquiries on the development and function of diverse OLs in healthy and disease brains (Gong et al., 2022), especially to uncover the relationship between their developmental origins and the functional and molecular heterogeneity.
Materials and methods
Mice
All mouse studies were carried out in strict accordance with the guidelines of the Institutional Animal Care and Use Committee of School of Basic Medical Sciences, Fudan University. All husbandry and experimental procedures were reviewed and approved by the same committee (Permit Number: 20210302-137). All applicable institutional and/or national guidelines for the care and use of animals were followed. The following transgenic mouse lines were used in this study: Nkx2.1Cre (Jax 008661) (Xu et al., 2008), Gsh2Cre (Jax 025806) (Kessaris et al., 2006), Emx1Cre (Jax 005628) (Gorski et al., 2002), Ai65 (Jax 021875) (Madisen et al., 2015), and RC::FLTG (Jax 026932) (Plummer et al., 2015). The tTA2-dependent tdTomato reporter (TRE-RFP) was derived from Ai62 (Jax 022731) (Madisen et al., 2015), by removing LoxP-STOP-LoxP with E2a-Cre (Jax 003724). The Flp-dependent H2B-GFP reporter (HG-FRT) was derived from HG-dual (Jax 028581) via removal of loxP flanking STOP cassette by CMV-Cre (He et al., 2016). The OpalinP2A-Flpo-T2A-tTA2 allele was generated by targeted insertion of the T2A-Flpo-P2A-tTA2 sequence immediately before the STOP codon of the endogenous Opalin gene using homologous recombination. Gene targeting vector was generated using PCR-based cloning approach as described before (He et al., 2016). More specifically, a 4.7-kb 5′ homology arm, a loxP flanking Neo-positive selection cassette, a T2A-Flpo-P2A-tTA2 cassette and a 2.7-kb 3′ homology arm were cloned into a building vector containing the DTA-negative selection cassette to generate the targeting vector. Targeting vector was linearized and transfected into a C57/black6 ES cell line. ES clones that survived through negative and positive selections were first screened by genomic PCR, then confirmed by Southern blotting using appropriate DIG-dUTP-labeled probes. One positive ES cell clone was used for blastocyst injection to obtain male chimera mice carrying the modified allele following standard procedures. Chimera males were bred with C57BL/6J females to confirm germline transmission by genomic PCR. The Neo selection cassette was self-excised during spermatogenesis of F0 chimeras. Heterozygous F1 siblings were bred with one another to establish the colony. Targeting vector construction, ES cell transfections and screening, blastocyst injections, and chimera breeding were performed by Cyagen.
Genomic PCR
Request a detailed protocolGenomic DNA was prepared from mouse tails. Tissue was lysed by incubation in tail lysis buffer (Viagen, 102-T) with 0.1 mg/ml proteinase K (Diamond, A100706) overnight at 55°C followed by 45 min at 90°C in an air bath to inactivate proteinase K. The lysate was cleared by centrifugation at maximum speed (21,130 G) for 15 min in a table-top centrifuge. Supernatant containing genomic DNA was used as the PCR template for amplifying DNA products. The following primers were used:
Opalin-F: 5′-GGCCTATGTTTGATTTCCAGCACTG-3′
Opalin-R: 5′-AGCACTTATGACTGCTGAGCCGTTC-3′
Immunohistochemistry and microscopy
Request a detailed protocolMice were anesthetized by intraperitoneal injection of 1.5% sodium pentobarbital (0.09 mg/g body weight) and then intracardially perfused with saline followed by 4% paraformaldehyde in 0.1 M phosphate buffer. Following post fixation at 4°C for 24 hr, brain samples were sectioned at 30 μm using a vibratome (Leica VT1000S), or transferred into 30% sucrose in 0.1 M PB for cryoprotection, embedded in optimal cutting temperature (OCT) compound, and sectioned using a cryostat (Leica CM1950). For CC1 immunostaining, antigen retrieval was performed prior to blocking by boiling for 3 min in 10 mM citrate buffer (pH 6.0). Sections were blocked in phosphate buffered saline (PBS) containing 0.05% Triton and 5% normal donkey serum and then incubated with the following primary antibodies in the blocking solution at 4°C overnight: RFP (goat polyclonal antibody, 1:2000, SICGEN AB0081-200; rabbit polyclonal antibody, 1:2000, Rockland 600-401-379), GFP (chicken polyclonal antibody, 1:1000, Aves Labs, GFP-1020), MBP (rat polyclonal antibody, 1:500, AbD Serotec, MCA409S), CC1 (rabbit polyclonal antibody, 1:500, Oasis Biofarm, OB-PRB070, mouse polyclonal antibody, 1:300, Millipore, OP80), ASPA (rat polyclonal antibody, 1:200, Oasis Biofarm, OB-PRT005), Sox9 (rabbit polyclonal antibody, 1:2000, Chemicon, AB5535), and NeuN (mouse monoclonal antibody, 1:500, Millipore, MAB377). Sections were then incubated with appropriate Alexa fluor dye-conjugated IgG secondary antibodies (1:500, Thermo Fisher Scientific) or CF dye-conjugated IgG secondary antibodies (1:250, Sigma) in blocking solution and mounted in Aqua-mount (Southern Biotech, 0100-01). Sections were counterstained with DAPI (4',6-diamidino-2-phenylindole). Sections were imaged with confocal microscopy (Olympus FV3000), fluorescence microscopy (Nikon Eclipse Ni; Olympus VS120; Olympus VS200), and fluorescent stereoscope (Nikon SMZ25). All quantifications were performed in 2-month-old adult mice from coronal sections between Bregma +1.94 and −2.80 mm. Anatomical regions were identified according to the Paxinos ‘The Mouse Brain’ Atlas and the Allen Reference Atlas, and their areas were measured in ImageJ for density calculations, whenever applicable. For cortical regions, every fourth section within the range of selection was analyzed. For whiter matter tracts, three consecutive sections at Bregma 0.14 were analyzed. At least three brains were analyzed for each genotype. To quantify density and co-localization, cells were identified and counted in Adobe Photoshop or ImageJ in conjugation with QuPath.
Statistical analysis
Request a detailed protocolGraphPad Prism version 8.0.1 was used for statistical calculations. No statistical methods were used to predetermine sample sizes, but our sample sizes are similar to those reported in previous publications. Data collection and analysis were performed blind to the conditions of the experiments whenever possible. No animals or data points were excluded from the analysis. Normalcy was assessed using Shapiro–Wilk test. Equal variances were assessed using F test or Bartlett’s test. Statistical significance was tested using two-tailed unpaired t-test, Welch’s t-test, one-way analysis of variance (ANOVA), and two-way ANOVA followed by Tukey’s or Bonferroni post hoc test, wherever appropriate. Data are presented as mean ± standard error of the mean. p < 0.05 was considered significant. Significance is marked as *p < 0.05, **p < 0.01, and ***p <0.001.
Data availability
All data generated or analyzed during this study are included in the manuscript. Source data have been provided for Figures 1 and 2.
References
-
Cortical excitatory neurons and glia, but not GABAergic neurons, are produced in the Emx1-expressing lineageThe Journal of Neuroscience 22:6309–6314.https://doi.org/10.1523/JNEUROSCI.22-15-06309.2002
-
Immunological markers for central nervous system gliaNeuroscience Bulletin 39:379–392.https://doi.org/10.1007/s12264-022-00938-2
-
Identification of Tmem10 as a novel late-stage oligodendrocytes marker for detecting hypomyelinationInternational Journal of Biological Sciences 10:33–42.https://doi.org/10.7150/ijbs.7526
-
Dorsally and ventrally derived oligodendrocytes have similar electrical properties but myelinate preferred tractsThe Journal of Neuroscience 31:6809–6819.https://doi.org/10.1523/JNEUROSCI.6474-10.2011
-
Fate mapping Nkx2.1-lineage cells in the mouse telencephalonThe Journal of Comparative Neurology 506:16–29.https://doi.org/10.1002/cne.21529
-
Opalin, a transmembrane sialylglycoprotein located in the central nervous system myelin paranodal loop membraneThe Journal of Biological Chemistry 283:20830–20840.https://doi.org/10.1074/jbc.M801314200
Article and author information
Author details
Funding
National Science and Technology Innovation 2030 Major Projects of China (STI2030-Major Projects-2022ZD0206500)
- Miao He
National Natural Science Foundation of China (32171087)
- Miao He
National Natural Science Foundation of China (32371145)
- Ling Gong
The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Dr. Yilin Tai for helpful discussion and Dr. Min Jiang from core facility of IOBS, Fudan University, for technical support on imaging experiments. This study was supported by funds from the National Science and Technology Innovation 2030 Major Projects of China (STI2030-Major Projects-2022ZD0206500), National Natural Science Foundation of China (32171087, 32371145).
Ethics
All mouse studies were carried out in strict accordance with the guidelines of the Institutional Animal Care and Use Committee of School of Basic Medical Sciences, Fudan University. All husbandry and experimental procedures were reviewed and approved by the same committee (Permit Number: 20210302-137). All applicable institutional and/or national guidelines for the care and use of animals were followed.
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.95406. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2024, Cai, Zhao, Shi et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 951
- views
-
- 100
- downloads
-
- 0
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
- Genetics and Genomics
Paternal obesity has been implicated in adult-onset metabolic disease in offspring. However, the molecular mechanisms driving these paternal effects and the developmental processes involved remain poorly understood. One underexplored possibility is the role of paternally-induced effects on placenta development and function. To address this, we investigated paternal high-fat diet-induced obesity in relation to sperm histone H3 lysine 4 tri-methylation signatures, the placenta transcriptome and cellular composition. C57BL6/J male mice were fed either a control or high-fat diet for 10 weeks beginning at 6 weeks of age. Males were timed-mated with control-fed C57BL6/J females to generate pregnancies, followed by collection of sperm, and placentas at embryonic day (E)14.5. Chromatin immunoprecipitation targeting histone H3 lysine 4 tri-methylation (H3K4me3) followed by sequencing (ChIP-seq) was performed on sperm to define obesity-associated changes in enrichment. Paternal obesity corresponded with altered sperm H3K4me3 at promoters of genes involved in metabolism and development. Notably, sperm altered H3K4me3 was also localized at placental enhancers. Bulk RNA-sequencing on placentas revealed paternal obesity-associated sex-specific changes in expression of genes involved in hypoxic processes such as angiogenesis, nutrient transport, and imprinted genes, with a subset of deregulated genes showing changes in H3K4me3 in sperm at corresponding promoters. Paternal obesity was also linked to impaired placenta development; specifically, a deconvolution analysis revealed altered trophoblast cell lineage specification. These findings implicate paternal obesity-effects on placenta development and function as one potential developmental route to offspring metabolic disease.
-
- Developmental Biology
- Stem Cells and Regenerative Medicine
Deficient Anterior pituitary with common Variable Immune Deficiency (DAVID) syndrome results from NFKB2 heterozygous mutations, causing adrenocorticotropic hormone deficiency (ACTHD) and primary hypogammaglobulinemia. While NFKB signaling plays a crucial role in the immune system, its connection to endocrine symptoms is unclear. We established a human disease model to investigate the role of NFKB2 in pituitary development by creating pituitary organoids from CRISPR/Cas9-edited human induced pluripotent stem cells (hiPSCs). Introducing homozygous TBX19K146R/K146R missense pathogenic variant in hiPSC, an allele found in congenital isolated ACTHD, led to a strong reduction of corticotrophs number in pituitary organoids. Then, we characterized the development of organoids harboring NFKB2D865G/D865G mutations found in DAVID patients. NFKB2D865G/D865G mutation acted at different levels of development with mutant organoids displaying changes in the expression of genes involved on pituitary progenitor generation (HESX1, PITX1, LHX3), hypothalamic secreted factors (BMP4, FGF8, FGF10), epithelial-to-mesenchymal transition, lineage precursors development (TBX19, POU1F1) and corticotrophs terminal differentiation (PCSK1, POMC), and showed drastic reduction in the number of corticotrophs. Our results provide strong evidence for the direct role of NFKB2 mutations in the endocrine phenotype observed in patients leading to a new classification of a NFKB2 variant of previously unknown clinical significance as pathogenic in pituitary development.