In utero gene expression in the Slc39a8(neo/neo) knockdown mouse

Slc39a8 encodes ZIP8, a divalent cation/bicarbonate symporter expressed in pluripotent mouse embryonic stem cells, and therefore ubiquitous in adult tissues; ZIP8 influxes Zn2+, Mn2+ and Fe2+. Slc39a8(neo/neo) knockdown mice exhibit 10–15% of wild-type ZIP8 mRNA and protein levels, and show pleiotropic phenotype of stunted growth, neonatal lethality, multi-organ dysmorphogenesis, and dysregulated hematopoiesis manifested as severe anemia. Herein we performed RNA-seq analysis of gestational day (GD)13.5 yolk sac and placenta, and GD16.5 liver, kidney, lung, heart and cerebellum, comparing Slc39a8(neo/neo) with Slc39a8(+/+) wild-type. Meta-data analysis of differentially-expressed genes revealed 29 unique genes from all tissues — having enriched GO categories associated with hematopoiesis and hypoxia and KEGG categories of complement, response to infection, and coagulation cascade — consistent with dysregulated hematopoietic stem cell fate. Based on transcription factor (TF) profiles in the JASPAR database, and searching for TF-binding sites enriched by Pscan, we identified numerous genes encoding zinc-finger and other TFs associated with hematopoietic stem cell functions. We conclude that, in this mouse model, deficient ZIP8-mediated divalent cation transport affects zinc-finger (e.g. GATA proteins) and other TFs interacting with GATA proteins (e.g. TAL1), predominantly in yolk sac. These data strongly support the phenotype of dysmorphogenesis and anemia seen in Slc39a8(neo/neo) mice in utero.

A Slc39a8 inducible-global knockout and a Slc39a8 hepatocyte-specific conditional knockout, created independently, exhibited striking decreases in Mn 2+ concentrations in multiple organs tested; authors also found lowered levels of two Mn 2+ -dependent enzymes -arginase and β-1,4-galactosyltransferase activities -and evidence that hepatic ZIP8 rescues Mn 2+ from bile and regulates whole-body Mn 2+ homeostasis 14 . The Slc39a8 global knockout shows striking cardiac extracellular matrix accumulation, suggesting that Slc39a8 expression is essential for development of ventricular compaction 15 .

Analysis of differentially-expressed genes in individual tissues. Differential-expression analysis
started with low-expression filtering, following which 6,249 genes were removed and 17,163 genes retained (having CPM >1 in at least three samples). The number of significant differentially-expressed genes ( Table 2) included 15 in yolk sac, two in placenta, one in liver, 13 in lung, but none in kidney, heart and cerebellum. Two genes (Adh7 and Zbtb8b) encode Zn 2+ -containing proteins. A heat map (Supplementary Data Figure S1) shows intensities of each gene expression.
In order to get sufficient numbers of differentially-expressed genes, a relaxed cut-off was used (significant at P < 0.05 and absolute fold-change >2) for all statistical functional enrichment analyses. KEGG pathways affected by differentially-expressed genes comprised: in yolk sac, two ("malaria" and "African trypanosomiasis") down-regulated in Slc39a8(neo/neo) compared with wild-type; in kidney one down-regulated ("tight junction"); and in lung two down-regulated ("complement" and "coagulation cascades") and one up-regulated ("Staphylococcus aureus infection"). Significant KEGG categories included "complement", "response to infection", and "coagulation cascade". All these pathways are consistent with ZIP8 deficiency causing dysregulated hematopoietic stem cell fate in Slc39a8(neo/neo).
Global Transcriptome Analysis. RNA-seq analysis confirmed that expression intensities are distinctly different across the seven tissues examined (Fig. 2). Subsequently, global analysis revealed that the number of significantly differentially-expressed genes totaled 695 for all tissues combined (Supplemental Table S1). There were 646 unique genes. Gbp1 (guanylate-binding protein-1) was found differentially-expressed in five tissues; there were six genes in three tissues (Fig. 3A), and 33 genes in two tissues (Fig. 3B). This resulted in 40 genes (including Slc39a8) differentially-expressed in more than one tissue (Supplemental Table S2). The remaining 606 genes were differentially-expressed in only one tissue. Notable among the 40 genes are four Slc genes (other than Slc39a8) and two hemoglobin genes.
More than 80 differentially-expressed genes (Supplemental Table S1) were associated with hematopoiesis-and hypoxia-related functions (red font); these include all myelogenesis functions (e.g. cytokines, interferons, innate immunity, and inflammatory functions), because all are derived from hematopoietic stem cells. During mouse in utero development, hematopoiesis transitions from the aorta-gonad-mesonephros (AGM) region and yolk sac (between GD8 and GD13), then to liver (GD11 to GD20), and also to spleen (after GD15.5), and finally to bone marrow after GD17.5. Mouse (but not human) placenta is also a hematopoietic organ 39 .
Many other classes of genes that are seen (Supplemental Table S1) -includes those encoding: Zn 2+ -finger transcription factors (TFs) and homeobox and other very-early-embryonic functions associated with cell cycle and cell division; other Zn 2+ -containing proteins; proteins that are posttranslationally glycosylated; enzymes that participate in lipid and glucose homeostasis; growth factors (including several oncogenes and proto-oncogenes); and G-protein-coupled receptors. Nine cytochrome P450 (Cyp) genes, and 27 solute carrier (Slc) genes other than Slc39a8, were observed; explanations and speculations about members of these two gene superfamilies are discussed in Supplemental Table S1.
Systematic meta-analysis. Assuming each tissue is independent from all other tissues, we searched for differentially-expressed genes that tend to change in the same direction in all tissues. A meta-analysis based on Fisher's combined probability test was applied to integrate empirical P-values across all tissues. Forty-seven genes tended to be consistently up-or down-regulated in all tissues (Supplemental Table S3). Not surprisingly, Slc39a8 was the most consistently down-regulated gene in all seven tissues with a combined P-value of 3.5e-04 after FDR adjustment. Six hemoglobin genes were also significantly down-regulated. Heat maps illustrate the 45 down-regulated genes in all tissues (Fig. 4A), and the two up-regulated in all tissues (Fig. 4B)

Analysis of transcription factor (TF)-binding sites in genes identified by meta-analysis.
Meta-analysis assumes there is a common impact of markedly decreased ZIP8 expression, genome-wide, across all tissues. Transcription is the first step in gene expression, and transcription is modulated by interaction of TFs, having their corresponding binding-sites within specific DNA modules of expressed genes.
After removing Slc39a8 from Pscan and searching for JASPAR TF motifs in promoter regions, and using a Bonferroni-corrected P-value of <0.1, among the "consistently down-regulated genes" discovered by meta-analysis (Supplementary Table S3), we found only two significant motifs -the Tal1::Gata1 heterodimer (JASPAR ID: MA0140.1) and Nfya (JASPAR ID: MA0060.1). TAL1 and GATA are of critical importance (vide infra). NFYA transcription factors participate in open chromatin of many types of stem cells -participating in differentiation, cell-cycle and transcriptional regulators, splicing and signaling factors, fatty acid and glucose homeostasis, and mitochondrial regulators 40 . Searching for JASPAR TF-binding sites, using the same criterion for the two "consistently up-regulated genes" (Supplemental Table S3), we found no significant motifs.

Analysis of TF-binding sites in individual tissues. Yolk Sac.
Based on TF profiles in the JASPAR database, we searched in yolk sac for TF-binding sites enriched by Pscan (Table 3). We found nine up-and down-regulated genes having TF-binding sites significantly enriched in promoters: Plag1, Arnt, Hif1a, Znf423, Klf13, Sp1, Gata1, Gata2, and Gata4. Note that several TF-binding sites are duplicates with different matrix ID

Organ or system
Properties/Phenotypes of mammalian SLC39A8 gene (and its encoded ZIP8 transporter) References

Central nervous system
Mouse Slc39a8(neo/neo) knockdown fetus shows trend a in smaller size of cerebrum and cerebellum; Human SLC39A8 variant associated with increased risk of schizophrenia; Human SLC39A8 variant associated with impaired intellectual ability and cerebellar atrophy 13,[18][19][20]33,34 Eye Human SLC39A8 variant associated with retinal iron accumulation 35 Spleen Mouse Slc39a8(neo) allele is associated with highly significantly decreased spleen size 13 Gastrointestinal tract Human SLC39A8 variant associated with Crohn disease and human gut microbiome composition 36 Musculoskeletal system Mouse Slc39a8(neo/neo) knockdown fetus exhibits shortened limbs, and deformed skull; Human SLC39A8 variant associated with pathogenesis of osteoarthritis; Human SLC39A8 variant associated with severe limb malformations 13,18,19,37,38 Reproductive system Mouse Slc39a8 variant is associated with Cd 2+ -induced testicular necrosis 1-3,16 Table 1. Properties and phenotypes reported, to date, for mammalian SLC39A8 gene. a The term "trend'" denotes P-value > 0.05 < 0.10. b ALT, alanine aminotransferase (common measurement used to assess damage largely in liver); AST, aspartate aminotransferase (common measurement used to assess damage in heart, skeletal muscle, kidney and brain, as well as liver. Scientific REPoRts | (2018) 8:10703 | DOI:10.1038/s41598-018-29109-y numbers. It is noteworthy that -except for Arnt and Hif1a, associated with the downstream hypoxic response [reviewed in 41 ] -the remaining seven are Zn 2+ -finger TFs. Pleomorphic adenoma gene-1 (Plag1) is one of three members in the PLAG family and is best known as an oncogene associated with certain cancers, most notably pleiomorphic tumors of salivary gland; PLAG1 participates in cell proliferation by directly regulating numerous target genes, including growth factors 42  Lung. We found 11 significant TF-binding sites in significantly enriched in promoters of up-and down-regulated genes (Table 3) : Hnf4a, Hnf4g, Foxa2, Hnf1b, Hnf1a, Foxa1, Cux1, Cux2, Cebpe, Bhlha15, and Rxrg. Of these 11 genes, the four Hnf genes encode Zn 2+ -finger TFs. The hepatocyte nuclear factor (Hnf) family comprises four members that contribute to innumerable critical-life processes. Forkhead genes, which previously had been termed Hnf genes, are now officially called Foxa1, Foxa2, Foxa3 and Foxm1 [https://www.genenames.org/]. Cux genes code for evolutionarily highly conserved homeobox-domain proteins 44 . With regard to Cebpe, Bhlha15 and Rxrg -CEBPE is a leucine-zipper TF associated with several myelogenous leukemias 45 , BHLHA15 is a member of the basic-helix-loop-helix (bHLH) family of TFs involved in numerous developmental programs 46 , and RXRG is a member of the nuclear receptor superfamily that also participates in many morphogenesis processes 47 .
Kidney. Prdm1 was the only TF-binding site significantly enriched in promoters of up-and down-regulated genes ( Table 3). Prdm1 acts as a tumor suppressor gene involved in B-cell and T-cell lymphomas 48 ; the PRDM family contains 15 members of epigenetic regulators that participate in autoimmune diseases and infections 49 .
Heart. Ebf1 was the only TF-binding site significantly enriched in promoters of up-and down-regulated genes ( Table 3). EBF1 is regarded as a master regulator during B-cell differentiation and also plays a role in B-cell malignancies 50 .
Placenta, liver, and cerebellum. We found no statistically significant TF-binding sites as significantly enriched in promoters of up-and down-regulated genes in these tissues.

Discussion
Because ZIP8 is a known uptake transporter of Zn 2+ , Mn 2+ and Fe 2+ -functions of these divalent cations are briefly reviewed.
Zinc. Intracellular Zn 2+ is critical in homeostasis-related signal transduction 51 , cell cycle and proliferation 52 , numerous processes that occur during development and differentiation, and maintenance of many Zn 2+ -requiring functions. In mammals there are >100 Zn 2+ -dependent enzymes 53 and >2,000 Zn 2+ -containing TFs 54 . An independent analysis estimated there are ~2,800 Zn 2+ -binding proteins, corresponding to ~10% of the human proteome 55 . Because these enzymes and TFs carry out vital critical-life functions throughout development -often exerting cell-specific effects on morphogenesis, growth, and differentiation -an embryo's ability to uphold Zn 2+ homeostasis is essential from the single-cell-zygote stage onward 56 . Thus, when Slc39a8(neo/neo) is compared with Slc39a8(+/+), it is likely that many genes listed in Supplemental Table S1 are differentially-expressed as a direct, or indirect, manifestation of deficient ZIP8-mediated Zn 2+ transport.
Manganese. Mn 2+ is also a potent substrate for ZIP8 7,16 . There are several known defects in ZIP-mediated Mn 2+ transport [reviewed in 57 ]. A human SLC39A8 nonsynonymous variant, resulting in decreased ZIP8 transport, impairs the function of Mn 2+ -dependent enzymes -most notably β1,4-galactosyltransferase, essential for biosynthesis of the carbohydrate moiety of glycoproteins; compromised galactosylation in children causes a disorder with deformed skull, severe seizures, short limbs, profound psychomotor retardation, and hearing loss 19 . Because numerous proteins are posttranslationally glycosylated, many genes listed in Supplemental  Table S1 might be differentially-expressed as a reflection of defective ZIP8-mediated Mn 2+ transport. There are no Mn 2+ -containing TFs.
Iron. ZIP8, along with its closest evolutionarily-related ZIP14 transporter, is a Fe 2+ -transport protein 5,58 . These data might be consistent with studies in the intact animal 13 -showing that Slc39a8(neo/neo) fetuses exhibit dysregulated hematopoiesis -suggestive of an iron-transport defect. As noted (vide supra), Supplemental Table S1, highlighted in red font, lists >80 differentially-expressed genes associated with dysregulated hematopoietic stem cell fate and subsequent downstream events including hypoxia, when Slc39a8(neo/neo) is compared with wild-type; how many of these genes might reflect the importance of Fe 2+ transport more so than Zn 2+ transport, will require further studies. Although there are a few Fe 2+ -containing TFs in prokaryotes 59 , there are none in eukaryotes.
Emergence of several critical-life genetic networks. By comparing transcriptomes in the present study, specific genetic pathways are noteworthy (vide infra). There are two possibilities to explain emergence of these genetic pathways in response to globally deficient ZIP8-mediated transport. One, perhaps ZIP8-mediated divalent cation transport might be tissue-or cell type-specific -depending upon ZIP8 levels and the abundance of other redundant transport systems in each particular tissue or cell-type. In other words, it is possible that pathways affected by deficient ZIP8-mediated Zn 2+ transport might be more significant in lung and kidney, those affected by diminished ZIP8-mediated Mn 2+ transport might be more prominent in liver, and those affected Scientific REPoRts | (2018) 8:10703 | DOI:10.1038/s41598-018-29109-y by deficient ZIP8-mediated Fe 2+ transport might be more important in hematopoietic tissues. Two, it is possible that defective Zn 2+ transport, primarily in early-embryo yolk sac, might override (i.e. occur upstream of) Mn 2+ -dependent and/or Fe 2+ -dependent gene products and their functions. These two possibilities will require further experiments, but the present study would argue in favor of the latter, i.e. defective Zn 2+ -finger TFs, primarily in yolk sac, cause irreparable impairment of hematopoietic stem cell fate, which is then manifesteddevelopmentally later -in the five fetal organs studied herein.
Fundamental Importance of TAL1. Tal1 was up-regulated in placenta and down-regulated in kidney (Supplementary Table S2). Differentially-expressed Gata genes were not seen among the 646 unique genes (Supplementary Table S1); this could reflect normal transcription, but then defective function due to its posttranslational Zn 2+ -and/or Mn 2+ -binding requirements. After searching JASPAR TF motifs in promoter regions, among "consistently down-regulated genes" by meta-analysis (Supplemental Table S3), we found the Tal1::Gata1 heterodimeric complex (JASPAR ID: MA0140.1) as an important TF-binding motif. This finding is strengthened when we saw that TF-binding sites enriched by Pscan included Gata1, Gata2 and Gata4 in yolk sac (Table 3).
TAL1 is a bHLH TF that plays an important role in hematopoietic stem cell development 60 . GATA TFs are Zn 2+ -finger DNA-binding proteins that participate in various biological processes -including hematopoiesis and T-cell development [reviewed in ref. 43 ]. Shifts in TAL1 occupancy of TF-binding sites during erythroid differentiation are associated with gene repression (i.e. dissociation), and activation (i.e. co-occupancy with GATA1); in fact, recruitment by GATA TFs appears to be a stronger determinant of TAL1-binding to chromatin than the canonical E-box binding-site motif 61 .
In order to validate the TAL1 result from PScan, we sought to perform an independent statistical test. From a list of manually-curated 80 TAL1-downstream targets 62 , we intersected the genes in this list with those in Supplemental Table S1 for all tissues (Supplemental Table S4). TAL1 targets were significantly enriched in yolk sac (Fisher's exact test P = 1.178e-05), placenta (P = 4.982e-06), and marginally enriched in lung (P = 0.01895). Twelve hematopoiesis-related genes were seen, predominantly in yolk sac and placenta. Egln3 was the most strongly up-regulated gene in yolk sac; EGLN3 is a bHLH TF that responds to hypoxia by decreasing the production of HIF-regulated angiogenic factors 63 . Cbfa2t3 was the most dramatically down-regulated gene in placenta; CBFA2T3 is clinically associated with several types of human leukemia. The only gene seen in heart was Zfpm1; ZFPM1 is a Zn 2+ -finger TF that (like TAL1) also forms heterodimers with TFs of the GATA family.
Moreover, TAL1 is known to activate Egln3, Hemgn and Ank1 62 , and all were up-regulated in Slc39a8(neo/ neo) yolk sac. TAL1 represses Hbb-b2 62 , which was down-regulated in Slc39a8(neo/neo) yolk sac (Supplemental Table S4). These data suggest that TAL1 protein activity might be increased in Slc39a8(neo/neo) yolk sac, whereas Tal1 gene differential expression was not seen. Similarly, in Slc39a8(neo/neo) placenta, Hemgn and Ank1 were up-regulated and Hbb-b2 down-regulated -suggesting that TAL1 protein might also be activated in that tissue.
None of the genes encoding Zn 2+ -finger TFs listed in Table 3 was differentially expressed to a level seen in Supplemental Table S1, i.e. these genes all appear to be transcribed somewhat normally. Thus, we suggest that, although transcription of these Zn 2+ -finger genes was within normal limits, only the protein functions of these Compared with Slc39a8(+/+) wild-type and Slc39a8(+/neo) heterozygotes that are pink in color and normal in size, the Slc39a8(neo/neo) littermates are extremely pale, show stunted growth, and deformed skulls and limbs. Slc39a8(neo/neo) liver, kidney, lung, spleen, cerebrum and cerebellum were all statistically significantly smaller in size than that in Slc39a8(+/+) or Slc39a8(+/neo) littermates 13 . For the RNA-seq analysis described herein, Slc39a8(+/neo) littermates were not studied.  Table 2. Genes differentially regulated with FDR < 0.1, when comparing Slc39a8(neo/neo) with Slc39a8(+/+) wild-type. a RPKM, "number of Reads-Per-Kilobase-of-transcript-per-Million mapped reads". b Genes that are erythropoiesis-and hypoxia-related -including those involved in iron transport; these also include all platelet and white-cell types and functions because all are derived from hematopoietic stem cells. c These two genes encode Zn 2+ -containing proteins. ADH7 is an enzyme, and ZBTB8B is a Zn 2+ -finger TF. d Although 31 genes are listed, these two are not statistically significant (i.e. P > 0.05 when P adj is taken into account). Fundamental Importance of GATA. The hematopoietically-expressed GATA family members of Zn 2+ -finger TFs are well known to function as key regulators of blood cell fate; for example, phenotypes of knockout mice -deficient in Gata1, Gata2, or Gata3 -suggest that these factors each play critical, but distinctly different, roles in hematopoiesis 64 . GATA2 mutations are clinically associated with familial myelodysplastic syndrome 65,66 , by affecting monocytes and/or B cells. In order to validate the GATA1 result from Pscan (Table 3), we used a list of 512 GATA1-downstream target genes 67 ; this list includes most of the well-established GATA1-mediated erythroid-specific targets, e.g. the Gata1 gene itself, Gata2, the β-globin locus (specifically the locus-control region), Epor, Nfe2, Slc4a1, Gypa, Tal1, Lrf, Klf1, Nrf2, Runx1 and Alas2. We intersected the genes in this list with those in Supplemental Table S1 for all Slc39a8(neo/neo) tissues (Supplemental Table S5). Yolk sac was the tissue with the most enriched number of GATA1 targets: excluding Slc39a8, there were 13 differentially-expressed GATA1-downstream target genes out of the 512 total (Fisher's exact test P = 0.001369). These findings provide strong evidence for key involvement of GATA1 in many of the differentially-expressed transcriptomic changes in Slc39a8(neo/neo) yolk sac.

NO SIGNIFICANT GENES FOUND
We found no enrichment of GATA1 targets among differentially-expressed genes in liver, kidney, heart or cerebellum. These data convincingly indicate that deficient ZIP8-mediated transport affects hematopoiesis-related processes primarily in yolk sac -as opposed to effects in placenta, or developmentally later, in early fetal liver.
Defective Hematopoiesis Leads to Hypoxia. Throughout our analysis presented herein, dysregulation of hematopoiesis-related stem cell functions is repeatedly described; decreases in hemoglobin content, red blood cell number, and iron uptake would likely elicit a "hypoxia signal" in many cell types -resulting in activation of hypoxic-response downstream-target genes. The striking anemia phenotype, seen in the earliest of visible Slc39a8(neo/neo) yolk sac, placenta, and embryos [ Fig. 1 and ref. 13 ], supports these findings.
A list of 87 manually-curated HIF-downstream targets 68 has been reported; examples include genes involved in fatty acid and glucose homeostasis, glycolysis pathway, cell proliferation, cell migration, and angiogenesis. We intersected the genes in this list with those in Supplemental Table S1 for all Slc39a8(neo/neo) tissues (Supplemental Table S6). Fifteen differentially-expressed genes were detected in yolk sac (Fisher's exact test P = 4.206e-15), one in kidney, two in lung, and none in placenta, liver, heart or cerebellum; as expected, all 15 of these were up-regulated (Supplemental Table S6). This independent statistical test confirms our Hif1a result from Pscan (Table 3).
In lung, Abcb1a was up-regulated; Cyp2s1 and Trf down-regulated. ABCB1A ("P-glycoprotein") is up-regulated in response to many types of stress, including hypoxia 69 . In the presence of hypoxia, HIF1A/ARNT dimers are known to bind to an aryl hydrocarbon receptor response element (AHRE) upstream of the Cyp2s1 gene 70 ; furthermore, the Arnt::Hif1a heterodimer TF-binding site enriched by Pscan was found to be up-regulated in yolk sac (Table 3). Trf encodes transferrin, which is important in Fe 2+ transport 71 ; defective hematopoiesis resulting in anemia and subsequent hypoxia is the likely cause of down-regulation of Trf in lung.
In kidney one HIF-downstream target gene (Met up-regulated) was found. Met, a proto-oncogene, participates in angiogenesis -which would be stimulated by hypoxia 72 . Hif1a expression itself was found not to be differentially-expressed in any tissue, whereas Hif3a was up-regulated in yolk sac (Supplementary Table S1). Hif1a and Hif3a in mice encode hypoxia-inducible TF subunits that respond to low oxygen (pO 2 ) levels. Either of these TFs can bind as a heterodimeric complex with constitutively-expressed ARNT (aryl hydrocarbon receptor nuclear translocator), product of the Arnt gene.
Members of the Egl-9 family of HIF-inducible TFs (Egln1, Egln2 and Egln3) were also differentially-expressed (Supplemental Table S1). All six of these above-mentioned genes are members of the bHLH per-Arnt-sim (bHLH/ PAS) family of TFs that detect endogenous and exogenous signals; the bHLH/PAS family contains at least 30 members in the human and mouse genomes 41 . In fact, two other bHLH/PAS differentially-expressed genes (Bhlhe40 & Bhlhe41) were found to be up-regulated in yolk sac (Supplemental Table S1) -again supporting the theme of the Slc39a8(neo/neo) response to hypoxic stress.

Ingenuity Pathway
Analysis. An independent Ingenuity Pathway Analysis (IPA) for the 165 differentiallyexpressed genes in yolk sac (Supplementary Table S7) revealed major pathways that were consistent with the data discussed throughout the text. Specifically, the top candidate from "Upstream Regulators Analysis" was the transcription factor HIF1A, and it is predicted to be activated. The top categories from the "Diseases and Bio Functions Analysis" included "Transport of molecule, " "Morphology of head, " "Glycolysis of cell, " "Morphology of nervous system, " and "Familial hemolytic anemia" -all of which confirms the data presented herein -as well as the associations of human SLC39A8 variants with most if not all of the clinical disorders (vide supra). With a FDR-adjusted P-value < 0.05 from "Canonical Pathway Analysis, " no significant canonical pathways were found.

Conclusions
A summary of many of these differentially-expressed genes, and their resultant effects, elicited by ZIP8 deficiency and described throughout this paper, is offered in Fig. 5.
To our knowledge, we know of no studies, reported to date, that dissect the mechanism behind why or how the knockdown of one critically important transporter gene can result in the up-and down-regulation of so many dozens of other transporter genes -and many of the transported moieties are not ions. These intriguing data represent many new avenues of possible experiments for future fascinating studies of whole-genome interactions.  Table 3. Significant transcription factor (TF)-binding sites in genes having differential-expression at P < 0.05 and absolute fold-change >2; TF-binding site enrichment by Pscan a . a Using this rigorous method of analysis, we found no significant TF-binding sites in placenta, liver or cerebellum. b TF name denotes the TF-binding site(s) enriched by Pscan. c Although 30 enriched TF-binding sites are listed, these five are not statistically significant (i.e. P > 0.05 when the Bonferroni-corrected P-value is taken into account).
All mouse experiments were conducted in accordance with the National Institutes of Health standards for the care and use of experimental animals. All experimental protocols were approved by the University of Cincinnati College of Medicine (IACUC) Institutional Animal Care and Use Committee [protocol #11-09-12-01].
Total RNA Extraction. RNeasy Micro kit (Qiagen) was used for total RNA extraction. Each sample, together with Buffer RLT, were added to ceramic-beads (1.4-mm) in tubes and homogenized by using the Precellys homogenizer (Cayman Chemical Company; Ann Arbor, MI). RNA extraction was then performed, following manufacturer's recommendations.
RNA-Seq for differential gene expression profiling. For each RNA-seq sample, 1 μg RNA was used (which required two or more mice per sample). RNA-seq was performed at the Genomics, Epigenomics and Sequencing Core at the University of Cincinnati. To prepare the RNA library for sequencing, we used TruSeq RNA Library Prep Kit (Illumina; San Diego, CA), according to manufacturer's instructions. In brief, poly(A+) RNA was purified from total RNA and then fragmented before cDNA synthesis. Double-strand cDNA fragments were then end-repaired and TA-ligated to the sequencing adapter; the specific index was added to the library during the PCR amplification step. After 12 cycles of PCR enrichment, library quality was evaluated by a Bioanalyzer High Sensitivity chip and quantified by Kapa Library Quantification kit (Kapa Biosystem; Wilmington, MA), using the ABI 9700HT real-time PCR system (Thermo Fisher; Waltham, MA). Next, six individually indexed cDNA libraries were pooled in equal amounts for clustering in the cBot system (Illumina). Libraries were clustered onto a flow cell at concentration of 15 pM, using Illumina's TruSeq SR Cluster Kit v3, and then sequenced for 50 cycles using TruSeq SBS kit on the Illumina HiSeq system. Each sample was expected to generate ~30 million reads.

Figure 5.
Illustration of critically affected genes and their downstream effects that most closely fit the data presented in the present study. ZIP8 deficiency (top), has a major impact on TAL1 and GATA transcription factors (TFs), which appear to function primarily in the hematopoietic stem cells of Slc39a8(neo/neo) GD13.5 yolk sac. This function then causes severe dysregulation of hematopoietic stem cell fate and striking anemia in yolk sac, which is visibly obvious [in Fig. S2 of ref. 13 ]. Moreover, hematopoiesis in yolk sac is well known to precede that in liver and then spleen and marrow [cf. Fig. S4 of ref. 13 ]. Downstream effects, as development proceeds, include severe anemia and defects in coagulation, innate immunity, and response to inflammation. The striking anemia leads to a hypoxia response which is seen in all tissues examined, but largely in yolk sac. TAL1, T-cell acute lymphocytic leukemia protein-1 TF. GATA, family of six zinc-finger TFs that regulate hematopoietic stem cell fate. Interactions between TAL1 and GATA exist 61 [see text]. ZIP8-deficiency, plus the result of all these downstream changes, alter the expression of nine Cyp genes and 27 Slc genes (excluding Slc39a8); these changes are mostly unique to one tissue, as detailed in Supplementary Table S1. Δ denotes "changes in". HIF, hypoxia-inducible factor.