Species differences in specific ligand-binding affinity and activation of AHR: The biological basis for calculation of relative effective potencies and toxic equivalence factors.

In 2022 the World Health Organization (WHO) published updated 'Toxic Equivalence Factors' (TEFs) for a wide variety of chlorinated dioxins, dibenzofurans and PCBs [collectively referred to as 'dioxin-like chemicals'; DLCs) that interact with the aryl hydrocarbon receptor (AHR)]. Their update used sophisticated statistical analysis of hundreds of published studies that reported estimation of 'Relative Effective Potency' (REP) values for individual DLC congeners. The weighting scheme used in their assessment of each study favored in vivo over in vitro studies and was based largely on rodent studies. In this Commentary, we highlight the large body of published studies that demonstrate large species differences in AHR-ligand activation and provide supporting evidence for our position that the WHO 2022 TEF values intended for use in human risk assessment of DLC mixtures will provide highly misleading overestimates of 'Toxic Equivalent Quotients' (TEQs), because of well-recognized striking differences in AHR ligand affinities between rodent (rat, mouse) and human. The data reviewed in our Commentary support the position that human tissue-derived estimates of REP/TEF values for individual DLC congeners, although uncertain, will provide much better, more realistic estimates of potential activation of the human AHR, when exposure to complex DLC mixtures occurs.

Gene-Environment Interactions: My Unique Journey.

I am deeply honored to be invited to write this scientific autobiography. As a physician-scientist, pediatrician, molecular biologist, and geneticist, I have authored/coauthored more than 600 publications in the fields of clinical medicine, biochemistry, biophysics, pharmacology, drug metabolism, toxicology, molecular biology, cancer, standardized gene nomenclature, developmental toxicology and teratogenesis, mouse genetics, human genetics, and evolutionary genomics. Looking back, I think my career can be divided into four distinct research areas, which I summarize mostly chronologically in this article: (a) discovery and characterization of the AHR/CYP1 axis, (b) pharmacogenomics and genetic prediction of response to drugs and other environmental toxicants, (c) standardized drug-metabolizing gene nomenclature based on evolutionary divergence, and (d) discovery and characterization of the SLC39A8 gene encoding the ZIP8 metal cation influx transporter. Collectively, all four topics embrace gene-environment interactions, hence the title of my autobiography.

Endocrine pancreas-specific Gclc gene deletion causes a severe diabetes phenotype.

Reduced glutathione (GSH) is an abundant antioxidant that regulates intracellular redox homeostasis by scavenging reactive oxygen species (ROS). Glutamate-cysteine ligase catalytic (GCLC) subunit is the rate-limiting step in GSH biosynthesis. Using the Pax6-Cre driver mouse line, we deleted expression of the Gclc gene in all pancreatic endocrine progenitor cells. Intriguingly, Gclc knockout (KO) mice, following weaning, exhibited an age-related, progressive diabetes phenotype, manifested as strikingly increased blood glucose and decreased plasma insulin levels. This severe diabetes trait is preceded by pathologic changes in islet of weanling mice. Gclc KO weanlings showed progressive abnormalities in pancreatic morphology including: islet-specific cellular vacuolization, decreased islet-cell mass, and alterations in islet hormone expression. Islets from newly-weaned mice displayed impaired glucose-stimulated insulin secretion, decreased insulin hormone gene expression, oxidative stress, and increased markers of cellular senescence. Our results suggest that GSH biosynthesis is essential for normal development of the mouse pancreatic islet, and that protection from oxidative stress-induced cellular senescence might prevent abnormal islet-cell damage during embryogenesis.

Fatal Epileptic Seizures in Mice Having Compromised Glutathione and Ascorbic Acid Biosynthesis.

Reduced glutathione (GSH) and ascorbic acid (AA) are the two most abundant low-molecular-weight antioxidants in mammalian tissues. GclmKO knockout mice lack the gene encoding the modifier subunit of the rate-limiting enzyme in GSH biosynthesis; GclmKO mice exhibit 10-40% of normal tissue GSH levels and show no overt phenotype. GuloKO knockout mice, lacking a functional Gulo gene encoding L-gulono-γ-lactone oxidase, cannot synthesize AA and depend on dietary ascorbic acid for survival. To elucidate functional crosstalk between GSH and AA in vivo, we generated the GclmKO/GuloKO double-knockout (DKO) mouse. DKO mice exhibited spontaneous epileptic seizures, proceeding to death between postnatal day (PND)14 and PND23. Histologically, DKO mice displayed neuronal loss and glial proliferation in the neocortex and hippocampus. Epileptic seizures and brain pathology in young DKO mice could be prevented with AA supplementation in drinking water (1 g/L). Remarkably, in AA-rescued adult DKO mice, the removal of AA supplementation for 2-3 weeks resulted in similar, but more severe, neocortex and hippocampal pathology and seizures, with death occurring between 12 and 21 days later. These results provide direct evidence for an indispensable, yet underappreciated, role for the interplay between GSH and AA in normal brain function and neuronal health. We speculate that the functional crosstalk between GSH and AA plays an important role in regulating glutamatergic neurotransmission and in protecting against excitotoxicity-induced brain damage.

Placing human gene families into their evolutionary context.

Following the draft sequence of the first human genome over 20 years ago, we have achieved unprecedented insights into the rules governing its evolution, often with direct translational relevance to specific diseases. However, staggering sequence complexity has also challenged the development of a more comprehensive understanding of human genome biology. In this context, interspecific genomic studies between humans and other animals have played a critical role in our efforts to decode human gene families. In this review, we focus on how the rapid surge of genome sequencing of both model and non-model organisms now provides a broader comparative framework poised to empower novel discoveries. We begin with a general overview of how comparative approaches are essential for understanding gene family evolution in the human genome, followed by a discussion of analyses of gene expression. We show how homology can provide insights into the genes and gene families associated with immune response, cancer biology, vision, chemosensation, and metabolism, by revealing similarity in processes among distant species. We then explain methodological tools that provide critical advances and show the limitations of common approaches. We conclude with a discussion of how these investigations position us to gain fundamental insights into the evolution of gene families among living organisms in general. We hope that our review catalyzes additional excitement and research on the emerging field of comparative genomics, while aiding the placement of the human genome into its existentially evolutionary context.

Update of the keratin gene family: evolution, tissue-specific expression patterns, and relevance to clinical disorders.

Intermediate filament (IntFil) genes arose during early metazoan evolution, to provide mechanical support for plasma membranes contacting/interacting with other cells and the extracellular matrix. Keratin genes comprise the largest subset of IntFil genes. Whereas the first keratin gene appeared in sponge, and three genes in arthropods, more rapid increases in keratin genes occurred in lungfish and amphibian genomes, concomitant with land animal-sea animal divergence (~ 440 to 410 million years ago). Human, mouse and zebrafish genomes contain 18, 17 and 24 non-keratin IntFil genes, respectively. Human has 27 of 28 type I "acidic" keratin genes clustered at chromosome (Chr) 17q21.2, and all 26 type II "basic" keratin genes clustered at Chr 12q13.13. Mouse has 27 of 28 type I keratin genes clustered on Chr 11, and all 26 type II clustered on Chr 15. Zebrafish has 18 type I keratin genes scattered on five chromosomes, and 3 type II keratin genes on two chromosomes. Types I and II keratin clusters-reflecting evolutionary blooms of keratin genes along one chromosomal segment-are found in all land animal genomes examined, but not fishes; such rapid gene expansions likely reflect sudden requirements for many novel paralogous proteins having divergent functions to enhance species survival following sea-to-land transition. Using data from the Genotype-Tissue Expression (GTEx) project, tissue-specific keratin expression throughout the human body was reconstructed. Clustering of gene expression patterns revealed similarities in tissue-specific expression patterns for previously described "keratin pairs" (i.e., KRT1/KRT10, KRT8/KRT18, KRT5/KRT14, KRT6/KRT16 and KRT6/KRT17 proteins). The ClinVar database currently lists 26 human disease-causing variants within the various domains of keratin proteins.

Overview of PAX gene family: analysis of human tissue-specific variant expression and involvement in human disease.

Paired-box (PAX) genes encode a family of highly conserved transcription factors found in vertebrates and invertebrates. PAX proteins are defined by the presence of a paired domain that is evolutionarily conserved across phylogenies. Inclusion of a homeodomain and/or an octapeptide linker subdivides PAX proteins into four groups. Often termed "master regulators", PAX proteins orchestrate tissue and organ development throughout cell differentiation and lineage determination, and are essential for tissue structure and function through maintenance of cell identity. Mutations in PAX genes are associated with myriad human diseases (e.g., microphthalmia, anophthalmia, coloboma, hypothyroidism, acute lymphoblastic leukemia). Transcriptional regulation by PAX proteins is, in part, modulated by expression of alternatively spliced transcripts. Herein, we provide a genomics update on the nine human PAX family members and PAX homologs in 16 additional species. We also present a comprehensive summary of human tissue-specific PAX transcript variant expression and describe potential functional significance of PAX isoforms. While the functional roles of PAX proteins in developmental diseases and cancer are well characterized, much remains to be understood regarding the functional roles of PAX isoforms in human health. We anticipate the analysis of tissue-specific PAX transcript variant expression presented herein can serve as a starting point for such research endeavors.

COVID-19 vulnerability: the potential impact of genetic susceptibility and airborne transmission.

The recent coronavirus disease (COVID-19), caused by SARS-CoV-2, is inarguably the most challenging coronavirus outbreak relative to the previous outbreaks involving SARS-CoV and MERS-CoV. With the number of COVID-19 cases now exceeding 2 million worldwide, it is apparent that (i) transmission of SARS-CoV-2 is very high and (ii) there are large variations in disease severity, one component of which may be genetic variability in the response to the virus. Controlling current rates of infection and combating future waves require a better understanding of the routes of exposure to SARS-CoV-2 and the underlying genomic susceptibility to this disease. In this mini-review, we highlight possible genetic determinants of COVID-19 and the contribution of aerosol exposure as a potentially important transmission route of SARS-CoV-2.

Genetics and functions of the retinoic acid pathway, with special emphasis on the eye.

Retinoic acid (RA) is a potent morphogen required for embryonic development. RA is formed in a multistep process from vitamin A (retinol); RA acts in a paracrine fashion to shape the developing eye and is essential for normal optic vesicle and anterior segment formation. Perturbation in RA-signaling can result in severe ocular developmental diseases-including microphthalmia, anophthalmia, and coloboma. RA-signaling is also essential for embryonic development and life, as indicated by the significant consequences of mutations in genes involved in RA-signaling. The requirement of RA-signaling for normal development is further supported by the manifestation of severe pathologies in animal models of RA deficiency-such as ventral lens rotation, failure of optic cup formation, and embryonic and postnatal lethality. In this review, we summarize RA-signaling, recent advances in our understanding of this pathway in eye development, and the requirement of RA-signaling for embryonic development (e.g., organogenesis and limb bud development) and life.

SLC39A8 gene encoding a metal ion transporter: discovery and bench to bedside.

SLC39A8 is an evolutionarily highly conserved gene that encodes the ZIP8 metal cation transporter in all vertebrates. SLC39A8 is ubiquitously expressed, including pluripotent embryonic stem cells; SLC39A8 expression occurs in every cell type examined. Uptake of ZIP8-mediated Mn2+, Zn2+, Fe2+, Se4+, and Co2+ represents endogenous functions-moving these cations into the cell. By way of mouse genetic differences, the phenotype of "subcutaneous cadmium-induced testicular necrosis" was assigned to the Cdm locus in the 1970s. This led to identification of the mouse Slc39a8 gene, its most closely related Slc39a14 gene, and creation of Slc39a8-overexpressing, Slc39a8(neo/neo) knockdown, and cell type-specific conditional knockout mouse lines; the Slc39a8(-/-) global knockout mouse is early-embryolethal. Slc39a8(neo/neo) hypomorphs die between gestational day 16.5 and postnatal day 1-exhibiting severe anemia, dysregulated hematopoiesis, hypoplastic spleen, dysorganogenesis, stunted growth, and hypomorphic limbs. Not surprisingly, genome-wide association studies subsequently revealed human SLC39A8-deficiency variants exhibiting striking pleiotropy-defects correlated with clinical disorders in virtually every organ, tissue, and cell-type: numerous developmental and congenital disorders, the immune system, cardiovascular system, kidney, lung, liver, coagulation system, central nervous system, musculoskeletal system, eye, and gastrointestinal tract. Traits with which SLC39A8-deficiency variants are currently associated include Mn2+-deficient hypoglycosylation; numerous birth defects; Leigh syndrome-like mitochondrial redox deficiency; decreased serum high-density lipoprotein-cholesterol levels; increased body mass index; greater risk of coronary artery disease, hypotension, cardiovascular death, allergy, ischemic stroke, schizophrenia, Parkinson disease, inflammatory bowel disease, Crohn disease, myopia, and adolescent idiopathic scoliosis; systemic lupus erythematosus with primary Sjögren syndrome; decreased height; and inadvertent participation in the inflammatory progression of osteoarthritis.

Update on the human and mouse lipocalin (LCN) gene family, including evidence the mouse Mup cluster is result of an "evolutionary bloom".

Lipocalins (LCNs) are members of a family of evolutionarily conserved genes present in all kingdoms of life. There are 19 LCN-like genes in the human genome, and 45 Lcn-like genes in the mouse genome, which include 22 major urinary protein (Mup) genes. The Mup genes, plus 29 of 30 Mup-ps pseudogenes, are all located together on chromosome (Chr) 4; evidence points to an "evolutionary bloom" that resulted in this Mup cluster in mouse, syntenic to the human Chr 9q32 locus at which a single MUPP pseudogene is located. LCNs play important roles in physiological processes by binding and transporting small hydrophobic molecules -such as steroid hormones, odorants, retinoids, and lipids-in plasma and other body fluids. LCNs are extensively used in clinical practice as biochemical markers. LCN-like proteins (18-40 kDa) have the characteristic eight ß-strands creating a barrel structure that houses the binding-site; LCNs are synthesized in the liver as well as various secretory tissues. In rodents, MUPs are involved in communication of information in urine-derived scent marks, serving as signatures of individual identity, or as kairomones (to elicit fear behavior). MUPs also participate in regulation of glucose and lipid metabolism via a mechanism not well understood. Although much has been learned about LCNs and MUPs in recent years, more research is necessary to allow better understanding of their physiological functions, as well as their involvement in clinical disorders.

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.

Hepatic ZIP8 deficiency is associated with disrupted selenium homeostasis, liver pathology, and tumor formation.

Zrt/Irt-like protein 8 (ZIP8) (encoded by Slc39a8) is a multifunctional membrane transporter that influxes essential metal cations Zn2+, Mn2+, Fe2+, and nonmetal inorganic selenite (HSeO3-). Physiological roles of ZIP8 in different cell types and tissues remain to be elucidated. We aimed to investigate ZIP8 functions in liver. Two mouse models were used in this study: 1) 13- to 21-mo-old Slc39a8(+/neo) hypomorphs having diminished ZIP8 levels and 2) a liver-specific ZIP8 acute knockdown mouse (Ad-shZip8). Histology, immunohistochemistry, and Western blotting were used to investigate ZIP8-deficiency effects on hepatic injury, inflammatory changes, and oxidative stress. Selenium (Se) and zinc (Zn) were quantified in tissues by inductively coupled plasma-mass spectrophotometry. We found that ZIP8 is required to maintain normal liver function; moderate or acute decreases in ZIP8 activity resulted in hepatic pathology. Spontaneous liver neoplastic nodules appeared in ~50% of Slc39a8(+/neo) between 13 and 21 mo of age, exhibiting features of inflammation, fibrosis, and liver injury. In Ad-shZip8 mice, significant hepatomegaly was observed; histology showed ZIP8 deficiency was associated with hepatocyte injury, inflammation, and proliferation. Significant decreases in Se, but not Zn, were found in Ad-shZip8 liver. Consistent with this Se deficit, liver expression of selenoproteins glutathione peroxidases 1 and 2 was downregulated, along with decreases in antioxidant superoxide dismutases 1 and 2, consistent with increased oxidative stress. Thus, ZIP8 plays an important role in maintaining normal hepatic function, likely through regulating Se homeostasis and redox balance. Hepatic ZIP8 deficiency is associated with liver pathology, including oxidative stress, inflammation, proliferation, and hepatocellular injury. NEW & NOTEWORTHY Zrt/Irt-like protein 8 (ZIP8) is a multifunctional membrane transporter that facilitates biometal and mineral uptake. The role of ZIP8 in liver physiology has not been previously investigated. Liu et al. discovered unique ZIP8 functions, i.e., regulation of hepatic selenium content and association of ZIP8 deficiency in mouse liver with liver defects.

Cytochrome P450 1A1 (CYP1A1) protects against nonalcoholic fatty liver disease caused by Western diet containing benzo[a]pyrene in mice.

The Western diet contributes to nonalcoholic fatty liver disease (NAFLD) pathogenesis. Benzo[a]pyrene (BaP), a prototypical environmental pollutant produced by combustion processes, is present in charcoal-grilled meat. Cytochrome P450 1A1 (CYP1A1) metabolizes BaP, resulting in either detoxication or metabolic activation in a context-dependent manner. To elucidate a role of CYP1A1-BaP in NAFLD pathogenesis, we compared the effects of a Western diet, with or without oral BaP treatment, on the development of NAFLD in Cyp1a1(-/-) mice versus wild-type mice. A Western diet plus BaP induced lipid-droplet accumulation in liver of Cyp1a1(-/-) mice, but not wild-type mice. The hepatic steatosis observed in Cyp1a1(-/-) mice was associated with increased cholesterol, triglyceride and bile acid levels. Cyp1a1(-/-) mice fed Western diet plus BaP had changes in expression of genes involved in bile acid and lipid metabolism, and showed no increase in Cyp1a2 expression but did exhibit enhanced Cyp1b1 mRNA expression, as well as hepatic inflammation. Enhanced BaP metabolic activation, oxidative stress and inflammation may exacerbate metabolic dysfunction in liver of Cyp1a1(-/-) mice. Thus, Western diet plus BaP induces NAFLD and hepatic inflammation in Cyp1a1(-/-) mice in comparison to wild-type mice, indicating a protective role of CYP1A1 against NAFLD pathogenesis.

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.

The basic-helix/loop/helix per-Arnt-sim (bHLH/PAS) family comprises many transcription factors, found throughout all three kingdoms of life; bHLH/PAS members "sense" innumerable intracellular and extracellular "signals" - including endogenous compounds, foreign chemicals, gas molecules, redox potential, photons (light), gravity, heat, and osmotic pressure. These signals then initiate downstream signaling pathways involved in responding to that signal. The term "PAS", abbreviation for "per-Arnt-sim" was first coined in 1991. Although the mouse Arnt gene was not identified until 1991, evidence of its co-transcriptional binding partner, aryl hydrocarbon receptor (AHR), was first reported in 1974 as a "sensor" of foreign chemicals, up-regulating cytochrome P450 family 1 (CYP1) and other enzyme activities that usually metabolize the signaling chemical. Within a few years, AHR was proposed also to participate in inflammation. The mouse [Ah] locus was shown (1973-1989) to be relevant to chemical carcinogenesis, mutagenesis, toxicity and teratogenesis, the mouse Ahr gene was cloned in 1992, and the first Ahr(-/-) knockout mouse line was reported in 1995. After thousands of studies from the early 1970s to present day, we now realize that AHR participates in dozens of signaling pathways involved in critical-life processes, affecting virtually every organ and cell-type in the animal, including many invertebrates.

Personalized medicine: Genetic risk prediction of drug response.

Pharmacogenomics (PGx), a substantial component of "personalized medicine", seeks to understand each individual's genetic composition to optimize drug therapy -- maximizing beneficial drug response, while minimizing adverse drug reactions (ADRs). Drug responses are highly variable because innumerable factors contribute to ultimate phenotypic outcomes. Recent genome-wide PGx studies have provided some insight into genetic basis of variability in drug response. These can be grouped into three categories. [a] Monogenic (Mendelian) traits include early examples mostly of inherited disorders, and some severe (idiosyncratic) ADRs typically influenced by single rare coding variants. [b] Predominantly oligogenic traits represent variation largely influenced by a small number of major pharmacokinetic or pharmacodynamic genes. [c] Complex PGx traits resemble most multifactorial quantitative traits -- influenced by numerous small-effect variants, together with epigenetic effects and environmental factors. Prediction of monogenic drug responses is relatively simple, involving detection of underlying mutations; due to rarity of these events and incomplete penetrance, however, prospective tests based on genotype will have high false-positive rates, plus pharmacoeconomics will require justification. Prediction of predominantly oligogenic traits is slowly improving. Although a substantial fraction of variation can be explained by limited numbers of large-effect genetic variants, uncertainty in successful predictions and overall cost-benefit ratios will make such tests elusive for everyday clinical use. Prediction of complex PGx traits is almost impossible in the foreseeable future. Genome-wide association studies of large cohorts will continue to discover relevant genetic variants; however, these small-effect variants, combined, explain only a small fraction of phenotypic variance -- thus having limited predictive power and clinical utility.