The iodide transporter Slc26a7 impacts thyroid function more strongly than Slc26a4 in mice

SLC26A4 is a known iodide transporter, and is localized at the apical membrane of thyrocytes. Previously, we reported that SLC26A7 is also involved in iodide transport and that Slc26a7 is a novel causative gene for congenital hypothyroidism. However, its detailed role in vivo remains to be elucidated. We generated mice that were deficient in Slc26a7 and Slc26a4 to delineate differences and associations in their roles in iodide transport. Slc26a7−/− mice showed goitrous congenital hypothyroidism and mild growth failure on a normal diet. Slc26a7−/− mice with a low iodine environment showed marked growth failure. In contrast, Slc26a4−/− mice showed no growth failure and hypothyroidism in the same low iodine environment. Double-deficient mice showed more severe growth failure than Slc26a7−/− mice. RNA-seq analysis revealed that the number of differentially expressed genes (DEGs) in Slc26a7−/− mice was significantly higher than that in Slc26a4−/− mice. These indicate that SLC26A7 is more strongly involved in iodide transport and the maintenance of thyroid function than SLC26A4.

and SLC26A4 in the mouse thyroid, we generated mice deficient in these proteins. The targeting constructs for the disruption of Slc26a7 and Slc26a4 were designed to create double-strand breaks in exons 5 and 2 of the Slc26a7 and Slc26a4 genes, respectively ( Fig. 1a-c). We generated heterozygous mutant mice with frameshift variants and crossed them to create deficient mice (Slc26a7 −/− and Slc26a4 −/− ). Genomic DNA and cDNA sequencing confirmed that the frameshift variants were homozygous ( Supplementary Fig. S1a,e). Furthermore, qRT-PCR and RNA-seq demonstrated that Slc26a7 −/− and Slc26a4 −/− mice showed significantly reduced mRNA expression of the corresponding gene compared to WT mice ( Supplementary Fig. S1b,c,f,g). In addition to Slc26a7 −/− mice, we generated another strain of Slc26a7-modified mice, which we designated as Slc26a7 del/del mice. This strain was found to harbor a 205-base pair deletion that contained exon 5 and the flanking intron in the gene, which resulted in the deletion of exon 5 during splicing ( Fig. 1b and Supplementary Fig. S1d). Exon 5 encodes transmembrane segment (TMS) 5, which is crucial for SLC26A7 function 7 . To generate double-deficient mice (Slc26a7 del/del Slc26a4 −/− ), Slc26a7 +/del Slc26a4 +/− female mice were crossed with Slc26a7 +/del Slc26a4 +/− male mice. Body weights, thyroid sections, and serum FT4, FT3, and TSH levels in WT, Slc26a7 −/− , and Slc26a4 −/− mice were analyzed by phenotypic analysis. Slc26a7 −/− mice had evident goiters, and their thyroid weights were significantly higher than those of WT mice (Fig. 1d,e). The characteristic follicular structure was lost in the thyroid glands of Slc26a7 −/− mice (Fig. 1f). The serum FT4 levels were significantly lower and the serum TSH levels were significantly higher in Slc26a7 −/− mice than in WT mice (Fig. 2a,b). The serum FT3 levels in Slc26a7 −/− mice were comparable to those in WT mice at day 45, but were significantly lower at day 90 ( Fig. 2a,b). Moreover, Slc26a7 −/− mice showed mild growth failure compared to WT mice (Slc26a7 −/− : 24.1 ± 0.71 g vs. WT: 27.5 ± 0.46 g; p = 0.0008) (Fig. 1g,h). In contrast, Slc26a4 −/− mice did not show hypothyroidism or growth failure (Figs. 1g, 2b), which is consistent with a previous report 6 . Slc26a4 −/− mice showed no phenotypic abnormality with reference to the thyroid gland (Figs. 1f, 2b), but showed characteristic degeneration of sensory hair cells in the cochlea (Supplementary Fig. S2a,b) 11 .
There were no differences in thyroid morphology or thyroid hormone levels between Slc26a7 −/− and Slc26a7 del/del mice ( Supplementary Fig. S3a,b).
Double-deficient mice showed severe growth failure. We also analyzed double-deficient mice to confirm the phenotype of mice with impairment of both the iodide transporters at the apical membrane of thyroid follicular cells. The size and histology of the thyroid glands of the double-deficient mice were almost identical to those of Slc26a7 −/− mice (Fig. 1f). Similar to Slc26a7 −/− mice, the double-deficient mice had significantly lower levels of FT4 and FT3 and significantly higher levels of TSH compared to WT mice (Fig. 2b). When assessed for body weights, the double-deficient mice showed severe growth failure compared to Slc26a7 −/− mice (double-deficient: 19.0 ± 0.57 g vs. Slc26a7 −/− : 24.1 ± 0.71 g; p < 0.0001) (Fig. 1g,h). However, there was no statistically significant difference between Slc26a7 −/− and double-deficient mice in thyroid function tests. Administering L-T4 to the double-deficient mice tended to improve their growth to some extent ( Supplementary Fig. S4a,b).
Iodine environment greatly influenced the phenotype of Slc26a7 −/− mice. To investigate the effect of iodine intake on the phenotypes of Slc26a7 −/− and Slc26a4 −/− mice, we fed Slc26a7 +/− and Slc26a4 +/− female mice a low iodine diet for 28 days. We found that urine iodine levels were significantly lower in the low iodine diet group (Fig. 3a). Subsequently, these mice were crossed with male mice of the same genotype (Slc26a7 +/− or Slc26a4 +/− ), and their offspring were raised in the same low iodine environment. The genotype ratios of the offspring were consistent with the expected single-locus Mendelian inheritance pattern (Slc26a7 +/+ : Slc26a7 +/− : Slc26a7 −/− = 34: 60: 27; p = 0.66, Chi-square test). Unlike WT and Slc26a7 +/− mice, Slc26a7 −/− mice showed marked growth failure, and exhibited a 100% mortality rate post-weaning (Fig. 3b). The serum FT4 levels of Slc26a7 −/− mice at day 21 (at weaning) were significantly lower than those of WT and Slc26a7 +/− mice (Fig. 3c). Serum FT4 levels in Slc26a7 −/− mice that were fed a low iodine diet tended to be lower than those in Slc26a7 −/− mice that were fed a normal diet (Fig. 3c). This suggests that hypothyroidism in Slc26a7 −/− mice was further aggravated by a low iodine diet. We were able to rescue Slc26a7 −/− mice that had been raised in a low iodine environment from the fetal stage by administering L-T4 and a high iodine diet after weaning (Fig. 3b). The gDNA target sequence is located in exon 5; black and orange underlining indicates the target sequence and PAM sequence, respectively. In the mutant allele, one base insertion (red letter) is introduced, resulting in a frame shift. (b) Strategy to generate Slc26a7 del/del mouse. The target sequence is the same as in Slc26a7 −/− mouse, and a 205 bases deletion here causes exon 5 to be skipped when splicing occurs. (c) Strategy to generate Slc26a4 −/− mouse. The target gDNA sequence is located in exon 2. In the mutant allele, one base duplication (red letter) is introduced, resulting in a frame shift. (d) Photographic image of the thyroid gland in male wild type (WT) and Slc26a7 −/− mice at day 90, demonstrating the presence of a goiter in the Slc26a7 −/− mouse. The area surrounded by dotted lines indicates the thyroid gland. (e) Relative weights of the thyroid gland in male WT, Slc26a7 +/− , and Slc26a7 −/− mice at day 90; n = 4-9 mice of each genotype. Oneway analysis of variance (ANOVA) showed significant differences among the three groups (F (2, 15) = 24.1, p < 0.0001). ***p < 0.001, ****p < 0.0001 determined by Tukey's test. (f) Hematoxylin and eosin-stained thyroid sections from male WT, Slc26a4 −/− , Slc26a7 −/− , and Slc26a7 del/del Slc26a4 −/− mice at day 90. (g,h) Growth curve (g) and body weight at day 90 (h) of male WT, Slc26a4 −/− , Slc26a7 −/− , and Slc26a7 del/del Slc26a4 −/− mice from day 21; n = 5-11 mice of each genotype. One-way ANOVA showed significant differences among the three groups (F (3, 33) = 34.3, p < 0.0001). ***p < 0.001, ****p < 0.0001 determined by Tukey's test. Slc26a7 has a stronger effect on the thyroid phenotype than Slc26a4. To elucidate the mechanisms of iodide transport in thyroid follicular cells and the cross-talk of the genes related to iodide transport, we compared the gene expression profiles of the thyroid glands of WT, Slc26a7 −/− , and Slc26a4 −/− mice at day 90 by RNA-seq analysis. First, we analyzed the differentially expressed genes (DEGs) in Slc26a7 −/− and Slc26a4 −/− mice in comparison with WT mice. The results showed that Slc26a7 −/− mice showed far more DEGs compared to Slc26a4 −/− mice ( Fig. 4a,b), suggesting that Slc26a7 has a stronger effect on the thyroid phenotype than Slc26a4. We then analyzed the differential expression of 26 genes involved in thyroid hormone synthesis and release ( Table 1). The results showed that significant expression changes occurred mainly in the genes involved in iodide transport, with Slc26a4 being the most upregulated (Fig. 4c). This suggests that SLC26A4 may compensate for the loss of function of SLC26A7. Finally, we quantified the expression levels of the genes encoding iodide transporters (Slc26a4, Slc26a7, Slc5a5, Slc5a8, Ano1, Cftr, and Clcn5) using the TPM value obtained by StringTie in WT mice. Results showed that Slc26a7 expression levels were higher than those of Slc26a4 (Fig. 4d).

Discussion
While SLC26A7 deficiency in humans causes goitrous CH, its precise role and the underlying mechanism of CH development are unclear 7,8,12 . In this study, we show that loss of SLC26A7 function in mice leads to goitrous CH (Figs. 1d-f, 2). The serological and morphological findings in Slc26a7 −/− mice were similar to those in Slc5a5deficient mice 13 . In contrast, Slc26a4 −/− mice showed characteristic abnormal inner ear findings ( Supplementary  Fig. S2b), though the thyroid phenotype remained normal. Although Slc5a5, Slc26a4, and Slc26a7 have been identified as iodide transporters and genetic abnormalities in these genes are known to cause CH in humans, phenotypic analysis of deficient mice suggests that the effect of Slc26a4 on the mouse thyroid is less prominent than that of Slc5a5 and Slc26a7 under normal conditions. Slc26a7 −/− mice fed a low iodine diet showed more severe growth failure than those fed a normal diet (Fig. 3b). This implies that iodine intake greatly influences the phenotype of Slc26a7 −/− mice. Moreover, all Slc26a7 −/− mice with a low iodine diet died when therapeutic intervention was not provided; this is presumably due to marked hypothyroidism based on the lower FT4 levels (Fig. 3c), but FT3 levels could not be evaluated at this time.
We have described the differences in the function of SLC26A7 and SLC26A4 in the mouse thyroid based on serological and growth data, and have also demonstrated such differences by RNA-seq analysis. We found that the number of DEGs in Slc26a7 −/− mice was clearly higher than that in Slc26a4 −/− mice (Fig. 4a,b). It is well known that the expression of Slc26a4, Slc5a5, and Tshr is positively regulated by TSH [14][15][16] , and the expression of these genes is upregulated in Slc26a7 −/− mice by elevated levels of TSH. Nevertheless, among the genes involved in thyroid hormone synthesis, Slc26a4 was identified as the most upregulated gene in Slc26a7 −/− mice (Fig. 4c), suggesting that SLC26A4 might have a compensatory relationship with the function of SLC26A7 at the apical membrane of thyroid follicular cells (Fig. 4e). However, we were unable to validate the RNA-seq results by qRT-PCR analysis. Furthermore, it would be more informative if a similar genetic analysis could be performed on   www.nature.com/scientificreports/ Slc26a7 −/− mice treated with L-T4 to normalize the TSH levels, as the results would reflect only the direct effects of SLC26A7 deficiency. These issues must be addressed in future studies. We also generated and analyzed double-deficient mice to verify the relationship and homology between SLC26A7 and SLC26A4. Interestingly, the double-deficient mice showed significant growth failure compared to Slc26a7 −/− mice (Fig. 1g,h). This finding also implies that despite its low expression level, SLC26A4 may functionally complement SLC26A7, at least partially. The growth failure of the double-deficient mice improved to some extent with the administration of L-T4 (Fig. 1g,h). Therefore, it is possible that hypothyroidism may have contributed to the growth failure. However, contrary to our expectations, there was no difference in the degree of hypothyroidism between Slc26a7 −/− and double-deficient mice at day 90 (Fig. 2b). Thus, it is possible that SLC26A4 was not able to complement SLC26A7 deficiency with reference to iodide uptake to an extent that could significantly improve thyroid function at day 90. The possible causes of the weight differences between Slc26a7 −/− and double-deficient mice might include non-thyroidal effects and a complementary role in thyroid function during the pre-and post-weaning periods. We were not able to examine the non-thyroid phenotype of Slc26a7 −/− and double-deficient mice in this study and cannot rule out the possibility that renal tubular acidosis, which has been reported in the past 10 , may affect growth retardation. Evaluation of thyroid function prior to day 90 and renal tubular function may provide more insight into the impact of SLC26A7 and SLC26A4 on phenotype.
In addition, this study has several limitations. First, although we analyzed Slc26a7 and Slc26a4 mRNA, we did not analyze the corresponding proteins due to the lack of specific antibodies. Additionally, the localization of SLC26A7 in the mouse thyroid remains unclear. Second, we could not investigate the pathophysiological mechanism by which functional Slc26a7 deficiency causes goitrous hypothyroidism. In our previous study, we demonstrated that SLC26A7 is involved in iodide transport in vitro 7 . In vivo data from the present study support the in vitro data, indicating that SLC26A7 is involved in iodide transport, and that there is likely to be cross-talk between SLC26A7 and SLC26A4. However, further mechanistic details could not be obtained. Third, because there is a paucity of clinical information on congenital hypothyroidism caused by loss of function of SLC26A7, we do not believe that these results in mice are entirely consistent with the human phenotype.
In conclusion, we revealed that the phenotypes of Slc26a7-and Slc26a4-deficient mice were distinct. SLC26A7 played a more important role in the mouse thyroid than SLC26A4, while cross-talk between SLC26A7 and SLC26A4 was present.

Methods
Animals and treatments. The mice were maintained under specific pathogen-free conditions in a laboratory animal facility at Nagoya City University. All experiments comply with the ARRIVE (Animal Research: Reporting of In Vivo Experiments) guidelines 2.0, were performed in compliance with the relevant Japanese and institutional laws and guidelines, and were approved by the Institutional Animal Care and Use Committee (IACUC) of Nagoya City University School of Medical Sciences (authorization number H30M-42). C57BL/6N mice were purchased from Japan SLC Inc. (Hamamatsu, Japan). Slc26a7 and Slc26a4 mutant mice were generated with a C57BL/6N background using the CRISPR/Cas9 system. For Cas9 guide RNA, the target sequences for Slc26a7 (5′-GCC ATA TAT ATC CGG ACC AC-3′) and SLC26A4 (5′-GTG TAC AGC GAG CTC GCC TT-3′) were selected in exon 5 and exon 2, respectively, and were purchased from Integrated DNA Technology (IDT; Coralville, IA, USA) ( Fig. 1a-c, Supplementary Fig. S1). RiboNucleoProtein (RNP) complexes with Alt-R S.p. Cas9 nuclease V3 (IDT) were transferred to C57BL/6N fertilized oocytes as previously described 17 . F0 mice with the desired mutations were crossed with wild-type (WT) mice to obtain heterozygous F1 mice. These heterozygous F1 mice were then crossed with each other to generate single-deficient mice (Slc26a7 −/− , Slc26a7 del/del , and Slc26a4 −/− ). To generate Slc26a7 and Slc26a4 double-deficient mice (Slc26a7 del/del Slc26a4 −/− ), double heterozygous mice with in-frame deletions in Slc26a7 (Slc26a7 +/del ) and frameshift insertions in Slc26a4 (Slc26a4 +/− ) were created and intercrossed. Genotyping was performed using the Sanger sequencing technique. Genomic DNA was extracted from tail biopsies using the KAPA MG kit (Kapa Biosystems, Woburn, MA), and analyzed using a SeqStudio genetic analyzer (ThermoFisher Scientific, Waltham, MA). The primer sequences used are listed in Supplementary Table 1.
Blood samples were collected from the retroorbital plexus of mice anesthetized with isoflurane; thyroid samples were collected after the mice had died from prolonged anesthesia with isoflurane. Mice were routinely fed a standard chow containing approximately 1.7 mg Iodine (I)/kg (CE-2; CLEA Japan, Tokyo, Japan). Female Slc26a7 +/− and Slc26a4 +/− mice were fed low iodine chow containing less than 500 µg I/kg (CLEA Japan) for analyzing the phenotypes associated with a low iodine diet. After 28 days, the mice were placed in metabolic cages and urine iodine was measured in 24-h urine samples. These mice were bred with Slc26a7 +/− or Slc26a4 +/− male mice. The same iodine environment was maintained during pregnancy and postpartum. Thyroid hormone replacement with albuminized drinking water containing 0.39 mg/L l-thyroxine (L-T4) (L-T4 sodium salt pentahydrate; Med Chem Express, NJ) was started after weaning on day 21. In the case of Slc26a7 −/− mice with a low iodine diet, in addition to providing water containing L-T4, subcutaneous injections of 40 ng/g body weight L-T4 were administered on day 21, 24, and 28; these mice were also fed a high iodine chow containing about 11.9 mg I/kg (CLEA Japan). www.nature.com/scientificreports/ tively. Additionally, the intra-assay variations for each assay were < 15%, < 15%, < 6% for FT4, FT3 and TSH, respectively.

Measurement of serum hormones.
Determination of urine iodine. Urine iodine concentrations in mice were determined by inductively coupled plasma mass spectrometry (ICP-MS) (ELAN DRC II, Perkin Elmer, Waltham, MA, USA). Sample and standard solutions were prepared by alkaline digestion with 25% tetramethyl ammonium hydroxide (TMAH). Iodine was detected at mass 127 I, with 130 Te as the internal standard.
Thyroid histology. For histological analysis, thyroid glands were fixed overnight in 10% neutralized buffered formalin (NBF), washed with 1 × PBS, and embedded in paraffin. The sections were then stained with hematoxylin and eosin using standard procedures.

RNA-Seq analysis.
qRT-PCR analysis. RNA was extracted from the thyroid tissue using the method described earlier, and 500 ng of total RNA was used for cDNA synthesis using SuperScript IV Reverse Transcriptase (Invitrogen, Carlsbad, CA). The cDNA was amplified using the FastStart Essential DNA Green Master Mix (Roche, Mannheim, Germany). All procedures were as per the manufacturer's protocols. A Light Cycler 96 system (Roche) was used for amplification and a large intron was selected. Primers were designed to be located in the exons before and after that intron (see Supplementary Table 2 for primer sequences). Each sample was analyzed in triplicate under the same conditions. The standard curve was prepared by making a stepwise dilution of the cDNA of the sample to be measured. The relative quantification was performed by taking into account the PCR efficiency obtained from the standard curve, and the data were normalized to the GAPDH internal control.