Abstract
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.
Similar content being viewed by others
Introduction
Starting with a Cd2+-induced testicular necrosis phenotype and genetically “sensitive” vs “resistant” inbred mouse strains, this laboratory identified Slc39a8 (encoding the ZIP8 transporter) as the major gene responsible for this trait1,2,3. ZIP8 functions as a Zn2+/(HCO3−)2, Mn2+/(HCO3−)2 or Cd2+/(HCO3−)2 symporter, in each case moving all three ions of the complex into the cell in an electroneutral manner4. ZIP8 has also been shown to import divalent iron and cobalt5. Moreover, ZIP8 transports selenium in the form of selenite into the cell — presumably as a Zn2+/(HCO3−)(HSeO3−) electroneutral complex6. In mammalian cell and Xenopus oocyte cultures, Mn2+ and Fe2+ are able to substitute for the Zn2+ cation5,7,8. It is possible that metal-ion specificity of ZIP8-mediated uptake in the intact animal might be dependent on organ and/or cell-type.
Developmentally, ZIP8 is expressed in mouse visceral endoderm at gestational day (GD)7.5 [ref.9], in the gastrula stage10, and in pluripotent embryonic stem cells11. These findings are compatible with the observation of early embryolethality, seen in Slc39a8(−/−) knockout mice when the gene is globally 100% ablated12.
During creation of a Slc39a8 conditional knockout construct12, the Frt-flanked neomycin-resistance mini-gene neo (transcribed in reverse orientation) was inserted into intron 3, combined with loxP sites in introns 3 and 6. Fortuitous retention of the neo mini-cassette produced a unique hypomorph displaying 10–15% of wild-type levels of ZIP8 mRNA and protein in all tissues examined1,2,3. The Slc39a8(neo/neo) mouse has been extensively characterized: pleiotropic effects include stunted growth, neonatal lethality, shortened limbs and deformed skull, severe anemia, dysregulation of hematopoiesis, and multi-organ dysmorphogenesis; consistent with the severe anemia — striking decreases in size of placenta, and size and number of hematopoietic islands are seen in Slc39a8(neo/neo) GD13.5 yolk sac and placenta, as well as in GD16.5 liver13.
A Slc39a8 inducible-global knockout and a Slc39a8 hepatocyte-specific conditional knockout, created independently, exhibited striking decreases in Mn2+ concentrations in multiple organs tested; authors also found lowered levels of two Mn2+-dependent enzymes – arginase and β-1,4-galactosyltransferase activities – and evidence that hepatic ZIP8 rescues Mn2+ from bile and regulates whole-body Mn2+ homeostasis14. The Slc39a8 global knockout shows striking cardiac extracellular matrix accumulation, suggesting that Slc39a8 expression is essential for development of ventricular compaction15.
These mouse data1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16, summarized in Table 1, underscore the likely fundamental importance of ZIP8-mediated physiological functions in a large variety of tissues and cell-types. The gene was originally discovered in human monocytes17. More recently, a growing number of clinical studies of SLC39A8 variants have shown the pleiotropic effect of deficient ZIP8 function: during early development18,19,20 and in liver14,18,19, kidney21, lung22,23,24,25, heart/cardiovascular system21,26,27,28,29,30,31, whole blood32, immune system17,22,23,24,25, brain20,33,34 including cerebellum18,19, eye35, gastrointestinal tract36, and musculoskeletal system18,19,37,38 (Table 1). In the present study, we chose to carry out RNA-seq analysis on seven tissues: GD13.5 yolk sac and placenta, plus GD16.5 fetal liver, kidney, lung, heart and cerebellum.
Results
Phenotype
Compared with Slc39a8(+/+) wild-type, the Slc39a8(neo/neo) fetuses and newborns were remarkably abnormal with dysmorphogenesis and severe anemia (Fig. 1); even placenta and yolk sac showed gross anemia13.
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 Zn2+-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 organ39.
Many other classes of genes that are seen (Supplemental Table S1) — includes those encoding: Zn2+-finger transcription factors (TFs) and homeobox and other very-early-embryonic functions associated with cell cycle and cell division; other Zn2+-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).
Among the statistically significant 45 genes being down-regulated, significantly enriched GO categories included: [a] “oxygen carrier activity” (GO:0005344) with Hbb-b1, Hbb-b2, and Hbb-y; [b] “homeostatic processes” (GO:0042592) with Ahsg, Apob, Bpgm, Cps1, Dcn, Gpr65, Hba-a1, Hba-a2, Hbb-b1, Hbb-b2, Otc, Slc25a37, Slc39a8, Slc4a1, Slc8a1, and Snca; and [c] “hemoglobin complexes” (GO:0005833) with Hba-a1, Hba-a2, Hbb-b1, Hbb-b2, Hbb-bs, Hbb-y. The down-regulated genes in GO categories GO:0005344 and GO:0005833— plus the four hemoglobin genes included in GO:0042592 category — provide further evidence indicating that globally diminished ZIP8-mediated function results in dysregulated essential hematopoiesis, leading to hypoxia.
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 regulators40. 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 numbers. It is noteworthy that – except for Arnt and Hif1a, associated with the downstream hypoxic response [reviewed in41] — the remaining seven are Zn2+-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 factors42. With regard to Znf423, Klf13 and Sp1 — there are about 515 Znf genes, 18 Klf genes, and nine Sp genes in mammalian genomes [https://www.genenames.org/]. With regard to Gata1, Gata2 and Gata4 — GATA factors (vide infra) are central to hematopoietic stem cell fate [reviewed in43].
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 Zn2+-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 proteins44. With regard to Cebpe, Bhlha15 and Rxrg — CEBPE is a leucine-zipper TF associated with several myelogenous leukemias45, BHLHA15 is a member of the basic-helix-loop-helix (bHLH) family of TFs involved in numerous developmental programs46, and RXRG is a member of the nuclear receptor superfamily that also participates in many morphogenesis processes47.
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 lymphomas48; the PRDM family contains 15 members of epigenetic regulators that participate in autoimmune diseases and infections49.
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 malignancies50.
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 Zn2+, Mn2+ and Fe2+ — functions of these divalent cations are briefly reviewed.
Zinc
Intracellular Zn2+ is critical in homeostasis-related signal transduction51, cell cycle and proliferation52, numerous processes that occur during development and differentiation, and maintenance of many Zn2+-requiring functions. In mammals there are >100 Zn2+-dependent enzymes53 and >2,000 Zn2+-containing TFs54. An independent analysis estimated there are ~2,800 Zn2+-binding proteins, corresponding to ~10% of the human proteome55. 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 Zn2+ homeostasis is essential from the single-cell-zygote stage onward56. 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 Zn2+ transport.
Manganese
Mn2+ is also a potent substrate for ZIP87,16. There are several known defects in ZIP-mediated Mn2+ transport [reviewed in57]. A human SLC39A8 nonsynonymous variant, resulting in decreased ZIP8 transport, impairs the function of Mn2+-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 loss19. Because numerous proteins are posttranslationally glycosylated, many genes listed in Supplemental Table S1 might be differentially-expressed as a reflection of defective ZIP8-mediated Mn2+ transport. There are no Mn2+-containing TFs.
Iron
ZIP8, along with its closest evolutionarily-related ZIP14 transporter, is a Fe2+-transport protein5,58. These data might be consistent with studies in the intact animal13 — 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 Fe2+ transport more so than Zn2+ transport, will require further studies. Although there are a few Fe2+-containing TFs in prokaryotes59, 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 Zn2+ transport might be more significant in lung and kidney, those affected by diminished ZIP8-mediated Mn2+ transport might be more prominent in liver, and those affected by deficient ZIP8-mediated Fe2+ transport might be more important in hematopoietic tissues. Two, it is possible that defective Zn2+ transport, primarily in early-embryo yolk sac, might override (i.e. occur upstream of) Mn2+-dependent and/or Fe2+-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 Zn2+-finger TFs, primarily in yolk sac, cause irreparable impairment of hematopoietic stem cell fate, which is then manifested — developmentally 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 Zn2+- and/or Mn2+-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 development60. GATA TFs are Zn2+-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 targets62, 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 factors63. 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 Zn2+-finger TF that (like TAL1) also forms heterodimers with TFs of the GATA family.
Moreover, TAL1 is known to activate Egln3, Hemgn and Ank162, and all were up-regulated in Slc39a8(neo/neo) yolk sac. TAL1 represses Hbb-b262, 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 Zn2+-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 Zn2+-finger genes was within normal limits, only the protein functions of these Zn2+-finger TFs are the result of ZIP8 deficiency; in other words, postttranslational incorporation of Zn2+, and/or deficient posttranslational glycosylation due to diminished Mn2+ levels, did not occur normally.
Fundamental Importance of GATA
The hematopoietically-expressed GATA family members of Zn2+-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 hematopoiesis64. GATA2 mutations are clinically associated with familial myelodysplastic syndrome65,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 genes67; 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.
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 targets68 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 hypoxia69. In the presence of hypoxia, HIF1A/ARNT dimers are known to bind to an aryl hydrocarbon receptor response element (AHRE) upstream of the Cyp2s1 gene70; 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 Fe2+ transport71; 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 hypoxia72.
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 (pO2) 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 genomes41. 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 differentially-expressed 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.
Material and Methods
Animals
Cloning details for generating the viable Slc39a8(+/neo) heterozygous mouse line were described previously12. At 6–10 weeks of age, Slc39a8(+/neo) heterozygous males and females were bred. The morning on which the vaginal plug was found – was considered GD0.5. Individual Slc39a8(+/+), Slc39a8(+/neo) and Slc39a8(neo/neo) yolk sacs and placentas at GD13.5 were collected and genotyped; heterozygotes were discarded. Fetal liver, kidney, lung, heart and cerebellum at GD16.5 were also collected and genotyped.
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.
Bioinformatics analysis
Sequence reads were aligned to the genome by using the standard Illumina sequence analysis pipeline, which was analyzed by the Statistical Genomics and Systems Biology Core at the University of Cincinnati. Analyses included determinations of: [a] RNA-seq data quality control (QC) results; [b] all gene expression levels for the RNA samples; and [c] significantly differentially-expressed genes between the two groups, Slc39a8(neo/neo) vs Slc39a8(+/+). Genes with counts per million (CPM) >1 in at least three samples were used for further analysis. Differential expression analysis was performed using DESeq.73. The Benjamini–Hochberg procedure74 was applied to adjust P-values for False Discovery Rate (FDR). Genes with FDR-adjusted P-values of <0.1 were considered as significantly differentially-expressed.
In addition to differential-expression analysis for individual tissues, meta-analysis was performed to identify genes that tend to change in the same direction across all tissues. The meta-analysis used Fisher’s combined probability test to integrate empirical P-values from all tissues. The empirical P-value for gene i in tissue j is defined as rij = rank(sij)/N, and sij = sign(fij) × [−log10(pij)], where fij is the log2 fold-change of gene i in tissue j, pij the P-value of differential expression of gene i in tissue j, and N the total number of genes. See footnote of Supplemental Table S2 for a detailed description. The empirical P-value, based on gene rank, was less prone to an extreme effect in any individual tissue; thus, it was more robust, when compared with the original P-value from differential analysis. FDR-adjusted P-values of <0.1 were used as the cut-off for genes that tend to be differentially-expressed in all tissues.
Functional analysis of differentially-expressed genes was further performed. Because too few genes were significant at FDR < 0.1, a less-stringent cut-off P-value of <0.05 and absolute fold-change of >2 was used, in order that a sufficient number of differentially-expressed genes could be used for statistical functional enrichment analysis. This would potentially result in a gene list having an increased number of false positives. Therefore, the purpose of the relaxed cut-off was to explore new hypotheses, rather than to derive individual candidate genes for validations. Specifically, MouseMine75 was used to identify significantly enriched gene ontology (GO) categories. Signaling-pathway impact analysis (SPIA) methodology76 – with number of differentially-expressed (NDE) genes >1; and FDR-corrected global P-value for each gene (pGFdr) <0.2 as cutoff] was used to identify significantly impacted Kyoto Encyclopedia of Genes and Genomes (KEGG)-signaling pathways. Ingenuity Pathway Analysis (IPA) was used to create the summary report and network figure.
In order to search for upstream regulators causing the observed changes in transcriptomes, the transcription-factor-binding site (TFBS) enrichment analysis was performed for promoters of differentially-expressed genes, as well as the significant genes from meta-analysis using Pscan77. For each TFBS, Pscan calculates a matching score for each promoter and compares the average matching score in a gene list vs those throughout the rest of the genome. If the average score is significantly higher in that gene list than the rest of the genome, then Pscan flags this TFBS as significant for this gene list. In this analysis, TF binding profiles in the JASPAR 2016 CORE database [http://jaspar.genered.net] was used, and promoter regions from −450 to +50 were considered. Significant TFBS genes were selected by Bonferroni corrected P-values of <0.1, as recommended by PScan.
All bioinformatics analyses, except for the tools mentioned above, were performed using the statistical computing platform R, and heat maps were plotted using the R package ComplexHeatmap78. We have uploaded our Slc39a8 RNA-seq data onto Gene Expression Omnibus (GEO). The link to our dataset is https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE111080.
References
Dalton, T. P. et al. Refining the mouse chromosomal location of Cdm, the major locus associated with susceptibility to cadmium-induced testicular necrosis. Pharmacogenetics 10, 141–151 (2000).
Dalton, T. P. et al. Identification of mouse SLC39A8 as the transporter responsible for cadmium-induced toxicity in the testis. Proc. Natl. Acad. Sci. USA 102, 3401–3406 (2005).
Wang, B. et al. Enhanced cadmium-induced testicular necrosis and renal proximal tubule damage caused by gene-dose increase in a Slc39a8-transgenic mouse line. Am J Physiol Cell Physiol 292, C1523–C1535 (2007).
Liu, Z. et al. Cd2+ versus Zn2+ uptake by the ZIP8 HCO3 −dependent symporter: kinetics, electrogenicity and trafficking. Biochem. Biophys. Res. Commun. 365, 814–820 (2008).
Wang, C. Y. et al. ZIP8 is an iron and zinc transporter whose cell-surface expression is up-regulated by cellular iron loading. J Biol Chem 287, 34032–34043 (2012).
McDermott, J. R. et al. Zinc- and bicarbonate-dependent ZIP8 transporter mediates selenite uptake. Oncotarget. 7, 35327–35340 (2016).
He, L. et al. ZIP8, member of the solute-carrier-39 (SLC39) metal-transporter family: characterization of transporter properties. Mol Pharmacol 70, 171–180 (2006).
Nebert, D. W. et al. ZIP14 and ZIP8 zinc/bicarbonate symporters in Xenopus oocytes: characterization of metal uptake and inhibition. Metallomics. 4, 1218–1225 (2012).
Moore-Scott, B. A., Opoka, R., Lin, S. C., Kordich, J. J. & Wells, J. M. Identification of molecular markers that are expressed in discrete anterior-posterior domains of the endoderm from the gastrula stage to mid-gestation. Dev. Dyn. 236, 1997–2003 (2007).
Harrison, S. M., Dunwoodie, S. L., Arkell, R. M., Lehrach, H. & Beddington, R. S. Isolation of novel tissue-specific genes from cDNA libraries representing the individual tissue constituents of the gastrulating mouse embryo. Development 121, 2479–2489 (1995).
Zhu, H., Yang, H. & Owen, M. R. Combined microarray analysis uncovers self-renewal-related signaling in mouse embryonic stem cells. Syst. Synth. Biol. 1, 171–181 (2007).
Wang, B., He, L., Dong, H., Dalton, T. P. & Nebert, D. W. Generation of a Slc39a8 hypomorph mouse: markedly decreased ZIP8 Zn2+/(HCO3 −)2 transporter expression. Biochem. Biophys. Res. Commun. 410, 289–294 (2011).
Gálvez-Peralta, M. et al. ZIP8 zinc transporter: indispensable role for both multiple-organ organogenesis and hematopoiesis in utero. PLoS ONE. 7, e36055 (2012).
Lin, W. et al. Hepatic metal ion transporter ZIP8 regulates manganese homeostasis and manganese-dependent enzyme activity. J Clin. Invest 127, 2407–2417 (2017).
Lin, W. et al. Zinc transporter Slc39a8 is essential for cardiac ventricular compaction. J Clin. Invest. 128, 826–833 (2018).
He, L., Wang, B., Hay, E. B. & Nebert, D. W. Discovery of ZIP transporters that participate in cadmium damage to testis and kidney. Toxicol. Appl. Pharmacol. 238, 250–257 (2009).
Begum, N. A. et al. Mycobacterium bovis BCG cell wall and lipopolysaccharide induce a novel gene, BIGM103, encoding a 7-TM protein: identification of a new protein family having Zn-transporter and Zn-metalloprotease signatures. Genomics 80, 630–645 (2002).
Boycott, K. M. et al. Autosomal-recessive intellectual disability with cerebellar atrophy syndrome caused by mutation of the manganese- and zinc-transporter gene SLC39A8. Am J Hum Genet 97, 886–893 (2015).
Park, J. H. et al. SLC39A8 deficiency: a disorder of manganese transport and glycosylation. Am J Hum Genet 97, 894–903 (2015).
Riley, L. G. et al. A SLC39A8 variant causes manganese deficiency, and glycosylation and mitochondrial disorders. J Inherit. Metab Dis. 40, 261–269 (2017).
Waterworth, D. M. et al. Genetic variants influencing circulating lipid levels and risk of coronary artery disease. Arterioscler. Thromb. Vasc. Biol. 30, 2264–2276 (2010).
Besecker, B. et al. The human zinc transporter SLC39A8 (ZIP8) is critical in zinc-mediated cytoprotection in lung epithelia. Am J Physiol Lung Cell Mol. Physiol 294, L1127–L1136 (2008).
Liu, M. J. et al. ZIP8 regulates host defense through zinc-mediated inhibition of NF-κB. Cell Rep. 3, 386–400 (2013).
Pyle, C. J. et al. Zinc modulates endotoxin-induced human macrophage inflammation through ZIP8 Induction and C/EBPβ Inhibition. PLoS ONE 12, e0169531 (2017).
Pyle,C.J. et al. Elemental ingredients in the macrophage cocktail: role of ZIP8 in host response to Mycobacterium tuberculosis. Int J Mol. Sci. 18 (2017).
Johansson, A. et al. Genome-wide association and Mendelian randomization study of NT-proBNP in patients with acute coronary syndrome. Hum Mol. Genet 25, 1447–1456 (2016).
Verdugo, R. A. et al. Graphical modeling of gene expression in monocytes suggests molecular mechanisms explaining increased atherosclerosis in smokers. PLoS One 8, e50888 (2013).
Ehret, G. B. et al. Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk. Nature 478, 103–109 (2011).
Ehret, G. B. et al. The genetics of blood pressure regulation and its target organs from association studies in 342, 415 individuals. Nat. Genet 48, 1171–1184 (2016).
Teslovich, T. M. et al. Biological, clinical and population relevance of 95 loci for blood lipids. Nature 466, 707–713 (2010).
Willer, C. J. et al. Discovery and refinement of loci associated with lipid levels. Nat. Genet 45, 1274–1283 (2013).
Ng, E. et al. Genome-wide association study of toxic metals and trace elements reveals novel associations. Hum Mol. Genet 24, 4739–4745 (2015).
Carrera, N. et al. Association study of nonsynonymous single nucleotide polymorphisms in schizophrenia. Biol Psychiatry 71, 169–177 (2012).
Schizophrenia Working Group of the Psychiatric Genomics Consortium Biological insights from 108 schizophrenia-associated genetic loci. Nature 511, 421-427 (2014).
Baumann, B. et al. Conditional Müller cell ablation leads to retinal iron accumulation. Invest Ophthalmol. Vis. Sci. 58, 4223–4234 (2017).
Li, D. et al. A pleiotropic missense variant in SLC39A8 is associated with Crohn disease and human gut microbiome composition. Gastroenterology 151, 724–732 (2016).
Song, J. et al. MicroRNA-488 regulates zinc transporter SLC39A8/ZIP8 during pathogenesis of osteoarthritis. J Biomed. Sci. 20, 31 (2013).
Kim, J. H. et al. Regulation of the catabolic cascade in osteoarthritis by the zinc-ZIP8-MTF1 axis. Cell 156, 730–743 (2014).
Alvarez-Silva, M., Belo-Diabangouaya, P., Salaun, J. & Dieterlen-Lievre, F. Mouse placenta is a major hematopoietic organ. Development 130, 5437–5444 (2003).
Guo, J. et al. Chromatin and single-cell RNA-seq profiling reveal dynamic signaling and metabolic transitions during human spermatogonial stem cell development. Cell Stem Cell 21, 533–546 (2017).
Nebert, D. W. Aryl hydrocarbon receptor (AHR): “pioneer member” of the basic-helix/loop/helix per-Arnt-sim (bHLH/PAS) family of “sensors” of foreign and endogenous signals. Progr. Lipid Res. 67, 38–57 (2017).
Juma, A. R., Damdimopoulou, P. E., Grommen, S. V., Van de Ven, W. J. & De, G. B. Emerging role of PLAG1 as a regulator of growth and reproduction. J Endocrinol. 228, R45–R56 (2016).
Fujiwara, T. GATA transcription factors: basic principles and related human disorders. Tohoku J Exp. Med. 242, 83–91 (2017).
Bürglin, T. R. & Affolter, M. Homeodomain proteins: an update. Chromosoma 125, 497–521 (2016).
Gutierrez-Camino, A., Martin-Guerrero, I. & Garcia-Orad, A. Genetic susceptibility in childhood acute lymphoblastic leukemia. Med Oncol. 34, 179 (2017).
Gyoja, F. Basic helix-loop-helix transcription factors in evolution: roles in development of mesoderm and neural tissues. Genesis. 55 https://doi.org/10.1002/dvg.23051 (2017).
Musser, J. M. & Arendt, D. Loss and gain of cone types in vertebrate ciliary photoreceptor evolution. Dev. Biol 431, 26–35 (2017).
Boi, M., Zucca, E., Inghirami, G. & Bertoni, F. PRDM1: a tumor suppressor gene in B-cell and T-cell lymphomas. Leuk. Lymphoma 56, 1223–1228 (2015).
Mzoughi, S., Tan, Y. X., Low, D. & Guccione, E. The role of PRDMs in cancer: one family, two sides. Curr. Opin. Genet Dev. 36, 83–91 (2016).
Hu, Y., Yoshida, T. & Georgopoulos, K. Transcriptional circuits in B-cell transformation. Curr. Opin. Hematol. 24, 345–352 (2017).
Hirano, T. et al. Roles of zinc and zinc signaling in immunity: zinc as an intracellular signaling molecule. Advanc. Immunol. 97, 149–176 (2008).
Li, Y. & Maret, W. Transient fluctuations of intracellular zinc ions in cell proliferation. Exp. Cell Res. 315, 2463–2470 (2009).
Food and Nutrition Board. Zinc: Dietary reference intakes for vitamin A, vitamin K, arsenic, boron, chromium, copper, iodine, iron, manganese, molybdenum, nickel, silicon, vanadium, and zinc. Institute of Medicine 1, 442–501 (2001).
Prasad, A. S. Zinc deficiency. BMJ 326, 409–410 (2003).
Andreini, C., Banci, L., Bertini, I. & Rosato, A. Counting the zinc-proteins encoded in the human genome. J Proteome. Res. 5, 196–201 (2006).
Kambe, T., Weaver, B. P. & Andrews, G. K. The genetics of essential metal homeostasis during development. Genesis. 46, 214–228 (2008).
Zogzas, C. E. & Mukhopadhyay, S. Inherited disorders of manganese metabolism. Advanc. Neurobiol. 18, 35–49 (2017).
Sterling, J. et al. Iron importers ZIP8 and ZIP14 are expressed in retina and regulated by retinal iron levels. Exp. Eye Res. 155, 15–23 (2017).
Fleischhacker, A. S. & Kiley, P. J. Iron-containing transcription factors and their roles as sensors. Curr. Opin. Chem Biol 15, 335–341 (2011).
Tan, S. H., Bertulfo, F. C. & Sanda, T. Leukemia-initiating cells in T-cell acute lymphoblastic leukemia. Front Oncol. 7, 218 (2017).
Wu, W. et al. Dynamic shifts in occupancy by TAL1 are guided by GATA factors and drive large-scale reprogramming of gene expression during hematopoiesis. Genome Res. 24, 1945–1962 (2014).
Kassouf, M. T. et al. Genome-wide identification of TAL1’s functional targets: insights into its mechanisms of action in primary erythroid cells. Genome Res. 20, 1064–1083 (2010).
Sciorra, V. A., Sanchez, M. A., Kunibe, A. & Wurmser, A. E. Suppression of glioma progression by EGLN3. PLoS One 7, e40053 (2012).
Ferreira, R. et al. Dynamic regulation of GATA factor levels is more important than their identity. Blood 109, 5481–5490 (2007).
Wlodarski, M. W. et al. Prevalence, clinical characteristics, and prognosis of GATA2-related myelodysplastic syndromes in children and adolescents. Blood 127, 1387–1397 (2016).
Hirabayashi, S., Wlodarski, M. W., Kozyra, E. & Niemeyer, C. M. Heterogeneity of GATA2-related myeloid neoplasms. Int J Hematol. 106, 175–182 (2017).
Papadopoulos, G. L. et al. GATA1 genome-wide occupancy associates with distinct epigenetic profiles in mouse fetal liver erythropoiesis. Nucleic Acids Res. 41, 4938–4948 (2013).
Benita, Y. et al. An integrative genomics approach identifies hypoxia-inducible factor-1 (HIF1)-target genes that form the core response to hypoxia. Nucleic Acids Res. 37, 4587–4602 (2009).
Dopp, J. M. et al. Influence of intermittent hypoxia on myocardial and hepatic P-glycoprotein (ABCB1A) expression in a rodent model. Pharmacotherapy 29, 365–372 (2009).
Rivera, S. P. et al. A novel promoter element containing multiple overlapping xenobiotic and hypoxia response elements mediates induction of cytochrome P450 2S1 by both dioxin and hypoxia. J Biol Chem 282, 10881–10893 (2007).
Frazer, D. M. & Anderson, G. J. The regulation of iron transport. Biofactors 40, 206–214 (2014).
Krock, B. L., Skuli, N. & Simon, M. C. Hypoxia-induced angiogenesis: good and evil. Genes Cancer 2, 1117–1133 (2011).
Anders, S. & Huber, W. Differential expression analysis for sequence count data. Genome Biol 11, R106 (2010).
Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Statist. Soc. Series B 57, 289–300 (1995).
Motenko, H., Neuhauser, S. B., O’Keefe, M. & Richardson, J. E. MouseMine: a new data warehouse for MGI. Mamm. Genome 26, 325–330 (2015).
Tarca, A. L. et al. A novel signaling-pathway impact analysis. Bioinformatics. 25, 75–82 (2009).
Zambelli, F., Pesole, G. & Pavesi, G. Pscan: Finding over-represented transcription-factor binding-site motifs in sequences from co-regulated or co-expressed genes. Nucleic Acids Res. 37, W247–W252 (2009).
Gu, Z., Eils, R. & Schlesner, M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics. 32, 2847–2849 (2016).
Speliotes, E. K. et al. Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index. Nat. Genet 42, 937–948 (2010).
Dickinson, M. E. et al. High-throughput discovery of novel developmental phenotypes. Nature 537, 508–514 (2016).
Acknowledgements
We appreciate the technical help of Saikumar Karyala, and the IPA assay account generously provided by Alvaro Puga. We thank our colleagues —especially S. Steven Potter — for valuable discussions during the course of this work and careful reading of this manuscript. This project was funded, in part, by out-of-pocket personal funds during the final 12 months (DWN) and a Center for Environmental Genetics (CEG) Trainee-Mentor award, 2011–12 (MG-P). Generation and characterization of the knockdown mouse line was funded earlier by the National Institute of Environmental Health Sciences of the National Institutes of Health [grant numbers R01 ES010416 (DWN), T32 ES016646 (MG-P), and P30 ES06096 (DWN)].
Author information
Authors and Affiliations
Contributions
D.W.N., M.G.-P. and X.Z. conceived and designed the study. M.G.-P. and J.D. performed most of the experiments. D.W.N. and J.C. analyzed the data and wrote the manuscript, with the help of X.Z., Z.L. and M.G.-P.J.C. and X.Z. provided the data analysis. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Chen, J., Gálvez-Peralta, M., Zhang, X. et al. In utero gene expression in the Slc39a8(neo/neo) knockdown mouse. Sci Rep 8, 10703 (2018). https://doi.org/10.1038/s41598-018-29109-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-018-29109-y
This article is cited by
-
Back-translating GWAS findings to animal models reveals a role for Hgfac and Slc39a8 in alcohol and nicotine consumption
Scientific Reports (2022)
-
The schizophrenia risk locus in SLC39A8 alters brain metal transport and plasma glycosylation
Scientific Reports (2020)
-
Hyperaccumulation of cadmium by scallop Chlamys farreri revealed by comparative transcriptome analysis
BioMetals (2020)
-
SLC39A8 gene encoding a metal ion transporter: discovery and bench to bedside
Human Genomics (2019)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.