Transcription factor map alignment of promoter regions

We address the problem of comparing and characterizing the promoter regions of genes with similar expression patterns. This remains a challenging problem in sequence analysis, because often the promoter regions of co-expressed genes do not show discernible sequence conservation. In our approach, thus, we have not directly compared the nucleotide sequence of promoters. Instead, we have obtained predictions of transcription factor binding sites, annotated the predicted sites with the labels of the corresponding binding factors, and aligned the resulting sequences of labels—to which we refer here as transcription factor maps (TF-maps). To obtain the global pairwise alignment of two TF-maps, we have adapted an algorithm initially developed to align restriction enzyme maps. We have optimized the parameters of the algorithm in a small, but well-curated, collection of human–mouse orthologous gene pairs. Results in this dataset, as well as in an independent much larger dataset from the CISRED database, indicate that TF-map alignments are able to uncover conserved regulatory elements, which cannot be detected by the typical sequence alignments.

Transcriptional regulation of MDR1, a gene involved in antifungal drug resistance in Candida albicans

ABSTRACT Upregulation of the Major Facilitator transporter gene MDR1 (Multi_drug Resistance 1) is one of the mechanisms observed in Candida albicans clinical isolates developing resistance to azole antifungal agents. To better understand this phenomenon, the cis-acting regulatory elements present in a modulatable reporter system under the control of the MDR1 promoter were characterized. In an azole-susceptible strain, transcription of this reporter is transiently upregulated in response to either benomyl or H2O2, whereas its expression is constitutively high in an azole-resistant strain (FR2). Two cis-acting regulatory elements, that are necessary and sufficient to convey the same transcriptional responses to a heterologous promoter (CDR2), were identified within the MDR1promoter. The first element, called BRE (for Benomyl Response Element, -296 to -260 with respect to the ATG start codon), is required for benomyl-dependent MDR1 upregulation and for constitutive high expression of MDR1 in FR2. The second element, termed HRE (for H2O2 Response Element, -561 to -520), is required for H2O2-dependent MDR1 upregulation, but is dispensable for constitutive high expression. Two potential binding sites (TTAG/CTAA) for the blip transcription factor Cap1p lie within the HRE. Moreover, inactivation of CAP1 abolished the transient response to H2O2 and diminished significantly the transient response to benomyl. Cap1p, which has been previously implicated in cellular responses to oxidative stress, may thus play a transacting and positive regulatory role in benomyl- and H2O2-dependent transcription of MDR1. However, it is not the only transcription factor involved in the response of MDR1 to benomyl. A minimal BRE element (-290 to -273) that is sufficient to detect in vitro sequence-specific binding of protein complexes in crude extracts prepared from C. albicans was also delimited. Genome-wide transcript profiling analyses undertaken with a matched pair of clinical isolates, one of which being azole-resistant and upregulating MDR1, and with an azole-susceptible strain exposed to benomyl, revealed that genes specifically upregulated by benomyl harbour in their promoters Cap1p binding site(s). This strengthened the idea that Cap1p plays a role in benomyl-dependent upregulation of MDR1. BRE-like sequences were also identified in several genes co-regulated with MDR1 in both conditions, which was consistent with the involvement of the BRE in both processes. A set of 147 mutants lacking a single transcription factor gene was next screened for loss of MDR1response to benomyl. Unfortunately, none of the tested mutants showed a loss of benomyl-dependent MDR1 upregulation. Nevertheless, a significant diminution of the response was observed in the mutants in which the MADS-box transcription factor Mcm1p and the C2H2 zinc finger transcription factor orf19.13374p were inactivated, suggesting that Mcm1p and orf19.13374p are involved in MDR1response to benomyl. Interestingly, the BRE contains a perfect match to the binding consensus of Mcm1p, raising the possibility that MDR1may be a direct target of this transcriptional activator. In conclusion, while the identity of the trans-acting factors that bind to the BRE and HRE remains to be confirmed, the tools we have developed during characterization of the cis-acting elements of the MDR1promoter should now serve to elucidate the nature of the components that modulate its activity. RESUME La surexpression du gène MDR1 (pour Résistance Multidrogue 1), qui code pour un transporteur de la famille des Major Facilitators, est l'un des mécanismes observés dans les isolats cliniques de la levure Candida albicans développant une résistance aux agents antifongiques appelés azoles. Pour mieux comprendre ce phénomène, les éléments de régulation agissant en cis dans un système rapporteur modulable sous le contrôle du promoteur MDR1 ont été caractérisés. Dans une souche sensible aux azoles, la transcription de ce rapporteur est transitoirement surélevée en réponse soit au bénomyl soit à l'agent oxydant H2O2, alors que son expression est constitutivement élevée dans une souche résistante aux azoles (souche FR2). Deux éléments de régulation agissant en cis, nécessaires et suffisants pour transmettre les mêmes réponses transcriptionnelles à un promoteur hétérologue (CDR2), ont été identifiés dans le promoteur MDR1. Le premier élément, appelé BRE (pour Elément de Réponse au Bénomyl, de -296 à -260 par rapport au codon d'initiation ATG) est requis pour la surexpression de MDR1dépendante du bénomyl et pour l'expression constitutive de MDR1 dans FR2. Le deuxième élément, appelé HRE (pour Elément de Réponse à l'H2O2, de -561 à -520), est requis pour la surexpression de MDR1 dépendante de l'H2O2, mais n'est pas impliqué dans l'expression constitutive du gène MDR1. Deux sites de fixation potentiels (TTAG/CTAA) pour le facteur de transcription Cap1p ont été identifiés dans l'élément HRE. De plus, l'inactivation de CAP1 abolit la réponse transitoire à l'H2O2 et diminua significativement la réponse transitoire au bénomyl. Cap1p, qui est impliqué dans les réponses de la cellule au stress oxydatif, doit donc jouer un rôle positif en trans dans la surexpression de MDR1 dépendante du bénomyl et de l'H2O2. Cependant, ce n'est pas le seul facteur de transcription impliqué dans la réponse au bénomyl. Un élément BRE d'une longueur minimale (de -290 à -273) a également été défini et est suffisant pour détecter une interaction spécifique in vitro avec des protéines provenant d'extraits bruts de C. albicans. L'analyse du profil de transcription d'une paire d'isolats cliniques comprenant une souche résistante aux azoles surexprimant MDR1, et d'une souche sensible aux azoles exposée au bénomyl, a révélé que les gènes spécifiquement surexprimés par le bénomyl contiennent dans leurs promoteurs un ou plusieurs sites de fixation pour Cap1p. Ceci renforce l'idée que Cap1p joue un rôle dans la surexpression de MDR1dépendante du bénomyl. Une ou deux séquences ressemblant à l'élément BRE ont également été identifiées dans la plupart des gènes corégulés avec MDR1 dans ces deux conditions, ce qui était attendu compte-tenu du rôle joué par cet élément dans les deux processus. Une collection de 147 mutants dans lesquels un seul facteur de transcription est inactivé a été testée pour la perte de réponse au bénomyl de MDR1. Malheureusement, la surexpression de MDR1 dépendante du bénomyl n'a été perdue dans aucun des mutants testés. Néanmoins, une diminution significative de la réponse a été observée chez des mutants dans lesquels le facteur de transcription à MADS-box Mcm1p et le facteur de transcription à doigts de zinc de type C2H2 orf19.13374p ont été inactivés, suggérant que Mcm1p et orf19.13374p sont impliqués dans la réponse de MDR1au bénomyl. Il est intéressant de noter que la BRE contient une séquence qui s'aligne parfaitement avec la séquence consensus du site de fixation de Mcm1p, ce qui soulève la possibilité que MDR1 pourrait être une cible directe de cet activateur transcriptionnel. En conclusion, alors que l'identité des facteurs agissant en trans en se fixant à la BRE et à la HRE reste à être confirmée, les outils que nous avons développés au cours de la caractérisation des éléments agissant en cis sur le promoteur MDR1 peut maintenant servir à élucider la nature des composants modulant son activité.

Evolutionary origin and functions of retrogene introns.

Retroposed genes (retrogenes) originate via the reverse transcription of mature messenger RNAs from parental source genes and are therefore usually devoid of introns. Here, we characterize a particular set of mammalian retrogenes that acquired introns upon their emergence and thus represent rare cases of intron gain in mammals. We find that although a few retrogenes evolved introns in their coding or 3' untranslated regions (untranslated region, UTR), most introns originated together with untranslated exons in the 5' flanking regions of the retrogene insertion site. They emerged either de novo or through fusions with 5' UTR exons of host genes into which the retrogenes inserted. Generally, retrogenes with introns display high transcription levels and show broader spatial expression patterns than other retrogenes. Our experimental expression analyses of individual intron-containing retrogenes show that 5' UTR introns may indeed promote higher expression levels, at least in part through encoded regulatory elements. By contrast, 3' UTR introns may lead to downregulation of expression levels via nonsense-mediated decay mechanisms. Notably, the majority of retrogenes with introns in their 5' flanks depend on distant, sometimes bidirectional CpG dinucleotide-enriched promoters for their expression that may be recruited from other genes in the genomic vicinity. We thus propose a scenario where the acquisition of new 5' exon-intron structures was directly linked to the recruitment of distant promoters by these retrogenes, a process potentially facilitated by the presence of proto-splice sites in the genomic vicinity of retrogene insertion sites. Thus, the primary role and selective benefit of new 5' introns (and UTR exons) was probably initially to span the often substantial distances to potent CpG promoters driving retrogene transcription. Later in evolution, these introns then obtained additional regulatory roles in fine tuning retrogene expression levels. Our study provides novel insights regarding mechanisms underlying the origin of new introns, the evolutionary relevance of intron gain, and the origin of new gene promoters.

Analysis of LEDGF/p75 expression regulation

Background To replicate, retroviruses must insert DNA copies of their RNA genomes into the host genome. This integration process is catalyzed by the viral integrase protein. The site of viral integration has been shown to be non-random and retrovirus-specific. LEDGF/p75, a splice variant encoded by PSIP1 gene and described as a general transcription coactivator, was identified as a tethering factor binding both to chromatin and to lentiviral integrases, thereby affecting integration efficiency as well as integration site selection. LEDGF/p75 is still a poorly characterized protein, and its cellular endogenous function has yet to be fully determined. In order to start unveiling the roles of LEDGF/p75 in the cell, we started to investigate the mechanisms involved in the regulation of LEDGF/p75. Materials and methods To identify PSIP1 minimal promoter and associated regulatory elements, we cloned a region starting 5 kb upstream the transcription start site (TSS, +1 reference position) to the ATG start codon (+816), as well as systematic truncations, in a plasmid containing the firefly luciferase reporter gene. These constructs were co-transfected into HEK293 cells with a plasmid encoding the Renilla luciferase under the pTK promoter as an internal control for transfection efficiency. Both luciferase activities were assessed by luminescence as an indicator of promoter activity. Results Luciferase assays identified regions -76 to +1 and +1 to +94 as two independent minimal promoters showing respectively a 3.7x and 2.3x increase in luciferase activity. These two independent minimal promoters worked synergistically increasing luciferase activity up to 16.3x as compared to background. Moreover, we identified five regulatory blocks which modulated luciferase activity depending on the DNA region tested, three enhancers (- 2007 to -1159, -284 to -171 and +94 to +644) and two silencers (-171 to -76 and +796 to +816). However, the silencing effect of the region -171 to -76 is dependent on the presence of the +94 to +644 region, ruling out the enhancer activity of the latter. Computational analysis of PSIP1 promoter revealed the absence of TATA box and initiator (INR) sequences, classifying this promoter as nonconventional. TATA-less and INR-less promoters are characterized by multiple Sp1 binding sites, involved in the recruitment of the RNA pol II complex. Consistent with this, PSIP1 promoter contains multiple putative Sp1 binding sequences in regions -76 to +1 and +1 to +94.

Disruption of gene expression in hybrids of the fire ants Solenopsis invicta and Solenopsis richteri.

Transcriptome analysis is a powerful tool for unveiling the distribution and magnitude of genetic incompatibilities between hybridizing taxa. The nature of such incompatibilities is closely associated with the evolutionary histories of the parental species and may differ across tissues and between the sexes. In eusocial insects, the presence of castes that experience divergent selection regimes may result in additional distinct patterns of caste-specific hybrid incompatibilities. We analysed levels of expression of >14 000 genes in two life stages of each caste in the fire ants Solenopsis invicta and Solenopsis richteri and in their hybrids. We found strong contributions of both developmental stage and caste to gene expression patterns. In contrast, variability in expression was only weakly associated with taxonomic identity, with hybrid scores falling between those of the two parental species. Hybrid incompatibilities were surprisingly modest, with only 32 genes being mis-expressed, indicating low levels of disruption in gene regulation in hybrids; males and workers each mis-expressed at least seven times as many genes as queens. Interestingly, homologues of many of the mis-expressed genes have been implicated in behavioural variation in Drosophila melanogaster. General expression profiles of hybrids consistently were more similar to those of S. richteri than S. invicta, presumably because S. richteri trans-regulatory elements tend to be dominant and/or because there is an overall bias in the genetic composition of the hybrids towards S. richteri. Altogether, our results suggest that selection acting on each caste may contribute differently to interspecific divergence and speciation in this group of ants.

Lentiviral vectors: a powerful tool to target astrocytes in vivo.

The morphological and functional diversity of astrocytes, and their essential contribution in physiological and pathological conditions, are starting to emerge. However, experimental systems to investigate neuron-glia interactions and develop innovative approaches for the treatment of central nervous system (CNS) disorders are still very limited. Fluorescent reporter genes have been used to visualize populations of astrocytes and produce an atlas of gene expression in the brain. Knock-down or knock-out of astrocytic proteins using transgenesis have also been developed, but these techniques remain complex and time-consuming. Viral vectors have been developed to overexpress or silence genes of interest as they can be used for both in vitro and in vivo studies in adult mammalian species. In most cases, high transduction efficiency and long-term transgene expression are observed in neurons but there is limited expression in astrocytes. Several strategies have been developed to shift the tropism of lentiviral vectors (LV) and allow local and controlled gene expression in glial cells. In this review, we describe how modifications of the interaction between the LV envelope glycoprotein and the surface receptor molecules on target cells, or the integration of cell-specific promoters and miRNA post-transcriptional regulatory elements have been used to selectively express transgenes in astrocytes.

Sense and antisense transcripts in the histone H1 (HIS-1) locus of Leishmania major.

Histone H1 in the parasitic protozoan Leishmania is a developmentally regulated protein encoded by two genes, HIS-1.1 and HIS-1.2. These genes are separated by approximately 20 kb of sequence and are located on the same DNA strand of chromosome 27. When Northern blots of parasite RNA were probed with HIS-1 strand-specific riboprobes, we detected sense and antisense transcripts that were polyadenylated and developmentally regulated. When the HIS-1.2 coding region was replaced with the coding region of the neomycin phosphotransferase gene, antisense transcription of this gene was unaffected, indicating that the regulatory elements controlling antisense transcription were located outside of the HIS-1.2 gene, and that transcription in Leishmania can occur from both DNA strands even in the presence of transcription of a selectable marker in the complementary strand. A search for other antisense transcripts within the HIS-1 locus identified an additional transcript (SC-1) within the intervening HIS-1 sequence, downstream of adenine and thymine-rich sequences. These results show that gene expression in Leishmania is not only regulated polycistronically from the sense strand of genomic DNA, but that the complementary strand of DNA also contains sequences that could drive expression of open reading frames from the antisense strand of DNA. These findings suggest that the parasite has evolved in such a way as to maximise the transcription of its genome, a mechanism that might be important for it to maintain virulence.

Comparison of human chromosome 21 conserved nongenic sequences (CNGs) with the mouse and dog genomes shows that their selective constraint is independent of their genic environment.

The analysis of conservation between the human and mouse genomes resulted in the identification of a large number of conserved nongenic sequences (CNGs). The functional significance of this nongenic conservation remains unknown, however. The availability of the sequence of a third mammalian genome, the dog, allows for a large-scale analysis of evolutionary attributes of CNGs in mammals. We have aligned 1638 previously identified CNGs and 976 conserved exons (CODs) from human chromosome 21 (Hsa21) with their orthologous sequences in mouse and dog. Attributes of selective constraint, such as sequence conservation, clustering, and direction of substitutions were compared between CNGs and CODs, showing a clear distinction between the two classes. We subsequently performed a chromosome-wide analysis of CNGs by correlating selective constraint metrics with their position on the chromosome and relative to their distance from genes. We found that CNGs appear to be randomly arranged in intergenic regions, with no bias to be closer or farther from genes. Moreover, conservation and clustering of substitutions of CNGs appear to be completely independent of their distance from genes. These results suggest that the majority of CNGs are not typical of previously described regulatory elements in terms of their location. We propose models for a global role of CNGs in genome function and regulation, through long-distance cis or trans chromosomal interactions.

KAP1 regulates gene networks controlling T-cell development and responsiveness.

Chromatin remodeling at specific genomic loci controls lymphoid differentiation. Here, we investigated the role played in this process by Kruppel-associated box (KRAB)-associated protein 1 (KAP1), the universal cofactor of KRAB-zinc finger proteins (ZFPs), a tetrapod-restricted family of transcriptional repressors. T-cell-specific Kap1-deleted mice displayed a significant expansion of immature thymocytes, imbalances in CD4(+)/CD8(+) cell ratios, and altered responses to TCR and TGFβ stimulation when compared to littermate KAP1 control mice. Transcriptome and chromatin studies revealed that KAP1 binds T-cell-specific cis-acting regulatory elements marked by the H3K9me3 repressive mark and enriched in Ikaros/NuRD complexes. Also, KAP1 directly controls the expression of several genes involved in TCR and cytokine signaling. Among these, regulation of FoxO1 seems to play a major role in this system. Likely responsible for tethering KAP1 to at least part of its genomic targets, a small number of KRAB-ZFPs are selectively expressed in T-lymphoid cells. These results reveal the so far unsuspected yet important role of KAP1-mediated epigenetic regulation in T-lymphocyte differentiation and activation.

RNA-based gene duplication: mechanistic and evolutionary insights.

Gene copies that stem from the mRNAs of parental source genes have long been viewed as evolutionary dead-ends with little biological relevance. Here we review a range of recent studies that have unveiled a significant number of functional retroposed gene copies in both mammalian and some non-mammalian genomes. These studies have not only revealed previously unknown mechanisms for the emergence of new genes and their functions but have also provided fascinating general insights into molecular and evolutionary processes that have shaped genomes. For example, analyses of chromosomal gene movement patterns via RNA-based gene duplication have shed fresh light on the evolutionary origin and biology of our sex chromosomes.

Detecting long-range chromatin interactions using the chromosome conformation capture sequencing (4C-seq) method.

Eukaryotic transcription is tightly regulated by transcriptional regulatory elements, even though these elements may be located far away from their target genes. It is now widely recognized that these regulatory elements can be brought in close proximity through the formation of chromatin loops, and that these loops are crucial for transcriptional regulation of their target genes. The chromosome conformation capture (3C) technique presents a snapshot of long-range interactions, by fixing physically interacting elements with formaldehyde, digestion of the DNA, and ligation to obtain a library of unique ligation products. Recently, several large-scale modifications to the 3C technique have been presented. Here, we describe chromosome conformation capture sequencing (4C-seq), a high-throughput version of the 3C technique that combines the 3C-on-chip (4C) protocol with next-generation Illumina sequencing. The method is presented for use in mammalian cell lines, but can be adapted to use in mammalian tissues and any other eukaryotic genome.

Neofunctionalization in vertebrates: the example of retinoic acid receptors.

Understanding the role of gene duplications in establishing vertebrate innovations is one of the main challenges of Evo-Devo (evolution of development) studies. Data on evolutionary changes in gene expression (i.e., evolution of transcription factor-cis-regulatory elements relationships) tell only part of the story; protein function, best studied by biochemical and functional assays, can also change. In this study, we have investigated how gene duplication has affected both the expression and the ligand-binding specificity of retinoic acid receptors (RARs), which play a major role in chordate embryonic development. Mammals have three paralogous RAR genes--RAR alpha, beta, and gamma--which resulted from genome duplications at the origin of vertebrates. By using pharmacological ligands selective for specific paralogues, we have studied the ligand-binding capacities of RARs from diverse chordates species. We have found that RAR beta-like binding selectivity is a synapomorphy of all chordate RARs, including a reconstructed synthetic RAR representing the receptor present in the ancestor of chordates. Moreover, comparison of expression patterns of the cephalochordate amphioxus and the vertebrates suggests that, of all the RARs, RAR beta expression has remained most similar to that of the ancestral RAR. On the basis of these results together, we suggest that while RAR beta kept the ancestral RAR role, RAR alpha and RAR gamma diverged both in ligand-binding capacity and in expression patterns. We thus suggest that neofunctionalization occurred at both the expression and the functional levels to shape RAR roles during development in vertebrates.

Manipulating the sensitivity of signal-induced repression: quantification and consequences of altered brinker gradients.

Traditionally, the analysis of gene regulatory regions suffered from the caveat that it was restricted to artificial contexts (e.g. reporter constructs of limited size). With the advent of the BAC recombineering technique, genomic constructs can now be generated to test regulatory elements in their endogenous environment. The expression of the transcriptional repressor brinker (brk) is negatively regulated by Dpp signaling. Repression is mediated by small sequence motifs, the silencer elements (SEs), that are present in multiple copies in the regulatory region of brk. In this work, we manipulated the SEs in the brk locus. We precisely quantified the effects of the individual SEs on the Brk gradient in the wing disc by employing a 1D data extraction method, followed by the quantification of the data with reference to an internal control. We found that mutating the SEs results in an expansion of the brk expression domain. However, even after mutating all predicted SEs, repression could still be observed in regions of maximal Dpp levels. Thus, our data point to the presence of additional, low affinity binding sites in the brk locus.

Restricted transgene expression in the brain with cell-type specific neuronal promoters.

Tissue-targeted expression is of major interest for studying the contribution of cellular subpopulations to neurodegenerative diseases. However, in vivo methods to investigate this issue are limited. Here, we report an analysis of the cell specificity of expression of fluorescent reporter genes driven by six neuronal promoters, with the ubiquitous phosphoglycerate kinase 1 (PGK) promoter used as a reference. Quantitative analysis of AcGFPnuc expression in the striatum and hippocampus of rodents showed that all lentiviral vectors (LV) exhibited a neuronal tropism; however, there was substantial diversity of transcriptional activity and cell-type specificity of expression. The promoters with the highest activity were those of the 67 kDa glutamic acid decarboxylase (GAD67), homeobox Dlx5/6, glutamate receptor 1 (GluR1), and preprotachykinin 1 (Tac1) genes. Neuron-specific enolase (NSE) and dopaminergic receptor 1 (Drd1a) promoters showed weak activity, but the integration of an amplification system into the LV overcame this limitation. In the striatum, the expression profiles of Tac1 and Drd1a were not limited to the striatonigral pathway, whereas in the hippocampus, Drd1a and Dlx5/6 showed the expected restricted pattern of expression. Regulation of the Dlx5/6 promoter was observed in a disease condition, whereas Tac1 activity was unaffected. These vectors provide safe tools that are more selective than others available, for the administration of therapeutic molecules in the central nervous system (CNS). Nevertheless, additional characterization of regulatory elements in neuronal promoters is still required.

Critical role for Ets, AP-1 and GATA-like transcription factors in regulating mouse Toll-like receptor 4 (Tlr4) gene expression.

TLR4 (Toll-like receptor 4) is essential for sensing the endotoxin of Gram-negative bacteria. Mutations or deletion of the TLR4 gene in humans or mice have been associated with altered predisposition to or outcome of Gram-negative sepsis. In the present work, we studied the expression and regulation of the Tlr4 gene of mouse. In vivo, TLR4 levels were higher in macrophages compared with B, T or natural killer cells. High basal TLR4 promoter activity was observed in RAW 264.7, J774 and P388D1 macrophages transfected with a TLR4 promoter reporter vector. Analysis of truncated and mutated promoter constructs identified several positive [two Ets (E twenty-six) and one AP-1 (activator protein-1) sites] and negative (a GATA-like site and an octamer site) regulatory elements within 350 bp upstream of the transcriptional start site. The myeloid and B-cell-specific transcription factor PU.1 bound to the proximal Ets site. In contrast, none among PU.1, Ets-1, Ets-2 and Elk-1, but possibly one member of the ESE (epithelium-specific Ets) subfamily of Ets transcription factors, bound to the distal Ets site, which was indispensable for Tlr4 gene transcription. Endotoxin did not affect macrophage TLR4 promoter activity, but it decreased TLR4 steady-state mRNA levels by increasing the turnover of TLR4 transcripts. TLR4 expression was modestly altered by other pro- and anti-inflammatory stimuli, except for PMA plus ionomycin which strongly increased promoter activity and TLR4 mRNA levels. The mouse and human TLR4 genes were highly conserved. Yet, notable differences exist with respect to the elements implicated in gene regulation, which may account for species differences in terms of tissue expression and modulation by microbial and inflammatory stimuli.