Mouse genetic studies unravel the implications of multiple main drivers underlying neuroanatomical phenotypes at the 16p11.2 locus in a sex-specific manner
To identify which gene(s) regulate(s) mammalian brain architecture at the de novo 16p11.2 deletion locus, we aim to assess the alterations in NeuroAnatomical Phenotypes (named NAPs by us [36]), independently in male and female heterozygous (het) knockout (KO) mice, in the autism-associated B6N Del(7Sult1a1-Spn)6Yah mouse model of the entire 16p11.2 deletion (hereafter named as Del/+) as well as individually for each of the 30 protein-coding genes within the interval.
We developed, or acquired through collaborations, 243 adult mutant and 221 colony-matched littermate wild-type (WT) mice, corresponding to the deletion of the entire locus as well as 20 unique gene deletions (highlighted in red in Fig. 1A), each studied with an average of four biological replicates, a number which we previously showed suitable to detect NAPs with effect size superior to 5% [36]. For the remaining 10 genes (in black in Fig. 1A), the germline transmission of the mutation failed despite multiple attempts or no mouse model was available during the course of the study. For one gene of interest (Kctd13), multiple allelic strategies were used. To ensure high comparability between the results, mouse mutants assessed in this study were all processed on identical genetic backgrounds (C57BL/6) and at the same age using the same standardized pipelines [37,38,39,40]. A detailed description of study samples and allelic constructions is provided in Additional file 2: Table S1 and Additional file 9.
Mouse neuroanatomical studies for the identification of NAP genes at the 16p11.2 autism-associated locus. A Schematic representation of the human 16p11.2 region showing gene content and order with genes that underwent mouse neuroanatomical studies in red. B At bregma + 0.98 mm and bregma − 1.34 mm, 22 and 45 parameters were measured, respectively. C These 67 parameters (listed in Additional file 2: Table S2) are grouped into five categories (brain size, ventricle, cortex, commissure, and subcortex) on the two coronal sections. D Bar plots of the percentage change of male Del/+ relative to male WT mice for each of the parameters measured. The colored region indicates the presence of a significant parameter at the 0.05 level. White indicates a p-value > 0.05. E Bar plots of percentage change of female Del/+ relative to female WT. F Left: representative brain image stained with Nissl-luxol of adult male WT and Del/+ mice showing the cortical layers of the somatosensory cortex at bregma − 1.34 mm. Right: histograms showing the thickness of the upper and lower layers of the somatosensory cortex in male Del/+ mice. G Histograms showing the thickness of the upper and lower layers of the cortex in female Del/+ mice. H Venn diagram illustrating NeuroAnatomical Phenotype (NAP) genes (genes whose mutations yield neuroanatomical phenotypes) in male mice positioned in each category. Phenotypic directionality is color-coded. Blue corresponds to a reduction in the size of the brain regions affected, red to increase, and yellow occurring of both. I, J Bar plots showing the number of NAPs per section (bregma + 0.98 mm in orange and bregma − 1.34 mm in blue) for each gene assessed for male (I) and female (J), respectively. Genes are listed on the x-axis and sorted according to the number of NAPs
Using a highly robust approach for the assessment of 67 coronal neuroanatomical parameters per brain sample (Fig. 1B and Additional file 1: Fig. S1), described in detail elsewhere [38], we systematically quantified the same two coronal brain sections at bregma + 0.98 mm and bregma − 1.34 mm and collected neuroanatomical measurements blind to the genotype (Additional file 3: Table S7). These parameters were grouped into five main categories: brain size, commissure, ventricle, cortex, and subcortex (Fig. 1C and Additional file 2: Table S2). After multiple quality control steps and critical evaluation of each phenotype, gene association was carried out within our internal database using a standardized statistical pipeline (Additional file 9). Heat maps of the assessed genes comprising percentage change relative to WTs and p-value are provided in Additional file 1: Fig. S2 for males, Additional file 1: Fig. S3 for females, and Additional file 4: Table S8.
Overall, the neuroanatomical profile associated with the male Del/+ mice indicated a decreased size of the majority of regions assessed, but only the somatosensory cortex reached a significance level and was reduced in size by 10% (p < 0.05) (Fig. 1D). We asked whether alterations in specific cortical layers may contribute to this phenotype and found that it stemmed from the upper II–IV cortical layers that were thinner by 10% (p < 0.0091) (Fig. 1F). We did not detect any significant NeuroAnatomical Phenotypes in the female Del/+ mice (Fig. 1E, G).
At the single-gene level, we identified 13 genes associated with NAPs (hereafter named NAP genes), when using a relaxed significance threshold of 0.05, associated with defects in commissure (ten genes), cortex (eight genes), subcortical structures (seven genes), brain size (five genes), and ventricle (three genes) in males (Fig. 1H). Eight genes (Bola2, Qprt, Maz, Mvp, Ppp4c, Slx1b, Taok2, and Zg16) gave significant results affecting two or more categories. Mvp was the only gene affecting all five main categories. The remaining five genes (Doc2a, Fam57b, Hirip3, Spn, and Tbx6) presented specific phenotypes in one brain category. In the 13 NAP genes, 38.5% decreased the size of the affected brain structures, while 38.5% increased their sizes, and 23% had bidirectional effects (genes in blue, red, and yellow, respectively, in Fig. 1H). When using a stringent significance threshold of 0.0001 to account for multiple comparisons, Mvp and Tbx6 remained significant NAP genes.
Figure 1I shows the number of NAPs for each of the 21 individually deleted alleles in males. While Mvp still stood out as the strongest candidate gene based on the number of affected parameters (n = 14; all decreased in size except the lateral ventricles that were enlarged in size when compared to littermate WT mice), Ppp4c showed ten NAPs all decreased in size, Zg16 implicated seven NAPs all increased in size, and Taok2 showed six NAPs with three increased and three decreased in size. A heat map is provided in Additional file 1: Fig. S2, and neuroanatomical phenotypes are described in detail for each gene in Additional file 10.
Next, sex differences were assessed to determine the impact of each mutation on the female brain. Overall, female mice displayed a reduced number of NAP genes (10 as opposed to 13 in males), and for each of the nine genes in common between males and females, there were less affected brain parameters in females (Fig. 1J). Directionality of the phenotypes was consistent between sexes, but NAPs were significantly different with an excess of parameters at bregma − 1.34 mm for females (p = 0.00006, Fisher test) (Fig. 1I, J). Gdpd3 was a female-only NAP gene, Bola2 showed more severe anomalies in females, and Mvp showed no anomalies in females (Additional file 1: Fig. S4A). Female NAPs are described for each gene in detail in Additional file 1: Fig. S3 and Additional file 10.
Four genes (Coro1a, Ino80e, Mapk3, and Kctd13) were categorized as non-NAP, both in males and females. Considering the discrepancies in the literature [20,21,22], heterozygous Kctd13+/− mice were assessed twice using two independent allelic constructions, which gave identical results with no NAPs (Additional file 9). We also studied the homozygous (hom) Kctd13−/− mice and found a reduction in the size of the hippocampus by 10% (p = 0.015) for males (Additional file 1: Fig. S2 and Additional file 1: Fig. S4C-D) and 6% (p = 0.014) for females (Additional file 1: Fig. S3). This reinforces the existing link between Kctd13 and hippocampal biology [21, 22].
To summarize, our systematic neuroanatomical screen demonstrates the implication of multiple main genes that drive NAPs at the 16p11.2 locus with the major vault protein, Mvp, gene being one of the top drivers. Our work also highlights the profound sex differences with NAPs being more prominent in males.
The major vault protein is expressed in the limbic system in both sexes
Focusing on Mvp, the one gene that gave the most widespread NAPs across multiple brain categories in males (no NAPs in females), we thought to establish its expression distribution independently in male and female WT mouse brains to determine whether Mvp sex-specific NAPs are correlated to sex-specific changes in MVP expression levels.
MVP transcripts, assessed at several developmental stages from embryonic day 16.5 (E16.5) to 30 weeks of age, were constant between males and females with a postnatal peak in line with an existing resource of expression profiles across multiple organs and developmental stages, which additionally revealed no expression of Mvp expression before E16.5 [41]. Expression was higher in peripheral tissues and the cerebellum than in the cortex and the hippocampus, both in males and females (Fig. 2A, B and Additional file 1: Fig. S5A-C).
MVP expression is highly specific to the limbic system. A, B Mvp qRT-PCR expression in the liver (LIV), cerebellum (CRB), hippocampus (HP), and cortex (CTX), in male (A) and female (B) WT mice. Normalization was done using the Gnas/Hprt ratio. Plots are shown as mean ± SEM. Immunofluorescence of MVP (green) in CRB (C), ventral HP (D), and deep layers of the cingulate gyrus (Cg) (E). Arrowheads point to MVP-positive cells. AV, arbor vitae; gr/mole, granular/molecular layer of CRB; or/py, oriens/pyramidal layer of HP. F–O Selection of MVP immuno-histo-fluorescence (IHF) images on the coronal brain sections, each of 30 µm in thickness, performed throughout the entire brain from bregma + 1.41 mm to bregma − 7.67 mm. MVP presence was revealed using an anti-vault antibody (N2-B15) and secondary antibodies fused to Alexa488 fluorochrome (green). Dashed lines delineate regions with MVP-positive cells (white), MVP-positive tract (yellow), and MVP-negative cells (gray). It is noteworthy to mention that the experiment was repeated on a new set of brain samples on sagittal orientation, ranging from lateral + 0.24 mm to lateral + 0.84mm. Quantitative analysis both on the coronal and sagittal sections was done, and MVP expression was scored as mild (+), moderate (++), or strong (+++) across all the positive brain regions. The complete analysis is available in Additional file 2: Tables S3-S4. Pir, piriform cortex; MS, medial septum; DBB, diagonal band of Broca; TS, triangular septal nucleus; ZI, zona incerta; RM, retromammillary nucleus; py/mol, pyramidal/molecular layer of the hippocampus (HP); DG, dentate gyrus of the hippocampus; 4V, fourth ventricle; CP, choroid plexus; MVe, medial vestibular nucleus; Sol, solitary nucleus; AP, area postrema; 10N, vagus nucleus; Raphe, raphe magnus nucleus. P, Q IHF of MVP, NeuN neurons, and Hoechst on the brain sections from male (P) and female (Q) WT mice. The dashed square delineates the region presented at a higher magnification on individual channels for MVP, NeuN, and Hoechst, and merge of the three. In all structures, the MVP signal came from the cytoplasm
The neuroanatomical spatial distribution of the major vault protein was then quantified by immunofluorescence using an anti-vault antibody, throughout four entire adult brains on consecutive histological sections by systematically counting the number of MVP/vault-positive cells and scoring as mild (+), moderate (++), or strong (+++). A total of 36 brain regions showed a positive signal (Additional file 2: Table S3) with equal scores between the sexes summarized in Additional file 2: Table S4. More specifically, MVP/vault signal was refined to the granular layer and the arbor vitae of the cerebellum (Fig. 2C), the oriens layer of the CA3 region of the ventral hippocampus (Fig. 2D), and the deep layers of the cingulate gyrus (Fig. 2E). Interestingly, we noticed MVP/vault signal pertained to specific nuclei of the limbic system (Fig. 2F–O), for example, the triangular septal nucleus (Fig. 2I), the zona incerta (Fig. 2J), the solitary nucleus (Fig. 2M), the vagal nucleus (Fig. 2O), the paraventricular hypothalamic nuclei, and the dorsal medial thalamic nucleus (Additional file 2: Table S4). At the subcellular level, MVP/vault signal was limited to the cytoplasm of neurons both in males and females (Fig. 2P, Q and Additional file 1: Fig. S5D-F). It was noteworthy that the patterns observed were consistent across sectioning planes (coronal and sagittal) and replicates.
To sum up, our data show no apparent sex differences in the pattern of MVP/vault expression that could explain the sex differences in NAPs. We cannot however exclude possible sex differences in the levels of MVP expression in other cell types of the brain that will require further investigations.
The major vault protein is implicated in the regulation of brain size and neuronal morphology in males only
The MVP mouse model used in this study was validated as loss-of-function (LoF) of Mvp, assessed at the transcript and protein levels (Fig. 3A, B, Additional file 1: Fig. S6A-D and Additional file 9). We also verified that the expression of neighboring genes was unaffected (Additional file 1: Fig. S6E). Among 617 successfully genotyped animals, we observed normal Mendelian ratio inheritance, indicating that the loss of Mvp has no effect on survival (Fig. 3C and Additional file 1: Fig. S6F).
Validation and phenotyping of Mvp-deficient mice. A Schematic representation of Mvp knockout model construction. The knockout mutants were obtained by a promoter gene trap design, resulting in the insertion of a β-galactosidase/neomycin phosphotransferase (β-GEO) cassette within the first intron of the Mvp gene. The polyA-tail (pA) serves to end the transcription and thus discard the transcription of Mvp gene. B Validation of Mvp knockout at the protein level using Western blot (WB) from liver extracts of Mvp+/+ (n = 2), Mvp+/− (n = 2), and Mvp−/− (n = 3) mice. Actin was used as a loading control for normalization. For quantification, see Additional file 1: Fig. S6D. C Viability screen of Mvp knockout mice. Left: schematic representation of heterozygous-by-heterozygous breeding scheme, with expected Mendelian ratio of segregation of the offspring. Right: graph represents the ratio of each genotype of successfully genotyped mice obtained from heterozygous-by-heterozygous breeding (n = 617). D, E Graph showing the total brain area (1_TBA) measured at four time points, embryonic age 18.5 (E18.5), postnatal day 10 (P10), P45, and P120, in male (D) and female (E) Mvp−/− mice (n = 5–7 in each group). F Half images of the coronal brain sections stained with Nissl-luxol at the four time points from male Mvp+/+ and Mvp−/− mice. Plots are shown as mean ± SEM. Two-tailed Student t test equal variance (D, E). *p < 0.05, **p < 0.01
We studied the dynamics of brain size across four time points: embryonic day 18.5 (E18.5), postnatal day 10 (P10), P45, and P120, using our systematic and standardized procedures [37,38,39,40]. At E18.5 and P10, the total brain area measurement in male Mvp−/− mice was normal but smaller at P45 (− 7%, p = 0.0061) and P120 (− 9%, p = 0.046), suggesting that the microcephaly is acquired between P10 and P45 (Fig. 3D, F). Interestingly, a recent study found that the protein level of MVP became higher at P35 [33], which could explain why the microcephaly phenotype became visible after this point at P45 in our study. Female Mvp−/− brain size was normal across the time points assessed (Fig. 3E). In addition to the total brain area, we quantified 40 other parameters at E18.5 (Additional file 2: Table S5) and did not find any phenotypes at this stage in both sexes.
At P120, het male Mvp+/− mice showed the same reduction in total brain size (− 9%, p = 0.016) when compared to hom male Mvp−/− mice (− 9%, p = 0.046) (Fig. 4A, B). Consistently, both Mvp+/− and Mvp−/− male mice showed small brain nuclei associated to the limbic system such as the cingulate gyrus (Mvp+/−: − 13%, p = 0.0016; Mvp−/−: − 8%, p = 0.007), the somatosensory cortex (Mvp+/−: − 12%, p = 0.0025; Mvp−/−: − 12%, p = 0.0016), and the hippocampus (Mvp+/−: − 20%, p = 0.023; Mvp−/−: − 24%, p = 0.038). In female Mvp−/− mice, no change was detected at any given time (Additional file 1: Fig. S7A). The number and effect size of NAPs being similar between male Mvp−/− and Mvp+/− mice (Fig. 4A, B: 17 NAPs versus 14, respectively), we performed subsequent studies specifically comparing Mvp−/− to Mvp+/+. NAPs in male Mvp−/− mice were also confirmed on sagittal planes using a previously described procedure [37] with 15 parameters significantly smaller in size at P120 including the cingulate gyrus (− 25%, p = 0.000065), the hippocampus (− 24%, p = 0.0069), the corpus callosum (− 26%, p = 0.000013), and the thalamus (− 24%, p = 0.00033) (Fig. 4E).
Involvement of the major vault protein in male-specific neuroanatomical phenotypes. A, B Left: schematic representation of the 24 brain regions assessed at bregma + 0.98 mm and − 1.34 mm in male Mvp+/− (A) and male Mvp−/− (B) mice at P120. Colored regions indicate the presence of at least one significant parameter within the brain region at the 0.05 level. White indicates a p-value > 0.05, and gray shows not enough data to calculate a p-value. Right: histograms of percentage change relative to Mvp+/+ for each of the 42 parameters (see Additional file 2: Table S2). C, D Left: coronal sections of Cg (C) and somatosensory cortex S2 (D) from Mvp+/+ and Mvp−/− at P120. Right: cell count and average cell area measures within Cg and S2. E Left: schematic representation of the 22 brain regions quantified at lateral + 0.60 mm on the parasagittal section from Mvp+/+ and Mvp−/−. Colored regions indicate the presence of at least one significant parameter within the brain region at the 0.05 level. White indicates a p-value > 0.05, and gray shows not enough data to calculate a p-value. Right: histograms of percentage change relative to Mvp+/+ (set as 0) for each of the measured parameters (listed in Additional file 2: Table S6). F Left: schematic representation of a hippocampal neuron in culture and the different measures taken. Middle: neurons stained with MAP2 (red) and SMI312 (green). Right: graphs showing measures of soma area, growth cone area, and axonal length from hippocampal neurons derived from multiple male Mvp+/+ and Mvp−/− embryos. G Graphs showing measures of soma area, growth cone area, and axonal length from hippocampal neurons derived from multiple female Mvp+/+ and Mvp−/− embryos (n > 3 embryos per group, Additional file 5: Table S9). Plots are shown as mean ± SEM (except for A, B, and E). Two-tailed Student t test equal variance (C, D, F, and G). *p < 0.05, **p < 0.01, ***p < 0.001
To investigate whether neuronal morphology may contribute to brain size phenotype in male Mvp−/− mice, we took advantage of the high resolution of our approach and developed a suite of automated tools to count the number of cells and calculate the average cell size within each brain region (Additional file 9). No significant change was detected in the number of cells in male Mvp−/− mice; however, cells were significantly smaller in size in all affected brain regions including the cingulate gyrus (− 15%, p = 0.0053) (Fig. 4C) and somatosensory cortex (− 15%, p = 0.001) (Fig. 4D). Cell size was normal in unaffected brain region such as the retrosplenial cortex (Additional file 1: Fig. S8F). The same set of studies were conducted in female Mvp−/− mice, but no defects in cell count and size were found (Additional file 1: Fig. S7B-D). To further characterize the cellular phenotype, we conducted hippocampal neuronal cultures, independently in males and females (Fig. 4F, Additional file 9 and Additional file 5: Table S9). Consistently, neurons derived from male Mvp−/− mice showed a reduction of the soma size by 6% (p = 0.0003). The growth cones were smaller by 23% (p = 0.0081), and no differences were detected for the axonal length (Fig. 4F). Neurons derived from female Mvp−/− showed no cell morphological defects (Fig. 4G).
Finally, to test if the anatomical changes relate to neural connectivity defects in male Mvp−/− mice, we measured miniature excitatory postsynaptic currents (mEPSCs) in pyramidal neurons of the anterior cingulate gyrus (Additional file 9). The mEPSC amplitude was smaller (p = 0.023) (Additional file 1: Fig. S9A), and the density of postsynaptic dendritic spines was reduced (− 5%, p = 0.02) (Additional file 1: Fig. S9D), indicating functional and morphological changes in synaptic connections. Neuronal ultrastructure was examined, but no anomalies were seen in the cingulate gyrus (Additional file 1: Fig. S9E-F). To identify potential biological pathways that might explain NAPs, we carried out transcriptomic analyses in the cingulate gyrus, and although the complete loss of Mvp was confirmed, no differentially expressed genes were found (Additional file 1: Fig. S9G), suggesting that Mvp has no major role on bulk transcriptional regulation that could explain the phenotypes at the time point sampled.
Altogether, these findings show that Mvp is not essential for survival and is not implicated in neuronal proliferation processes but instead in the regulation of neuronal size, morphology, and function, as well as in the maintenance of brain homeostasis after birth, specifically in males.
The double hemideletion of Mvp and Mapk3 alters behavioral performances in male mice
To determine if NAPs in male Mvp−/− mice lead to specific behavioral phenotypes, we assessed a broad range of sixteen behavioral tests evaluating eleven core behaviors (anxiety, depression, anhedonia, memory, locomotion, coordination, motricity, sociability, schizophrenia, autism, and epilepsy). The raw data is available in Additional file 6: Tables S10, Additional file 7: Tables S11, and Additional file 8: Tables S12, and the behavioral pipeline is described in detail in Additional file 9.
Intriguingly, we found no major behavioral anomalies in male Mvp−/− mice (Fig. 5 and Additional file 1: Fig. S10). We thought to test whether introducing an additional stressor might underscore behavioral implications. To explore this idea, we generated and characterized the behavior of double male KOs of Mvp and Mapk3 (Additional file 1: Fig. S12 and Additional file 9) since studies have consistently shown MVP-mediated regulation of ERK signaling pathway [31,32,33], which is known to be one of the main pathways disrupted at the 16p11.2 locus [23].
Behavioral analysis of Mvp and Mapk3 genes. Selection of behavioral paradigms among 17 tested in 12 Mvp+/+;Mapk3+/+ and 12 Mvp+/−;Mapk3+/− mice, from 11 to 25 weeks of age. Elevated plus maze (A). Open field (B). Fear conditioning (C). Tail suspension (D). PTZ test (E). Forced swim (F). G Summary table of core functions assessed in various cohorts. Green indicates significance, gray no difference, and white indicates not done (not d.). Arrows indicate the directionality of the effect, and the equal signs indicate no changes. A comprehensive description of behavioral tests, raw datasets, and additional findings are provided in Additional file 9, Additional file 6: Tables S10, Additional file 7: Tables S11, Additional file 8: Tables S12, and Additional file 1: Figs. S10-S12, respectively. The mean and standard error of the mean are shown in the graphs. Two-tailed Student t test equal variance (A, D–F) and two-way ANOVA with Sidak post hoc (B, C). *p < 0.05, **p < 0.01, ***p < 0.001
Double Mvp+/−;Mapk3+/− male mice showed normal body weight (Additional file 1: Fig. S12A) as well as normal motricity, coordination, reward-seeking behaviors, working and short-term memory, associative learning, circadian activity, and no signs of autism and schizophrenic-like behaviors (Additional file 1: Fig. S12B-K). However, double Mvp+/−;Mapk3+/− male mice spent more time in the open arm of the elevated plus maze (Fig. 5A). While the open-field test indicated no increase in traveled distance suggesting normal basic activity, it showed increased time spent in the center of the apparatus (Fig. 5B) and decreased latency to enter the center (Additional file 1: Fig. S12L), suggesting resistance to anxiogenic behavior. Fear conditioning in contextual and cued testing showed decreased freezing time (Fig. 5C). The tail suspension test showed a trend of increased latency to immobility (Fig. 5D), while the forced swim test did not show any apparent differences (Fig. 4F). The duration of clonic seizures was dramatically reduced in the pentylenetetrazol (PTZ)-induced seizure paradigm (Fig. 5E), supporting the idea of a resistance to epilepsy also. Additionally, we studied the same series of eleven core functions and assessed the phenotypes of male Mapk3+/− mice. Like male Mvp−/− mice, we found no behavioral anomalies in male Mapk3+/− mice (Fig. 5 and Additional file 1: Fig. S11), suggesting that Mvp or Mapk3 alone is not sufficient to potentiate the behaviors seen in the double Mvp+/−;Mapk3+/− male mice. A summary of behavioral findings is provided in Fig. 5G. Taken together, these results confirm the in vivo interaction between MVP and ERK1. Finally, preliminary findings, where we examined ERK activity in the cortex of our various mouse models in males and females, indicate that ERK activity is not associated with NAPs or behavioral phenotypes presented in this study (Additional file 1: Fig. S12M-P and Additional file 10).
Add Comment