Sage Journals: Discover world-class research

Abstract

Background:

A finely tuned balance between excitation and inhibition is essential for proper brain function. Disruptions in the GABAergic system, which alter this equilibrium, are a common feature in various types of neurological disorders. Understanding GABAergic neuron maturation processes is thus currently a major challenge in basic neuroscience. Thyroid hormones (THs) are required for the proper maturation of parvalbumin (PV)-expressing GABAergic interneurons in the mouse neocortex. However, the timeline of this TH action has yet to be elucidated. The aim of the present study was to define better the time window during which THs promote the postnatal maturation of PV neurons in the mouse neocortex.

Methods:

We used genetically engineered mouse models expressing dominant-negative mutations of the TH nuclear receptor α1 (TRα1). The consequences of blocking TH signaling at different times in development were assessed in PV neurons of the somatosensory cortex, in terms of histology and gene expression.

Results:

Histological observations in mice revealed that the action of THs during the first three postnatal weeks was necessary to initiate the expression of PV and the elaboration of a specialized extracellular matrix called the perineuronal net (PNN). By contrast, after the third postnatal week, TH action on PV neuron maturation appeared to be somewhat dispensable. Transcriptome analysis of neocortical GABAergic neurons two weeks after birth identified a small set of putative target genes for TRα1. Several of these genes are involved in the postnatal remodeling of the repertoire of ion channels within PV neurons and in the elaboration of PNNs.

Conclusions:

These data suggest that THs act as a timer to define the temporal boundaries of the critical period of heightened cortical plasticity, which plays a fundamental role in the development of neuronal circuits. Unveiling the molecular underpinnings of TH action in PV neurons may help understand better neurological disorders associated with alterations of TH signaling, such as hypothyroidism, resistance to THs, or Allan–Herndon–Dudley syndrome, but also more widely, neurological disorders associated with an imbalance in the excitation/inhibition ratio in the brain, including attention-deficit/hyperactivity disorder, autism, and epilepsy.

Keywords

thyroid hormones TRα1 PV neurons GABAergic neurons neocortex

Introduction

Parvalbumin-expressing neurons (PV neurons) are GABAergic inhibitory interneurons that play a central role in the organization of neuronal circuits in the neocortex.¹ They are fast-spiking, synchronize network oscillations, and play an important role in information processing. In mice, PV neurons begin to express parvalbumin (PV) during the first postnatal weeks.^2,3 During this time, the electrical activity of these neurons progressively increases,⁴ and synaptogenesis is very active. This opens the so-called critical period of heightened synaptic plasticity⁵ during which sensory input and experience shape neuronal networks. The critical period ends when circuit rewiring becomes actively dampened.⁶ Notably, the elaboration of the perineuronal net (PNN) around PV neurons, a densely packed extracellular matrix enriched in proteoglycans, stabilizes synapses and restricts plasticity.⁷ Defects in PV neuron maturation during the critical period cause a persistent imbalance between neuronal excitation and inhibition, which is at the origin of several pathological conditions, including attention-deficit/hyperactivity disorder, autism, and epilepsy.⁸

Thyroid hormones (THs; including triiodothyronine [T3] or 3,3′,5-triiodo-l-thyronine, and the less active precursor T4 or thyroxine) are critically required for the proper maturation of PV neurons.^9–11 Mutations of TH nuclear receptor α1 (TRα1), the most ubiquitous nuclear receptor of T3, also have dramatic and irreversible consequences on neurodevelopment, similar to those of TH deficiency in both humans^12,13 and mice.^14–16 Furthermore, the restricted expression of TRα1^L400R, a dominant-negative TRα1 receptor, in the GABAergic lineage is sufficient to dramatically inhibit the expression of PV and the elaboration of PNNs.¹⁷

To gain more insight into the timeline of TH action on PV neuron maturation, we generated novel mouse models with restricted or inducible expression of mutant TRα1 receptors. We evidenced a lifelong persistence of defects in PV neuron maturation when TH signaling was blocked in GABAergic neurons from early stages of development. Transcriptome analysis identified a set of genes whose regulation is expected to be involved in TH-mediated maturation of PV neurons and in the closure of the critical period.

Materials and Methods

Mouse models and treatments

Experiments involving the use of live animals for the current project were approved by local ethics committees (C2EA015 and C2EA017) and subsequently authorized by the French Ministry of Research (Projects #33279-2021082516194165 and #25496-2020032811174506).

The Thra^AMI allele allows the expression of the dominant-negative TRα1^L400R after Cre/loxP-mediated deletion of a cassette with polyadenylation signals.¹⁸ Thra^Slox is identical, except for the presence of a C-terminus frameshift mutation, which leads to the expression of the dominant-negative TRα1^E395fs401X. The ROSA-tdTomato^lox reporter transgene (also known as Ai9, MGI Cat# 4436851, RRID: MGI: 4436851) drives the Cre-dependent expression of a red fluorescent protein, mainly localized in cytoplasm.¹⁹ The ROSA-GSGFPL10^lox transgene encodes the GS-EGFP-L10a protein, which is a green fluorescent protein fused to the N-terminus part of the large subunit ribosomal protein L10a, mainly localized in cell nuclei.²⁰ All transgenes were in the C57Bl6/J genetic background. Gad2-Cre and Gad2-CreER^T2 are knock-in alleles in which the Cre recombinase reading frame is inserted into the Gad2 gene in order to restrict the expression of the recombinase in GABAergic neurons.²¹ In Gad2-Cre mice, Cre induces DNA recombination from around E12.5, whereas in Gad2-CreER^T2 mice, CreER^T2 is translocated to the cell nucleus only after mice are given tamoxifen (T5648, Sigma, 10 mg/mL in sunflower oil), which allows for control of the time at which a mutated TRα1 starts to be expressed. PND7 mice received a single 50 µL subcutaneous tamoxifen injection, whereas 3-week-old mice received five daily intraperitoneal injections (daily dose: 7.5 µL/g).

Gestating mice were fed for 2 weeks with iodine-deficient food supplemented with 0.15% propyl-thiouracyl (PTU Envigo ref TD.95125) to cause deep hypothyroidism. TH levels were restored by a single intraperitoneal injection in half of the pups (10 µg T4 + 1 µg T3 dissolved in 50 µL of phosphate-buffered saline [PBS], all chemicals from Sigma Aldrich, France). The other PTU-treated animals were injected with 50 µL of PBS.

RNA sequencing analysis from sorted nuclei

Nuclei were isolated from whole individual cortex frozen in liquid nitrogen, and nuclear RNA was extracted as previously described.²⁰ Nuclear RNA was extracted from sorted nuclei (RNeasy Micro Kit, Qiagen ref 74004) and quantified using Tapestation4150 (Agilent). One nanogram of RNA was reverse-transcribed using the SMART-Seq v4 Low Input RNA kit (Takara). Complementary DNA (cDNA) was quantified and qualified using Qubit (Invitrogen) and Tapestation4150 (Agilent). Libraries were then prepared from 1 ng cDNA using the Nextera XT DNA Library Kit (Illumina). Libraries were sequenced (>2.10⁷ reads/library) on a NextSeq500 DNA sequencer (Illumina). Reads were aligned on the mouse genome (mm10 GRCm38 release) with Bowtie2 (Galaxy Version 2.4.2 + galaxy 0).²² Count table (Supplementary Table S1) was prepared using htseq-count (Galaxy Version 0.9.1 + galaxy 1, mode union, feature type: gene).²³ Differential gene expression analysis was performed with DEseq2 (Galaxy Version 2.1.8.3)²⁴ (Supplementary Table S1) using the following thresholds: false discovery rate <0.05; p-adjusted value <0.05; maximum adjusted basecount >30; and fold change >2.0. Hierarchical clustering was performed using Euclidean distance and complete distance with ClustVis.²⁵ Gene ontology analysis was performed with Gorilla (https://cbl-gorilla.cs.technion.ac.il/). Raw RNA sequencing (RNA-seq) data are available at NCBI Gene Expression Omnibus (GSE290072 and GSE290207).

RT-qPCR analysis

RNAs were extracted using TRIzol (Invitrogen, Carlsbad, CA). Total RNA was converted to cDNA using M-MLV reverse transcriptase (Promega, Madison, WI). RT-qPCR was performed using SYBRGreen mix (BioRad iQ supermix). The results were analyzed according to the ΔΔCT method.²⁶ The housekeeping gene Hprt was used as the reference.

Brain slice histology

Brain slice collection was performed as described.¹⁷ Coronal sections (50 μm) were cut with a vibrating microtome (Integraslice 7550 SPDS; Campden Instruments, Loughborough, UK), in PBS at room temperature. Brain sections were stored at −20°C in cryoprotectant (30% ethylene glycol and 20% glycerol in 10 mM low-salt PBS). Immunohistochemistry was performed on free-floating brain sections. For PV labeling, a mouse anti-PV primary antibody (PARV19, Sigma P3088, 1:2000) and a secondary antibody made in donkey (anti-mouse DyLight 633, ThermoFisher Scientific, 1:1000) were used. For the concurrent PNN labeling, a biotinylated Wisteria floribunda lectin (WFL/WFA, Vector Laboratories B-1355, 20 µg/mL) and streptavidin coupled to DyLight 488 (Vector Laboratories, SA-5488, 1:1000) were used. Nonspecific binding sites were blocked by incubating sections for 1 hour in PBS with 1% bovine serum albumin (BSA) and 0.2% Triton X-100. Brain sections were incubated overnight at 4°C with PV antibody and biotinylated WFA, which were diluted in PBS with 1% BSA, 0.2% Triton X-100, and 1% dimethyl sulfoxide. Sections were washed in PBS and further incubated for 15 minutes at room temperature with 4’,6-diamidino-2-phenylindole (1:5000, Sigma). Incubation with the secondary antibody lasted for three hours at room temperature. Sections were mounted in Fluoroshield™ (Sigma), coverslipped, and imaged using an inverted confocal microscope (Zeiss LSM 780 and Zeiss LSM 800).

Image analysis

Image analysis was performed using ImageJ. Final images resulted from the z projection of three optical sections (maximum intensity tool) that were 2 µm apart (the resulting images thus reflected the fluorescence collected over a thickness of 4 µm). PV+ WFA+ Tomato+ cells were detected from final images using the Cellpose wrapper for Fiji.²⁷ Numbers of mono-, double-, or triple-labeled cells per slice were measured automatically from 3 to 4 slices per animal. The code of the original macro may be found at https://github.com/jbrocardplatim/MacroPV.

The densities of Tomato+, PV+, WFA+ cells, and their colocalizations were counted in the whole area of one image taken in the somatosensory cortex. Three to four images in each hemisphere were analyzed for each mouse. Quantitative data were analyzed using a two-tailed unpaired T-test. Data are expressed as mean ± standard error of the mean. Data analyses were performed with Excel 2016 and GraphPad Prism software (v.8). The level of significance was set at p < 0.05.

Results

TRα1^E395fs401X expression in GABAergic neurons from an early developmental stage permanently inhibits PV neuron maturation

As the expression of TRα1^L400R in the GABAergic lineage, starting at embryonic stage, is lethal after birth,¹⁷ we used TRα1^E395fs401X hypomorphic mutation, whose expression can be triggered by Cre/loxP recombination from the Thra^Slox allele²⁸ to address the long-term consequences of preventing T3 response in GABAergic neurons. One-third of Thra^Slox/+ ROSA-tdTomato^lox/+ Gad2-Cre mice survived to adulthood, while all the Thra⁺^/⁺ ROSA-tdTomato^lox/+ Gad2-Cre littermates survived and were used as controls.

Red fluorescence allowed one to trace the entire GABAergic lineage and did not indicate any cell loss (Supplementary Fig. S1). On the contrary, the density of PV+ cells was drastically reduced in mutant mice (Fig. 1a, b). The density of WFA+ cells (Fig. 1a, c) and the ratio of PV+ cells labeled with WFA were reduced (Fig. 1a, d). Some recovery of PV+ cell density was observed at the adult stage (Fig. 1e, f), while the PNN defect persisted (Fig. 1e, g, h). Therefore, early TRα1^E395fs401X expression permanently affects the maturation of PV neurons.

FIG. 1.

Persistent defect in PV neuron maturation in mice expressing TRα1^E395fs401X. Thra⁺^/⁺ ROSA-tdTomato^lox/+ Gad2-Cre littermates (controls) and Thra^Slox/+ ROSA-tdTomato^lox/+ Gad2-Cre mice (mutants) express TRα1^E395fs401X and a red fluorescent protein only in GABAergic neurons and in their progenitors from embryonic day 12.5. The neocortex of mutant mice did not display any major anatomical alteration. Histological analysis of the somatosensory neocortex at PND21 (a–d, controls: n = 6, mutants: n = 6) and at adult stage (e–h, controls: n = 6, mutants: n = 5). (a, e) Immunohistochemistry of parvalbumin and WFA-labeled PNNs in the somatosensory cortex (Bregma 1.0–0.0). WFA binds N-acetyl-galactosamine-β1 residues of PNN glycoproteins. (b, f) The density of PV+ cells was lower in mutants at both ages. As PV neurons normally represent 40% of the GABAergic neuronal population in the neocortex,²⁹ and as the overall density of GABAergic neurons is unchanged (as evidenced by the density of tdTomato+ cells, see Supplementary Fig. S1), the decreased density of PV+ cells likely reflects a defect in their maturation. A similar phenotype has previously been reported in Thra^AMI/+ ROSA-tdTomato^loxI/+ Gad2-Cre mice.¹⁷ (c, g) The density of WFA+ cells was lower in mutants at both ages. (d, h) The proportion of PV+ cells surrounded by PNNs was also lower in mutants at both ages. Data are presented as mean ± SEM and analyzed using an unpaired two-tailed T-test (*p < 0.05, **p < 0.01, ***p < 0.001). PNN, perineuronal net; PV, parvalbumin; SEM, standard error of the mean; TRα1, TH nuclear receptor α1; WFA, Wisteria Floribunda Agglutinin.

Postnatal expression of TRα1^L400R in GABAergic neurons is sufficient to inhibit PV neuron maturation

We replaced the Gad2Cre driver with the Gad2-CreERT2 inducible version. This allows for the induction of TRα1^L400R expression by tamoxifen treatment at a chosen time and addresses a possible temporal restriction on T3 action in GABAergic neurons. Treatment with tamoxifen at PND7 was sufficient to significantly reduce PV+ cell density at PND14 (Fig. 2a, b). The density of WFA+ cells (Fig. 2a, c) and the fraction of PV+ cells labeled with WFA (Fig. 2a, d) remained unchanged, perhaps because the elaboration of PNNs was not completed at this early stage. By contrast, a treatment with tamoxifen between PND21 and PND25 did not have visible histological consequences in the adult stage (Fig. 2e–h). Therefore, THs exert their influence during a postnatal time window, which includes the second postnatal week and precedes weaning (∼PND25).

FIG. 2.

A time window for TH action on PV neuron maturation. Histological analysis of Thra^+/+ ROSA-tdTomato^lox/+ Gad2-CreER^T2 or Gad2-CreER^T2 negative littermates (controls) and Thra^AMI^/⁺ ROSA-tdTomato^lox/+ Gad2-CreER^T2 mice (mutants) that express TRα1^L400R in the GABAergic lineage from PND7 to PND14 (a–d, controls: n = 8, mutants: n = 6) and from PND21 to adult (e–h, controls: n = 11, mutants: n = 7). The neocortex of mutant mice did not display any major anatomical alteration. (a, e) Immunohistochemistry of parvalbumin and WFA-labeled PNNs in the somatosensory cortex (Bregma 1.0–0.0). (b, f) The density of PV+ cells was lower in mice expressing TRα1^L400R from PND7 but not from PND21. (c, g) The density of WFA+ cells was not significantly altered in either case. Note that WFA staining is weak at PND14 as PNN elaboration is not terminated. (d, h) The proportion of PV neurons surrounded by PNNs was not significantly altered in either case. Data are presented as mean ± SEM and analyzed using unpaired two-tailed T-test (*p < 0.05, **p < 0.01, ***p < 0.001). PNN, perineuronal net; PV, parvalbumin; SEM, standard error of the mean; TH, thyroid hormone; TRα1, TH nuclear receptor α1; WFA, Wisteria Floribunda Agglutinin.

Consequences of TRα1^L400R expression on gene expression in neocortical GABAergic neurons

We used the ROSA-GSGFPL10^lox reporter transgene, which labels cell nuclei, to extract nuclear RNA from neocortical GABAergic neurons at PND15 and perform a bulk RNA-seq analysis. A total of 1306 genes were sensitive to either TRα1^L400R expression, hypothyroidism, or TH stimulation (Deseq2 adjusted p-value >0.05; fold change >2.0; and basemean >30; Fig. 3a, b). Among them, a large set of genes was only deregulated in hypothyroid mice. In the end, only 61 genes behaved as bona fide TRα1 target genes, with a robust negative response to both TRα1^L400R expression and hypothyroidism, and clear induction after short TH treatment of hypothyroid mice (Fig. 3a). Clustering analysis revealed a general trend toward recovery of hypothyroid mice after TH treatment (Fig. 3c, d). Five genes appeared to be negatively regulated by T3-bound TRα1 (Fig. 3b). qPCR analyses carried out on the cortex of mice with different circulating TH levels confirmed part of the results of the RNA-seq analysis (Supplementary Fig. S2). Notably, many of the genes identified as regulated by TH/TRα1 signaling are expressed in PV neurons (Supplementary Table S2).

FIG. 3.

RNA sequencing analysis of gene expression in sorted GABAergic neuron nuclei. Taking advantage of the selective expression of a fluorescent protein present in the nuclei, RNA was extracted from nuclei of GABAergic neurons sorted from the neocortex of four groups of PND14 mice. Thra^+/+ ROSA-GSGFPL10^lox Gad2-Cre mice (Ctrl group) are control transgenic mice that express only fluorescent protein. Thra^AMI/+ ROSA-GSGFPL10^lox Gad2-Cre mice (Mut group) express TRα1^L400R in GABAergic neurons. Thra^+/+ ROSA-GSGFPL10^lox Gad2-Cre mice were made hypothyroid by a PTU treatment started during gestation (Hypoth. group). Some of these hypothyroid Thra^+/+ ROSA-GSGFPL10^lox Gad2-Cre mice received a single TH injection at PND13 (hypoth. +TH group). The rapid response of genes after TH treatment is a hallmark of genes whose transcription is directly regulated by the liganded TRα1 nuclear receptor. (a) Venn diagram of genes that were downregulated by TRα1^L400R, downregulated by hypothyroidism, or induced after TH treatment (thresholds: fold change >2 or <0.5, maximum adjusted base counts >30, adjusted p-value >0.05). (b) Venn diagram of genes that were either upregulated by TRα1^L400R, upregulated by hypothyroidism, or repressed after TH treatment (thresholds: fold change >2 or <0.5, maximal adjusted base counts >30, adjusted p-value >0.05). (c) Heatmap showing genes that were downregulated by both TRα1^L400R and hypothyroidism. Brown is for high expression, and blue is for low expression. Note that although TH treatment tends to restore the initial expression level for most genes, a full recovery is only observed for the upper gene cluster. (d) Heatmap showing genes that were upregulated by both TRα1^L400R and hypothyroidism. Recovery after TH treatment is less visible than in (c), probably due to the stability and slow turnover of nuclear RNA in neurons. (e) Heatmap of markers of PV neuron maturation. The marks on the right side are for genes whose expression normally increases (red) or decreases (green) between PND7 and PND21. Note the reduced expression of genes which are normally upregulated over time, in PND14 hypothyroid and TRα1^L400R-expressing mice. These include a number of genes encoding potassium channels (Kcna1, Kcna2, Kcnab2, Kcnab3, Kcnb1, Kcnc1, Kcnd2, Kcng4, Kcnj9, and Kcnq2), sodium channels (Scn1, Scn4b, and Scn9a) and proteins playing a key role in the elaboration of the PNN (Bcan, Hapln4, and Kank2). The opposite is generally observed for genes that are normally downregulated over time, several of which encode components of calcium voltage-gated channels (Cacna1b, Cacna1g, Cacnb1, Cacng4, and Cacng5). These alterations of gene expression reveal a blockade of the PV neuron maturation program, caused either by TRα1^L400R expression or by hypothyroidism. A 24-hour TH treatment is not sufficient to restore the expression of these marker genes. PNN, perineuronal net; PV, parvalbumin; TH, thyroid hormone; TRα1, TH nuclear receptor α1.

Gene ontology analysis indicated that many of the genes deregulated in hypothyroid and mutant mice encode components of the extracellular matrix (GO:0031012; enrichment 4.4, p-value 3.10⁻¹³). This includes genes encoding collagens or proteins binding to the glycosaminoglycans of PNNs (GO:0005539; enrichment 4.2, p-value 5.10⁻⁶). The expression pattern of two of these genes (Sema7a, encoding semaphorin 7A and Hapnl4, encoding hyaluronan and proteoglycan link protein 4) identifies them as genuine TRα1 target genes. This observation is consistent with the hypothesis that TH/TRα1 signaling stimulates the elaboration of PNNs around PV neurons.

The large group of genes that are deregulated in hypothyroid and mutant mice but are not responsive to TH stimulation is also enriched in voltage-gated potassium channels (GO:0005249; enrichment 6.1, p-value 10⁻⁶). We extracted the genes that mainly encode ion channels and whose expression changes during postnatal maturation of GABAergic neurons.⁴ This analysis (Fig. 3e) revealed that the immature status of PV neurons persists in hypothyroid and mutant mice at PND14.

Discussion

In the present study, we extend our previous observations showing that a cell-autonomous action of TH/TRα1 in GABAergic neurons is required for the maturation of neocortical PV neurons in mice. We identified a time window that spans from birth to the end of the third postnatal week, during which the action of THs is critical for the proper maturation of PV neurons. Analysis of gene expression also reveals that at this stage, TH/TRα1 normally activates the transcription of a small group of genes necessary for the elaboration of the PNN around PV neurons.

PV neurons show near-normal maturation when the blockade of TH/TRα1 signaling occurs at a late stage, or if it occurs only after the onset of the expression of Pvalb, the gene that encodes PV.²⁰ This confirms a previous whole-cortex RNA-seq analysis, which showed that gene expression changes rapidly in GABAergic neurons during the second postnatal week under the control of TH/TRα1 signaling.¹⁷ Our novel transcriptome analysis indicates that blocking TH/TRα1 signaling compromises the expression of a number of genes involved in a well-documented transition, during which the repertoire of ion channels is extensively modified.⁴ It is remarkable that the second postnatal week, during which THs seem to play a crucial role in PV neuron maturation, coincides with a peak in circulating levels of THs.³⁰ This strongly suggests that THs represent a temporal cue to synchronize the postnatal maturation of PV neurons.

We provide evidence that the postnatal maturation of PV neurons is promoted in a cell-autonomous manner. However, we did not notice any obvious change in the radial distribution of GABAergic neurons. Therefore, the reported defect in prenatal migration of GABAergic neuron progenitors might not be a cell-autonomous consequence of early TH deficiency.³¹ A complete recovery of PV neuron density in adults has been reported for mice expressing TRα1^R384C,14 which we did not observe for TRα1^E395fs401. This discrepancy is probably explained by the residual capacity of TRα1^R384C to bind T3 and activate transcription.³² Furthermore, the late onset of TRβ1 expression in PV neurons³³ could be sufficient to counteract the negative influence of TRα1^R384C.

The mode of action of THs in PV neurons during the second postnatal week is clarified by RNA-seq data. We found that PV neurons were much more sensitive to variations in TH signaling at PND14 than in adulthood, where the same method only showed minor changes in gene expression.²⁰ Another finding is that genes that are consistently deregulated by blocking TH/TRα1 signaling encode extracellular matrix proteins, specific components of or proteins interacting with PNN glycosaminoglycans in PV neurons. Rapid recovery of gene expression after TH treatment in hypothyroid mice provides a clear indication of direct transcriptional control by TRα1 for a small subset of genes. Deregulation of these genes is likely to represent the initial event preventing PV neuron maturation. Among these were genes encoding transcription factors or cofactors (Dbp, Hr, Klf9), also regulated by TH in other cell types³⁴ and again genes encoding PNN components. Thus, it appears that liganded TRα1 notably activates the transcription of genes encoding collagens (Col24a1 and Col9a2), ephrin A4 (Epha4), laminin alpha 3 subunit (Lama3), semaphorin 7A (Sema7a), hyaluronan, and proteoglycan link protein 4 (Hapln4), an essential component of PNNs.^35–37

The novel data presented in this article provide strong support for the hypothesis that TH is a major temporal signal for the opening and closing of the critical period of heightened plasticity, as previously suggested.³⁸ The critical period occurs at slightly different times in different neocortical areas, as does the onset of TH signaling.³⁹ It is characterized by a series of cellular events that play a crucial role in the elaboration of functional neuronal circuits.⁵ The beginning of PV synthesis marks the opening of the critical period, while the elaboration of PNNs around PV neurons stabilizes synapses and decides their closure. Thus, our findings suggest that both the opening and the closure of the critical period are directly stimulated by liganded TRα1. Therefore, THs appear to act in the postnatal neocortex as it does during amphibian metamorphosis, as an external cue that represents an important temporal signal. Of note, oligodendrocyte maturation, which is also a TH-dependent process,^40,41 has also been shown to play a role in PV neuron maturation and in the closure of the critical period.⁴²

Patients with early TH deficiency caused by congenital hypothyroidism⁴³ or carrying TRα1 mutations¹³ are currently treated with TH replacement therapy. Treatment can have a positive influence on cognitive development only if it is initiated soon after birth. This is an indication that the action of TH on the maturation of the human neocortex is restricted to an early temporal window, as we have found in mice. An alteration in the timing of the critical period is also suspected to cause autism spectrum disorders.⁸ In that regard, it is noticeable that some of the rare patients in whom THRA mutations have been discovered were first diagnosed for autism spectrum disorders.^44,45 Within the framework of our hypothesis, the irreversible cognitive and neurological defects that are seen both in patients with TH-related diseases and in patients with autism spectrum disorders might be similarly regarded as consequences of a disruption of the critical period of heightened cortical plasticity. Delving deeper into the detailed molecular events occurring in the cortex during the critical period will be key to finding therapeutic approaches to these debilitating pathologies.

Study limitation

Our study identifies TH-regulated genes in PV neurons but does not address chromatin binding of TRα1 to the neighboring regulatory sequences. This is not currently feasible with the small number of PV neurons sorted from the mouse neocortex. Also, the histological methods that we have used only indirectly inform on the timing of the critical period, which can be defined in a number of different ways.^46,47

Conclusions

THs appear as major players of the critical period of heightened cortical plasticity, which is a decisive step in neurodevelopment. Unveiling the underpinnings of TH action in PV neurons should help understand neurological disorders deriving from alterations in TH signaling, as well as those associated with an imbalance in the excitation/inhibition ratio in the brain.

Authors’ Contributions

J.R.: Investigation, writing—original draft, and writing—review and editing. S.M., J.B., R.G., and D.A.: Investigation and writing—review and editing. J.W.: Supervision. F.F.: Conceptualization, funding acquisition, supervision, writing—original draft, and writing—review and editing. S.R.: Conceptualization, investigation, supervision, and writing—review and editing.

Footnotes

Acknowledgments

The authors acknowledge the contribution of SFR Biosciences (Université Claude Bernard Lyon 1, CNRS UAR3444, Inserm US8, ENS de Lyon): Nadine Aguilera and the Plateau de Biologie Expérimentale de la Souris (ANIRA-PBES) for mouse breeding; ANIRA-CYTOMETRIE and particularly Sébastien Dussurgey for nuclei sorting; and ANIRA-AGC for mouse genotyping. The authors thank Benjamin Gillet and Sandrine Hughes of the deep sequencing facility (PSI IGFL, Lyon). They also thank Elodie Martel for carrying out histology experiments as part of her internships and Isabelle Dusart (Institut de Biologie Paris Seine, France) for helpful discussions.

Author Disclosure Statement

The authors declare no conflicts of interest.

Funding Information

Funding for this research was provided by the China Scholarship Council (supporting Juan REN’s PhD project at ENS de Lyon), the European Union’s Horizon 2020 program (grant agreement No. 825753, ERGO), the Agence Nationale de la Recherche (Thyromut2 program, ANR-15-CE14-0011-01), and the Joint Research Institute for Science and Society (JORISS 2022, ).

Supplemental Material

References

Hijazi

, Smit

, van Kesteren

. Fast-spiking parvalbumin-positive interneurons in brain physiology and Alzheimer’s disease. Mol Psychiatry, 2023; 28(12):4954–4967; doi: 10.1038/s41380-023-02168-y

del Rio

, de Lecea

, Ferrer

, et al. The development of parvalbumin-immunoreactivity in the neocortex of the mouse. Brain Res Dev Brain Res, 1994; 81(2):247–259; doi: 10.1016/0165-3806(94)90311-5

Huang

, Kirkwood

, Pizzorusso

, et al. BDNF regulates the maturation of inhibition and the critical period of plasticity in mouse visual cortex. Cell, 1999; 98(6):739–755.

Okaty

, Miller

, Sugino

, et al. Transcriptional and electrophysiological maturation of neocortical fast-spiking GABAergic interneurons. J Neurosci, 2009; 29(21):7040–7052; doi: 10.1523/JNEUROSCI.0105-09.2009

Hensch

. Critical period plasticity in local cortical circuits. Nat Rev Neurosci, 2005; 6(11):877–888; doi: 10.1038/nrn1787

Takesian

, Hensch

. Balancing plasticity/stability across brain development. Prog Brain Res, 2013; 207:3–34; doi: 10.1016/B978-0-444-63327-9.00001-1

Wingert

, Sorg

. Impact of perineuronal nets on electrophysiology of parvalbumin interneurons, principal neurons, and brain oscillations: A review. Front Synaptic Neurosci, 2021; 13:673210; doi: 10.3389/fnsyn.2021.673210

Sohal

, Rubenstein

JLR

. Excitation-inhibition balance as a framework for investigating mechanisms in neuropsychiatric disorders. Mol Psychiatry, 2019; 24(9):1248–1257; doi: 10.1038/s41380-019-0426-0

Bernal

, Morte

, Diez

. Thyroid hormone regulators in human cerebral cortex development. J Endocrinol, 2022; 255(3):R27–R36; doi: 10.1530/Joe-22-0189

10.

Ren

, Flamant

. Thyroid hormone as a temporal switch in mouse development. Eur Thyroid J, 2023; 12(2):e220225; doi: 10.1530/ETJ-22-0225

11.

Richard

, Ren

, Flamant

. Thyroid hormone action during GABAergic neuron maturation: The quest for mechanisms. Front Endocrinol (Lausanne), 2023; 14:1256877; doi: 10.3389/fendo.2023.1256877

12.

Berbel

, Navarro

, Roman

. An evo-devo approach to thyroid hormones in cerebral and cerebellar cortical development: Etiological implications for autism. Front Endocrinol (Lausanne), 2014; 5:146; doi: 10.3389/fendo.2014.00146

13.

van Gucht

ALM

, Moran

, Meima

, et al. Resistance to thyroid hormone due to heterozygous mutations in thyroid hormone receptor alpha. Curr Top Dev Biol, 2017; 125:337–355; doi: 10.1016/bs.ctdb.2017.02.001

14.

Wallis

, Sjogren

, van Hogerlinden

, et al. Locomotor deficiencies and aberrant development of subtype-specific GABAergic interneurons caused by an unliganded thyroid hormone receptor alpha1. J Neurosci, 2008; 28(8):1904–1915; doi: 10.1523/JNEUROSCI.5163-07.2008

15.

Uchida

, Hasuoka

, Fuse

, et al. Thyroid hormone insufficiency alters the expression of psychiatric disorder-related molecules in the hypothyroid mouse brain during the early postnatal period. Sci Rep, 2021; 11(1):6723; doi: 10.1038/s41598-021-86237-8

16.

Mayerl

, Heuer

. Thyroid hormone transporter Mct8/Oatp1c1 deficiency compromises proper oligodendrocyte maturation in the mouse CNS. Neurobiol Dis, 2023; 184:106195; doi: 10.1016/j.nbd.2023.106195

17.

Richard

, Guyot

, Rey-Millet

, et al. A pivotal genetic program controlled by thyroid hormone during the maturation of GABAergic neurons. iScience, 2020; 23(3):100899; doi: 10.1016/j.isci.2020.100899

18.

Quignodon

, Vincent

, Winter

, et al. A point mutation in the activation function 2 domain of thyroid hormone receptor alpha1 expressed after CRE-mediated recombination partially recapitulates hypothyroidism. Mol Endocrinol, 2007; 21(10):2350–2360; doi: 10.1210/me.2007-0176

19.

Madisen

, Zwingman

, Sunkin

, et al. A robust and high-throughput Cre reporting and characterization system for the whole mouse brain. Nat Neurosci, 2010; 13(1):133–140; doi: 10.1038/nn.2467

20.

Ren

, Markossian

, Guyot

, et al. Thyroid hormones maintain parvalbumin neuron functions in mouse neocortex. bioRxiv, 2024; 2024(07.16.603713); doi: 10.1101/2024.07.16.603713

21.

Taniguchi

, He

, Wu

, et al. A resource of Cre driver lines for genetic targeting of GABAergic neurons in cerebral cortex. Neuron, 2011; 71(6):995–1013; doi: 10.1016/j.neuron.2011.07.026

22.

Langmead

, Trapnell

, Pop

, et al. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol, 2009; 10(3):R25; doi: 10.1186/gb-2009-10-3-r25

23.

Anders

, Pyl

, Huber

. HTSeq—A Python framework to work with high-throughput sequencing data. Bioinformatics, 2015; 31(2):166–169; doi: 10.1093/bioinformatics/btu638

24.

Love

, Huber

, Anders

. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol, 2014; 15(12):550; doi: 10.1186/s13059-014-0550-8

25.

Metsalu

, Vilo

. ClustVis: A web tool for visualizing clustering of multivariate data using Principal Component Analysis and heatmap. Nucleic Acids Res, 2015; 43(W1):W566–W70; doi: 10.1093/nar/gkv468

26.

Bookout

, Cummins

, Mangelsdorf

, et al. High-throughput real-time quantitative reverse transcription PCR. Curr Protoc Mol Biol, 2006;Chapter 15(Unit 15 8):Unit 15.8; doi: 10.1002/0471142727.mb1508s73

27.

Stringer

, Wang

, Michaelos

, et al. Cellpose: A generalist algorithm for cellular segmentation. Nat Methods, 2021; 18(1):100–106; doi: 10.1038/s41592-020-01018-x

28.

Markossian

, Guyot

, Richard

, et al. CRISPR/Cas9 editing of the mouse Thra gene produces models with variable resistance to thyroid hormone. Thyroid, 2018; 28(1):139–150; doi: 10.1089/thy.2017.0389

29.

Tasic

, Menon

, Nguyen

, et al. Adult mouse cortical cell taxonomy revealed by single cell transcriptomics. Nat Neurosci, 2016; 19(2):335–346; doi: 10.1038/nn.4216

30.

Friedrichsen

, Christ

, Heuer

, et al. Regulation of iodothyronine deiodinases in the Pax8-/- mouse model of congenital hypothyroidism. Endocrinology, 2003; 144(3):777–784.

31.

da Cunha Menezes

, de Abreu

, Davis

, et al. Effects of gestational hypothyroidism on mouse brain development: Gabaergic systems and oxidative stress. Dev Biol, 2024; 515:112–120; doi: 10.1016/j.ydbio.2024.07.010

32.

Tinnikov

, Nordstrom

, Thoren

, et al. Retardation of post-natal development caused by a negatively acting thyroid hormone receptor alpha1. Embo J, 2002; 21(19):5079–5087.

33.

Hochbaum

, Hulshof

, Urke

, et al. Thyroid hormone remodels cortex to coordinate body-wide metabolism and exploration. Cell, 2024; 187(20):5679–5697.e23; doi: 10.1016/j.cell.2024.07.041

34.

Zekri

, Guyot

, Flamant

. An atlas of thyroid hormone receptors’ target genes in mouse tissues. Int J Mol Sci, 2022; 23(19):11444; doi: 10.3390/ijms231911444

35.

Nojima

, Miyazaki

, Hori

, et al. Assessment of possible contributions of hyaluronan and proteoglycan binding link protein 4 to differential Perineuronal net formation at the calyx of held. Front Cell Dev Biol, 2021; 9:730550; doi: 10.3389/fcell.2021.730550

36.

Edamatsu

, Miyano

, Fujikawa

, et al. Hapln4/Bral2 is a selective regulator for formation and transmission of GABAergic synapses between Purkinje and deep cerebellar nuclei neurons. J Neurochem, 2018; 147(6):748–763; doi: 10.1111/jnc.14571

37.

Popelář

, Díaz Gómez

, Lindovský

, et al. The absence of brain-specific link protein Bral2 in perineuronal nets hampers auditory temporal resolution and neural adaptation in mice. Physiol Res, 2017; 66(5):867–880; doi: 10.33549/physiolres.933605

38.

Batista

, Hensch

. Critical period regulation by thyroid hormones: Potential mechanisms and sex-specific aspects. Front Mol Neurosci, 2019; 12:77; doi: 10.3389/fnmol.2019.00077

39.

Quignodon

, Legrand

, Allioli

, et al. Thyroid hormone signaling is highly heterogeneous during pre- and postnatal brain development. J Mol Endocrinol, 2004; 33(2):467–476; doi: 10.1677/jme.1.01570

40.

Durand

, Raff

. A cell-intrinsic timer that operates during oligodendrocyte development. Bioessays, 2000; 22(1):64–71.

41.

Picou

, Fauquier

, Chatonnet

, et al. Deciphering direct and indirect influence of thyroid hormone with mouse genetics. Mol Endocrinol, 2014; 28(4):429–441; doi: 10.1210/me.2013-1414

42.

Kalish

, Barkat

, Diel

, et al. Single-nucleus RNA sequencing of mouse auditory cortex reveals critical period triggers and brakes. Proc Natl Acad Sci U S A, 2020; 117(21):11744–11752; doi: 10.1073/pnas.1920433117

43.

Van Vliet

, Deladoey

. Diagnosis, treatment and outcome of congenital hypothyroidism. Endocr Dev, 2014; 26:50–59; doi: 10.1159/000363155

44.

Yuen

, Thiruvahindrapuram

, Merico

, et al. Whole-genome sequencing of quartet families with autism spectrum disorder. Nat Med, 2015; 21(2):185–191; doi: 10.1038/nm.3792

45.

Kalikiri

, Mamidala

, Rao

, et al. Analysis and functional characterization of sequence variations in ligand binding domain of thyroid hormone receptors in autism spectrum disorder (ASD) patients. Autism Res, 2017; 10(12):1919–1928; doi: 10.1002/aur.1838

46.

Quast

, Reh

, Caiati

, et al. Rapid synaptic and gamma rhythm signature of mouse critical period plasticity. Proc Natl Acad Sci U S A, 2023; 120(2):e2123182120; doi: 10.1073/pnas.2123182120

47.

Reh

, Dias

, Nelson

3rd , et al. Critical period regulation across multiple timescales. Proc Natl Acad Sci U S A, 2020; 117(38):23242–23251; doi: 10.1073/pnas.1820836117

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

3.74 MB

Thyroid Hormones Act as a Timer for the Postnatal Maturation of Parvalbumin Neurons in Mouse Neocortex

Abstract

Background:

Methods:

Results:

Conclusions:

Keywords

Introduction

Materials and Methods

Mouse models and treatments

RNA sequencing analysis from sorted nuclei

RT-qPCR analysis

Brain slice histology

Image analysis

Results

TRα1E395fs401X expression in GABAergic neurons from an early developmental stage permanently inhibits PV neuron maturation

Postnatal expression of TRα1L400R in GABAergic neurons is sufficient to inhibit PV neuron maturation

Consequences of TRα1L400R expression on gene expression in neocortical GABAergic neurons

Discussion

Study limitation

Conclusions

Authors’ Contributions

Footnotes

Acknowledgments

Author Disclosure Statement

Funding Information

Supplemental Material

References

Supplementary Material

TRα1^E395fs401X expression in GABAergic neurons from an early developmental stage permanently inhibits PV neuron maturation

Postnatal expression of TRα1^L400R in GABAergic neurons is sufficient to inhibit PV neuron maturation

Consequences of TRα1^L400R expression on gene expression in neocortical GABAergic neurons