This Article Abstract Full Text (PDF) Supplemental Research Data Alert me when this article is cited Alert me if a correction is posted Citation Map Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Zhang, Z. Articles by Wheeler, D. A. Search for Related Content PubMed PubMed Citation Articles by Zhang, Z. Articles by Wheeler, D. A. Social Bookmarking                   What's this?

Letter

Genomic Analysis of the Nuclear Receptor Family: New Insights Into Structure, Regulation, and Evolution From the Rat Genome

¹ Human Genome Sequencing Center, Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas 77030, USA ² Department of Molecular and Cellular Biology, Baylor College of Medicine, Houston, Texas 77030, USA ³ Huffington Center on Aging, Department of Otolaryngology, Baylor College of Medicine, Houston, Texas 77030, USA ⁴ Department of Microbiology and Molecular Genetics, University of Texas Medical School, Houston, Texas 77225, USA

	ABSTRACT

Top ABSTRACT RESULTS AND DISCUSSION METHODS REFERENCES WEB SITE REFERENCES

 
Completion of the Rattus norvegicus genome sequence enabled a global inventory and analysis of the nuclear receptors (NRs) in three mammalian species. Forty-nine NR members were found in mouse, 48 in human. Forty-seven were found in the rat, with gaps at the locations expected for the other two. Pairwise comparisons of their distribution in rat, mouse, and human identified 11 syntenic NR gene blocks, including three small clusters of two or three closely related genes, each spanning 40 kb to 1700 kb. The exon structure of the ligand-binding domain suggests that exon shuffling has played a role in the evolution of this family. An invariant splice junction in all members of the NR family except LXR suggests a functional role for the intron. The ligand-binding domains of PXR and CAR are among the most divergent in the family. Their higher nucleotide substitution rates may be related to the central role played by these two NRs in the metabolism of the foreign compounds and may have resulted from limited positive selection.

Nuclear receptors (NRs) are transcription factors capable of exerting regulation of gene expression in the nucleus in response to various extracellular and intracellular signals (Tsai and O'Malley 1994; Mangelsdorf et al. 1995). They are activated by binding of small hydrophobic compounds, such as steroids, retinoids, and thyroid hormones. Ligand binding triggers a conformational change in the receptor proteins, which enables an interaction with cofactors and specific cis-regulatory DNA sequences called hormone response elements (HREs) to subsequently modify gene expression. Cognate ligands are not identified for all nuclear receptors. Those that currently lack identified ligand molecules are termed "orphan" NRs (Giguere 1999). Because NRs bind small molecules which can be easily modified by drug design, and regulate a group of diverse and crucial biological functions such as metabolism, homeostasis, development, and disease, they have become promising pharmacological targets.

NRs share a similar modular domain structure, which includes, from N-terminus to C-terminus, the variable modulatory A/B domain, the DNA-binding domain (DBD), the hinge D-region, the ligand-binding domain (LBD), and an F-domain that is not found in all NRs. The DBD contains two zinc fingers in tandem that encompass 80 amino acid residues in total and are directly involved in recognition of the cognate HRE. The LBD harbors a hydrophobic ligand-binding pocket, deep within its core, that is specific to and thus variable among different receptors. The DBD and LBD are the two most conserved domains of NRs and, as a result, are regarded as dual signatures of this protein family.

NRs constitute one of the largest groups of transcription factors in animals. Twenty-one NR genes are identified in the complete sequence of the Drosophila melanogaster genome (Adams et al. 2000), and over 270 are found in Caenorhabditis elegans (Sluder and Maina 2001). The latest estimate of the number in the human genome sequence, based on sequence alignment and phylogenetic analysis, is 49 NR genes and three NR pseudogenes (Robinson-Rechavi et al. 2001). In a detailed study of the evolutionary relationship among NRs (Laudet 1997), the majority of them were assigned to six well defined subfamilies whose inter-relationships remain unresolved. As a result of the work of Laudet, a systematic naming convention was proposed (Nuclear Receptors Committee 1999) including the creation of a new subfamily 0, which consists of the nuclear receptors lacking either the DBD or the LBD.

With the draft rat genome sequence available (Rat Genome Sequencing Project Consortium 2004), it is now possible to conduct a three-way study of the NR genes comparing the human, mouse, and rat. To gain new insights into the structure, regulation, and evolution of this fascinating family we sought to determine their genomic location and their gene structure, and re-evaluate their phylogenetic relationships in Homo sapiens and the two most medically important model systems.

	RESULTS AND DISCUSSION

Top ABSTRACT RESULTS AND DISCUSSION METHODS REFERENCES WEB SITE REFERENCES

 
Nuclear Receptor Inventory in Rat, Mouse, and Human Genomes
The presence of six NR domains was examined in the rat, mouse, and human genomic sequences using GENEWISEDB (see Methods). The numbers of NR domains identified in the three genomes are summarized in Table 1 (see Supplemental Table 1 for a detailed inventory and genomic coordinates). Grouping and subsequent assignment of these domains to different NRs by BLAST revealed that most of the known mammalian NR genes are present in the current three genome sequences (Suppl. Table 1, Fig. 1); however, the sequences encoding several receptors are partially or completely missing in the rat and mouse genomes. The absence of the sequences encoding Rev-erb (NR1D2) and PNR (NR2E3) and the LBD of TLX (NR2E1) in the rat genome and the DBD of LXR (NR1H2) in the mouse genome can be explained by gaps in these two assemblies at the expected syntenic locations. The final tally of complete or partially identified NR genes was 48 for human, 49 for mouse, and 47 for rat.

View this table: [in this window] [in a new window]   Table 1. Numbers of Sequences Encoding Nuclear Receptors Found in Rat, Mouse, and Human Genomesa

View larger version (29K): [in this window] [in a new window]   Figure 1 The chromosomal landscape of rat nuclear receptor genes. NR genes on the forward strand were placed on the right of the chromosomes, and NR genes on the reverse strand were placed on the left. NR1D2, NR2E3, and the sequences encoding the LBDs of NR1B2 and NR2E1 are missing due to sequence gaps in the current rat genome assembly. Their genomic locations are indicated in the square brackets (L for LBD). The syntenic blocks containing NR genes are highlighted in green (see also Table 3).

Although some domain singletons might be a result of false positive identifications, others defy so quick a dismissal and remain puzzling. For example, a 522-bp sequence identified on human chromosome 16 encodes a partial A/B domain of the glucocorticoid receptor (GR, NR3C1), and is 95% identical to a portion of the first coding exon of GR. The observation that the intron downstream of the first coding exon of GR harbors a potentially active family-Y Alu element (Batzer et al. 1990) and that the 522-bp partial copy is immediately surrounded by two Alu elements from the families Y and C, suggested that the creation of this GR domain singleton may be related to the retrotransposition activity of the nearby Alu-Y element.

NR pseudogenes () were identified in each species (Table 2). Our results confirmed the existence of the three known pseudogenes in the human genome including FXR, the only unprocessed pseudogene (Maglich et al. 2001; Robinson-Rechavi et al. 2001). Because the mouse and rat orthologs of FXR are expressed (identified and experimentally proven to be active by one of the authors, J.W., and Otte et al. 2003) and because FXR may share some functions with FXR in cholesterol metabolism, it remains unclear under what circumstances FXR was silenced and how its loss was tolerated and fixed in the ancestral primate population.

View this table: [in this window] [in a new window]   Table 2. Human and Rodent Nuclear Receptor Pseudogenes

Genomic Distribution of Nuclear Receptors
The genomic locations of NRs were mapped onto the rat karyogram (Fig. 1). NR genes were distributed throughout the rat genome except for chromosomes 9, 12, 14, and 17. Although rat chromosome Y was unavailable, no NR genes are found on the human Y chromosome, and none were expected there for rat or mouse. The Poisson test rejected the random distribution (P < 0.001) of NRs in the rat genome. We identified 11 syntenic blocks common to all three genomes; that is, in each block, the same set of NR genes locate on a single chromosome in all three genomes (Table 3). The sizes of these 11 blocks vary from 0.21 Mb to 54.33 Mb. Except for the blocks I, II, and IV, all syntenic blocks have similar sizes in all three genomes. ROR (NR1F3) and FXR (NR1H5) in block IV are less than 9 Mb apart in the rodent genomes; however, they are separated by a 34-Mb interval that includes the centromere in human.

View this table: [in this window] [in a new window]   Table 3. Syntenic Blocks Containing NR Genes in Rat, Mouse, and Human Genomes

The genome sequences bring details of this organization into focus. The gene order, spacing, and orientations are different in the extant clusters i and ii (Fig. 2). Although TR and Rev-erb maintain the same tail-to-tail orientation, the pair is inverted relative to RAR in the two clusters. Among these six genes, only TR has two splice variants with downstream extended 3' coding exons, that is, the ones overlapping Rev-erb (see inset, Fig. 2A). This would suggest that the TR gene structure reflects the ancestral state of TR and therefore recruitment of the terminal exon occurred as a result of the juxtaposition of the two NR genes. It will be interesting to determine whether this is a mammalian invention, as suggested by the negative findings in chicken, or is a general feature of vertebrates.

View larger version (11K): [in this window] [in a new window]   Figure 2 Chromosomal location of related NR gene clusters i and ii. Genes are labeled on each figure, and closely related paralogs are similarly shaded. Coordinate positions for each chromosome are indicated below the number lines. Gene arrow lengths are proportional to the size of each gene. (A) Cluster i spans 270 kb. The inset gives a scale drawing of the relationship of the NR1A1 and 1D1 genes. The three known variants of NR1A1 are shown. Coding exons are shaded boxes, 3' UTRs are open. A filled inverted triangle marks the splice acceptor of the invariant LBD splice junction (see also Fig. 3B). (B) Cluster ii, 1.4 Mb. The rat gene for NR1D2 is only presumed to exist at the indicated position. Sequences for this gene are absent from the assembly, and a gap exists at this position (indicated by the broken line). The rat NR1B2 is a partial gene, containing a DBD but not an LBD, most likely as a result of incomplete assembly of this draft genome. Note that the 1A and 1D genes are on opposite strands in each cluster, in the same orientation relative to each other; their order changes relative to the 1B gene.

Phylogenetic Analysis
The NR DBDs and LBDs were tested separately to reinvestigate the possibility that recombination between the two domains accounted for some of the diversity in the NR family (see Laudet et al. 1992), and to enable investigation of the relationships with the NR0B group. The overall topologies of the two trees gave the expected subfamily and group clades, upon which their systematic nomenclature is based (Laudet 1997). They differed in small details in a way that could possibly be consistent with one or more exchange events between these two domains early in NR history. For example, subfamily NR4 was closer to NR1 in the LBD tree (Fig. 3A) but was closer to NR5 in the DBD tree (Fig. 3C). However, 68% bootstrap support was marginal for the DBD configuration. This issue warrants further investigation.

View larger version (27K): [in this window] [in a new window]   Figure 3 Unrooted phylogenetic trees of the NR family. The same color scheme for NR subfamilies is used as in Fig. 1. Group-level designations (e.g., 0B, 1A, 1B,..., 6A) label the interior branches, but common gene names label the terminal branches. Bootstrap values expressed in percentage are indicated at the nodes (branch bifurcations). (A) A complete tree constructed from the multiple sequence alignment of the LBDs of all NRs found in rat, mouse, and human. Shading highlights groups exhibiting rapid evolution. (B) NR2 subfamily clade (orphan receptors) taken from the DBD tree. Shading highlights a group exhibiting increased conservation. (C) Portion of DBD tree showing the relationship between subfamilies NR4 and NR5.

The LBDs of most NR members have changed little since the divergence of humans and rodents. This is manifested in the tree as extremely short terminal branch lengths, that is, those branches representing the last common ancestor of the three species. However, three groups, (NR1I2-3, NR1H5, and NR0B1-2, see Fig. 3A, shaded groups) were significantly more divergent among the three species. Nucleotide substitution analysis revealed that the synonymous rates in the LBDs of CAR (NR1I3) and PXR (NR1I2) are average for the family, whereas the nonsynonymous rates were 6.4 and 3.7 times higher than the average (Suppl. Fig. 1).

Evaluation of the terminal branch lengths of all NR members revealed cases where the rat sequence was closer to human than the mouse was: RAR (NR1B1), GR (NR3C1), and LXR (NR1H3). For many others members, the human, rat, and mouse were virtually indistinguishable. These observations may be of practical value in choosing model systems for pharmacological studies.

There was too little variation in the 80-aa DBD to form a well resolved tree, so DNA sequences were used for this domain. The terminal branches were again of most interest to the inter-species comparison. Long terminal branches were observed for all but two NR members: COUP-TFI (NR2F1) and COUP-TFII (NR2F2; Fig. 3B, shaded portion). The relative absence of inter-species variation in the DNA encoding these two NR domains suggested the possibility of selection operating on the DNA sequence itself. Conserved regulatory sequences could be one explanation for this observation. It may therefore be significant that the DBD is uninterrupted by introns in these two NRs (see below).

The K_A/K_S ratios of the LBD domains indicate that the NRs are subject to strong purifying selection. No positive selection was detected by Student's t-test. However, because the K_A/K_S ratios of the LBDs of PXR and CAR were 4.0 and 5.6 times greater than the averages, respectively, these two domains may have experienced limited positive selection in the context of the NR evolution. For PXR and CAR, the increased K_A/K_S ratios in the LBDs could be more readily explained by their biological functions. PXR, an orphan NR preferentially expressed in the liver and intestine, responds to potentially harmful chemicals by activating the expression of cytochrome P-450 genes crucial for the detoxification of a wide variety of structurally diverse xenobiotics and endobiotics (Kliewer et al. 1998; Lehmann et al. 1998). The K_A/K_S ratios of the remaining orphans were much more conserved and thus their ligands, if any, are not likely to be species-specific.

We investigated the structural implications of the LBD sequence variation in the PXR group. Thirty-three variable sites in the multiple sequence alignment of the LBDs of PXR from human, mouse, rat, rhesus, pig, rabbit, dog, chicken, and zebrafish were mapped on the tertiary structure of the LBD of the human PXR (Watkins et al. 2003). Seven sites line the inner surface of the ligand-binding pocket (Watkins et al. 2001); eight variable sites were distributed along -helix 9 (9), which is involved in protein–protein interactions; the remaining sites were distributed uniformly throughout the LBD (Fig. 4). The set lining the ligand-binding pocket was in a position that could possibly form direct contacts with the bound ligand and may therefore contribute to the difference between the ligand-binding properties of the human and rodent PXRs. The set distributed along 9 was outwardly oriented.

View larger version (51K): [in this window] [in a new window]   Figure 4 Variable sites in the LBD of PXR. (A) The variable sites, highlighted in gray, in the protein sequence alignment of the LBDs of the human, mouse, and rat PXRs. The corresponding secondary structure is indicated below the sequence alignment: the -helix is represented by the cylinder, and the -sheet by the parallelogram. The seven variable sites at which the amino acid residues line the ligand-binding pocket and the ones found in the -helix 9 and its vicinity are boxed with the solid line and broken line, respectively. Other available LBD sequences of the rhesus, pig, rabbit, dog, chicken, and zebrafish PXRs are omitted from the sequence alignment presented here, because the inclusion of them does not introduce changes to the general variation pattern. (B) The same sites, highlighted in yellow, in the tertiary structure of the LBD of the human PXR (the blue solid ribbon). The agonist is shown as the small red molecular structure bound in the receptor's ligand-binding pocket, and the coactivator in fragment is shown as the green solid ribbon. Details of the variable sites, with the side chains of the amino acid residues at these sites shown in yellow, in the ligand-binding pocket (C) and in the -helix 9 and its vicinity (D) are shown.

Exon Structures of DBD and LBD
Information on the splice junctions, derived from BLAT alignments of amino acid sequences of NR mRNAs to the genome, was used to further characterize the family. All splice junctions were conserved among orthologs of the various NR family members. When comparing paralogous members, informative patterns of conservation emerge within these two domains (Fig. 5).

View larger version (38K): [in this window] [in a new window]   Figure 5 The gene structure encoding the DBD and LBD domains of the NR genes. Open bars are exons, drawn to scale; line segments, drawn at fixed length, give intron locations. (A) DBD splice junctions. Sequences are 75–78 aa in length. The shaded boxes indicate the location of the two C4 zinc finger motifs within this highly conserved domain. Introns may be found at seven different locations in the DBD across the entire family, or may be absent. Vertical hash marks indicate the location of junctions that were shared in the following groupings: a, NR2B, NR2C1,2; b, NR1A,1B,1C,1D,1F,1H4-5, NR5A; c, NR1I, NR4A; d, NR3A,3B,3C; e, NR2E3; f, NR2A, NR2F6; g, NR2E1; and not shown are group h, NR1H, NR2F, and NR6A, which have no intron in the DBD. (B) LBD splice junctions. Sequences are 170 (NR1D2) to 208 (NR0B1) aa in length. Each row is a schematic drawing giving the relative location of the splice junction and the group of NRs sharing the splice junction pattern. The position of splice junctions in orthologs was always the same, and thus species designations are omitted. Two conserved motifs (I and II, see text) in the LBDs are shown as the hatched areas. The location of a highly conserved negatively charged amino acid residue (aspartic acid or glutamic acid) in motif II is marked by an inverted triangle. The four regions within which introns were found are indicated by slash marks: "\" in motif 1, "|" intermotif region, no slash in motif II, and "/" after motif II (see text). (C) The consensus sequences of motifs I and II. The secondary structure of the corresponding part of the LBD, derived from crystallographic studies, is indicated below the sequence. Letters in bold correspond to the residues of the NR signature, involved in stabilizing the canonical fold of the NR LBDs (see Wurtz et al. 1996).

The splice junction was lost from NR1H2&3, NR2F1&2, and NR6A1. Because these do not form a monophyletic group in the tree (Fig. 3A), the intron was probably lost in three separate events. Members of subfamilies NR1 and NR3 show little variation in junction location, whereas subfamily NR2 has several variants. In two cases members of different subfamilies shared the same splice junction: NR1 and NR5, and NR1I and NR4. This result, taken together with the phylogenetic results described above, may suggest a complex evolutionary relationship between the subfamilies NR1, 4, and 5 (see Phylogenetic Analysis above). Alternatively, there could be preferred sites for acquiring introns. Elucidation of the principles governing the dynamics of intron acquisition and change over long evolutionary timescales is needed to understand these relationships.

The LBD was less conserved, overall, than the DBD compared across the whole family (alignments in Suppl. Fig. 2). Within it, three sequence motifs were identified (see Methods), although none of those were as conserved as the zinc finger motifs in the DBD. Motif I, spanning -helices 3 and 4, was previously described (see Wang et al. 1989 and Wurtz et al. 1996 for definition of helices in LBD). Motif II, -helices 7–9, was described in part (Wang et al. 1989 identified conservation in helices 7 and 8). Steroid receptor groups NR3A and NR3C differed from all others in that motif II was not detected. However, sequence similarity to the left half of the consensus pattern (Fig. 5C) was easily observable (note crosshatching of NR3A1-3 and NR3C1-4 in Fig. 5B, and see Suppl. Fig. 2; see also Wang et al. 1989). Motif III spans most of -helices 10 and 11. Motifs I and II are indicated schematically in Figure 5B (motif III is altered or deleted in some NR isoforms, so we have omitted it from the figure pending further analysis of all of the individual splice variants).

Up to four splice junctions were found in the peptide sequences of the region of the LBDs to which our analysis was confined (see the Pfam profiles used to identify LBD, and Methods). The locations of the four introns were confined to distinct regions of the LBD as defined by the aforementioned structural motifs. The first is within motif I; second, between motifs I and II; third, within motif II; and fourth, downstream from motif II. Introns were lost multiple times at regions 1, 2, and 4. Moreover, the precise location of introns 1, 2, and 4 was variable. In distinct contrast, intron 3 was invariant in that it was present in all of the NRs except LXR (NR1H2), and it was always a phase-1 intron at the same amino acid position in motif II. Although conserved in position and phase, this intron varied in size from 123 bp in mouse NR2A2 to 53,000 bp in human NR5A2. Except in TLX and the steroid hormone receptor genes, on the splice acceptor side of the intron, there was a highly conserved aspartic acid (occasionally replaced by glutamic acid) which contributes to the polar interactions involved in the NR dimerization of those members in which it is present.

The conserved LBD splice junction was likely to have originated early in the family and was subsequently conserved in evolution: it was also observed in the LBD of the Danio rerio SVP46 (NR2F5, data not shown). The selective pressure maintaining the splice junction could arise from conservation of amino acid sequence. The aforementioned aspartic acid codon is split by this phase-1 splice junction. However, motif II is much less conserved than the zinc finger motifs or motif I, and some NR sub-families have neither the aspartic acid nor a glutamic acid at the splice junction. Thus, some sequence or structural motif in the NR mRNA involved in its regulation, processing, or stability may be the determinant of the conservation of this splice junction. Because LXR (NR1H2)—as the single exception—lost this splice junction, the comparison of its expression to other NR genes may shed light on this phenomenon.

Comparison of the 26 different splicing patterns of the sequence encoding the LBD (Fig. 5B) conveys the sense that large-scale sequence changes, intron loss or gain, and exon addition or substitution played an important role in shaping the evolution of this family. The loss or gain of the first, second, or fourth introns in the LBD occurred within many NRs. Large-scale innovative changes in the coding sequence of the LBD region may have contributed significantly to the rise of some new NR genes. FXR (NR1H5) appears to have added one exon between the conserved motifs I and II. In the steroid receptors (ER of the group NR3A; GR, MR, PR, and AR of the group NR3C), the helix 9 homology of motif II abruptly disappears beginning with the loss of the conserved aspartic acid on the splice acceptor side of the exon (see above). Although the downstream C-terminal boundary of the loss of homology is difficult to determine, this observation could be neatly explained by substitution of the exon downstream from the conserved third splice junction, an event leading to specificity for steroid ligands in groups NR3A and 3C.

Further variation in LBD splice junction patterns may exist in other isoforms, and thus a full accounting of all isoforms, in these and other species, will be important.

Conclusion
The genomic comparison of the NR families from three related mammals affords new insight and raises new questions about the structure, function, and evolution of this important family of transcription factors. Despite the high degree of conservation among the NR sequences, there was clearly distinguishable species-specific variation in three groups. Among them, PXR and CAR, in group NR1I, share some ligands (Moore et al. 2000) and regulate overlapping but distinct sets of genes involved in xenobiotic detoxification (Maglich et al. 2002). Given the central role of CAR and PXR in the xenobiotic metabolism and ingestion, these two NRs may have evolved faster in response to different sets of environmental challenges encountered by humans, mice, and rats. The nature of the species-specific adaptations represented by the other two rapidly evolving groups, NR1H5 and NR0B1-2, awaits improved understanding of their functional roles.

Paralogous NR family members exhibit a variety of different exon structures in both their DBDs and LBDs. Among the variation, the conserved location of the splice junction in the second motif of the LBD stands out as a peculiar phenomenon. It may prove to be a more reliable signature for the NR genes than the C4 zinc finger. Very similar findings are reported in other gene families, for example, chemoreceptor superfamily (Robertson et al. 2003) and DEAD helicase genes (Boudet et al. 2001). An understanding of the selective constraints that preserve such ancient introns may lead to new understanding of protein or mRNA structure and processing.

	METHODS

Top ABSTRACT RESULTS AND DISCUSSION METHODS REFERENCES WEB SITE REFERENCES

 
Identification of Nuclear Receptor Genes in Human, Mouse, and Rat Genomes
Six structural and functional domains specific for members of the NR family were obtained from Pfam (Bateman et al. 2002). They are the ligand-binding domain (Pfam database entry name: hormone_rec), found in all members of the family, the C4-type zinc finger DNA-binding domain (zf-C4), found in all but two members, and the four modulator A/B domains, each specific for a given steroid receptor: androgen receptor (Androgen_recep), glucocorticoid receptor (GCR), estrogen receptor (Oest_recep), and progesterone receptor (Prog_receptor). The DBD sequence corresponded to a 75–78-residue segment, starting at the location two amino acid residues before the first conserved cysteine, and encompassing both C4 zinc fingers. The LBD began at the twelfth residue of -helix 3 and extended through -helix 10 (Wurtz et al. 1996; Greschik et al. 1999).

The mRNA and protein sequences of 62 representative NRs (Robinson-Rechavi et al. 2001) were downloaded from GenBank. If the human gene sequence of an NR was available, the mouse and rat gene sequences of the same NR were also retrieved, if available. The Pfam domains present in these 62 NRs were identified using HMMPFAM (HMMER 2.3.1; Eddy 1998). Because the E-values of the identification of the NR domains are 10²⁰ 10⁵⁰ times less than those of other domains identified, their identification and presence in NRs were unambiguous.

The human, mouse, and rat genomic sequences used in this study were human genome build 34 of June 2003, mouse genome build of February 2003, and rat genome build of April 2003. To take advantage of parallel computing, each of these three genomes was partitioned into 750-kb segments with 2-kb overlaps. Only domains of the NRs identified at the previous step with stringently high E-values were searched in the genomic sequences using GENEWISEDB (Wise 2.2.0; Birney and Durbin 2000). Usually GENEWISEDB can predict the presence of a domain in a genome based on the domain profile in Pfam without any modifications to the genomic sequence, but occasionally it introduced one or more frameshifts to make sensible prediction alignments. Although GENEWISEDB labeled such predictions as pseudogenes, we treated them with extra care because the necessity of introducing frameshifts may well result from sequencing errors in the genome.

Domains identified in each genome were grouped together based on their orientation and the coordinates of their genomic locations, and were compared to the 62 NR protein sequences using the best BLASTP hit as the identity. The GENEWISEDB search results were also parsed to create custom annotation tracks in the UCSC genome browser (http://genome.ucsc.edu/) to depict the exon-intron structure of the predicted domains and to enable cross-examination with mRNA/EST evidence, synteny, and genomic sequence conservation across species.

Pseudogenes were identified among NR genes which had more than one copy in a genome and when the sequence of the mRNA transcript of this gene or its orthologs was available. The mRNA sequence was aligned using TBLASTN and BLAT (Kent 2002) to the genomic sequences at the locations where the different copies of multiple NR genes were found. A copy of an NR gene in the genome was considered to be a pseudogene if frameshifts or nonsense mutations were found in its sequence which could not be credibly attributed to the sequencing errors.

Statistical Test for Clustering of Nuclear Receptor Genes in the Rat Genome
The spatial distribution of the NRs in the rat genome was tested for clustering by ² (Zar 1984). The rat genome was divided into nonoverlapping 2.25-Mb segments, and the number (X) of NRs was tallied in each segment. The observed frequency (f_o) of X was tallied, and the corresponding expected frequency (f_e) was calculated from the Poisson probability P(X). The ² value was ² = 17.702, degrees of freedom 1. Because ²_{0.001, 1} = 10.828, the random distribution is rejected (P < 0.001).

Sequence Analyses
The peptide sequences of the DBD and the LBD were identified and extracted from genomic sequence using GENWISEDB. Sequences in each set were aligned using CLUSTALW, and the multiple sequence alignments were then inspected and manually refined in BioEdit. The C-terminal 5–10 residues were incorrect in about half of all sequences extracted from the genomes by GENEWISEDB. They were corrected to match the corresponding GenBank sequence. Nucleotide sequences of each domain were aligned in accordance with their corresponding amino acid sequence alignment.

Corrected but unaligned LBD peptide sequences were searched for conserved sequence motifs (http://blocks.fhcrc.org/). BLOCKMAKER returns two sets of motifs generated by complimentary methods of detecting ungapped regions of similarity (Henikoff et al. 1995). Positive identifications required that both methods agreed on the presence and location of the given motif within the sequence. Three motifs were obtained for all sequences except the steroid receptors and TLX (data not shown). Motif I included the "canonical LBD signature" spanning -helices 3–5 (see Fig.4B,C; Wang et al. 1989; Wurtz et al. 1996) of 12 helices in the LBD. Motif II corresponded to a central portion of the LBD spanning helices 7–9, involved in dimerization (see Fig. 4B,C). Conservation in the first half of this domain, up to the aspartic acid (Fig. 4C) was observed by Wang et al. (1989), but others have noted extended conservation (Laudet et al. 1992). Motif III spanned helices 10–11. This region is subject to alternate splicing in some NR genes, so it was set aside pending complete description of the family isoforms.

CLUSTALW correctly aligned the residues corresponding to motifs I and III but not motif II. In particular, the subfamily NR0B alignment was greatly improved using motif II as a guide; minor adjustments were required in some other subfamilies. Phylogenetic tree reconstruction of both the protein and DNA alignments was performed using an implementation of the neighbor-joining method in the PAUP^*4.0 software package (Swofford 2003) together with a bootstrap of 1000 replicates.

K_A/K_S of every orthologous gene pair was calculated as the measure of sequence evolution (Li et al. 1985). Student's t-test was used to detect positive Darwinian selection.

Splice Junction Analysis
Splice junctions in the coding sequences were located using BLAT to match all protein sequences (62 representative members from GenBank described above) to the corresponding genome. The BLAT exon segments were manually aligned in a manner that brought into register the DBD and LBD of each protein from the three genomes using EXCEL. This enabled rapid curation of exons found by BLAT, which included elimination of false positive exons due to such things as single-residue indels, missing small N-terminal exons, and other splice site ambiguities that may have tricked BLAT.

	Acknowledgements

 
We thank Hugh M. Robertson for helpful discussions. A.J.C., R.B.L., and F.A.P. were supported in part by the NIH NURSA orphan receptor program, grant #U19DK62434-01. This work was supported by grants 5-U54 HG02051 (Human/MB) from the NHGRI and 5-U54 HG02345 (Rat) from the NHGRI/NHLBI to R.A.G.

The publication costs of this article were defrayed in part by payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 USC section 1734 solely to indicate this fact.

	Footnotes

 
Article and publication are at http://www.genome.org/cgi/doi/10.1101/gr.2160004.

⁵ Corresponding author.
E-MAIL wheeler{at}bcm.tmc.edu; FAX (713)798-6977.

[Supplemental material is available online at www.genome.org.]

	REFERENCES

Top ABSTRACT RESULTS AND DISCUSSION METHODS REFERENCES WEB SITE REFERENCES

 

Adams, M.D., Celniker, S.E., Holt, R.A., Evans, C.A., Gocayne, J.D., Amanatides, P.G., Scherer, S.E., Li, P.W., Hoskins, R.A., Galle, R.F., et al. 2000. The genome sequence of Drosophila melanogaster. Science 287: 2185-2195.[Abstract/Free Full Text]

Bailey, J.A., Gu, Z., Clark, R.A., Reinert, K., Samonte, R.V., Schwartz, S., Adams, M.D., Myers, E.W., Li, P.W., and Eichler, E.E. 2002. Recent segmental duplications in the human genome. Science 297: 1003-1007.[Abstract/Free Full Text]

Bateman, A., Birney, E., Cerruti, L., Durbin, R., Etwiller, L., Eddy, S.R, Griffiths-Jones, S., Howe, K.L., Marshall, M., and Sonnhammer, E.L. 2002. The Pfam protein families database. Nucleic Acids Res. 30: 276-280.[Abstract/Free Full Text]

Batzer, M.A., Kilroy, G.E., Richard, P.E., Shaikh, T.H., Desselle, T.D., Hoppens, C.L., and Deininger, P.L. 1990. Structure and variability of recently inserted Alu family members. Nucleic Acids Res. 18: 6793-6798.[Abstract/Free Full Text]

Birney, E. and Durbin, R. 2000. Using GeneWise in the Drosophila annotation experiment. Genome Res. 10: 547-548.[Abstract/Free Full Text]

Bonnelye, E., Vanacker, J.M., Desbiens, X., Begue, A., Stehelin, D., and Laudet, V. 1994. Rev-erb , a new member of the nuclear receptor superfamily, is expressed in the nervous system during chicken development. Cell Growth Differ. 5: 1357-1365.[Abstract]

Boudet, N., Aubourg, S., Toffano-Nioche, C., Kreis, M., and Lecharny, A. 2001. Evolution of intron/exon structure of DEAD helicase family genes in Arabidopsis, Caenorhabditis, and Drosophila. Genome Res. 11: 2101-2114.[Abstract/Free Full Text]

Eddy, S. 1998. Profile hidden Markov models. Bioinformatics 14: 755-763.[Abstract/Free Full Text]

Forrest, D., Sjoberg, M., and Vennstrom, B. 1990. Contrasting developmental and tissue-specific expression of and thyroid hormone receptor genes. EMBO J. 9: 1519-1528.[Medline]

Gampe Jr., R.T., Montana, V.G., Lambert, M.H., Miller, A.B., Bledsoe, R.K., Milburn, M.V., Kliewer, S.A., Willson, T.M., and Xu, H.E. 2000. Asymmetry in the PPAR/RXR crystal structure reveals the molecular basis of heterodimerization among nuclear receptors. Mol. Cell 5: 545-555.[CrossRef][Medline]

Giguere, V. 1999. Orphan nuclear receptors: From gene to function. Endocr. Rev. 20: 689-725.[Abstract/Free Full Text]

Greschik, H., Wurtz, J.-M., Hublitz, P., Kohler, F., Moras, D., and Schule, R. 1999. Characterization of the DNA-binding and dimerization properties of the nuclear orphan receptor germ cell nuclear factor. Mol. Cell. Biol. 19: 690-703.[Abstract/Free Full Text]

Guo, W., Burris, T.P., Zhang, Y.H., Huang, B.L., Mason, J., Copeland, K.C., Kupfer, S.R., Pagon, R.A., and McCabe, E.R. 1996. Genomic sequence of the DAX1 gene: An orphan receptor responsible for X-linked adrenal hypoplasia congenital and hypergonadotropic hypogonadism. J. Clin. Endocrinol. Metab. 81: 2481-2486.[Abstract]

Gurates, B., Amsterdam, A., Tamura, M., Yang, S., Zhou, J., Fang, Z., Amin, S., Sebastian, S., and Bulun, S.E. 2003. WT1 and DAX-1 regulate SF-1-mediated human P450arom gene expression in gonadal cells. Mol. Cell Endocrinol. 208: 61-75.[CrossRef][Medline]

Henikoff, S., Henikoff, J.G., Alford, W.J., and Pietrokovski, S. 1995. Automated construction and graphical presentation of protein blocks from unaligned sequences. Gene 163: GC17-26.[CrossRef][Medline]

Johansson, L., Thomsen, J.S., Damdimopoulos, A.E., Spyrou, G., Gustafsson, J.A., and Treuter, E. 2000. The orphan nuclear receptor SHP utilizes conserved LXXLL-related motifs for interactions with ligand-activated estrogen receptors. Mol. Cell. Biol. 20: 1124-1133.[Abstract/Free Full Text]

Kent, W.J. 2002. BLAT: The BLAST-like Alignment Tool. Genome Res. 4: 656-664.

Kliewer, S.A., Moore, J.T., Wade, L., Staudinger, J.L., Watson, M.A., Jones, S.A., McKee, D.D., Oliver, B.B., Willson, T.M., Zetterstrom, R.H., et al. 1998. An orphan nuclear receptor activated by pregnanes defines a novel steroid signaling pathway. Cell 92: 73-82.[CrossRef][Medline]

Koh, Y.S. and Moore, D.S. 1999. Linkage of the nuclear hormone receptor genes NR1D2, THRB, and RARB: Evidence for an ancient, large-scale duplication. Genomics 57: 289-292.[CrossRef][Medline]

Laudet, V. 1997. Evolution of the nuclear receptor superfamily early diversification from an ancestral orphan receptor. J. Mol. Endocrinol. 19: 207-226.[Abstract]

Laudet, V., Hanni, C., Coll, J., Catzeflis, F., and Stehelin, D. 1992. Evolution of the nuclear receptor gene superfamily. EMBO J. 11: 1003-1013.[Medline]

Lazar, M.A., Hodin, R.A., Darling, D.S., and Chin, W.W. 1989. A novel member of the thyroid/steroid hormone receptor family is encoded by the opposite strand of the rat c-erbA transcriptional unit. Mol. Cell. Biol. 9: 1128-1136.[Abstract/Free Full Text]

Lehmann, J.M., McKee, D.D., Watson, M.A., Willson, T.M., Moore, J.T., and Kliewer, S.A. 1998. The human orphan nuclear receptor PXR is activated by compounds that regulate CYP3A4 gene expression and cause drug interactions. J. Clin. Invest. 102: 1016-1023.[Medline]

Li, W.-H., Wu, C.I., and Luo, C.C. 1985. A new method for estimating synonymous and nonsynonymous rates of nucleotide substitution considering the relative likelihood of nucleotide and codon changes. Mol. Biol. Evol. 2: 150-174.[Abstract]

Maglich, J.M., Sluder, A., Guan, X., Shi, Y., McKee, D.D., Carrick, K., Kamdar, K., Willson, T.M., and Moore, J.T. 2001. Comparison of complete nuclear receptor sets from the human, Caenorhabditis elegans and Drosophila genomes. Genome Biol. 2: RESEARCH0029.

Maglich, J.M., Stoltz, C.M., Goodwin, B., Hawkins-Brown, D., Moore, J.T., and Kliewer, S.A. 2002. Nuclear pregnane X receptor and constitutive androstane receptor regulate overlapping but distinct sets of genes involved in xenobiotic detoxification. Mol. Pharmacol. 62: 638-646.[Abstract/Free Full Text]

Mangelsdorf, D.J., Thummel, C., Beato, M., Herrlich, P., Schutz, G., Umesono, K., Blumberg, B., Kastner, P., Mark, M., Chambon, P., et al. 1995. The nuclear receptor superfamily: The second decade. Cell 83: 835-839.[CrossRef][Medline]

Moore, L.B., Parks, D.J., Jones, S.A., Bledsoe, R.K., Consler, T.G., Stimmel, J.B., Goodwin, B., Liddle, C., Blanchard, S.G., Willson, T.M., et al. 2000. Orphan nuclear receptors constitutive androstane receptor and pregnane X receptor share xenobiotic and steroid ligands. J. Biol. Chem. 275: 15122-15127.[Abstract/Free Full Text]

Nuclear Receptors Committee. 1999. A unified nomenclature system for the nuclear receptor subfamily. Cell 97: 1-20.[CrossRef][Medline]

Otte, K., Kranz, H., Kober, I., Thompson, P., Hoefer, M., Haubold, B., Remmel, B., Voss, H., Kaiser, C., Albers, M., et al. 2003. Identification of farnesoid X receptor as a novel mammalian nuclear receptor sensing lanosterol. Mol. Cell Biol. 23: 864-872.[Abstract/Free Full Text]

Rat Genome Sequencing Project Consortium. 2004. Genome sequence of the Brown Norway Rat yields insights into mammalian evolution. Nature (in press).

Robertson, H.M., Warr, C.G., and Carlson, J.R. 2003. Molecular evolution of the insect chemoreceptor superfamily in Drosophila melanogaster. Proc. Natl. Acad. Sci. 100: 14537-14542.[Abstract/Free Full Text]

Robinson-Rechavi, M., Carpentier, A.-S., Duffraisse, M., and Laudet, V. 2001. How many nuclear hormone receptors in the human genome? Trends Genet. 17: 554-556.[CrossRef][Medline]

Sluder, A.E. and Maina, C.V. 2001. Nuclear receptors in nematodes: Themes and variations. Trends Genet. 17: 206-213.[CrossRef][Medline]

Swofford, D.L. 2003. PAUP^*. Phylogenetic Analysis Using Parsimony (^*and Other Methods). Version 4. Sinauer Associates, Sunderland, MA.

Tsai, M.J. and O'Malley, B.W. 1994. Molecular mechanisms of action of steroid/thyroid receptor superfamily members. Annu. Rev. Biochem. 63: 451-486.[CrossRef][Medline]

Tuzun, E., Bailey, J.A., and Eichler, E.E. 2004. Recent segmental duplications in the working draft assembly of the brown Norway rat. Genome Res. (this issue).

Wang, L.H., Tsai, S.Y., Cook, R.G., Beattie, W.G., Tsai, M.J., and O'Malley, B.W. 1989. COUP transcription factor is a member of the steroid receptor superfamily. Nature 340: 163-166.[CrossRef][Medline]

Watkins, R.E., Wisely, G.B., Moore, L.B., Collins, J.L., Lambert, M.H., Williams, S.P., Willson, T.M., Kliewer, S.A., and Redinbo, M.R. 2001. The human nuclear xenobiotic receptor PXR: Structural determinants of directed promiscuity. Science 292: 2329-2333.[Abstract/Free Full Text]

Watkins, R.E., Davis-Searles, P.R., Lambert, M.H., and Redinbo, M.R. 2003. Coactivator binding promotes the specific interaction between ligand and the pregnane X receptor. J. Mol. Biol. 331: 815-828.[CrossRef][Medline]

Wurtz, J.M., Bourguet, W., Renaud, J.P., Vivat, V., Chambon, P., Moras, D., and Gronemeyer, H. 1996. A canonical structure for the ligand-binding domain of nuclear receptors. Nat. Struct. Biol. 3: 87-94.[CrossRef][Medline]

Zar, J.H. 1984. The Poisson distribution and randomness. Biostatistical analysis, 2nd ed. Prentice-Hall, Inc., Englewood Cliffs, NJ.

Zhang, M. and Chiang, J.Y. 2001. Transcriptional regulation of the human sterol 12-hydroxylase gene (CYP8B1): Roles of heaptocyte nuclear factor 4 in mediating bile acid repression. J. Biol. Chem. 276: 41690-41699.[Abstract/Free Full Text]

	WEB SITE REFERENCES

Top ABSTRACT RESULTS AND DISCUSSION METHODS REFERENCES WEB SITE REFERENCES

 

http://genome.ucsc.edu/; UCSC genome bioinformatics.

http://blocks.fhcrc.org/; blocks WWW server.

Received November 12, 2003; accepted in revised format December 10, 2003.
CiteULike    Connotea    Del.icio.us    Digg    Reddit    Technorati    What's this?

This article has been cited by other articles:

				D. J. Johnson, A. Owen, N. Plant, P. G. Bray, and S. A. Ward Drug-Regulated Expression of Plasmodium falciparum P-Glycoprotein Homologue 1: a Putative Role for Nuclear Receptors Antimicrob. Agents Chemother., April 1, 2008; 52(4): 1438 - 1445. [Abstract] [Full Text] [PDF]

				Z. D. Zhang, P. Cayting, G. Weinstock, and M. Gerstein Analysis of Nuclear Receptor Pseudogenes in Vertebrates: How the Silent Tell Their Stories Mol. Biol. Evol., January 1, 2008; 25(1): 131 - 143. [Abstract] [Full Text] [PDF]

				C. D. Fisher, L. M. Augustine, J. M. Maher, D. M. Nelson, A. L. Slitt, C. D. Klaassen, L. D. Lehman-McKeeman, and N. J. Cherrington Induction of Drug-Metabolizing Enzymes by Garlic and Allyl Sulfide Compounds via Activation of Constitutive Androstane Receptor and Nuclear Factor E2-Related Factor 2 Drug Metab. Dispos., June 1, 2007; 35(6): 995 - 1000. [Abstract] [Full Text] [PDF]

				G. Benoit, A. Cooney, V. Giguere, H. Ingraham, M. Lazar, G. Muscat, T. Perlmann, J.-P. Renaud, J. Schwabe, F. Sladek, et al. International Union of Pharmacology. LXVI. Orphan Nuclear Receptors Pharmacol. Rev., December 1, 2006; 58(4): 798 - 836. [Abstract] [Full Text] [PDF]

				Y. Zhao, Z. Yang, J. K. Phelan, D. A. Wheeler, S. Lin, and E. R. B. McCabe Zebrafish dax1 Is Required for Development of the Interrenal Organ, the Adrenal Cortex Equivalent Mol. Endocrinol., November 1, 2006; 20(11): 2630 - 2640. [Abstract] [Full Text] [PDF]

				A. K. Iyer, Y.-H. Zhang, and E. R. B. McCabe Dosage-Sensitive Sex Reversal Adrenal Hypoplasia Congenita Critical Region on the X Chromosome, Gene 1 (DAX1) (NR0B1) and Small Heterodimer Partner (SHP) (NR0B2) Form Homodimers Individually, as Well as DAX1-SHP Heterodimers Mol. Endocrinol., October 1, 2006; 20(10): 2326 - 2342. [Abstract] [Full Text] [PDF]

				A. M. Nedelcu and R. E. Michod The Evolutionary Origin of an Altruistic Gene Mol. Biol. Evol., August 1, 2006; 23(8): 1460 - 1464. [Abstract] [Full Text] [PDF]

				X. Hu and J. W. Funder The Evolution of Mineralocorticoid Receptors Mol. Endocrinol., July 1, 2006; 20(7): 1471 - 1478. [Abstract] [Full Text] [PDF]