Harry Newmark
Widad Dantoft
Peter Ghazal- Division of Infection and Pathway Medicine, School of Biomedical Sciences, University of Edinburgh, Edinburgh, UK
In vertebrate animals, the sterol metabolic network is emerging as a central player in immunity and inflammation. Upon infection, flux in the network is acutely moderated by the interferon (IFN) response through direct molecular and bi-directional communications. How sterol metabolism became linked to IFN control and for what purpose is not obvious. Here, we deliberate on the origins of these connections based on a systematic review of the literature. A narrative synthesis of publications that met eligibility criteria allowed us to trace an evolutionary path and functional connections between cholesterol metabolism and immunity. The synthesis supports an ancestral link between toxic levels of cholesterol-like products and the vitamin D receptor (VDR). VDR is an ancient nuclear hormone receptor that was originally involved in the recognition and detoxification of xenobiotic marine biotoxins exhibiting planar sterol ring scaffolds present in aquatic environments. Coadaptation of this receptor with the acquisition of sterol biosynthesis and IFNs in vertebrate animals set a stage for repurposing and linking a preexisting host-protection mechanism of harmful xenobiotics to become an important regulator in three key interlinked biological processes: bone development, immunity, and calcium homeostasis. We put forward the hypothesis that sterol metabolites, especially oxysterols, have acted as evolutionary drivers in immunity and may represent the first example of small-molecule metabolites linked to the adaptive coevolution and diversification of host metabolic and immune regulatory pathways.
Introduction
Host-protection pathways against foreign harmful exogenous agents, inclusive of biotoxins and pathogens, exist in all branches of life. Pathways that allow the removal of biotoxins and metabolic by-products are considered to be distinct from those that neutralize and eliminate pathogens. For instance, it is understood that the P450 enzymes, which represent an ancient detoxification system, and interferon (IFN) pathways, that are central for immunity against infection in animals, are biologically unrelated. However, could specific metabolic pathways and metabolites provide an interconnection?
The substrates for P450 enzymes, while highly diverse, are lipophilic molecules often containing multiple planar ring structures. Notably, the most highly related P450s across the different kingdoms are involved in the metabolism of sterols (constituting multiple planar ring lipophilic molecules) and which further contribute an essential enzymatic role in the production of endogenous lipid metabolites, in particular as part of the sterol biosynthesis pathway (1–3). It has been debated whether the adaptation of P450 enzymes to the biosynthesis of sterols became firmly established in early eukaryotic (or late-stage prokaryotic) evolution with the arrival of atmospheric oxygen leading to the production of cholesterol in animals, ergo-sterol in fungi, and phyto-sterols in plants (4–7). The primary driver for sterol biosynthesis evolution was likely the selective advantage imparted by cholesterol toward modulation of membrane properties. However, too much cholesterol in membranes of cells, especially in the endoplasmic reticulum, the site of biosynthesis, can be highly toxic and accordingly sterol production, storage, and elimination is under stringent homeostatic regulation.
Sterols are not only required for membranes but also for the synthesis of steroid hormones, which regulate diverse physiological functions ranging from reproduction to stress and immunity. Outside the well-known functions of steroids, sterols, and in particular oxidized cholesterol and sterol metabolites, oxysterols, have been more recently found to have key roles in immunity (8, 9). Most importantly, the regulation of metabolic flux in cholesterol biosynthesis is directly linked to immune control through coupling to IFN signaling (10–13). Also see Robertson and Ghazal (14) for a review of our most current understanding of how IFN regulation is molecularly wired to sterol biosynthesis. We, therefore, posit that natural selection may have coadapted sterol metabolism and secondary metabolites as a link between functionally unrelated host-protection pathways in countering harmful chemical and biological agents. This proposition evokes the question of whether there is evidence for an ancestral gene that supports a link between these distinct host-protection pathways?
To address this question we sought to systematically review and provide a narrative synthesis of the literature based on investigating the ancestral connections between sterol metabolites, immunity, and xenobiotics. We find evidence supporting an evolutionary course for co-opting the ancestral, xenobiotic binding, vitamin D receptor (VDR) to adaptively recognize a specific non-typical oxysterol molecule, 1,25-dihydroxyvitamin D3, that in present day mammals governs prominent functions in calcium homeostasis, and immunity. It is important to clarify that vitamin D3 is a ring-opened version of 7-dehydrocholesterol and hence of the general class of sterols and steroids. For this reason and although vitamin D3 metabolites are not derived from cholesterol, we consider 25-hydroxyvitamin D3 (the inactive form found in serum) and 1,25-dihydroxyvitamin D3 (the active ligand to VDR) as non-typical oxysterols; as they are oxidized forms of the ring-opened cholesterol precursor 7-dehydrocholesterol.
On the basis of evidence presented, we further hypothesize that sterols and their oxidized metabolites have contributed as key evolutionary drivers for repurposing ancestral nuclear hormone receptors, in particular VDR, from protecting against harmful lipids to become important regulators of immunity.
The Nuclear Hormone Receptor Family Connection
Central to the recognition of sterol-like molecules and activation of detoxification systems and immunity are the nuclear hormone receptors (15, 16). In particular, the subfamily known as NR1I that includes the pregnane X receptor (PXR), constitutive androstane receptor (CAR), and the VDR (15, 17). Although each have important individual functions in humans, these receptors act as important regulators of P450 enzymes and have strong genetic evidence suggesting they originated from a single ancestral nuclear receptor (18).
Notably, VDR that is activated by a specific ligand, 1,25-dihydroxyvitamin D3, generated from vitamin D that is derived from a precursor of cholesterol, 7-dehydrocholesterol from the sterol biosynthesis pathway and synthesis in humans begins in the skin upon exposure to ultraviolet B (UVB) light emitted from the sun. The vitamin D synthesis pathway is summarized in Figure 1 (for notation see Table 1), and involves the skin, liver, and kidneys (19). Interestingly, animals with fur and feathers are still able to synthesize vitamin D from sunlight despite UVB not reaching the skin (20). Here, vitamin D synthesis occurs through the sebaceous glands producing oily secretions (containing 7-dehydrocholesterol) that cover fur or feathers and ingested after grooming (21–23).
Figure 1. The vitamin D synthesis pathway. Vitamin D3 can be acquired both endogenously from cholesterol in the skin and exogenously through diet (vitamin D2 and vitamin D3). In the skin, 7-dehydrocholesterol, a cholesterol precursor, is converted to previtamin D3 upon ultraviolet B (UVB) exposure. Through a series of cytochrome P450 enzyme-mediated reactions, previtamin D3 is converted first into 25(OH)2D3 in liver hepatocytes and then activated in the kidney by 1α-hydroxylation, to form 1,25(OH)2D3. The degradation of 1,25(OH)2D3 and intermediate metabolites is mediated by negative feedback mechanisms (24). Figure created using the SBGN format on VANTED (25). See Table 1 for notation.
Table 1. Systems biology graphical notation legend.
Vitamin D deficiency can result in clinical disorders, the most notable being the characteristic bow-legged musculoskeletal manifestation known as rickets. Additional studies have also linked deficiency to cardiovascular disease, cancer, autoimmune conditions, and decreased antimicrobial protection (26–28). The VDR is known to heterodimerize with retinoid X receptor (RXR) and exerts biological effects as a ligand activated transcription factor by binding to specific vitamin D response element (VDRE) in gene promoters of over 200 genes (29). The active ligand of VDR is 1,25-dihydroxyvitamin D3 (1,25(OH)2D3), also known as calcitriol, although certain bile acids are also capable of inducing transactivation to a lesser degree (Figure 1). 1,25(OH)2D3 is the active form of vitamin D, produced by enzymatic hydroxylation of the circulating 25-hydroxyvitamin D3 by cytochrome P450 (CYP) enzyme 27B1 (see Figure 1).
While vitamin D and its receptor have been long regarded as mediators of calcium and phosphate homeostasis, VDR has additional roles in innate and acquired immunity and xenobiotics (30). Vitamin D-mediated calcium homeostasis has been around since the first terrestrial vertebrates, including amphibians, which have also been observed to suffer from calcium deficient ailments such as rickets (31). Before the calcium endocrine system, ancient VDR functioned as a xenobiotic receptor, mediating the degradation of marine biotoxins (32). It still retains this ability, and in humans VDR is important in detoxifying the toxic secondary bile acid lithocholic acid (LCA) in the colon by activating the CYP3A4 P450 enzyme (30).
Vitamin D is itself an ancient sterol–steroid, present in phytoplankton and zooplankton (33). VDR orthologs have likewise been observed in ancient vertebrates (34) and invertebrates (35). Accordingly, we next examine the origin and early evolutionary progression of VDR and its role in detoxification.
Primordial Nuclear Hormone Receptor Family of the VDR and Biological Roles in Detoxification
The VDR is descended from a group of xenobiotic nuclear hormone receptors known as NR1I, as shown in Figure 2. The NR1J nuclear receptor subfamily has important roles for xenobiotic detoxification in arthropods and nematodes (15, 16, 36). For instance, the related DHR96 receptor in Drosophila melanogaster can bind and detoxify a phenobarbital insult through a CYP transcriptional response (37), and NHR-8, in the Caenorhabditis elegans gut, senses colchicine and targets the activation of its cognate detoxification pathway (38). Many of these related nuclear receptors have also shown the ability to bind small lipophilic molecules such as cholesterol or steroid hormones (15, 16, 36). In this way, conservation of function can be seen throughout the evolution of these nuclear receptors. We will first discuss the ancestral xenobiotic role of VDR before considering how evolutionary pressures may have promoted the functional repurposing of this receptor with the acquisition of new roles including detoxification of endogenous compounds (e.g., vitamin D metabolites and bile acids), lipid regulation, and immunity.
Figure 2. Phylogenetic relationship of vertebrate vitamin D receptor (VDR) with invertebrate orthologs. Both vertebrate and invertebrate VDR, belonging to the closely related NR1I and NR1J subgroups, respectively, regulate the expression of genes involved in xenobiotic metabolism (detoxification). Evolutionary change in the vertebrate VDR has resulted in the acquisition of new immunological functions. Distances between nodes are not to scale.
The tunicate, Ciona intestinalis, represents the closest extant invertebrate relative of vertebrates possessing an ancestral VDR gene (39). Fidler et al. (35) investigated the potential ligand-binding properties of this receptor, named as CiVDR/PXRα for its homology with both the VDR and PXR. Interestingly, vitamin D, or indeed any bile salts, were unable to produce any transactivation potential (40, 41). Despite speculation that the CiVDR/PXRα ortholog may be used for calcium homeostasis, there is evidence for closer functional similarity to PXRs current role in xenobiotics. This possibility is supported by the argument that the ocean is a plentiful source of calcium, making any need for homeostasis redundant. Present day PXR function in humans is to detoxify foreign toxic compounds by sensing and then activating the enzyme CYP3A4. There is good experimental evidence to support a functional role for the C. intestinalis CiVDR/PXRα ortholog to be ligand activated by microalgal biotoxins, including okadaic acid and pectenotoxin-2 (35) (Figure 3). Filter feeding tunicates like C. intestinalis accumulate these biotoxins through the large quantity of microalgae in their diet. As high concentrations of these chemicals are able to kill cells, CiVDR/PXRα’s ability to bind and detoxify them would be appropriate and consistent with PXRs current role in humans. Indeed, it has been shown that orthologous NR1Jβ receptors in mollusks likewise respond to xenobiotic insult from okadaic acid by activating detoxification pathways (42).
Figure 3. Molecular structures of calcitriol, okadaic acid, and pectenotoxin-2. (A) Molecular structure of 1,25(OH)2D3, otherwise known as calcitriol (PubChem CID = 5280453), a vitamin D receptor (VDR) agonist. (B,C) Atomic structure of two natural Ciona intestinalis CiVDR/PXRα analogs, okadaic acid (446512) (B) and pectenotoxin-2 (6437385) (C), from left to right. Figures produced using MarvinSketch (http://www.chemaxon.com).
Similar to CiVDR/PXRα in C. intestinalis, NHR-8 and DHR96, members of the NR1J subfamily, in C. elegans and D. melanogaster, respectively, have also been shown to be essential for mediating xenobiotic resistance by promoting the expression of genes involved in metabolism of endo- and xenobiotics (37, 43, 44). Increased expression of genes involved in xenobiotic metabolism, together with resistance to xenobiotics, are frequently correlated with lifespan extension in C. elegans, D. melanogaster, and mice, suggesting detoxification of diet-acquired toxins is a host-protection mechanism against accumulation of specifically lipophilic toxins that negatively impact health during aging (The Green Theory of Aging) (16, 45). However, a recent study by Afschar et al. showed that DHR96 is indeed essential for mediating resistance to xenobiotics but not for increasing lifespan of insulin-mutant flies (44), indicating that xenobiotic resistance and longevity may not be causally connected. It has been suggested that the co-occurrence of xenobiotic resistance and lifespan extension may have co-evolved because lowered insulin/insulin-like growth factor signaling (IIS) can also signal the presence of pathogens (44). In line with this concept, in C. elegans and D. melanogaster, genes involved in xenobiotic metabolism have also been shown to be indirectly activated by toxic microbial by-products that directly cause dysfunction in cellular processes such as an altered metabolism, and decreased host translation and IIS (43, 46, 47).
Evolutionary Repurposing of Primordial VDR from Exogenous to Endogenous Detoxification Pathways
Pharmacophore modeling of CiVDR/PXRα ligands revealed specific chemical scaffolds were required for receptor binding, comprising two hydrophobic features (in particular aromatic rings) and one hydrogen bond acceptor in a planar arrangement (35). The structure of activated vitamin D exhibits resemblance to this scaffold presenting a planar conformation with aromatic rings and hydrophobic features (48).
Further sequence identification by Ekins et al. (41) demonstrated a 67.6 and 17.1% similarity between the DNA binding domain (DBD) and the ligand-binding domain (LBD), respectively, between CiVDR/PXR and hVDR. Lower conservation of the LBD suggests evolutionary adaptive changes in ligand affinity. It is perhaps this promiscuous ligand-binding quality that allowed the ancestral VDR to function as a xenobiotic receptor, by binding and detoxifying new toxic chemicals on exposure. The similarities in structure between 1,25(OH)2D3 and exogenous marine biotoxins, alongside the genetic variability of VDRs’ LBD, allowed for the eventual binding and regulation of 1,25(OH)2D3 levels, a sterol-derived metabolite. As described below, observations in basal vertebrates provide insight toward understanding how this evolutionary pressure may have been applied.
The lamprey (Petromyzon marinus) is the most basal extant vertebrate, and therefore, provides valuable information regarding VDR functional evolution. The lamprey (lampVDR) has an 87 and 60% homology with hVDR DBD and LBD, respectively (18, 34). This differential increase in homology of the LBD in comparison with the tunicate CiVDR/PXR ancestral gene suggests VDRs’ ligand activated role may have changed. There have been three rounds of whole genome duplications (WGDs), represented as 1R, 2R, and 3R, since C. intestinalis. Humans diverged after the second round while Teleost fish were subjected to a third round and therefore, have an extra copy of the VDR gene (18). It is thought that after 1R, the combined VDR/PXR gene split giving rise to separate VDR and PXR genes with different but also overlapping functions. The lamprey was first to diverge after the 1R WGD event, as shown in Figure 2. While, a PXR homolog has yet to be identified in the lamprey, lampVDR has been shown to have high affinity binding and transactivation by 1,25(OH)2D3, which functionally activates CYP3A4, and possibly CYP24A1, in order to detoxify high levels of this sterol (34) and providing an opportunity for an extended regulatory role for the VDR in lipid metabolism. Although the lamprey lacks a calcium endocrine system, lampVDR may contribute to the regulation of other processes such as skin differentiation, as an observed increase in VDR presence in mucous glands and keratinized teeth has been reported (34, 49).
Thus, in addition to the metabolism of exogenous xenobiotic compounds, VDR further acquired an ability to detoxify certain lipophilic endogenous molecules such as bile acids. Bile salt pathways are important vertebrate mechanisms by which cholesterol can be removed from the body. There are at least three evolutionary classified bile salt pathways, referred to as early, later, and recent pathways (41) (Figure 4). The lamprey uses the early fish pathway, while chondrichthyes, such as the Little skate, use the later pathway. In mammals, the “recent” pathway converts cholesterol to 24-carbon (C24) bile acids, which can subsequently be converted to toxic secondary bile acids in the intestine by resident microorganisms (50). Indeed, one of the major roles of the hVDR is its ability to detoxify the secondary bile acid LCA by transactivation of CYP3A4 (30). As LCA is a product of the most recent C24 bile acid pathway, basal vertebrates such as the lamprey, which employ the “early” or “later” fish pathways, are unable to bind this molecule with their VDR (51). In mammals, LCA affinity is almost certainly a more recent evolutionary adaption to changes in the gut microbiome resulting in the production of toxic secondary bile acids. This ability of VDR to subsequently transactivate the CYP3A4 gene in response to the binding of a toxic chemical traces an evolutionary path from VDRs’ ancestral xenobiotic function. In this regard, it is noteworthy that the NHR-8 and DAF-12 nuclear receptors in C. elegans, homologs of VDR and part of the NR1J subgroup, have convergently evolved to control and bind dafachronic acids, a bile acid look-alike, important in the life-cycle of this species (52, 53).
Figure 4. Schematic overview of the three major bile salt pathways throughout vertebrate evolution. All bile acids are derived from cholesterol, a 27-carbon molecule. In early fish, e.g., hagfish and sea lamprey, a 7α-hydrocholesterol is converted to 5α-bile alcohols, followed by conversion to 5α-bile alcohol sulfate. The production of bile salts in mammals, birds, cartilaginous fish, and in some teleost and amphibians is dependent on the conversion of 7α-hydrocholesterol into 5β-bile alcohol. In cartilaginous fish, mainly, but also in some teleost and amphibians, 5β-bile alcohol is further converted to 5β-bile alcohol sulfates (“Later pathway”). In mammals and birds, as depicted in the “Recent pathway,” bile salt production involves the conversion of 5β-bile alcohol to 24-carbon bile acids in the liver (cholic acid and chenodeoxycholic acid) by cytochrome P450-mediated oxidation. When these bile acids are secreted into the lumen of the intestine, cholic acid and chenodeoxycholic acid are converted, by colonic bacteria, to the secondary bile acids deoxycholic acid and lithocholic acid, respectively. While secondary bile acids in higher concentrations are potentially toxic to cells, they can, together with primary bile acids, be taken up into the blood stream and liver for re-secretion. In all pathways, the production of 7α-hydroxylation is the rate-limiting step in these reactions. Figure based on Ekins et al. (41).
Functional Diversification of VDR in Lipid Metabolism and Immunity
Vitamin D receptor’s ancestral roots clearly stem from its ability to recognize and regulate detoxification pathways for environmental toxic chemicals, as observed in both NR1I and NR1J subgroups. From this stemmed, its ability to bind and detoxify endogenous lipophilic molecules such as 1,25(OH)2D3 and bile acids using P450 enzymes. Sterol metabolites are important biological molecules requiring careful regulation for incorporation into cell membranes and steroid hormones. Accordingly, this opens a new opportunity for VDR to adopt a regulatory role through negative feedback mechanisms conferred by its ligand binding. The next section will discuss how sterol metabolites may have further coadapted VDR for driving the diversification of VDR functionality into lipid regulation, immunity, and bone prior to a role for VDR in calcium homeostasis (34).
As described above, VDR origin is based on recognizing and regulating the levels of lipophilic exogenous and endogenous molecules through the transcriptional regulation of P450 enzymes. To more fully understand the functional relationship by which VDR interacts with various lipid metabolism and immune signaling pathways, we constructed a pathway biology diagram (Figure 5, for notation see Table 1) from the research synthesis of literature mapping all known interactions. Figure 5 shows that VDR is deeply embedded in a network of signaling pathways including PPAR-α and PPAR-γ, nuclear factor-kappa B (NF-κB), p38 mitogen-activated protein kinase (MAPK), transforming growth factor-beta (TGF-β), and eicosanoid synthesis. Of note, all of these pathways have instrumental roles in immunity, bone regulation, cell proliferation, and lipid metabolism.
Figure 5. Interactions of the vitamin D receptor (VDR) with PPAR-α, nuclear factor-kappa B (NF-κB), p38 mitogen-activated protein kinase (MAPK), transforming growth factor-beta (TGF)-β, β-oxidation, eicosanoid production, and nuclear factor of activated T-cells (NFAT) pathways. From left to right, the pathways represented are PPAR-α/γ (blue), β-oxidation (light purple), NF-κB (red), p38 MAPK (green), TGF-β (turquoise), NFAT (dark purple), and eicosanoid synthesis (yellow). Figure created using the SBGN format on VANTED with edges representing high and low confidence interactions (25). See Table 1 for notation.
A range of biological functions associated with VDR and targeted signaling pathways are shown in Table 2. Notably, most of these pathways not only regulate lipid metabolism, but have direct roles in immunity and bone homeostasis too. All these processes are interconnected. Lipids have a fundamental role in the immune system and directly affect immune cell function. They alter membrane fluidity, lipid peroxidation, gene expression, and eicosanoid production (54). Bone and immunity have a very close relationship, with bone marrow being a “high fat” primary hematopoietic tissue, controlling the production of B cells and other important innate and adaptive immune responses (55). Sterol biosynthesis is at the heart of this control network, using sterol-based interactions with the VDR to effectively regulate lipid metabolism and its associated features within bone and immunity. These pathways provide a link between the seemingly unconnected and multiple divergences of the VDR before additional more specific gene targets became evolutionary fixed and conserved, such as calcium homeostasis. In addition, it was also found that VDR has the ability to directly regulate fatty acid beta-oxidation by interacting with the enzymes hexokinase, CBT1, and CBT2, possibly leading to a role for energy metabolism in adipose tissue (56, 57).
Table 2. The function of various signaling pathways and their corresponding interaction with the vitamin D receptor (VDR).
It is worth noting that all these interactions are compartmentalized in terms of tissue specific, time dependent, and multifactorial control levels. This view is consistent with the possibility that cholesterol, as a toxic molecule and precursor sterols to vitamin D, repurposed VDRs’ ancestral function to provide a wider regulation over lipid metabolism and immune pathways.
Vitamin D receptor has long been known to promote immune tolerance in the acquired immune system while providing protective innate mechanisms against pathogen infection. The acquired immune functions of the VDR are complicated, involving the regulation of multifaceted signaling pathways such as PPAR-γ and NF-κB. The net outcome of these cross-regulatory responses results in attenuation of the immune adaptive response, involving stimulation of interleukin (IL)-10 and downregulation of IL-12 (75). The most notable case for innate immunity is the VDR function to induce expression of the cathelicidin antimicrobial peptide (CAMP) gene, an important host defense protein. Gombart et al. (76) demonstrated that exaptation of an AluSx short interspersed element provided a perfect VDRE in the Camp promoter. Moreover, VDR is upregulated during infection in a toll-like receptor 2/1 (TLR2/1)-dependent manner (77–79). This shows a direct connection to innate immunity and is where the historic treatment of tuberculosis with cod liver oil (high in vitamin D) may have a possible molecular explanation (80). Unsurprisingly, VDR-mediated innate immune responses have become targets of pathogen evasion techniques (81).
The Connection with IFN Signaling
Interferons are a group of signaling proteins required for antiviral defense. Released by virally infected cells or leukocytes, they mediate a variety of innate and adaptive immunological responses by upregulating over 300 interferon-stimulated genes (ISGs) (82). IFNs can be classified into three types depending upon the receptor to which they bind and the signal transduction pathways they activate (83). Type I IFNs are split into multiple subtypes including -α, -β, -ω, -ε, -τ, -δ, and -κ, each with independent and redundant functions. Type II IFNs are conserved to just higher mammals and have only one member, IFNγ. Lastly, type III, containing IFN-λ genes (IL-29, IL-28A, and IL-28B), have similar biological properties as type I, but their genetic sequence contains non-coding intron sequences (84–86). The IFN system has displayed remarkable conservation throughout vertebrate evolution, demonstrating its importance for immunological defense (87). It is also inherently linked to cholesterol metabolism, as described elsewhere (8–10, 14).
Interferons appeared to have originated soon after the evolution of vertebrates as IFN homologous genes and their transcription factors have not yet been observed in primitive chordates such as the tunicate and sea urchin (87) or even closely related basal vertebrates such as the jawless fish, lamprey (88). We briefly discuss the evolution of type I and III IFN genes, as their ancestral homologs have coexisted with the evolving metabolic interactions between VDR, cholesterol, and immunity.
Interferon genes are present in many different varieties of fish, including the teleost clade of ray-finned fish that diverged from our evolutionary line 450 million years ago (87). Fish can possess singular or multiple IFN genes depending on the species and it is likely that vertebrate groups have independently evolved a vast array of structurally similar IFN molecules that perform different host protective functions (83, 89). Thus, the expanding role of IFNs coincided with the evolutionary changing roles of sterols and VDR.
Interferons act through the Janus kinase/signal transducer and activator of transcription (JAK/STAT) signal transduction pathways inducing and suppressing hundreds of genes. In mammals, IFNs are activated through a variety of pathogen pattern recognition systems, notably downstream of toll-like receptor (TLR) activation and by stimulator of interferon genes (STING)-activation that are now known to also target the sterol biosynthesis pathway. IFN signaling via STAT1 and IRF1 induces a cholesterol hydroxylase gene, CH25H and its cognate metabolite 25-hydroxycholesterol as well as a microRNA (miR342-5p) that dramatically suppresses the flux in the sterol biosynthesis pathway. A change in the flux of sterol biosynthesis in turn activates STING that further re-enforces the IFN response (14). It is worth mentioning that VDR can interact with STAT1 and curtail the nuclear translocation of ISGF3 and which may contribute to VDRs inherent immunosuppressive attributes (90). Furthermore, as mentioned above, numerous studies have reported on the importance of both IFN-expression and activation of the vitamin D-pathway on the expression of downstream effector molecules [e.g., antimicrobial peptides (AMPs)] that subsequently influence infection and inflammation (78, 91–99). Increased expression of type II IFN (IFNγ) have been correlated with macrophage activation, macrophage-dependent AMP gene expression, as well as with controlled growth of pathogenic intracellular microbes and better disease outcome (96, 97, 99). Fabri et al. reported that IFNγ, released by T cells induce in a vitamin D-pathway-dependent in human macrophages, autophagy, phagosomal maturation, and antimicrobial activity against Mycobacterium tuberculosis (97).
The role of IFN responses and the vitamin D-pathway has also been investigated in human leprosy. Teles et al. have revealed an inverse correlation between IFNβ, IL-10, and IFNγ, where IFNβ, in an IL-10-pathway-dependent manner, inhibited the IFNγ-induced and vitamin D-dependent, AMP response in disseminated and progressive lepromatous lesions (99). By contrast, IFNγ-specific genes were enriched in self-healing tuberculoid lesions. Both studies underscore the importance of adequate amounts of vitamin D in human populations for sustaining both innate and acquired immunity against infection. The close connection between IFN and the vitamin D-pathway have also been reported in experimental autoimmune encephalitis (EAE), a model for multiple sclerosis (MS) (100, 101), and diabetes (102), two non-infectious disorders characterized by excessive and uncontrolled inflammation and macrophage foam cell formation (characterized by accumulation of esterified cholesterol) (103, 104). Early studies have shown that IFNγ plays a crucial role in the induction of 1,25(OH)2D3 (105, 106), the active vitamin D metabolite that bind to VDR, in initializing VDR dimerization with RXR and, in VDR–RXR activation of VDRE-containing target genes (107–109). Adams et al. showed that IFNγ induces production of 1,25(OH)2D3 in macrophages and that the effect was abolished by addition of anti-IFNγ to the culture medium (110). The tissue availability of 1,25(OH)2D3 in immune cells is dependent on the expression of the activating enzyme 1α-hydroxylase (Cyp27b1) and its catabolic counterpart 25-OHD3-24-hydroxylase (Cyp24a1) (111). In addition to TLR signaling, expression of 1α-hydroxylase can be induced by IFNγ stimulation (112–114). Activation by IFNγ stimulation require, however, the cells to be differentiated, as IFNγ stimulation of undifferentiated monocyte THP1 cells failed to induce Cyp27b1 expression (1α-hydroxylase) in the absence of a second stimulus [lipopolysaccharide (LPS)] (112). On the other hand, Cyp24a1 (25-OHD3-24-hydroxylase) expression is induced by the type 2 T helper cell cytokine IL-4, in toll-like receptor 2/1 ligand-activated monocytes, but not by IFNγ (114). Stoffel et al. revealed that, in cultured monocytes, synergistic induction of Cyp27b1 gene expression by IFNγ and LPS, required, not only activation of the JAK/STAT pathway and NF-κB binding but also binding of phosphorylated C/EBPβ (by the p38 MAPK pathway) (113).
In addition to regulation of VDR through the transcriptional regulation of Cyp27b1, VDR gene expression was shown to be dependent on the IFNG gene (101). Ifng knockout (GKO) mice exhibited very low Vdr gene expression in the central nervous system. Correlating with the low Vdr expression, GKO mice also demonstrated an increased pathogenic T cells burden as well as a more severe EAE phenotype, suggesting that the aggressive autoimmune CD4+ T cell phenotype may be a consequence of inadequate Vdr gene expression (101). Furthermore, treatment of TLR2/1L-activated monocytes with IFNγ has been shown to not only stimulate Cyp27b1 expression but also Vdr gene expression (114). Collectively, these studies suggest that IFN-regulation of VDR activity is complex and regulated indirectly via regulation of 1α-hydroxylase, and possibly directly at the gene level by IFNγ-induced STAT binding.
Present Day IFN-Sterol Metabolic Link
Primitive NR1I and NR1J receptors have inherent affinities for lipophilic molecules, making them useful xenobiotic sensors in vertebrate and invertebrate organisms. Figure 6 depicts the adaption of this receptor from the detoxification of endogenous compounds including bile acids and vitamin D in vertebrate organisms to acquiring new biological roles associated with the increasing importance of sterols and especially oxysterols for immune cell function. From this foundation, VDR continued to develop further direct and indirect interactions with lipid metabolism pathways, placing VDR in an opportunistic position to influence the development and regulation of other important biological systems including, but not limited to, immunity, bone development, and cell proliferation. The immune roles have been further consolidated with the acquisition of fixed mutations, such as the VDR regulated production of innate AMPs.
Figure 6. Vitamin D receptor (VDR) is an ancient nuclear hormone receptor that was originally involved in the recognition and detoxification of xenobiotic marine biotoxins exhibiting sterol ring scaffolds present in aquatic environments. Coadaptation of this receptor with the acquisition of sterol biosynthesis and interferons in vertebrate animals set a stage for repurposing and linking a preexisting host-protection mechanism of harmful xenobiotics to become an important regulator in immunity.
In conclusion, our findings support an evolutionary basis to the IFN–sterol–immune metabolic link, arising from the xenobiotic origins of a nuclear hormone receptor, VDR (Figure 6). Xenobiotics and immunity are considered separate pathways facilitating the removal of toxic molecules and pathogens, respectively. However, our narrative synthesis suggests that sterol metabolites acted as an evolutionary driver integrating a complex metabolic network under bi-directional IFN and VDR control. We believe this is likely to be a more general evolutionary mechanism for other nuclear hormone receptor functions in immunity such as glucocorticoid receptor, liver X receptor (LXR), and farnesoid X receptor (FXR) and possibly for other distinct immune directed metabolic pathways and associated ligand activated receptor systems.
Concluding Remarks on the General Hypothesis of Microbial–Host Small-Molecule Metabolites as Functional and Evolutionary Drivers for Recognizing and Defending Against Non-Self
The question of whether a functional–evolutionary link between planar sterol-like molecules, immunity, and detoxification is specifically unique to VDR would implicate a limited selective role rather than one pertaining to a more central evolutionary principle with broader biological significance. Notably, in this context other evolutionary related subfamily 1 nuclear receptors, including LXR, FXR, PXR, CAR, retinoic-acid receptor (RAR), and RAR-related orphan receptor (32), similarly exhibit, at a number of levels, cross talk between mevalonate-sterol metabolism and immunity, as well as in xenobiotic detoxification [for examples see Ref. (115–118)]. Interestingly, RXRs partner with many of these receptors and which recruits corepressor or coactivator molecules to regulate transcriptional responses. Homologs and orthologs of this highly conserved nuclear receptor have been identified in marine and terrestrial invertebrates (119–123), so it is likely that RXR coevolved with VDR and the many other nuclear receptors it is associated with. RAR and RXR are activated by the vitamin A derived ligand, retinoic acid and its 9-cis conformer, respectively. While this lipophilic molecule is not a steroid-derived molecule, the heterodimerization of RXR provides an important but insufficient role of vitamin A metabolites in integrating with these receptor systems and sterol metabolism (124–127).
Of the related subfamily 1 nuclear receptors, strong conservation across vertebrate species can be found for the LXRs, with approximately 75% sequence identity in the LBD between human and non-mammalian LXRs (128). Consistent with the high degree of sequence conservation, ligand specificities between mammalian and non-mammalian vertebrates homologs, as well as between vertebrate and non-vertebrate LXR orthologs, are very similar (128). While vertebrate LXR agonists do not activate Ciona LXR, it is activated by a number of oxysterols as well as some pregnane and androstane steroids. In insects on the other hand, the ecdysone receptor (EcR) has been identified as the ortholog for both LXR and FXR combined (120). Together with its interaction partner, ultraspiracle (ortholog of RXR), EcR play an essential role in insect development and reproduction, as well as in basic metabolism and immunity (129, 130). In mammals, FXRs, with their high expression in liver, adrenal glands, intestine, and kidney, are activated by farnysol and its metabolites, part of the mevalonate–sterol biosynthesis pathway (131) and by primary bile acids such as chenodeoxycholic acid (132–134). Together with PXRs and VDR, FXR serves as one of the major transcriptional regulators of bile salt synthesis, partly by regulating the expression of CYP7A1 and CYP8B1 (135). Bile salts have to date not been detected in invertebrates, suggesting that regulation of bile salt synthesis is an evolutionary acquired trait specific for the vertebrate nuclear receptors. Two major evolutionary shifts in bile salt structure have been identified (50, 136) and hypothesized that the bile alcohols found in jawless fish (Agnatha) represents the “ancestral” bile salt phenotype from which the more “recent” vertebrate bile acids are derived (51). Secondary bile acids produced as a consequence of microbial metabolism in the gut are further detoxified through recognition by these nuclear receptors. In this regard, it is also notable that the related NRI family member, PXR that is well known for its xenobiotic detoxification role has also been shown to regulate intestinal inflammation by sensing bacterial metabolites (137). Furthermore, another unrelated xenobiotic receptor, the aryl hydrocarbon receptor, has been linked to the antibacterial response through sensing bacterial pigments and in enhancing IL-22 barrier immunity (138, 139).
While a comprehensive evolutionary history of LXR, FXR, and the other nuclear hormone receptors is beyond the present scope of this article; we believe the evolutionary path exemplified by VDR, in selecting dual functional roles in detoxification of harmful planar lipids and immune recognition, represents a central driver for the evolution of these receptors and evokes a general hypothesis for the coevolution of microbial–host metabolism underlying host-protection pathways.
Author Contributions
HN, WD, and PG conceived, performed analysis, and wrote the paper.
Conflict of Interest Statement
The authors declare that this document was written in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors would like to thank Bill Griffiths for comments on oxysterol chemistry and helpful critiquing by the reviewers.
Funding
The Centre for Systems Biology at Edinburgh is a Centre for Integrative Systems Biology (CISB) supported by the BBSRC and EPSRC. This work was supported by the BBSRC (BB/K019112/1), the BBSRC/EPSRC (BB/D019621/1), and the Wellcome Trust (WT066784/Z/02/Z) to PG. The funders had no role in the preparation of the manuscript.
Abbreviations
AluSx, Alu Sx subfamily; AHR, aryl hydrocarbon receptor; AMP, antimicrobial peptide; C24, 24-carbon; CAMP, cathelicidin antimicrobial peptide; CAR, constitutive androstane receptor; CBT, cytochrome b termination protein; CH25H, cholesterol 25-hydroxylase; CNS, central nervous system; CYP, cytochrome P450; DBD, DNA binding domain; DHR96, Drosophila hormone receptor-like in 96; EAE, experimental autoimmune encephalitis; EcR, ecdysone receptor; FXR, farnesoid X receptor; GKO, interferon-gamma knock out; GR, glucocorticoid receptor; HK, hexokinase; JAK/STAT, janus kinase/signal transducer and activator of transcription; IIS, insulin/insulin-like growth factor signaling; IL, interleukin; IFN, interferon; IRF, interferon response factor; ISG, interferon-stimulated gene; ISGF, interferon-stimulated gene factor; LBD, ligand-binding domain; LCA, lithocholic acid; LPS, lipopolysaccharide; LXR, liver X receptor; MAPK, mitogen-activated protein kinase; miR342-5p, microRNA 342-5p; MS, multiple sclerosis; NF-κB, nuclear factor-kappa B; NHR-8, nuclear hormone receptor family member nhr-8; NFAT, nuclear factor of activated T-cells; PPAR, peroxisome proliferator-activated receptor; PXR, pregnane X receptor; RAR, retinoic-acid receptor; ROR, RAR-related orphan receptor; RXR, retinoid X receptor; RSMAD, receptor-regulated SMAD; SMAD, mothers against decapentaplegic homolog; STING, stimulator of interferon genes; TGF-β, transforming growth factor-beta; Th2, type 2 T helper cell; TLR, toll-like receptor; TLR2/1L, toll-like receptor 2/1 ligand; USP, ultraspiracle; UVB, ultraviolet B; VDR, vitamin D receptor; VDRE, vitamin D response element; WGD, whole genome duplication.
References
2. Nebert DW, Wikvall K, Miller WL. Human cytochromes P450 in health and disease. Philos Trans R Soc B (2013) 368(1612):20120431. doi: 10.1098/rstb.2012.0431
3. Parvez M, Qhanya LB, Mthakathi NT, Kgosiemang IKR, Bamal HD, Pagadala NS, et al. Molecular evolutionary dynamics of cytochrome P450 monooxygenases across kingdoms: special focus on mycobacterial P450s. Sci Rep (2016) 6:33099. doi:10.1038/srep33099
4. Yoshida Y, Noshiro M, Aoyama Y, Kawamoto T, Horiuchi T, Gotoh O. Structural and evolutionary studies on sterol 14-demethylase P450 (CYP51), the most conserved P450 monooxygenase: II. Evolutionary analysis of protein and gene structures. J Biochem (1997) 122(6):1122–8. doi:10.1093/oxfordjournals.jbchem.a021870
5. Režen T, Debeljak N, Kordiš D, Rozman D. New aspects on lanosterol 14α-demethylase and cytochrome P450 evolution: lanosterol/cycloartenol diversification and lateral transfer. J Mol Evol (2004) 59(1):51–8. doi:10.1007/s00239-004-2603-1
6. Galea AM, Brown AJ. Special relationship between sterols and oxygen: were sterols an adaptation to aerobic life? Free Radic Biol Med (2009) 47(6):880–9. doi:10.1016/j.freeradbiomed.2009.06.027
7. Brown AJ, Galea AM. Cholesterol as an evolutionary response to living with oxygen. Evolution (2010) 64(7):2179–83. doi:10.1111/j.1558-5646.2010.01011.x
8. Fessler MB. The intracellular cholesterol landscape: dynamic integrator of the immune response. Trends Immunol (2016) 37(12):819–30. doi:10.1016/j.it.2016.09.001
9. Tall AR, Yvan-Charvet L. Cholesterol, inflammation and innate immunity. Nat Rev Immunol (2015) 15(2):104–16. doi:10.1038/nri3793
10. Blanc M, Hsieh WY, Robertson KA, Watterson S, Shui G, Lacaze P, et al. Host defense against viral infection involves interferon mediated down-regulation of sterol biosynthesis. PLoS Biol (2011) 9(3):e1000598. doi:10.1371/journal.pbio.1000598
11. Blanc M, Hsieh WY, Robertson KA, Kropp KA, Forster T, Shui G, et al. The transcription factor STAT-1 couples macrophage synthesis of 25-hydroxycholesterol to the interferon antiviral response. Immunity (2013) 38(1):106–18. doi:10.1016/j.immuni.2012.11.004
12. Watterson S, Guerriero ML, Blanc M, Mazein A, Loewe L, Robertson KA, et al. A model of flux regulation in the cholesterol biosynthesis pathway: immune mediated graduated flux reduction versus statin-like led stepped flux reduction. Biochimie (2013) 95(3):613–21. doi:10.1016/j.biochi.2012.05.024
13. York AG, Williams KJ, Argus JP, Zhou QD, Brar G, Vergnes L, et al. Limiting cholesterol biosynthetic flux spontaneously engages type I IFN signaling. Cell (2015) 163(7):1716–29. doi:10.1016/j.cell.2015.11.045
14. Robertson KA, Ghazal P. Interferon control of the sterol metabolic network: bidirectional molecular circuitry-mediating host protection. Front Immunol (2016) 7:634. doi:10.3389/fimmu.2016.00634
15. Timsit YE, Negishi M. CAR and PXR: the xenobiotic-sensing receptors. Steroids (2007) 72(3):231–46. doi:10.1016/j.steroids.2006.12.006
16. Hoffmann JM, Partridge L. Nuclear hormone receptors: roles of xenobiotic detoxification and sterol homeostasis in healthy aging. Crit Rev Biochem Mol Biol (2015) 50(5):380–92. doi:10.3109/10409238.2015.1067186
17. Reschly E, Krasowski MD. Evolution and function of the NR1I nuclear hormone receptor subfamily (VDR, PXR, and CAR) with respect to metabolism of xenobiotics and endogenous compounds. Curr Drug Metab (2006) 7(4):349. doi:10.2174/138920006776873526
18. Kollitz EM, Zhang G, Hawkins MB, Whitfield GK, Reif DM, Kullman SW. Molecular cloning, functional characterization, and evolutionary analysis of vitamin D receptors isolated from basal vertebrates. PLoS One (2015) 10(4):e0122853. doi:10.1371/journal.pone.0122853
19. Deeb KK, Trump DL, Johnson CS. Vitamin D signalling pathways in cancer: potential for anticancer therapeutics. Nat Rev Cancer (2007) 7(9):684–700. doi:10.1038/nrc2196
20. Hymøller L, Jensen SK. Vitamin D3 synthesis in the entire skin surface of dairy cows despite hair coverage. J Dairy Sci (2010) 93(5):2025–9. doi:10.3168/jds.2009-2991
21. Gaylor J, Sault FM. Localization and biosynthesis of 7-dehydrocholesterol in rat skin. J Lipid Res (1964) 5(3):422–31.
22. Carpenter KJ, Zhao L. Forgotten mysteries in the early history of vitamin D. J Nutr (1999) 129(5):923–7.
23. Agarwal SC, Stout SD. Bone Loss and Osteoporosis: An Anthropological Perspective. New York, NY: Springer (2003).
24. Christakos S, Dhawan P, Verstuyf A, Verlinden L, Carmeliet G. Vitamin D: metabolism, molecular mechanism of action, and pleiotropic effects. Physiol Rev (2016) 96(1):365–408. doi:10.1152/physrev.00014.2015
25. Rohn H, Junker A, Hartmann A, Grafahrend-Belau E, Treutler H, Klapperstück M, et al. VANTED v2: a framework for systems biology applications. BMC Syst Biol (2012) 6:139. doi:10.1186/1752-0509-6-139
26. Garland CF, Garland FC, Gorham ED, Lipkin M, Newmark H, Mohr SB, et al. The role of vitamin D in cancer prevention. Am J Public Health (2006) 96(2):252–61. doi:10.2105/AJPH.2004.045260
27. Holick MF, Chen TC. Vitamin D deficiency: a worldwide problem with health consequences. Am J Clin Nutr (2008) 87(4):1080S–6S.
28. Wang TJ, Pencina MJ, Booth SL, Jacques PF, Ingelsson E, Lanier K, et al. Vitamin D deficiency and risk of cardiovascular disease. Circulation (2008) 117(4):503–11. doi:10.1161/CIRCULATIONAHA.107.706127
29. Evans RM, Mangelsdorf DJ. Nuclear receptors, RXR, and the big bang. Cell (2014) 157(1):255–66. doi:10.1186/1471-2148-8-103
30. Haussler MR, Haussler CA, Bartik L, Whitfield GK, Hsieh J-C, Slater S, et al. Vitamin D receptor: molecular signaling and actions of nutritional ligands in disease prevention. Nutr Rev (2008) 66(Suppl 2):S98–112. doi:10.1111/j.1753-4887.2008.00093.x
31. Bouillon R, Suda T. Vitamin D: calcium and bone homeostasis during evolution. Bonekey Rep (2014) 3:1–10. doi:10.1038/bonekey.2013.214
32. Krasowski MD, Ni A, Hagey LR, Ekins S. Evolution of promiscuous nuclear hormone receptors: LXR, FXR, VDR, PXR, and CAR. Mol Cell Endocrinol (2011) 334(1):39–48. doi:10.1016/j.mce.2010.06.016
34. Whitfield GK, Dang HT, Schluter SF, Bernstein RM, Bunag T, Manzon LA, et al. Cloning of a functional vitamin D receptor from the lamprey (Petromyzon marinus), an ancient vertebrate lacking a calcified skeleton and teeth. Endocrinology (2003) 144(6):2704–16. doi:10.1210/en.2002-221101
35. Fidler AE, Holland PT, Reschly EJ, Ekins S, Krasowski MD. Activation of a tunicate (Ciona intestinalis) xenobiotic receptor orthologue by both natural toxins and synthetic toxicants. Toxicon (2012) 59(2):365–72. doi:10.1016/j.toxicon.2011.12.008
36. Markov GV, Laudet V. Origin and evolution of the ligand-binding ability of nuclear receptors. Mol Cell Endocrinol (2011) 334(1):21–30. doi:10.1016/j.mce.2010.10.017
37. King-Jones K, Horner MA, Lam G, Thummel CS. The DHR96 nuclear receptor regulates xenobiotic responses in Drosophila. Cell Metab (2006) 4(1):37–48. doi:10.1016/j.cmet.2006.06.006
38. Lindblom TH, Pierce GJ, Sluder AE. A C. elegans orphan nuclear receptor contributes to xenobiotic resistance. Curr Biol (2001) 11(11):864–8. doi:10.1016/S0960-9822(01)00236-6
39. Delsuc F, Brinkmann H, Chourrout D, Philippe H. Tunicates and not cephalochordates are the closest living relatives of vertebrates. Nature (2006) 439(7079):965–8. doi:10.1038/nature04336
40. Reschly EJ, Bainy CDA, Mattos JJ, Hagey LR, Bahary N, Mada SR, et al. Functional evolution of the vitamin D and pregnane X receptors. BMC Evol Biol (2007) 7(1):1. doi:10.1186/1471-2148-7-222
41. Ekins S, Reschly EJ, Hagey LR, Krasowski MD. Evolution of pharmacologic specificity in the pregnane X receptor. BMC Evol Biol (2008) 8(1):1. doi:10.1186/1471-2148-8-103
42. Cruzeiro C, Lopes-Marques M, Ruivo R, Rodrigues-Oliveira N, Santos MM, Rocha MJ, et al. A mollusk VDR/PXR/CAR-like (NR1J) nuclear receptor provides insight into ancient detoxification mechanisms. Aquat Toxicol (2016) 174:61–9. doi:10.1016/j.aquatox.2016.02.007
43. Hang S, Purdy AE, Robins WP, Wang Z, Mandal M, Chang S, et al. The acetate switch of an intestinal pathogen disrupts host insulin signaling and lipid metabolism. Cell Host Microbe (2014) 16(5):592–604. doi:10.1016/j.chom.2014.10.006
44. Afschar S, Toivonen JM, Hoffmann JM, Tain LS, Wieser D, Finlayson AJ, et al. Nuclear hormone receptor DHR96 mediates the resistance to xenobiotics but not the increased lifespan of insulin-mutant Drosophila. Proc Natl Acad Sci U S A (2016) 113(5):1321–6. doi:10.1073/pnas.1515137113
45. McElwee JJ, Schuster E, Blanc E, Piper MD, Thomas JH, Patel DS, et al. Evolutionary conservation of regulated longevity assurance mechanisms. Genome Biol (2007) 8(7):1. doi:10.1186/gb-2007-8-7-r132
46. McEwan DL, Kirienko NV, Ausubel FM. Host translational inhibition by Pseudomonas aeruginosa exotoxin A triggers an immune response in Caenorhabditis elegans. Cell Host Microbe (2012) 11(4):364–74. doi:10.1016/j.chom.2012.02.007
47. Melo JA, Ruvkun G. Inactivation of conserved C. elegans genes engages pathogen- and xenobiotic-associated defenses. Cell (2012) 149(2):452–66. doi:10.1016/j.cell.2012.02.050
48. Okamura WH, Midland MM, Hammond MW, Rahman NA, Dormanen MC, Nemere I, et al. Chemistry and conformation of vitamin D molecules. J Steroid Biochem Mol Biol (1995) 53(1):603–13. doi:10.1016/0960-0760(95)00107-B
49. Zaccone G, Howie A, Mauceri A, Fasulo S, Lo CP, Youson J. Distribution patterns of cytokeratins in epidermis and horny teeth of the adult sea lamprey, Petromyzon marinus. Folia Histochem Cytobiol (1994) 33(2):69–75.
50. Hagey LR, Møller PR, Hofmann AF, Krasowski MD. Diversity of bile salts in fish and amphibians: evolution of a complex biochemical pathway. Physiol Biochem Zool (2010) 83(2):308. doi:10.1086/649966
51. Krasowski MD, Ai N, Hagey LR, Kollitz EM, Kullman SW, Reschly EJ, et al. The evolution of farnesoid X, vitamin D, and pregnane X receptors: insights from the green-spotted pufferfish (Tetraodon nigriviridis) and other non-mammalian species. BMC Biochem (2011) 12(1):1. doi:10.1186/1471-2091-12-5
52. Groen AK, Kuipers F. Bile acid look-alike controls life span in C. elegans. Cell Metab (2013) 18(2):151–2. doi:10.1016/j.cmet.2013.07.009
53. Magner DB, Wollam J, Shen Y, Hoppe C, Li D, Latza C, et al. The NHR-8 nuclear receptor regulates cholesterol and bile acid homeostasis in C. elegans. Cell Metab (2013) 18(2):212–24. doi:10.1016/j.cmet.2013.07.007
54. De Pablo MA, De Cienfuegos GÁ. Modulatory effects of dietary lipids on immune system functions. Immunol Cell Biol (2000) 78(1):31–9. doi:10.1046/j.1440-1711.2000.00875.x
55. Zhao E, Xu H, Wang L, Kryczek I, Wu K, Hu Y, et al. Bone marrow and the control of immunity. Cell Mol Immunol (2012) 9(1):11–9. doi:10.1038/cmi.2011.47
56. Wong KE, Szeto FL, Zhang W, Ye H, Kong J, Zhang Z, et al. Involvement of the vitamin D receptor in energy metabolism: regulation of uncoupling proteins. Am J Physiol Endocrinol Metab (2009) 296(4):E820–8. doi:10.1152/ajpendo.90763.2008
57. Wong KE, Kong J, Zhang W, Szeto FL, Ye H, Deb DK, et al. Targeted expression of human vitamin D receptor in adipocytes decreases energy expenditure and induces obesity in mice. J Biol Chem (2011) 286(39):33804–10. doi:10.1074/jbc.M111.257568
58. Craig TA, Zhang Y, McNulty MS, Middha S, Ketha H, Singh RJ, et al. Research resource: whole transcriptome RNA sequencing detects multiple 1α,25-dihydroxyvitamin D3-sensitive metabolic pathways in developing zebrafish. Mol Endocrinol (2012) 26(9):1630–42. doi:10.1210/me.2012-1113
59. Sakuma T, Miyamoto T, Jiang W, Kakizawa T, Nishio S-I, Suzuki S, et al. Inhibition of peroxisome proliferator-activated receptor α signaling by vitamin D receptor. Biochem Biophys Res Commun (2003) 312(2):513–9. doi:10.1016/j.bbrc.2003.10.131
60. Narvaez C, Simmons K, Brunton J, Salinero A, Chittur S, Welsh J. Induction of STEAP4 correlates with 1,25-dihydroxyvitamin D3 stimulation of adipogenesis in mesenchymal progenitor cells derived from human adipose tissue. J Cell Physiol (2013) 228(10):2024–36. doi:10.1002/jcp.24371
61. Sertznig P, Seifert M, Tilgen W, Reichrath J. Activation of vitamin D receptor (VDR)-and peroxisome proliferator-activated receptor (PPAR)-signaling pathways through 1,25(OH)2D3 in melanoma cell lines and other skin-derived cell lines. Dermatoendocrinology (2009) 1(4):232–8. doi:10.4161/derm.1.4.9629
62. Wu S, Liao AP, Xia Y, Li YC, Li J-D, Sartor RB, et al. Vitamin D receptor negatively regulates bacterial-stimulated NF-κB activity in intestine. Am J Pathol (2010) 177(2):686–97. doi:10.2353/ajpath.2010.090998
63. Chen Y, Zhang J, Ge X, Du J, Deb DK, Li YC. Vitamin D receptor inhibits nuclear factor κB activation by interacting with IκB kinase β protein. J Biol Chem (2013) 288(27):19450–8. doi:10.1074/jbc.M113.467670
64. Tse AK-W, Zhu G-Y, Wan C-K, Shen X-L, Yu Z-L, Fong W-F. 1α,25-Dihydroxyvitamin D3 inhibits transcriptional potential of nuclear factor kappa B in breast cancer cells. Mol Immunol (2010) 47(9):1728–38. doi:10.1016/j.molimm.2010.03.004
65. Wu S, Xia Y, Liu X, Sun J. Vitamin D receptor deletion leads to reduced level of IκBα protein through protein translation, protein–protein interaction, and post-translational modification. Int J Biochem Cell Biol (2010) 42(2):329–36. doi:10.1016/j.biocel.2009.11.012
66. Zhang Y, Leung DY, Richers BN, Liu Y, Remigio LK, Riches DW, et al. Vitamin D inhibits monocyte/macrophage proinflammatory cytokine production by targeting MAPK phosphatase-1. J Immunol (2012) 188(5):2127–35. doi:10.4049/jimmunol.1102412
67. Boland R, Buitrago C, De Boland AR. Modulation of tyrosine phosphorylation signalling pathways by 1α,25(OH)2-vitamin D3. Trends Endocrinol Metab (2005) 16(6):280–7. doi:10.1016/j.tem.2005.06.003
68. Wöbke TK, Sorg BL, Steinhilber D. Vitamin D in inflammatory diseases. Front Physiol (2014) 5:244. doi:10.3389/fphys.2014.00244
69. Kim CS, Joo SY, Lee KE, Choi JS, Bae EH, Ma SK, et al. Paricalcitol attenuates 4-hydroxy-2-hexenal-induced inflammation and epithelial-mesenchymal transition in human renal proximal tubular epithelial cells. PLoS One (2013) 8(5):e63186. doi:10.1371/journal.pone.0063186
70. Ito I, Waku T, Aoki M, Abe R, Nagai Y, Watanabe T, et al. A nonclassical vitamin D receptor pathway suppresses renal fibrosis. J Clin Invest (2013) 123(11):4579–94. doi:10.1172/JCI67804
71. Halder SK, Goodwin JS, Al-Hendy A. 1,25-Dihydroxyvitamin D3 reduces TGF-β3-induced fibrosis-related gene expression in human uterine leiomyoma cells. J Clin Endocrinol Metab (2011) 96(4):E754–62. doi:10.1210/jc.2010-2131
72. Moreno J, Krishnan AV, Swami S, Nonn L, Peehl DM, Feldman D. Regulation of prostaglandin metabolism by calcitriol attenuates growth stimulation in prostate cancer cells. Cancer Res (2005) 65(17):7917–25. doi:10.1158/0008-5472
73. Brungs M, Radmark O, Samuelsson B, Steinhilber D. On the induction of 5-lipoxygenase expression and activity in HL-60 cells: effects of vitamin D3, retinoic acid, DMSO and TGFβ. Biochem Biophys Res Commun (1994) 205(3):1572–80. doi:10.1006/bbrc.1994.2846
74. Takeuchi A, Reddy GS, Kobayashi T, Okano T, Park J, Sharma S. Nuclear factor of activated T cells (NFAT) as a molecular target for 1α,25-dihydroxyvitamin D3-mediated effects. J Immunol (1998) 160(1):209–18.
75. Aranow C. Vitamin D and the immune system. J Investig Med (2011) 59(6):881–6. doi:10.2310/JIM.0b013e31821b8755
76. Gombart AF, Saito T, Koeffler HP. Exaptation of an ancient Alu short interspersed element provides a highly conserved vitamin D-mediated innate immune response in humans and primates. BMC Genomics (2009) 10(1):321. doi:10.1186/1471-2164-10-321
77. Gombart AF, Borregaard N, Koeffler HP. Human cathelicidin antimicrobial peptide (CAMP) gene is a direct target of the vitamin D receptor and is strongly up-regulated in myeloid cells by 1,25-dihydroxyvitamin D3. FASEB J (2005) 19(9):1067–77. doi:10.1096/fj.04-3284com
78. Liu PT, Schenk M, Walker VP, Dempsey PW, Kanchanapoomi M, Wheelwright M, et al. Convergence of IL-1β and VDR activation pathways in human TLR2/1-induced antimicrobial responses. PLoS One (2009) 4(6):e5810. doi:10.1371/journal.pone.0005810
79. Park K, Elias PM, Oda Y, Mackenzie D, Mauro T, Holleran WM, et al. Regulation of cathelicidin antimicrobial peptide expression by an endoplasmic reticulum (ER) stress signaling, vitamin D receptor-independent pathway. J Biol Chem (2011) 286(39):34121–30. doi:10.1074/jbc.M111.250431
80. Martineau AR, Honecker FU, Wilkinson RJ, Griffiths CJ. Vitamin D in the treatment of pulmonary tuberculosis. J Steroid Biochem Mol Biol (2007) 103(3):793–8. doi:10.1016/j.jsbmb.2006.12.052
81. Yenamandra SP, Hellman U, Kempkes B, Darekar SD, Petermann S, Sculley T, et al. Epstein–Barr virus encoded EBNA-3 binds to vitamin D receptor and blocks activation of its target genes. Cell Mol Life Sci (2010) 67(24):4249–56. doi:10.1007/s00018-010-0441-4
82. Hermant P, Michiels T. Interferon-λ in the context of viral infections: production, response and therapeutic implications. J Innate Immun (2014) 6(5):563–74. doi:10.1159/000360084
83. Xu L, Yang L, Liu W. Distinct evolution process among type I interferon in mammals. Protein Cell (2013) 4(5):383–92. doi:10.1007/s13238-013-3021-1
84. Robertsen B. The interferon system of teleost fish. Fish Shellfish Immunol (2006) 20(2):172–91. doi:10.1016/j.fsi.2005.01.010
85. Randall RE, Goodbourn S. Interferons and viruses: an interplay between induction, signalling, antiviral responses and virus countermeasures. J Gen Virol (2008) 89(1):1–47. doi:10.1099/vir.0.83391-0
86. Sang Y, Liu Q, Lee J, Ma W, McVey DS, Blecha F. Expansion of amphibian intronless interferons revises the paradigm for interferon evolution and functional diversity. Sci Rep (2016) 6:29072. doi:10.1038/srep29072
87. Langevin C, Aleksejeva E, Passoni G, Palha N, Levraud J-P, Boudinot P. The antiviral innate immune response in fish: evolution and conservation of the IFN system. J Mol Biol (2013) 425(24):4904–20. doi:10.1016/j.jmb.2013.09.033
88. Kasamatsu J, Oshiumi H, Matsumoto M, Kasahara M, Seya T. Phylogenetic and expression analysis of lamprey toll-like receptors. Dev Comp Immunol (2010) 34(8):855–65. doi:10.1016/j.dci.2010.03.004
89. Zou J, Tafalla C, Truckle J, Secombes CJ. Identification of a second group of type I IFNs in fish sheds light on IFN evolution in vertebrates. J Immunol (2007) 179(6):3859–71. doi:10.4049/jimmunol.179.6.3859
90. Lange CM, Gouttenoire J, Duong FH, Morikawa K, Heim MH, Moradpour D. Vitamin D receptor and Jak–STAT signaling crosstalk results in calcitriol-mediated increase of hepatocellular response to IFN-α. J Immunol (2014) 192(12):6037–44. doi:10.4049/jimmunol.1302296
91. Herrera G. Vitamin D in massive doses as an adjuvant to the sulfones in the treatment of tuberculoid leprosy. Int J Lepr (1949) 17(1/2):35–42.
92. Roy S, Frodsham A, Saha B, Hazra S, Mascie-Taylor C, Hill A. Association of vitamin D receptor genotype with leprosy type. J Infect Dis (1999) 179(1):187–91. doi:10.1086/314536
93. Liu PT, Stenger S, Li H, Wenzel L, Tan BH, Krutzik SR, et al. Toll-like receptor triggering of a vitamin D-mediated human antimicrobial response. Science (2006) 311(5768):1770–3. doi:10.1126/science.1123933
94. Krutzik SR, Hewison M, Liu PT, Robles JA, Stenger S, Adams JS, et al. IL-15 links TLR2/1-induced macrophage differentiation to the vitamin D-dependent antimicrobial pathway. J Immunol (2008) 181(10):7115–20. doi:10.4049/jimmunol.181.10.7115
95. Montoya D, Cruz D, Teles RM, Lee DJ, Ochoa MT, Krutzik SR, et al. Divergence of macrophage phagocytic and antimicrobial programs in leprosy. Cell Host Microbe (2009) 6(4):343–53. doi:10.1016/j.chom.2009.09.002
96. Rayamajhi M, Humann J, Penheiter K, Andreasen K, Lenz LL. Induction of IFN-αβ enables Listeria monocytogenes to suppress macrophage activation by IFN-γ. J Exp Med (2010) 207(2):327–37. doi:10.1084/jem.20091746
97. Fabri M, Stenger S, Shin D-M, Yuk J-M, Liu PT, Realegeno S, et al. Vitamin D is required for IFN-γ-mediated antimicrobial activity of human macrophages. Sci Transl Med (2011) 3(104):104ra102. doi:10.1126/scitranslmed.3003045
98. Liu PT, Wheelwright M, Teles R, Komisopoulou E, Edfeldt K, Ferguson B, et al. MicroRNA-21 targets the vitamin D-dependent antimicrobial pathway in leprosy. Nat Med (2012) 18(2):267–73. doi:10.1038/nm.2584
99. Teles RM, Graeber TG, Krutzik SR, Montoya D, Schenk M, Lee DJ, et al. Type I interferon suppresses type II interferon-triggered human anti-mycobacterial responses. Science (2013) 339(6126):1448–53. doi:10.1126/science.1233665
100. Mayne CG, Spanier JA, Relland LM, Williams CB, Hayes CE. 1,25-Dihydroxyvitamin D3 acts directly on the T lymphocyte vitamin D receptor to inhibit experimental autoimmune encephalomyelitis. Eur J Immunol (2011) 41(3):822–32. doi:10.1002/eji.201040632
101. Spanier JA, Nashold FE, Olson JK, Hayes CE. The Ifng gene is essential for Vdr gene expression and vitamin D3-mediated reduction of the pathogenic T cell burden in the central nervous system in experimental autoimmune encephalomyelitis, a multiple sclerosis model. J Immunol (2012) 189(6):3188–97. doi:10.4049/jimmunol.1102925
102. Wang Q, Li H, Xie H, Fu M, Guo B, Ding Y, et al. 25-Hydroxyvitamin D3 attenuates experimental periodontitis through downregulation of TLR4 and JAK1/STAT3 signaling in diabetic mice. J Steroid Biochem Mol Biol (2013) 135:43–50. doi:10.1016/j.jsbmb.2013.01.008
103. Kaplan M, Aviram M, Hayek T. Oxidative stress and macrophage foam cell formation during diabetes mellitus-induced atherogenesis: role of insulin therapy. Pharmacol Ther (2012) 136(2):175–85. doi:10.1016/j.pharmthera.2012.08.002
104. Vogel DY, Vereyken EJ, Glim JE, Heijnen PD, Moeton M, van der Valk P, et al. Macrophages in inflammatory multiple sclerosis lesions have an intermediate activation status. J Neuroinflammation (2013) 10(1):1. doi:10.1186/1742-2094-10-35
105. Dusso A, Kamimura S, Gallieni M, Zhong M, Negrea L, Shapiro S, et al. γ-Interferon-induced resistance to 1, 25-(OH) 2 D3 in human monocytes and macrophages: a mechanism for the hypercalcemia of various granulomatoses 1. J Clin Endocrinol Metab (1997) 82(7):2222–32. doi:10.1210/jcem.82.7.4074
106. Sharma OP. Hypercalcemia in granulomatous disorders: a clinical review. Curr Opin Pulm Med (2000) 6(5):442–7. doi:10.1097/00063198-200009000-00010
107. Chawla A, Repa JJ, Evans RM, Mangelsdorf DJ. Nuclear receptors and lipid physiology: opening the X-files. Science (2001) 294(5548):1866–70. doi:10.1126/science.294.5548.1866
108. Lin R, White JH. The pleiotropic actions of vitamin D. Bioessays (2004) 26(1):21–8. doi:10.1002/bies.10368
109. Yasmin R, Williams RM, Xu M, Noy N. Nuclear import of the retinoid X receptor, the vitamin D receptor, and their mutual heterodimer. J Biol Chem (2005) 280(48):40152–60. doi:10.1074/jbc.M507708200
110. Adams JS, Modlin RL, Diz MM, Barnes PF. Potentiation of the macrophage 25-hydroxyvitamin D1-hydroxylation reaction by human tuberculous pleural effusion fluid*. J Clin Endocrinol Metab (1989) 69(2):457–60. doi:10.1210/jcem-69-2-457
111. Adams JS, Chen H, Chun R, Ren S, Wu S, Gacad M, et al. Substrate and enzyme trafficking as a means of regulating 1,25-dihydroxyvitamin D synthesis and action: the human innate immune response. J Bone Miner Res (2007) 22(S2):V20–4. doi:10.1359/jbmr.07s214
112. Overbergh L, Stoffels K, Waer M, Verstuyf A, Bouillon R, Mathieu C. Immune regulation of 25-hydroxyvitamin D-1α-hydroxylase in human monocytic THP1 cells: mechanisms of interferon-γ-mediated induction. J Clin Endocrinol Metab (2006) 91(9):3566–74. doi:10.1210/jc.2006-0678
113. Stoffels K, Overbergh L, Giulietti A, Verlinden L, Bouillon R, Mathieu C. Immune regulation of 25-hydroxyvitamin-D3-1α-hydroxylase in human monocytes. J Bone Miner Res (2006) 21(1):37–47. doi:10.1359/JBMR.050908
114. Edfeldt K, Liu PT, Chun R, Fabri M, Schenk M, Wheelwright M, et al. T-cell cytokines differentially control human monocyte antimicrobial responses by regulating vitamin D metabolism. Proc Natl Acad Sci U S A (2010) 107(52):22593–8. doi:10.1073/pnas.1011624108
115. Stojancevic M, Stankov K, Mikov M. The impact of farnesoid X receptor activation on intestinal permeability in inflammatory bowel disease. Can J Gastroenterol Hepatol (2012) 26(9):631–7. doi:10.1155/2012/538452
116. Griffiths WJ, Abdel-Khalik J, Hearn T, Yutuc E, Morgan AH, Wang Y. Current trends in oxysterol research. Biochem Soc Trans (2016) 44(2):652–8. doi:10.1042/BST20150255
117. Luu W, Sharpe LJ, Capell-Hattam I, Gelissen IC, Brown AJ. Oxysterols: old tale, new twists. Annu Rev Pharmacol Toxicol (2016) 56:447–67. doi:10.1146/annurev-pharmtox-010715-103233
118. Pavek P. Pregnane X receptor (PXR)-mediated gene repression and cross-talk of PXR with other nuclear receptors via coactivator interactions. Front Pharmacol (2016) 7:456. doi:10.3389/fphar.2016.00456
119. Shea C, Hough D, Xiao J, Tzertzinis G, Maina CV. An rxr/usp homolog from the parasitic nematode, Dirofilaria immitis. Gene (2004) 324:171–82. doi:10.1016/j.gene.2003.09.032
120. King-Jones K, Thummel CS. Nuclear receptors—a perspective from Drosophila. Nat Rev Genet (2005) 6(4):311–23. doi:10.1038/nrg1581
121. Maeng S, Lee JH, Choi S-C, Kim MA, Shin YK, Sohn YC. The retinoid X receptor in a marine invertebrate chordate: evolutionary insights from urochordates. Gen Comp Endocrinol (2012) 178(2):380–90. doi:10.1016/j.ygcen.2012.06.019
122. Li Z, Wang W, Zhang E, Qiu G. Identification of spliced mRNA isoforms of retinoid X receptor (RXR) in the Oriental freshwater prawn Macrobrachium nipponense. Genet Mol Res (2014) 13(2):3914–26. doi:10.4238/2014.May.23.2
123. Gong J, Huang C, Shu L, Bao C, Huang H, Ye H, et al. The retinoid X receptor from mud crab: new insights into its roles in ovarian development and related signaling pathway. Sci Rep (2016) 6:23654. doi:10.1038/srep23654
124. Allenby G, Bocquel M-T, Saunders M, Kazmer S, Speck J, Rosenberger M, et al. Retinoic acid receptors and retinoid X receptors: interactions with endogenous retinoic acids. Proc Natl Acad Sci U S A (1993) 90(1):30–4. doi:10.1073/pnas.90.1.30
125. De Lera ÁR, Krezel W, Rühl R. An endogenous mammalian retinoid X receptor ligand, at last! ChemMedChem (2016) 11(10):1027–37. doi:10.1002/cmdc.201600105
126. Castillo AI, Sánchez-Martínez R, Moreno JL, Martínez-Iglesias OA, Palacios D, Aranda A. A permissive retinoid X receptor/thyroid hormone receptor heterodimer allows stimulation of prolactin gene transcription by thyroid hormone and 9-cis-retinoic acid. Mol Cell Biol (2004) 24(2):502–13. doi:10.1128/MCB.24.2.502-513.2004
127. Pérez E, Bourguet W, Gronemeyer H, de Lera AR. Modulation of RXR function through ligand design. Biochim Biophys Acta (2012) 1821(1):57–69. doi:10.1016/j.bbalip.2011.04.003
128. Reschly EJ, Ai N, Welsh WJ, Ekins S, Hagey LR, Krasowski MD. Ligand specificity and evolution of liver X receptors. J Steroid Biochem Mol Biol (2008) 110(1):83–94. doi:10.1016/j.jsbmb.2008.02.007
129. Beckstead RB, Lam G, Thummel CS. The genomic response to 20-hydroxyecdysone at the onset of Drosophila metamorphosis. Genome Biol (2005) 6(12):1. doi:10.1186/gb-2005-6-12-r99
130. Wang S, Liu S, Liu H, Wang J, Zhou S, Jiang R-J, et al. 20-hydroxyecdysone reduces insect food consumption resulting in fat body lipolysis during molting and pupation. J Mol Cell Biol (2010) 2(3):128–38. doi:10.1093/jmcb/mjq006
131. Forman BM, Goode E, Chen J, Oro AE, Bradley DJ, Perlmann T, et al. Identification of a nuclear receptor that is activated by farnesol metabolites. Cell (1995) 81(5):687–93. doi:10.1016/0092-8674(95)90530-8
132. Makishima M, Okamoto AY, Repa JJ, Tu H, Learned RM, Luk A, et al. Identification of a nuclear receptor for bile acids. Science (1999) 284(5418):1362–5. doi:10.1126/science.284.5418.1362
133. Parks DJ, Blanchard SG, Bledsoe RK, Chandra G, Consler TG, Kliewer SA, et al. Bile acids: natural ligands for an orphan nuclear receptor. Science (1999) 284(5418):1365–8. doi:10.1126/science.284.5418.1365
134. Wang H, Chen J, Hollister K, Sowers LC, Forman BM. Endogenous bile acids are ligands for the nuclear receptor FXR/BAR. Mol Cell (1999) 3(5):543–53. doi:10.1016/S1097-2765(00)80348-2
135. Kalaany NY, Mangelsdorf DJ. LXRS and FXR: the yin and yang of cholesterol and fat metabolism. Annu Rev Physiol (2006) 68:159–91. doi:10.1146/annurev.physiol.68.033104.152158
136. Hofmann AF, Hagey LR, Krasowski MD. Bile salts of vertebrates: structural variation and possible evolutionary significance. J Lipid Res (2010) 51(2):226–46. doi:10.1194/jlr.R000042
137. Venkatesh M, Mukherjee S, Wang H, Li H, Sun K, Benechet AP, et al. Symbiotic bacterial metabolites regulate gastrointestinal barrier function via the xenobiotic sensor PXR and toll-like receptor 4. Immunity (2014) 41(2):296–310. doi:10.1016/j.immuni.2014.06.014
138. Moura-Alves P, Faé K, Houthuys E, Dorhoi A, Kreuchwig A, Furkert J, et al. AhR sensing of bacterial pigments regulates antibacterial defence. Nature (2014) 512(7515):387–92. doi:10.1038/nature13684
Keywords: sterol, metabolism, vitamin D receptor, cholesterol, xenobiotics, immunity
Citation: Newmark H, Dantoft W and Ghazal P (2017) Evolutionary Origin of the Interferon–Immune Metabolic Axis: The Sterol–Vitamin D Link. Front. Immunol. 8:62. doi: 10.3389/fimmu.2017.00062
Received: 31 October 2016; Accepted: 16 January 2017;
Published: 09 February 2017
Edited by:
Jorg Hermann Fritz, McGill University, CanadaReviewed by:
Andrew John Brown, The University of New South Wales, AustraliaCaroline Pot, Lausanne University Hospital, Switzerland
Copyright: © 2017 Newmark, Dantoft and Ghazal. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Peter Ghazal, p.ghazal@ed.ac.uk