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.

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.

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.

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.

Feedback control of AHR signalling regulates intestinal immunity.

The aryl hydrocarbon receptor (AHR) recognizes xenobiotics as well as natural compounds such as tryptophan metabolites, dietary components and microbiota-derived factors, and it is important for maintenance of homeostasis at mucosal surfaces. AHR activation induces cytochrome P4501 (CYP1) enzymes, which oxygenate AHR ligands, leading to their metabolic clearance and detoxification. Thus, CYP1 enzymes have an important feedback role that curtails the duration of AHR signalling, but it remains unclear whether they also regulate AHR ligand availability in vivo. Here we show that dysregulated expression of Cyp1a1 in mice depletes the reservoir of natural AHR ligands, generating a quasi AHR-deficient state. Constitutive expression of Cyp1a1 throughout the body or restricted specifically to intestinal epithelial cells resulted in loss of AHR-dependent type 3 innate lymphoid cells and T helper 17 cells and increased susceptibility to enteric infection. The deleterious effects of excessive AHR ligand degradation on intestinal immune functions could be counter-balanced by increasing the intake of AHR ligands in the diet. Thus, our data indicate that intestinal epithelial cells serve as gatekeepers for the supply of AHR ligands to the host and emphasize the importance of feedback control in modulating AHR pathway activation.

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.

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.

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.

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.

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.

SLC39A8 Deficiency: A Disorder of Manganese Transport and Glycosylation.

SLC39A8 is a membrane transporter responsible for manganese uptake into the cell. Via whole-exome sequencing, we studied a child that presented with cranial asymmetry, severe infantile spasms with hypsarrhythmia, and dysproportionate dwarfism. Analysis of transferrin glycosylation revealed severe dysglycosylation corresponding to a type II congenital disorder of glycosylation (CDG) and the blood manganese levels were below the detection limit. The variants c.112G>C (p.Gly38Arg) and c.1019T>A (p.Ile340Asn) were identified in SLC39A8. A second individual with the variants c.97G>A (p.Val33Met) and c.1004G>C (p.Ser335Thr) on the paternal allele and c.610G>T (p.Gly204Cys) on the maternal allele was identified among a group of unresolved case subjects with CDG. These data demonstrate that variants in SLC39A8 impair the function of manganese-dependent enzymes, most notably ß-1,4-galactosyltransferase, a Golgi enzyme essential for biosynthesis of the carbohydrate part of glycoproteins. Impaired galactosylation leads to a severe disorder with deformed skull, severe seizures, short limbs, profound psychomotor retardation, and hearing loss. Oral galactose supplementation is a treatment option and results in complete normalization of glycosylation. SLC39A8 deficiency links a trace element deficiency with inherited glycosylation disorders.

Autosomal-Recessive Intellectual Disability with Cerebellar Atrophy Syndrome Caused by Mutation of the Manganese and Zinc Transporter Gene SLC39A8.

Manganese (Mn) and zinc (Zn) are essential divalent cations used by cells as protein cofactors; various human studies and animal models have demonstrated the importance of Mn and Zn for development. Here we describe an autosomal-recessive disorder in six individuals from the Hutterite community and in an unrelated Egyptian sibpair; the disorder is characterized by intellectual disability, developmental delay, hypotonia, strabismus, cerebellar atrophy, and variable short stature. Exome sequencing in one affected Hutterite individual and the Egyptian family identified the same homozygous variant, c.112G>C (p.Gly38Arg), affecting a conserved residue of SLC39A8. The affected Hutterite and Egyptian individuals did not share an extended common haplotype, suggesting that the mutation arose independently. SLC39A8 is a member of the solute carrier gene family known to import Mn, Zn, and other divalent cations across the plasma membrane. Evaluation of these two metal ions in the affected individuals revealed variably low levels of Mn and Zn in blood and elevated levels in urine, indicating renal wasting. Our findings identify a human Mn and Zn transporter deficiency syndrome linked to SLC39A8, providing insight into the roles of Mn and Zn homeostasis in human health and development.

Head-and-neck squamous cell carcinoma risk in smokers: no association detected between phenotype and AHR, CYP1A1, CYP1A2, or CYP1B1 genotype.

BACKGROUND: Head-and-neck squamous cell carcinoma (HNSCC) differs between smokers and nonsmokers in etiology and clinical presentation. Because of demonstrated unequivocal involvement in smoking-induced cancer in laboratory animals, four candidate genes--AHR, CYP1A1, CYP1A2, and CYP1B1--were selected for a clinical genotype-phenotype association study of HNSCC risk in smokers. Thirty-six single-nucleotide variants (mostly tag-SNPs) within and near these four genes [16 (AHR), 4 (CYP1A1), 4 (CYP1A2), and 12 (CYP1B1)] were chosen. METHODS: Extreme discordant phenotype (EDP) method of analysis was used to increase statistical power. HNSCC patients--having smoked 1-40 cigarette pack-years--represented the "highly-sensitive" (HS) population; heavy smokers having smoked ≥80 cigarette-pack-years without any type of cancer comprised the "highly-resistant" (HR) group. The vast majority of smokers were intermediate and discarded from consideration. Statistical tests were performed on N = 112 HS and N = 99 HR DNA samples from whole blood. CONCLUSIONS: Among the four genes and flanking regions--one haploblock, ACTTGATC in the 5' portion of CYP1B1, retained statistical significance after 100,000 permutations (P = 0.0042); among our study population, this haploblock was found in 36.4% of African-American, but only 1.49% of Caucasian, HNSCC chromosomes. Interestingly, in the 1000 Genomes Project database, frequency of this haplotype (in 1322 African and 1006 Caucasian chromosomes) is 0.356 and 0.003, respectively. This study represents an excellent example of "spurious association by population stratification". Considering the cohort size, we therefore conclude that the variant alleles chosen for these four genes, alone or in combinations, are not statistically significantly associated with risk of cigarette-smoking-induced HNSCC.

Genetic risk prediction: individualized variability in susceptibility to toxicants.

Genetic risk prediction uses genetic data to individualize prediction of outcome or effect from a known harmful toxicant. Several examples of toxicogenetics (usually binary traits) are discussed, reflecting largely Mendelian traits before the Human Genome Project began in 1990. Numerous complexities of the genome and what constitutes "a gene" have emerged during these past two decades. Examples of toxicogenomics (continuous outcomes, gradients) are examined. Most xenobiotic-induced environmental diseases resemble human complex diseases or other multifactorial traits such as height; these traits result from hundreds of low-effect genes. Consequently, uncovering an association between a trait and a genetic variant in a large cohort can provide important information about underlying biology; however, screening for a specific variant in an individual worker or patient has poor predictive value and little clinical utility. Individualized risk assessment for toxicants that cause environmental diseases, although a lofty goal, remains to be achieved.

Update of the human and mouse Fanconi anemia genes.

Fanconi anemia (FA) is a recessively inherited disease manifesting developmental abnormalities, bone marrow failure, and increased risk of malignancies. Whereas FA has been studied for nearly 90 years, only in the last 20 years have increasing numbers of genes been implicated in the pathogenesis associated with this genetic disease. To date, 19 genes have been identified that encode Fanconi anemia complementation group proteins, all of which are named or aliased, using the root symbol "FANC." Fanconi anemia subtype (FANC) proteins function in a common DNA repair pathway called "the FA pathway," which is essential for maintaining genomic integrity. The various FANC mutant proteins contribute to distinct steps associated with FA pathogenesis. Herein, we provide a review update of the 19 human FANC and their mouse orthologs, an evolutionary perspective on the FANC genes, and the functional significance of the FA DNA repair pathway in association with clinical disorders. This is an example of a set of genes--known to exist in vertebrates, invertebrates, plants, and yeast--that are grouped together on the basis of shared biochemical and physiological functions, rather than evolutionary phylogeny, and have been named on this basis by the HUGO Gene Nomenclature Committee (HGNC).

Smoke carcinogens cause bone loss through the aryl hydrocarbon receptor and induction of Cyp1 enzymes.

Smoking is a major risk factor for osteoporosis and fracture, but the mechanism through which smoke causes bone loss remains unclear. Here, we show that the smoke toxins benzo(a)pyrene (BaP) and 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) interact with the aryl hydrocarbon receptor (Ahr) to induce osteoclastic bone resorption through the activation of cytochrome P450 1a/1b (Cyp1) enzymes. BaP and TCDD enhanced osteoclast formation in bone marrow cell cultures and gavage with BaP stimulated bone resorption and osteoclastogenesis in vivo. The osteoclastogenesis triggered by BaP or RANK-L was reduced in Ahr(-/-) cells, consistent with the high bone mass noted in Ahr(-/-) male mice. The receptor activator of NF-κB ligand (RANK-L) also failed to induce the expression of Cyp1 enzymes in Ahr(-/-) cells. Furthermore, the osteoclastogenesis induced by TCDD was lower in Cyp1a1/1a2(-/-) and Cyp1a1/1a2/1b1(-/-) cultures, indicating that Ahr was upstream of the Cyp enzymes. Likewise, the pharmacological inhibition of the Cyp1 enzymes with tetramethylsilane or proadifen reduced osteoclastogenesis. Finally, deletion of the Cyp1a1, Cyp1a2, and Cyp1b1 in triple knockout mice resulted in reduced bone resorption and recapitulated the high bone mass phenotype of Ahr(-/-) mice. Overall, the data identify the Ahr and Cyp1 enzymes not only in the pathophysiology of smoke-induced osteoporosis, but also as potential targets for selective modulation by new therapeutics.