Genetic mechanisms control the linear scaling between related cortical primary and higher order sensory areas
Abstract
In mammals, the neocortical layout consists of few modality-specific primary sensory areas and a multitude of higher order ones. Abnormal layout of cortical areas may disrupt sensory function and behavior. Developmental genetic mechanisms specify primary areas, but mechanisms influencing higher order area properties are unknown. By exploiting gain-of and loss-of function mouse models of the transcription factor Emx2, we have generated bi-directional changes in primary visual cortex size in vivo and have used it as a model to show a novel and prominent function for genetic mechanisms regulating primary visual area size and also proportionally dictating the sizes of surrounding higher order visual areas. This finding redefines the role for intrinsic genetic mechanisms to concomitantly specify and scale primary and related higher order sensory areas in a linear fashion.
https://doi.org/10.7554/eLife.11416.001eLife digest
The neocortex is the most recently evolved part of the human brain. It is associated with higher thought processes, including language and the processing of information from our senses. Anatomically, the neocortex is organised into different regions called ‘primary areas’ and ‘higher order areas’, and perturbations to this organisation are associated with disorders such as autism.
There are many more higher order areas than primary areas in a mammalian brain. But, while primary areas are known to be specified by developmental genes in the embryo, little is known about how the development of higher order areas is controlled. Recent findings suggested that primary areas might themselves influence the emergence of higher order areas via a series of developmental events.
Now, Zembrzycki, Stocker et al. have investigated the developmental mechanisms that organise the mouse neocortex into its different regions. The experiments involved mouse mutants that produce either too much or too little of a protein called Emx2. This protein is known to determine the size and position of the primary visual area (commonly called V1) during embryonic development. In the mutant mice with too much Emx2, the primary visual area was about 150% larger than it should be, even though the neocortex was a normal size. Zembrzycki, Stocker et al. then went on to show that the higher order areas associated with the primary visual area also grew proportionally larger in these mutant mice. The opposite was true for mice that didn’t produce Emx2 in their brains, and these mice had a much smaller primary visual area than normal mice.
Together, these findings reveal a previously unknown linear relationship between the size of the primary visual area and higher order visual areas that is controlled by the genes that pattern the neocortex during development. This and other new insights will inform future studies of the development and organization of the neocortex and our understanding of how it evolved.
https://doi.org/10.7554/eLife.11416.002Introduction
The mouse neocortex is patterned into functionally distinct fields that include the primary sensory areas (visual, somatosensory and auditory), which receive modality-specific sensory inputs from thalamocortical axons (TCAs) originating from nuclei of the dorsal thalamus (O'Leary et al., 2013). In the cortex, the connections of TCAs establish precise topographic representations (or maps) of the sensory periphery (Krubitzer and Kaas, 2005; O'Leary et al., 2013). Primary areas are flanked by higher order sensory areas (HO), which are interconnected with them and also contain topographic maps (Felleman and Van Essen, 1991). In mammals, this evolutionarily conserved general layout of the intra-areal neural circuits is responsible for the orderly progression of sensory information, sensory perception and the integration of higher cortical functions (Felleman and Van Essen, 1991; Geschwind and Rakic, 2013; Krubitzer and Kaas, 2005; Laramée and Boire, 2014; O'Leary et al., 2013). Disrupted layouts of cortical area layouts appear to be associated with neurodevelopmental disorders including autism (Courchesne et al., 2011; Voineagu et al., 2011). Studies of cortical arealization, the mechanisms that pattern the neocortex into areas, have focused almost exclusively on the primary areas and have led to the prevailing model that genetic mechanisms intrinsic to the neocortex control arealization during early cortical development (Greig et al., 2013; Krubitzer and Kaas, 2005; O'Leary et al., 2013). For example, the graded expression of the homeodomain transcription factor Emx2 in neocortical progenitors determines the size and position of the primary visual area (V1) in mice (Bishop et al., 2000; Hamasaki et al., 2004). Although higher order areas outnumber primary areas by roughly 10-fold (Marshel et al., 2011; Wang and Burkhalter, 2007), mechanisms that specify them and define their proportions relative to primary areas have yet to be explored.
Results
To investigate the impact of altered primary area size on higher order areas, we have used the cortical visual area V1 as a model. Previous studies have shown that genetic manipulation of patterning genes, including Fgf17 and Emx2, results in altered V1 size (Cholfin and Rubenstein, 2007; Hamasaki et al., 2004). Here we have analyzed transgenic mice that overexpress Emx2 (ne-Emx2) and show area patterning defects including a V1 that is ~150% of the normal size, while retaining overall normal neocortex size (Hamasaki et al., 2004 Leingärtner et al., 2007). By revealing the targeting patterns of TCAs projecting from thalamic sensory nuclei into the cortex (Fujimiya et al., 1986), the perimeters of primary sensory areas and the border between the neocortex and entorhinal cortex (ECT) can be visualized by serotonin (5HT) staining using a single postnatal day (P) 7 tangential section of the flattened cortical hemisphere (Figure 1A). The staining shows that, in addition to the previously reported enlarged V1 (Hamasaki et al., 2004,; Leingärtner et al., 2007), the cortical tissue that is nested between V1 and the surrounding primary areas (primary somatosensory cortex: S1, auditory areas: Aud) and the ECT laterally appears qualitatively larger in ne-Emx2 brains, when compared to wildtype (wt) sections (Figure 1A). We have defined this caudal cortical territory, which lies outside of V1, S1, and the auditory areas and shows no or weak 5HT staining as a joint higher order cortical area complex and have termed it HO-5HT. The 5HTstaining revealed that this region contains the higher order visual areas surrounding V1 (Wang and Burkhalter, 2007), the retrosplenial cortex (RSC) medially, and the ventral posterior temporal cortex laterally. The accurate distribution of staining across cortical layers can only be estimated using tangential sections. However, using P7 sagittal section, we confirmed in layer 4 that the caudal 5HT-positive cortical area (corresponding to V1) and the anteriorly adjacent 5HT-negative area between V1 and S1 (corresponding to HO-5HT) appears larger in ne-Emx2 brains than in wt ones (Figure 1B). Next, we labeled TCAs projecting to V1 by filling the dLG with crystals of the lipophilic neuronal tracer DiI. On the P7 sagittal sections that were derived from five different medial to lateral levels, anterograde DiI labeling in the cortex revealed that TCAs from the dLG terminate in a smaller region in wt than in ne-Emx2 brains (Figure 1—figure supplement 1). Across genotypes, the DiI staining revealed a sharp border with adjacent cortical tissues that did not receive TCAs input from the dLG (Figure 1—figure supplement 1). This finding is consistent with the 5HT staining and indicates a well-defined border between V1 neighboring higher order areas that is anteriorly shifted in ne-Emx2 brains.
Cortical areas can also be distinguished by area-specific gene expression patterns, which overlap with anatomical area borders and shift similarly when area patterning is disrupted (O'Leary et al., 2013). For example, Rorb expression is strongly induced by thalamic input to primary areas (Jabaudon et al., 2012) like S1 and V1 but is low in areas that do not receive their major inputs from the principal thalamic sensory nuclei, such as cortical higher order areas surrounding V1 (Chou et al., 2013; Wang and Burkhalter, 2007). In situ hybridization (ISH) on sagittal sections adjacent to 5HT-stained ones revealed sharp Rorb gene expression borders between areas in layer 4. Notably in ne-Emx2 brains, the high-to-low Rorb expression border is located more anteriorly, and the area showing low Rorb expression and resembling HO-5HT (Chou et al., 2013) is larger than in wt sections (Figure 1B). This reveals that characteristic molecular markers that delineate the borders between V1 and surrounding higher order areas remain expressed at normal levels, but their sharp expression borders shift anteriorly in ne-Emx2.
Projection neurons in layer 5, which extend axons into subcortical targets, are similarly determined by a molecular code (Greig et al., 2013). We therefore predicted that the areal shifts in ne-Emx2 brains would be accompanied by corresponding changes in layer 5 output projections. We labeled two distinct types of layer 5 subcerebral projection neurons by inserting DiI crystals either into the superior colliculus, which labels corticotectal projections from V1 and HO, or else into the pyramidal decussation, which labels corticospinal projections from the frontal cortex and S1 (Greig et al., 2013; Zembrzycki et al., 2015). We found that layer 5 corticotectal projections extended more anteriorly in ne-Emx2 sagittal sections. Vice versa, the layer 5 corticospinal projections extended less posteriorly in ne-Emx2 brains. These staining patterns are consistent with an altered balance of projection neuron identity in layer 5 (Greig et al., 2013; Zembrzycki et al., 2015) and an overall expansion of visual areas, demonstrating that areal patterning changes in ne-Emx2 brains are not limited to the cortical layers that receive thalamic input. These findings complement previous reports describing Emx2 patterning functions (Bishop et al., 2000; Hamasaki et al., 2004; Leingärtner et al., 2003; Leingärtner et al., 2007) and indicate for the first time that V1 and higher order area sizes are altered concomitantly in ne-Emx2 brains at the level of area-specific connectivity and gene expression in multiple cortical layers. Taken together, our results suggest that changes in primary area size are paralleled by similar changes in higher order area size.
It is commonly assumed that areal patterning changes also alter area-specific functional neuronal properties and topographic sensory maps, but this has never been demonstrated conclusively. Therefore, to compare functional neuronal properties of an enlarged visual cortex to a normal-sized one, we have used Fourier intrinsic signal optical imaging to construct topographic visual response maps to light bars that were moved across the visual field of the retina (up and down: elevation maps; left to right: azimuth maps) (Kalatsky and Stryker, 2003). Visual responses in V1 of wt and ne-Emx2 mice produced intrinsic signal maps that were indistinguishable in strength, and the axes of azimuth and elevation were organized in the same way in all tested brains (Figure 2A), revealing that functional topographic organization of the visual cortex was intact. However, the representations of elevation and azimuth were expanded in ne-Emx2 animals, and their retinotopic maps were overall larger (elevation: 138% ± 8.7% of wt; azimuth: 143% ± 8.2% of wt). For example, the green region in the response maps representing ~20 to ~30 degrees of elevation/azimuth is clearly enlarged in ne-Emx2 brains, compared to wt brains (Figure 2A). To investigate the relationship between the location and size of the V1 functional response area and the histochemically delineated V1, as indicated by 5HT staining, multiple injections of DiI were placed lining the border of the V1 optical response map after the imaging procedure. On 5HT-stained, flattened tangential cortical sections, the DiI injection sites were found in all cases to be located near the border of the 5HT staining in V1, confirming the overall enlarged V1 perimeters in ne-Emx2 brains compared with the wt brains (Figure 2B). This shows that the 5HT-stained V1 area accurately corresponds to the intrinsic functional V1 map, suggesting that enlarged HO in ne-Emx2 brains have not acquired ectopic V1-like functional properties.
To further characterize the shifted border between visual areas in ne-Emx2 animals, we used additional neuronal tracing approaches. Stereotypically, V1 is connected with the dLG, whereas the HO areas are wired to the posterior thalamic nucleus (PO) (Leyva-Díaz and López-Bendito, 2013; López-Bendito and Molnár, 2003). To first label axonal connections between the cortex and the thalamus in wt brains, we administered dual tracer injections into locations that approximate to HO (DiD: green dye) and another injection around the approximated border area between V1 and HO (DiI: red dye). After diffusion of the tracers, we performed 5HT staining on flattened cortex sections to identify the areas in which the injections were administered and analyzed the patterns of retrograde dye labeling on coronal sections of the thalamus. In representative cases (Figure 2C) where 5HT staining confirmed that DiI was injected at the border between V1 and HO and DiD was injected into HO, retrogradely labeled green DiD cells were apparent in the dLG and the PO. Conversely, red DiI cells were only labeled in the dLG. In ne-Emx2 brains, we administered similar dual tracer injections: DiI was targeted to V1 and a DiD injection was administered around the approximate border between V1 and HO. Due to their enlarged V1, all ne-Emx2 injections were administered at more anterior coordinates than in wt brains (compare dashed lines in Figure 2C). In representative cases (Figure 2C) where 5HT staining confirmed that the DiI injection was administered into V1 and DiD was injected around the V1/HO border, red cells were found in the dLG, whereas green-labeled cells were apparent in the PO and the dLG. The dual tracings in wt and ne-Emx2 brains are consistent with the predicted connectivity of cortical neurons around the injection sites (Leyva-Díaz and López-Bendito, 2013; López-Bendito and Molnár, 2003). Although located more anteriorly in ne-Emx2 brains, the subcortical connectivity patterns around the V1/HO borders were similar, demonstrating that these neurons show connectivity patterns that are consistent with their intrinsic areal identity and not their topographic location on the cortical sheet. Taken together, our results indicate that increased V1 size in ne-Emx2 brains is accompanied by a concomitant enlargement and anterior shift of HO.
To define the individual magnitudes of the V1 and HO size increases in ne-Emx2 brains, we next used gene expression domains as molecular markers delineating visual areas and quantified them (Figure 3). An accurate assessment of area sizes using flattened and/or sectioned cortical tissues could potentially be hampered by imperfect flattening of the tissues or by cutting artifacts. Therefore, we have used RNA in situ hybridization on intact whole brains (whole mount in situ hybridization: WMISH) at P7 for quantification purposes, which has the advantage that quantifications can be made using single images without sectioning and artifacts that may arise from such tissue processing. We first used a set of two marker genes, Unc5d and Igfbp5, whose expression delineates V1 at P7 (Chou et al., 2013). The gene expression domains on WMISH-stained brains were outlined and their sizes quantified as a measure of V1 area size. The mean value of wt brains was defined as 100% and the area size percentages of ne-Emx2 brains displayed accordingly as ‘percent of wt’ (Figure 3A). V1 gene expression domains of both markers were larger in ne-Emx2 brains (Unc5d-V1: 148% ± 6.1%, p = 0.0003; Igfbp5-V1: 148 ± 4.5%, p < 0.0001). The magnitude of the increased in V1 size labeled genetically in ne-Emx2 brains is comparable to V1 area measurements derived from 5HT-stained P7 flattened cortical sections (Figure 3A: 5HT-V1: 142.1% ± 3.1%, p < 0.0001), indicating that molecular markers on whole brains reliably delineate V1 and can therefore be used to quantify and compare area sizes between samples and mouse lines.
We next have used additional markers to quantify higher order area sizes (Figure 3B). Previous studies have parsed higher order visual areas using neuroanatomical tracers (Wang and Burkhalter, 2007) (see also schematic in Figure 3) and have revealed genes that are expressed at different levels in V1 and higher order visual areas (Chou et al., 2013). For example, Cdh8 and Lmo4 expression is higher in the area surrounding V1, where higher order visual areas are located (Chou et al., 2013; Marshel et al., 2011; Wang and Burkhalter, 2007). The domains of high Cdh8 and Lmo4 expression appear to label higher order visual areas uniformly (Chou et al., 2013), without revealing subdivisions between them (compare schematic in Figure 3 showing approximate location and outline of higher order visual areas as identified by Wang and Burkhalter, 2007 to Cdh8 and Lmo4 gene expression domains around V1). On P7 WMISH-stained wt brains, we quantified the V1 and HO sizes in the medial cortex on the basis of low and high gene expression domains (see dotted lines in Figure 3, Figure 3—figure supplement 2). Anatomically, these gene expression domains surrounding V1, which show much stronger staining compared with V1, extend anteriorly up to the S1 border, laterally to the border of the auditory areas and the ECT and medially up to the border to the RSC, respectively (Figure 3—figure supplements 1,2). Hence, compared with the above-mentioned HO complex that was identified using 5HT staining (Figure 1), the higher order area complex labeled by Cdh8 and Lmo4 relates to a smaller cortical region that more closely relates to higher order visual areas, but excludes the RSC. The size (Cdh8-V1: 145.3 ± 7.4%, p < 0.0001; Lmo4-V1: 146.6% ± 6.7%, p = 0.0005) and shape of the gene expression domains in V1 were similar in Cdh8- and Lmo4-stained brains. Similarly, the gene expression domains nested around V1 largely overlapped between the two probes. The only apparent difference between them is around the anteromedial edge of the higher order visual areas (Wang and Burkhalter, 2007), where Cdh8 is expressed at much lower levels compared to more lateral regions around V1 across genotypes (arrowheads in Figure 3B). The wt values of the measurements were again defined as 100%. The overall shapes of the two HO marker gene domains were similar and the sizes larger in ne-Emx2 brains compared with wt brains (Cdh8- HO: 145.7 ± 6.4%, p = 0.0015; Lmo4- HO: 144.9 ± 3.8%, p = 0.00157). The analysis of different area-specific sets of marker genes, either showing unique expression in V1, or discernable expression levels between visual areas, revealed an increase in visually-related HO in ne-Emx2 brains that was proportionate to the V1 size increase. The extrastriate areas that we have identified on the basis of 5HT staining (Figure 1) included the RSC, which is not a higher order visual area (Garrett et al., 2014; Marshel et al., 2011; Vann et al., 2009; Wang and Burkhalter, 2007) raising the possibility that only related primary and higher order areas (e.g. vision) could scale proportionately. To test this possibility, we have used WMISH of Lypd1 on P7 wt and ne-Emx2 brains as a specific marker labeling the caudomedial cortex, where the RSC is located (Figure 3—figure supplement 1). We found that the specific Lypd1 gene expression domain in the caudomedial cortex is significantly enlarged in ne-Emx2 brains (114.3 ± 5.2%, p = 0.0225, n = 4), compared with wt brains. This size increase is not proportionate to the size increases of V1 and the higher order visual area complex labeled by Cdh8 and Lmo4 in ne-Emx2 brains (Figure 3) suggesting that increased V1 size is specifically accompanied by a proportionate size increase of related higher order visual areas.
To test if related HO size matches V1 size only when it is larger than normal, or if primary area size bi-directionally is accompanied by according scaling of related higher order areas, we next analyzed HO sizes excluding the RSC in brains with a smaller than normal V1 (Figure 4). Constitutive Emx2 mutant mice have an overall smaller brain and visual cortex, but homozygous mutants die perinatally (Bishop et al., 2000), preventing the analysis of cortical areas, which arise at later stages. To overcome this limitation, we generated a novel mouse line with floxed Emx2 alleles (Figure 4—figure supplement 1), allowing conditional inactivation of Emx2. We crossed Emx2 floxed mice with Emx1-IRES-Cre expressing mice (Gorski et al., 2002) to generate conditional, cortex-specific Emx2 mutant mice. These cKO mice are viable, fertile and have an anatomically normal neocortex (Figure 4—figure supplement 2). Confirming the prediction that reduced Emx2 expression levels in cortical progenitors would lead to smaller visual areas, cKO brains show areal patterning changes (Figure 4—figure supplement 3) that are similar to those previously reported in heterozygous Emx2 mutant brains (e.g. larger frontal cortex) (Hamasaki et al., 2004), but are opposite to those apparent in ne-Emx2 brains (e.g. smaller frontal cortex) (Hamasaki et al., 2004; Leingärtner et al., 2007). As in ne-Emx2 brains, V1 in cKO brains was characterized using 5HT staining and DiI injections into the dLG (Figure 4—figure supplement 4), revealing that V1 in cKO is greatly reduced in size.
To complement the quantification of V1 and its related HO sizes in ne-Emx2 brains, we next have used cKO brains to perform WMISH with both sets of marker genes noted above (Figure 3) and measured their sizes (Figure 4). Measurements of the molecular V1 marker domains (Unc5d-V1: 67.1 ± 2.6%, p = 0.00029; Igfbp5-V1: 68 ± 3.3%, p = 0.00017), as well as the V1 expression domains of Cdh8 (67.7 ± 4.5%, p < 0.0001) and Lmo4 (68.6 ± 0.8%, p < 0.0001) revealed that the molecularly defined V1 in cKO was smaller than in wt brains. These reductions matched the reduced V1 in 5HT-stained flattened cortex sections (Figure 4A: 5HT-V1: 68.5 ± 3%, p < 0.0001). Subsequent quantification of the gene expression domains of Cdh8 and Lmo4 surrounding V1 (Figure 4B: Cdh8- HO: 66.8 ± 2.9%, p < 0.0001; Lmo4- HO: 69.3 ± 2.4%, p < 0.0001) revealed that the cortical region that contains visually-related HO was also reduced in cKO brains to a degree proportional to the reduction in V1 size. These data demonstrate that when V1 size is reduced, related HO size is reduced to a similar extent.
In order to reveal a correlation between primary and related higher order area size between brains with larger and smaller than normal visual areas, we calculated the ratios between the genetically defined V1 and related HO sizes (Figure 4C). Ranging over an ~80% variation of the normal V1 size, Figure 4C reveals that related HO size is bi-directionally altered in a linear fashion (Cdh8 regression: y = −0.0071x + 1.7463; Lmo4 regression: y = 0.0208 + 1.7463). Taken together, our results are consistent with a proportional scaling relationship between the size of primary and related higher order visual areas: The size of V1 is determined by the activity of transcription factors including Emx2 during development, and this mechanism likewise controls the linear matching of the proportions of higher order visual areas in the mouse neocortex.
Discussion
The present findings address the mechanisms that specify and regulate the size of higher order sensory areas, an issue that has been largely neglected. They reveal a novel, prominent role for intrinsic genetic information in this process. Genetically altering the size of V1 over a range of ~80% of its normal size using a Emx2 gain-of function mouse line and a novel conditional Emx2 loss-of function mouse line showed that the specification of both primary and related higher order cortical areas during development was linearly scaled by driving the unique properties that characterize both, V1 and higher order visual areas.
Regardless of whether V1 was larger or smaller than in wt mice, related HO exhibited normal cytoarchitecture, genetic profiles, functional properties, and characteristic patterns of connectivity that resulted in an overall uniformly altered ‘visual cortical field’ in the occipital cortex that remained accurately and proportionally subdivided into V1 and higher order visual areas. This demonstrates that Emx2 (and perhaps additional intrinsic area patterning regulators) specify a ‘sensory cortical field’ that includes primary and higher order areas and a defined border between them. This model of cortical area patterning is not consistent with the possibility that the core properties of primary and higher order areas are specified sequentially or through parallel genetic mechanisms.
Our results further reveal that higher order areas do not have a fixed size. Rather their relative size is flexible. By using mouse models with bi-directional changes of V1 size as a model, our study revealed that higher order areas scale linearly together with their related primary sensory areas, This observation is important for at least two reasons: (i) it re-emphasizes a sequential model of primary sensory area formation that likewise influences the properties of related higher order areas. In this model, cortical intrinsic mechanisms specify all generic primary and higher order visual cortex properties during early development. Much later during postnatal development, geniculocortical input is needed to terminally differentiate the genetic distinctions between V1 and HO (Chou et al., 2013; Vue et al., 2013). (ii) It contradicts the hypothesis that cortical structure/function evolution mainly is driven by a disproportionate increase in the size of related higher order areas relative to primary areas. To the contrary, our results show that primary and related higher order areas remain proportionate when primary area size is altered through genetic mechanisms, suggesting that an increase in the complexity of connections and micro-circuits among higher order cortical processing centers likely accounts for gains in cortical functions that are characteristic for gyrencephalic mammals with larger cortical surface areas, compared to simpler lissencephalic mammals. In summary, the newly discovered linear scaling relationship between primary and related higher order areas has major implications for the basic understanding of the development and organization of the neocortical bauplan and its evolution and variability in normal and affected conditions.
Materials and methods
Mouse lines and conditional Emx2 gene targeting
Request a detailed protocolAll experiments were approved and conducted following the guidelines of the Institutional Animal Care and Use Committee at the Salk Institute and were in full compliance with the guidelines of the National Institutes of Health for the care and use of laboratory animals. When mice were mated, the morning of the identified vaginal plug was designated as E0.5. The morning on which pups were born was designated P 0.5. Transgenic mice overexpressing Emx2 under the Nestin promoter (ne-Emx2) were previously described (Hamasaki et al., 2004). For generating Emx2 floxed mice (Emx2fl/fl), gene targeting was carried out using homologous recombination in embryonic stem cells. A targeting construct was designed in which the 5’ loxP site was upstream of the Emx2 transcriptional start site and the 3’ loxP site downstream of Exon 2, followed by a FRT-site-flanked PGK-Neo cassette, Figure 4—figure supplement 1). Targeted embryonic stem cell clones were screened by Southern blot with probes A, B, and C and by PCR to identify Emx2floxed-neo/+ clones (Figure 4—figure supplement 1). Positive clones were injected into C57BL/6J blastocysts at the Salk Transgenic Core Facility and chimeras were mated to C57BL/6J females to obtain germline transmission. Heterozygous mice were mated with mice expressing FLPe (Rodríguez et al., 2000) to remove the neo cassette and then mated to obtain homozygous Emx2fl/fl mice. Cortex specific deletion of Emx2 (cKO) was obtained by crossing Emx2fl/fl mice with Emx1-IRES-Cre mice (Gorski et al., 2002). Specificity of Emx1-IRES-Cre-mediated deletion of Emx2 floxed alleles was analyzed by WMISH (described below) staining using a full-length Emx2 antisense RNA probe on E11 embryos. Genotyping was performed using primers for Emx2 floxed alleles (Emx2 forward: GAC-TCC-TTT-CCC-AAA-TAA-CCC-C, Emx2 reverse: GTA-AGC-GGG-CGG-GGA-CTG-GTT-C) and for the Cre recombinase (cre forward: GCT-AAA-CAT-GCT-TCA-TCG-TCG-G, cre reverse: GAT-CTC-CGG-TAT-TGA-AAC-TCC-AGC), and the ne-Emx2 transgene (nestin forward: TCA-ACC-CCT-AAA-AGC-TCC, Emx2 reverse: GGA-CGG-AGA-GAA-GGC-GGT).
In situ hybridization, immunostainings, and tangential cortical sections
Request a detailed protocolTissues were dissected, washed in phosphate-buffered saline (PBS), fixed overnight in 4% phosphate-buffered paraformaldehyde (PFA), washed in PBS, and cryopreserved in 30% sucrose in PBS. Postnatal brains were perfused with PFA, postfixed overnight in PFA, washed with PBS, and cryopreserved in 30% sucrose in PBS. Tissues were embedded in Tissue-Tek OCT (Sakura Finetek , Japan) and sectioned on a cryostat (Leica, Germany). Antisense RNA probes were labeled using a DIG-RNA labeling kit (Roche, Switzerland). ISH on 18-μm cryostat sections and WMISH using P7 brains were carried out as previously described (Chou et al., 2013; Hamasaki et al., 2004; Zembrzycki et al., 2007; Zembrzycki et al., 2015). For tangential cortical sections, cortical hemispheres were dissected, flattened, postfixed between slide glasses, and then cryoprotected. Tangential sections were cut into 40-μm slices from flattened cortical hemispheres with a sliding microtome and then they were immunostained for Serotonin (5HT, ImmunoStar, Hudson, WI). Immunostaining was developed using the diaminobenzidine colorimetric reaction and the Vectastain kit (Vector, Burlingame, CA). For Nissl staining, sections were stained with 0.5% cresyl violet and then dehydrated with graded alcohols.
Axon tracings
Request a detailed protocolLipophilic tracers DiI and DiD (all from Molecular Probes, Eugene, OR) were used to label corticothalamic-, thalamocortical-, corticotectal-, and corticospinal projections. For each experiment 4–6 brains with comparable tracer injection sites were cut and used for further data analysis, representative example images are shown in the figures. Analysis of thalamocortical axons by thalamic DiI injections (Figure 1B): P7 brains were fixed in 4% PFA, hemisected, and a coronal cut between the diencephalon and mesencephalon was made in order to expose thalamic nuclei at the section surface and DiI crystals were implanted to cover the dorsal lateral geniculate nucleus (dLG). After incubation for 1 to 2 months at 30°C to 60°C, preparations were sectioned sagittally on a vibratome (Leica). Sections were counterstained with DAPI (Vector) and analyzed under a fluorescence microscope to determine the tangential distribution of labeled thalamocortical axons in the neocortex.
Analysis of layer 5 subcerebral projection neurons (Figure 1C): Corticospinal neurons in cortical layer 5 were retrogradely labeled by inserting DiI crystals into the pyramidal decussation in 4% PFA fixed brains. Layer 5 corticotectal neurons were labeled in 4% PFA fixed brains by implantation of small DiI crystals into the upper layers of superior colliculus. Brains were incubated at 37°C for 2–3 months before 100 μm sagittal vibratome sections were cut and analyzed under fluorescent light. Analysis of area-specific thalamocortical and corticothalamic connectivity of caudal cortex (Figure 2): P7 pups were anesthetized by hypothermia and a small area of skull was removed to expose the cortical surface. DiI crystals and a small piece of DiD were implanted into cortical locations around V1 and the V1/HO border. After 1 day of survival, brains were removed after 4% PFA perfusion and their cortices and thalami dissected. Cortices were then flattened and stained for 5-HT to reveal primary sensory areas relative to the dye injection sites. The thalami were preserved sectioned coronally, stained with DAPI and DiI/DiD labeled cells analyzed under a fluorescence microscope.
Statistical analysis, area measurements, and intrinsic signal optical imaging
Request a detailed protocolData collection and analyses were performed blind to genotype and the conditions of the experiments, data were collected and processed randomly, and no data points were excluded. No statistical methods were used to predetermine sample sizes, but our sample sizes were similar to those reported in previous publications (for example, (Chou et al., 2013; Zembrzycki et al., 2013). Data met the assumptions of the statistical tests used, and the data distribution was assumed to be normal but was not formally tested. Statistics were calculated with Microsoft Excel. Quantifications show mean values of the tested groups and are displayed as a percentage of the wt group. Quantified sample sizes (number of brains: n) are indicated in the figure legends. The examples shown in each figure are representative and were reproducible for each set of experiments. Individual experiments were successfully repeated at least three times using different litters.
Area size measurements on 5HT stained sections and statistical analysis was performed as previously described (Leingärtner et al., 2007; Zembrzycki et al., 2015; 2013). To quantify molecular V1 and HO sizes, gene expression domains were quantified on single images of WMISH-stained brains (examples of measured area outlines are shown as dashed lines in Figures 3,4) using ImageJ (Rasband 1997−2013). Derived wt mean values were defined as 100% and values of the other mouse lines calculated accordingly. Statistical significance was determined using unpaired two-tailed t test, p values < 0.05 (indicated as *) were considered as statistically significant. Variance is indicated in the main text sections reflecting standard error of the mean. Intrinsic signal optical imaging was performed as previously described (Kalatsky and Stryker, 2003). To determine the spatial relationship between V1 functional maps and V1, defined histochemically by 5-HT staining, animals were imaged to determine the V1 map, and after completion, small DiI injections were made outside of the functional map perimeters. Animals were then perfused with 4% PFA, cortices dissected, flattened, sectioned tangentially, and stained for 5-HT.
References
-
Patterning of frontal cortex subdivisions by Fgf17Proceedings of the National Academy of Sciences of the United States of America 104:7652–7657.https://doi.org/10.1073/pnas.0702225104
-
Postnatal development of serotonin nerve fibers in the somatosensory cortex of mice studied by immunohistochemistryThe Journal of Comparative Neurology 246:191–201.https://doi.org/10.1002/cne.902460205
-
Topography and areal organization of mouse visual cortexJournal of Neuroscience 34:12587–12600.https://doi.org/10.1523/JNEUROSCI.1124-14.2014
-
Cortical excitatory neurons and glia, but not GABAergic neurons, are produced in the Emx1-expressing lineageThe Journal of Neuroscience : The Official Journal of the Society for Neuroscience 22:6309–6314.
-
Molecular logic of neocortical projection neuron specification, development and diversityNature Reviews. Neuroscience 14:755–769.https://doi.org/10.1038/nrn3586
-
ROR induces barrel-like neuronal clusters in the developing neocortexCerebral Cortex 22:996–1006.https://doi.org/10.1093/cercor/bhr182
-
New paradigm for optical imaging: temporally encoded maps of intrinsic signalNeuron 38:529–545.
-
The evolution of the neocortex in mammals: how is phenotypic diversity generated?Current Opinion in Neurobiology 15:444–453.https://doi.org/10.1016/j.conb.2005.07.003
-
Visual cortical areas of the mouse: comparison of parcellation and network structure with primatesFrontiers in Neural Circuits 8:149.https://doi.org/10.3389/fncir.2014.00149
-
Cortical area size dictates performance at modality-specific behaviorsProceedings of the National Academy of Sciences of the United States of America 104:4153–4158.https://doi.org/10.1073/pnas.0611723104
-
Thalamocortical development: how are we going to get there?Nature Reviews Neuroscience 4:276–289.https://doi.org/10.1038/nrn1075
-
BookArea Patterning of the Mammalian CortexIn: O'Leary DDM, Stocker AM, Zembrzycki A, John Rubenstein, Pasko Rakic, editors. Patterning and Cell Type Specification in the Developing Cns and Pns. Oxford: Elsevier. pp. 61–85.https://doi.org/10.1016/B978-0-12-397265-1.00021-6
-
High-efficiency deleter mice show that FLPe is an alternative to cre-loxPNature Genetics 25:139–140.https://doi.org/10.1038/75973
-
What does the retrosplenial cortex do?Nature Reviews Neuroscience 10:792–802.https://doi.org/10.1038/nrn2733
-
Thalamic control of neocortical area formation in miceJournal of Neuroscience 33:8442–8453.https://doi.org/10.1523/JNEUROSCI.5786-12.2013
-
Area map of mouse visual cortexThe Journal of Comparative Neurology 502:339–357.https://doi.org/10.1002/cne.21286
-
Sensory cortex limits cortical maps and drives top-down plasticity in thalamocortical circuitsNature Neuroscience 16:1060–1067.https://doi.org/10.1038/nn.3454
-
Postmitotic regulation of sensory area patterning in the mammalian neocortex by Lhx2Proceedings of the National Academy of Sciences of the United States of America 112:6736–6741.https://doi.org/10.1073/pnas.1424440112
Article and author information
Author details
Funding
National Institutes of Health (NS31558)
- Dennis DM O'Leary
National Institutes of Health (MH086147)
- Dennis DM O'Leary
National Institutes of Health (EY02874)
- Michael P Stryker
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Kevin Jones for providing Emx1-IRES-CRE mice. This work was supported by the Vincent J Coates Chair of Molecular Neurobiology at the Salk Institute for Biological Sciences (DDMO).
Ethics
Animal experimentation: All experiments were approved under Protocol #09-012 and conducted following the guidelines of the Institutional Animal Care and Use Committee at the Salk Institute and were in full compliance with the guidelines of the National Institutes of Health for the care and use of laboratory animals.
Copyright
© 2015, Zembrzycki 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
-
- 2,210
- views
-
- 484
- downloads
-
- 13
- 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
- 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.
-
- Developmental Biology
- Genetics and Genomics
We present evidence implicating the BAF (BRG1/BRM Associated Factor) chromatin remodeler in meiotic sex chromosome inactivation (MSCI). By immunofluorescence (IF), the putative BAF DNA binding subunit, ARID1A (AT-rich Interaction Domain 1 a), appeared enriched on the male sex chromosomes during diplonema of meiosis I. Germ cells showing a Cre-induced loss of ARID1A arrested in pachynema and failed to repress sex-linked genes, indicating a defective MSCI. Mutant sex chromosomes displayed an abnormal presence of elongating RNA polymerase II coupled with an overall increase in chromatin accessibility detectable by ATAC-seq. We identified a role for ARID1A in promoting the preferential enrichment of the histone variant, H3.3, on the sex chromosomes, a known hallmark of MSCI. Without ARID1A, the sex chromosomes appeared depleted of H3.3 at levels resembling autosomes. Higher resolution analyses by CUT&RUN revealed shifts in sex-linked H3.3 associations from discrete intergenic sites and broader gene-body domains to promoters in response to the loss of ARID1A. Several sex-linked sites displayed ectopic H3.3 occupancy that did not co-localize with DMC1 (DNA meiotic recombinase 1). This observation suggests a requirement for ARID1A in DMC1 localization to the asynapsed sex chromatids. We conclude that ARID1A-directed H3.3 localization influences meiotic sex chromosome gene regulation and DNA repair.