<title>Introduction</title> <p>Class I major histocompatibility complex (MHC) molecules are well known for their pivotal role in the recognition of altered self cells (e.g. virus-infected cells) by T cytotoxic (T<sub>C</sub>) cells. In this process, short peptide fragments derived from pathogens within the host cell are bound to class I receptor structures and transported to the cell surface. Here, the paired receptor/antigen complex is recognized by CD8<sup>+</sup> T<sub>C</sub> cells, initiating a sequence of events that ultimately leads to the lysis of the infected host cell ##UREF##0##[1]##, ##REF##9523112##[2]##.</p> <p>Since class I antigen presentation is essential for the cell-mediated clearance of intracellular pathogens, it is not surprising that these molecules are expressed on most somatic cells and in the majority of host tissue types ##REF##2115762##[3]##. An additional characteristic of the genes encoding class I receptors is a high level of polymorphism, which enables the host to recognize and bind a wide array of foreign peptides ##REF##21232350##[4]##. However, as the number of descriptions of vertebrate MHC loci has increased, it has become clear that many class I genes do not possess the features described above.</p> <p>The classification of vertebrate class I loci has been divided into two general categories: classical (or class Ia) and nonclassical (or class Ib). Classical genes are those found within the MHC region, possessing high polymorphism, strong and wide expression, and are involved in the presentation of endogenous antigens to T<sub>C</sub> cells. Classical loci have been well characterized in humans (HLA-A, -B, and -C) and mice (H-2K, D, and L) but are also well described in some other mammals and in fish, and to a lesser to degree in amphibians, birds, and non-avian reptiles.</p> <p>Nonclassical genes can be located inside or outside of the MHC and typically possess little or no polymorphism and weak expression that is often limited to specific tissue types. Class Ib molecules are known to perform a diverse array of functions which include the recognition of antigenic lipids (human and mouse CD1) ##REF##15032598##[5]##, the binding and transportation of classical class I molecules within the host cell (human HLA-E and mouse Qa-1<sup>b</sup>) ##REF##10047540##[6]##, and the targeting of evolutionarily conserved protein epitopes of pathogens (mouse H2-M3) ##REF##17514210##[7]##. There are even nonclassical loci whose functions are not related to immunity. For example, the human HFE molecule is known to play an important role in iron metabolism ##REF##10047540##[6]##.</p> <p>Phylogenetic studies depict an extremely complex evolutionary history of class I genes, where both classical and nonclassical loci exhibit little orthology across mammals as well as among vertebrates in general ##REF##2484936##[8]##. Since particular species, or closely related taxa, often possess exclusive sets of paralogous genes, it appears that class I lineages have undergone repeated, independent expansion and diversification events over the course of vertebrate evolution. This pattern has been shown to be congruent with a birth-and-death model of evolution, where loci are frequently duplicated and lost, even over short timescales ##REF##9223266##[9]##. However, concerted evolution is also thought to contribute to the close relationship of class I genes within species by the homogenization of even divergent lineages through inter-locus gene conversion ##REF##2315309##[10]##, ##UREF##1##[11]##. Although the function of most nonclassical genes is not well understood, their presence within these independently expanded class I clades is a testament to their importance in vertebrate immunity.</p> <p>Since class I evolution is known to be highly erratic, large taxonomic gaps in MHC characterization prevent the examination of the events and processes which have led to class I diversification. Therefore, a comparative phylogenetic approach may serve as an important tool for understanding class I history. This was recently demonstrated in the gray short-tailed opossum (<italic>Monodelphis domestica</italic>) where 11 class I loci were shown to have diverged after the split of marsupials from other mammals, yielding insight into the timing and breadth of class I differentiation within this lineage ##UREF##2##[12]##.</p> <p>Among vertebrates, non-avian reptiles are the most poorly represented taxon in terms of MHC data. For example, the order Squamata, which includes lizards and snakes, contains over 7,000 species that have been independently evolving for over 200 million years ##REF##13678684##[13]##–##REF##12856236##[15]##. Yet, only a single study characterizing full-length class I genes has been published in the last 15 years ##REF##18488213##[16]##. Squamates, together with the only two extant species in the order Sphenodontia, make up the superorder Lepidosauria, which is the other main lineage of amniotes in addition to archosaurs (birds and crocodilians) and mammals. Given the vast differences in the organization of mammalian and avian MHC regions ##UREF##1##[11]##, ##UREF##3##[17]##, the description of squamate class I genes can therefore provide a firmer phylogenetic basis upon which to reconstruct the characteristics of the ancestral amniote MHC, an important launching point for understanding how the mode of MHC evolution differs among major vertebrate lineages.</p> <p>In this study, we use both cDNA and genomic sequence data to characterize a class I gene from the Galápagos marine iguana (<italic>Amblyrhynchus cristatus</italic>), a member of the squamate subfamily Iguaninae. We also present partial fragments of a similar gene from other iguanine species. The marine iguana sequence possesses several hallmarks of a nonclassical locus and can serve as the basis for more detailed functional studies of class I genes in this group. We compare this sequence to the few others available from squamates, including those from a recent study of class Ia genes from marine iguanas and two other iguanines, the Galápagos land iguana (<italic>Conolophus subcristatus</italic>) and the common green iguana (<italic>Iguana iguana</italic>) ##REF##18488213##[16]##. The information provided here will hopefully aid in the collection of additional class I data from squamates and fill in the vast phylogenetic gap that is missing in our understanding of amniote MHC evolution.</p>
<title>Materials and Methods</title> <title>RNA isolation and cDNA synthesis</title> <p>Total RNA was isolated from the blood of a single marine iguana from the island of Santa Cruz, Galápagos, using the Tri Reagent BD kit (Molecular Research Center, Cincinnati, OH, USA) and the protocol provided. Complementary DNA (cDNA) was synthesized using the SuperScript III RT enzyme and reagents provided in the GeneRacer kit (Invitrogen, Carlsbad, CA, USA). The complete kit protocol was followed for generating first-strand cDNA pools with intact 5′ ends except that both random hexamers and the oligo dT primer provided in the kit were used in a single reverse transcription reaction to obtain full-length cDNA fragments.</p> <title>5′ and 3′ RACE of marine iguana cDNA</title> <p>In order to design gene-specific primers for 3′ RACE, a small fragment (201 base pairs [bp]) in the class I α-2 domain was amplified by PCR using a forward primer described in Radtkey et al. ##UREF##4##[18]## and a degenerate primer (MHCI-R3) designed from an alignment of a range of vertebrate class I sequences (see ##TAB##0##Table 1## for primer list and sequences; see ##FIG##0##Figure 1## for primer locations). Sequences derived from the small, amplified fragment were used to design two forward-facing gene-specific primers for use in 3′ RACE (IgNC-F1 and IgNC-F2).</p> <p>For the 3′ RACE, nested PCR was performed on marine iguana first-strand cDNA. The first round of PCR utilized the oligo dT adaptor-specific primer provided in the GeneRacer kit and the primer Ig-NC-F1. The PCR was performed with the following reagents and concentrations: 1U Phusion High Fidelity DNA Polymerase (Finnzymes, Espoo, Finland), 200 µM of each dNTP, 0.5 µM gene-specific primer, 1× Phusion HF Buffer with MgCl<sub>2</sub>. Because the melting temperature of the gene-specific primer was &gt;72°C, a two-step PCR was carried out with an initial denaturation (denat) of 98°C for 30 s followed by 30 cycles of 98°C for 10 s and 72°C for 60 s, and a final extension (ext) of 72°C for 10 min. Nested 3′ RACE was carried out under the same conditions as the initial PCR using the gene-specific primer Ig-NC-F2, the nested oligo dT adaptor primer provided, and 1 µl of PCR product from the initial 3′ RACE.</p> <p>Sequence data derived from 3′ RACE products were then used to design two gene-specific reverse primers (MHC17H3-5R-R1 and MHC17H3-5R-R2) in the 3′ UTR for use in 5′ RACE. In the first round of PCR, an adaptor-specific primer provided by the kit was used in combination with the gene-specific primer MHC17H3-5R-R1. PCR was performed in a 50 µl reaction with the following reagents and concentrations: 2.5 U Platinum <italic>Taq</italic> DNA Polymerase High Fidelity (Invitrogen), 200 µM of each dNTP, 0.2 µM gene-specific primer, 2 mM MgSO<sub>4</sub>, 1× High Fidelity PCR Buffer. Cycling conditions were as follows: initial denat at 94°C for 2 min followed by 5 cycles of 94°C for 30 s and 72°C for 2 min; 5 cycles of 94°C for 30 s and 70°C for 2 min; 23 cycles of 94°C for 30 s, 59°C for 30 s, and 68°C for 2 min; and a final ext for 10 min at 68°C. The nested 5′ RACE step was carried out with the nested 5′ adaptor primer and the primer MHC17H3-5R-R2 under the following conditions: initial denat at 94°C for 2 min; 25 cycles with 30 s of denat at 94°C, 30 s of annealing at 61°C, 2 min ext at 68°C; and a final ext of 68°C for 20 min.</p> <p>The nested RACE products from both the 5′ and 3′ procedures were ligated into the pCR4-TOPO vector, transformed, and grown overnight on LB plates following the manufacturer's protocol for the TOPO TA Cloning Kit for Sequencing (Invitrogen). Approximately ten clones each from 5′ and 3′ RACE were added to 25 µl of water and heated at 94°C for 10 min to lyse cells. 1 µl of this solution was added to a 25 µl PCR reaction containing 1.5 U AmpliTaq Gold DNA Polymerase (Applied Biosystems, Foster City, CA, USA), 200 µM of each dNTP, 1 µM standard T3 and T7 primers, 1.5 mM MgCl<sub>2</sub>, and 1× PCR Buffer without MgCl<sub>2</sub>. Cycling conditions were as follows: initial denat at 94°C for 10 min, followed by 30 cycles with 30 s of denat at 94°C, 30 s of annealing at 55°C, and 90 s ext at 72°C; and a final ext of 72°C for 10 min. PCR products were purified by filtration using either the QIAquick PCR purification kit (Qiagen; Valencia, CA, USA) or Millipore (Billerica, MA, USA) filter plate and sequenced using a 3730 DNA Analyzer (Applied Biosystems).</p> <title>PCR amplification of Amcr-UA from genomic DNA</title> <p>RACE products yielded a complete coding class I sequence designated <italic>Amcr-UA</italic>. (see <xref ref-type="sec" rid="s3">Results</xref>). In order to amplify this gene from genomic DNA, total DNA was extracted from the whole blood of a single marine iguana using standard phenol-chloroform extraction ##UREF##5##[19]##. Primers designed in the 5′ and 3′ UTR of the marine iguana cDNA sequence were used to amplify the corresponding genomic fragment from the total DNA (##TAB##0##Table 1##). This individual was not the same as the one used in RNA extraction and RACE but did come from the same population (Santa Cruz) and is likely to be closely related. PCR conditions were similar to those used in 3′ RACE with the Phusion enzyme, except that the two-step ext time was increased to 4 min. Due to the large size of the resultant fragment (∼8 kb), numerous internal primers (available upon request) had to be designed in addition to the original PCR primers in order to obtain the complete sequence.</p> <title>Population-level data of Amcr-UA</title> <p>PCR targeting a portion of exon 3, which corresponds to the α-2 domain, and intron 3 was performed on genomic DNA from 65 marine iguanas representing three different islands in the Galápagos archipelago: Santa Cruz, N = 31; Isabela, N = 18; Fernandina, N = 16. Exon 3 was targeted because it is known to be highly variable in class Ia genes and contains sites that are involved in antigen recognition ##REF##9928486##[20]##. Therefore, we felt that the presence or absence of nucleotide variation here would be satisfactory for determining whether this gene is polymorphic or not. The forward primer described in Radtkey et al. ##UREF##4##[18]## was used in combination with the primer MHCI-Int-R1 in a 25 µl PCR reaction with the following reagents and conditions: 2.5 U AmpliTaq DNA Polymerase (Applied Biosystems), 200 µM of each dNTP, 1 µM of each primer, 1.5 mM MgCl<sub>2</sub>, and 1× PCR Buffer without MgCl<sub>2</sub>. Cycling conditions were as follows: initial denat at 94°C for 5 min, followed by 33 cycles with 30 s of denat at 94°C, 30 s of annealing at 59°C, and 60 s ext at 72°C; and a final ext of 72°C for 10 min. PCR products were purified and sequenced as described above.</p> <title>Iguaninae species data collection</title> <p>A genomic fragment homologous to the <italic>Amcr-UA</italic> sequence was targeted in six other species in the same subfamily as <italic>Amblyrhynchus</italic> (Iguaninae): <italic>Conolophus subcristatus</italic>, <italic>Cyclura carinata</italic>, <italic>Cyclura cornuta</italic>, <italic>Cyclura rileyi</italic>, <italic>Ctenosaura defensor</italic>, and <italic>Ctenosaura clarki</italic> (for all table and figures, taxa are labeled with the following abbreviations: <italic>Amcr = Amblyrhynchus cristatus; Cosu = Conolophus subcristatus; Cyco = Cyclura cornuta; Cyca = Cyclura carinata; Cyri = Cyclura rileyi; Ctde = Ctenosaura defensor; Ctcl = Ctenosaura clarki</italic>). This fragment corresponds to the majority of exon 3, the entire intron 3, and a small portion of exon 4 in the <italic>Amcr-UA</italic> sequence, and was amplified in two separate but overlapping fragments. The first fragment, which spans exon 3 and the first part of intron 3, is identical to the one amplified in the population sample of <italic>Amblyrhynchus</italic>, and was generated using the same PCR protocol described above. The second fragment covers the second part of intron 3 and the very beginning of exon 4. The forward primer used to amplify this fragment was the same for all taxa (MHC1-Int-F2), but a different reverse primer (MHC1-R4) was used to amplify <italic>Amblyrhynchus</italic> and <italic>Conolophus</italic> specimens than for the other species (MHC1-R6). The PCR reagent concentrations for both primer pairs were the same as for the exon 3/intron 3 fragment, but the cycling conditions for the MHC1-Int-F2/MHC1-R4 primer combination were as follows: initial denat at 94°C for 5 min, followed by 35 cycles with 40 s of denat at 94°C, 40 s of annealing at 60°C, and 2 min 30 s ext at 72°C; and a final ext of 72°C for 15 min. The PCR cycling conditions for the MHC1-Int-F2/MHC1-R6 primer combination were 94°C for 5 min for initial denat, followed by 35 cycles with 45 s of denat at 94°C, 45 s of annealing at 55°C, and 60 s ext at 72°C; and a final ext of 72°C for 10 min.</p> <title>Data analysis</title> <p>Sequences from the <italic>Amcr-UA</italic> genomic fragment were aligned in the program SEQUENCHER 4.2.2 (Gene Codes Corporation, Ann Arbor, MI, USA). Comparison of the genomic and cDNA sequences was used to identify the exon/intron structure of <italic>Amcr-UA</italic>.</p> <p>The program MUSCLE v3.6 ##REF##15034147##[21]## was used to produce full-length amino acid and exon 4 (α-3 domain) nucleotide alignments for <italic>Amcr-UA</italic> and the following vertebrate class I sequences from GenBank: Galápagos marine iguana (<italic>A. cristatus</italic>), <italic>Amcr-UB*01</italic> EU604308, <italic>Amcr-UB*02</italic> EU604309, <italic>Amcr-UB*03</italic> EU604310, <italic>Amcr-UB*0401</italic> EU604311, <italic>Amcr-UB*0402</italic> EU604312; Galápagos land iguana (<italic>C. subcristatus</italic>), <italic>Cosu-UB*0101</italic> EU604313, <italic>Cosu-UB*0102</italic> EU604314, <italic>Cosu-UB*02</italic> EU604315, <italic>Cosu-UB*03</italic> EU604316; Green iguana (<italic>Iguana iguana</italic>), <italic>Igig-UB*0101</italic> EU604317, <italic>Igig-UB*0102</italic> EU604318, <italic>Igig-UB*02</italic> EU604319; Ameiva lizard, <italic>LC5</italic> M81095, <italic>LC25</italic> M91097; Northern water snake (<italic>Nerodia sipedon</italic>), <italic>SC1</italic> M81099; Chinese soft-shelled turtle (<italic>Pelodiscus sinensis</italic>), AB185243; Chicken (<italic>Gallus gallus</italic>), <italic>B-F10</italic> X12780; Mallard (<italic>Anas platyrhynchos</italic>), <italic>Du2</italic> AB115242; Great reed warbler (<italic>Acrocephalus arundinaceus</italic>), <italic>cN3</italic> AJ005503; Axolotl (<italic>Ambystoma mexicanum</italic>), <italic>Amme-3</italic> U83137; African clawed frog (<italic>Xenopus laevis</italic>), <italic>UAA-1f</italic> L20733; Mouse (<italic>Mus musculus</italic>), <italic>H2K</italic> L36312, <italic>H2-D1</italic> NM_010380, <italic>H2-Q1</italic> NM_010390, <italic>H2-Q10</italic> NM_010391; Wallaby (<italic>Macropus rufogriseus</italic>), <italic>Maru-UB*01</italic> L04952; Platypus (<italic>Ornithorhynchus anatinus</italic>), <italic>Oran2-1</italic> AY112715; Possum (<italic>Trichosurus vulpecula</italic>), <italic>Trvu-UB</italic> AF359509; Rainbow trout (<italic>Oncorhynchus mykiss</italic>), <italic>Onmy-UBA</italic> AF287487; Zebrafish (<italic>Danio rerio</italic>), <italic>Dare-UBA</italic> NM131471; Human, <italic>HLA-B7</italic> U29057, <italic>HLA-Cw</italic> D50852. This set of sequences is similar to the one used in a study of tuatara class I genes by Miller et al. ##REF##16434436##[22]## as well as the recent study of iguanine class Ia loci ##REF##18488213##[16]##, and was chosen for consistency. For the protein data, percent identity between <italic>Amcr-UA</italic> and other vertebrate class I genes was derived from <italic>p</italic>-distance values calculated separately for each structural domain in the program MEGA 4.0 ##REF##17488738##[23]##. Conserved vertebrate class I amino acid positions were identified following Kaufman et al. ##REF##7654997##[24]##.</p> <p>The exon 4 alignment provided the basis for Bayesian, maximum likelihood (ML), and neighbor-joining (NJ) phylogenetic reconstruction. The program MRMODELTEST v2 ##UREF##6##[25]##, which is based on code from the MODELTEST software ##REF##9918953##[26]##, was used to compare the fit of different nucleotide substitution models to the vertebrate dataset. The general time reversible model (GTR) with additional parameters for gamma distribution and fraction of invariable sites provided the best fit to the data according to both the hierarchical likelihood ratio test and the Akaike information criterion. This model was implemented in a Bayesian framework using the program MRBAYES ##REF##12912839##[27]## as well as in ML reconstruction using the TREEFINDER software ##REF##15222900##[28]##. For Bayesian analysis, the default software settings were used, and the search was run for 2,000,000 generations with the first 10% of parameter samples discarded as burn-in. In the ML analysis, 1,000 bootstrap replicates were run to assess support for specific nodes. NJ search and bootstrap analysis were carried out in PAUP using the GTR substitution model with 500 replicates.</p> <p>Sequence fragments spanning the majority of exon 3 and all of intron 3 in the seven iguanine species were aligned in MEGA, and <italic>p</italic>-distance values were calculated separately for the exon and intron.</p> <p>Sequences generated in this study were deposited in GenBank under the following accession numbers.: EU839663 (full-length <italic>Amcr-UA</italic> cDNA); EU839664 (<italic>Amcr-UA</italic> genomic fragment); EU839665-EU839670 (Iguaninae exon 3/intron 3).</p>
<title>Results</title> <title>Characterization of cDNA sequences from <italic>Amblyrhynchus cristatus</italic>
</title> <p>Multiple clones were sequenced from PCR products generated by both 5′and 3′ RACE of marine iguana cDNA. For 3′ RACE, resulting clones carried two unique class I-like sequences. One of these became the subjected of another study ##REF##18488213##[16]##, while the other was the focus of this paper, and was used to design specific 5′ RACE primers. The sequence of this 3′ RACE clone was identical in the area of overlap with the single class I fragment obtained from 5′ RACE. When aligned, these sequences comprised a 1,294 bp fragment which spanned the complete coding sequence (CDS) as well as the 5′ and 3′ UTRs of a single class I sequence type (##FIG##0##Figure 1##). This gene was labeled <italic>Amcr-UA</italic> based on established nomenclature rules ##REF##2329006##[29]##. Although the 3′ UTR sequence reaches the site of polyadenylation, neither of the canonical polyadenylation signals (AATAAA or ATTAAA) were found. However, there are several alternative motifs at a range of positions in the 3′ UTR that are known to be involved in polyadenylation, including AATACA, AATATA, GATAAA or ATAAA, and AAGAAA, which are located 64, 101, 131, and 150 bp upstream of the poly(A) region respectively (##FIG##1##Figure 2##; see ##REF##10899149##[30]##–##REF##15946807##[34]## for characterization of alternative polyadenylation sites).</p> <p>One of the most striking features of the <italic>Amcr-UA</italic> sequence is the absence of the transmembrane (Tm) and cytoplasmic (Cyt) domains that are characteristic of most classical and nonclassical class I MHC genes. The CDS extends only 11 amino acid residues downstream from the expected end of the α-3 domain. The first three of these residues (EEP) are identical to the first three Tm domain residues of the classical <italic>UB</italic> genes from the three iguanine species as well as the Ameiva lizard (<italic>LC1</italic>) (##FIG##2##Figure 3##); but the other eight residues have no clear homology to any of the other loci. There was no evidence of a longer 3′ RACE fragment that might contain a full-length version of the truncated <italic>Amcr-UA</italic> sequence.</p> <p>While the structure and composition of class I genes are quite variable within and between vertebrate groups, the <italic>Amcr</italic>-<italic>UA</italic> sequence possesses many of the conserved amino acid residues that are characteristic of class I genes, and even some residues that are common in classical class I genes (indicated by asterisks in ##FIG##2##Figure 3##). For example, there are nine class I α chain residues which are highly conserved across vertebrates and are known to be involved in the binding of N- and C-termini of antigenic peptides ##REF##7654997##[24]##, ##REF##7612235##[35]##, ##REF##10380691##[36]##. Seven of these residues are shared between the <italic>Amcr-UA</italic> sequence (positions in <italic>Amcr-UA</italic>: Y9, Y62, Y124, T144, W148, Y160, and Y175) and the common reference type found in the classical Human HLA-A, -B, and -C, molecules. Moreover, one of the two remaining positions (R87) matches the predominant character state in non-mammalian vertebrates, including lizard, fish, birds, tuatara, and amphibians. Several studies have used the presence of these conserved peptide-binding residues in different species to distinguish classical antigen-presenting class I genes from nonclassical class I genes with alternative functions ##REF##10380691##[36]##–##REF##9396864##[38]##.</p> <p>Other features that are shared between the marine iguana sequence and nearly all class I genes of other vertebrates include four cysteine residues at positions 103, 165, 204, and 260, which are involved in intradomain disulfide bridge formation, two salt bridges at positions H5-D31 and H95-D120, and the highly conserved FYP motif at positions 209–211 in the α-3 domain. Class I genes of most species contain a NQS or NQT glycosylation site near the end of the α1 domain; however, the <italic>Amcr-UA</italic> sequence contains a histidine (H) rather than a glutamine (Q) at position 90, which is not the case in other vertebrates, including the classical sequences from the three iguanine species as well as the Ameiva lizard. But this residue is not known to be important for N-linked glycan formation. The amino acid residues which correspond to the CD8 binding site of the HLA-A locus are at positions 224–230 in the <italic>Amcr-UA</italic> fragment, but three of the seven residues are non-polar, which would be unexpected at these sites for a classical class I gene; however this region is known to be highly variable among species, and the coevolution of the CD8 ligand and its class I binding site is not well understood ##REF##7654997##[24]##.</p> <title>Genomic sequence of <italic>Amcr-UA</italic>
</title> <p>Amplification of the <italic>Amcr-UA</italic> gene was carried out on marine iguana genomic DNA using specific primers designed in the 5′ and 3′ UTR of the cDNA sequence. The resulting PCR product was approximately 8 kb and was identical to the cDNA in all areas of overlap, except for a single nucleotide position in the expected signal peptide. Some degree of mismatch between the cDNA and genomic sequences was anticipated since these fragments were derived from different individuals. A small portion (&lt;100 bp) of intron 1 in the amplified genomic fragment could not be sequenced from either direction after several attempts. Secondary structure analysis revealed that the ∼80 bp before and after the unsequenced region are highly complementary (∼80%) and could form a secondary loop. However, even after addition of DMSO, which is known to relax secondary structure during cycle sequencing, there was still a consistent and steep drop-off of chromatogram peaks in this region.</p> <p>The structure of the genomic fragment is shown in ##FIG##0##Figure 1##. This sequence revealed that the 11 amino acid residues (plus stop codon) that extend from the expected end of the α-3 domain are encoded in the same exon as the α-3 domain itself (##FIG##0##Figure 1##). In other vertebrates, separate exons usually code for the Tm/Cyt and α-3 domains. Yet the “EEP” amino acid motif, which marks the beginning of the Tm/Cyt domain in the classical class I sequences from the three iguanine species and the Ameiva lizard, is part of the α-3 encoding exon in <italic>Amcr-UA</italic>, and might be the result of some type of genomic rearrangement.</p> <p>Besides the missing Tm/Cyt domains, the position of exons and introns are the same as in the majority of class I genes, where the signal peptide and each of the three extracellular domains (α-1, α-2, and α-3) are encoded by separate exons with intervening introns. The missing fragment from intron 1 was estimated to be ∼100 bp by comparing the band size to a DNA ladder on an agarose gel.</p> <title>Relationship of <italic>Amcr-UA</italic> to Class I sequences of other vertebrates</title> <p>Over the three extracellular domains, the amino acid identity between the <italic>Amcr-UA</italic> sequence and other published vertebrate class I genes ranges from 33% in the Rainbow trout (<italic>Oncorhynchus mykiss</italic>) to between 56.2% and 60.0% for the iguanine <italic>UB</italic> sequences (##TAB##1##Table 2##). Considered separately, the α-2 and α-3 domains are most similar to the classical iguanine sequences, but the α-1 domain is approximately as close to the tuatara (45.5%–50.0%) and Ameiva lizard (46.7%) as it is to the iguanine <italic>UB</italic> sequences (43.5%–50%).</p> <p>Phylogenetic reconstruction was carried out on exon 4 nucleotide sequences from a similar set of taxa as in Miller et al. ##REF##16434436##[22]##. In this previous study, bootstrap values supporting many of the class I gene clusters among vertebrates were low, as is the case in our study when the <italic>Amcr-UA</italic> and iguanine <italic>UB</italic> fragments are included (##FIG##3##Figure 4##). However, there is strong nodal support for the monophyly of all squamate class I sequences (<italic>A. cristatus</italic>, <italic>C. subcristatus</italic>, <italic>I. iguana</italic>, <italic>A. ameiva</italic>, and <italic>N. sipedon</italic>) included in the analysis, suggesting that these genes derive from a common ancestral locus whose descendant lineages have not been identified in any other vertebrate group including <italic>Sphenodon</italic>, which is also a member of the Lepidosauria superorder. However, within squamates, there is moderate support for the basal position of <italic>Amcr-UA</italic> relative to all other published class I sequences.</p> <title>Polymorphism of <italic>Amcr-UA</italic> in <italic>Amblyrhynchus cristatus</italic> populations</title> <p>The majority of exon 3 (248 bp) and the first 235 bp of intron 3 were amplified from 65 marine iguanas from three separate islands representing distinct parts of the Galápagos archipelago. Direct sequencing of PCR products for 62 of these individuals generated a single sequence with no double peaks, which is identical to the full-length genomic fragment. Due to the limited polymorphism exhibited at the population level, these data are not displayed graphically. But it is worth noting that two individuals possess a double peak at a single nucleotide position in exon 3, where one base matches the main <italic>Amcr-UA</italic> type and the other results in a nonsynonymous amino acid change (from threonine to serine at amino acid position 133). Another individual possesses double peaks at two nucleotide positions, one site in exon 3 and another in intron 3; here also, the change in the exon is nonsynonymous (from isoleucine to valine at amino acid position 168)</p> <title>Patterns of variation between Iguana species</title> <p>Primers designed on the <italic>Amcr</italic>-<italic>UA</italic> sequence successfully amplified the majority of exon 3 (248 bp) and the entire intron 3 (∼1,030 bp) of six other species from the Iguaninae subfamily, including a representative of <italic>Conolophus</italic> (Galápagos land iguana), which is thought to be the sister genus of the monospecific <italic>Amblyhrynchus</italic>
##UREF##7##[39]##, ##REF##12116477##[40]##. No stop codons were found in the amino acid sequence of any of the species examined. Uncorrected pairwise genetic distance (<italic>p-</italic>distance) between species was calculated separately for the exon and intron (##TAB##2##Table 3##). With the exception of <italic>Conolophus</italic> and <italic>Ctenosaura clarki</italic>, the <italic>p</italic>-distance was lower in exon 3 than in intron 3 for all pairwise comparisons, even though the α-2 domain is known to be polymorphic and under strong diversifying selection in classical class I genes ##REF##9928486##[20]##.</p> <p>There were 19 polymorphic nucleotide positions in the alignment of iguanine <italic>UA</italic> exon 3 sequences, but only five of these variable sites contributed to nonsynonymous differences between one or more species (##FIG##4##Figure 5##). As a result, only three of the 82 amino acid positions were polymorphic. While the low divergence of exon 3 compared to intron 3 suggests that the structure of this gene may be conserved, there was insufficient exon variation to carry out <italic>d</italic>
<sub>N</sub>/<italic>d</italic>
<sub>S</sub>-based analyses as a means of detecting purifying selection acting on this locus.</p>
<title>Discussion</title> <p>The basic structure of class I loci is highly conserved across vertebrates. Typically, a single gene encodes a large α chain that contains separate exons for the two membrane-distal peptide-binding domains (α-1 and α-2), a third exon for the membrane-proximal Ig-like domain (α-3), and a number of additional, smaller exons that are responsible for the transmembrane (Tm) and cytoplasmic (Cyt) domains ##UREF##0##[1]##. The Tm/Cyt domains are essential for the standard antigen presentation function of class I molecules because they anchor the receptor structure into the cell membrane of altered self-cells where they can be recognized by CD8<sup>+</sup> T<sub>C</sub> cells ##UREF##0##[1]##, ##REF##9523112##[2]##. However, in the <italic>Amcr-UA</italic> cDNA sequence, the Tm/Cyt domains are largely missing, except for an 11 amino acid extension of the α-3 domain that is followed by a stop codon. The first three residues of this extension are identical to the start of the Tm domain identified in iguanine <italic>UB</italic> and Ameiva lizard class I sequences, but the other eight residues show no apparent match between these sequences.</p> <p>Interestingly, the genomic data did not reveal the presence of a Tm/Cyt coding region in between the unexpected early stop codon and the downstream 3′ UTR sequence, and there is no evidence of a cryptic, hidden, or vestigial splice site in the 3′ end of the gene. Therefore, the lack of these domains in the cDNA is not likely to be a result of incomplete transcription or a splicing event. Rather, it seems apparent that the Tm/Cyt region is simply absent at this locus. One possible explanation for the match between the <italic>Amcr-UA</italic> α-3 extension and the start of the iguanine <italic>UB</italic> and Ameiva lizard Tm sequences (the “EEP” motif) is that the ancestor of <italic>Amcr-UA</italic> underwent a genomic rearrangement in which a small portion of the Tm-containing exon was transferred to the end of the α-3 encoding exon. Characterization and comparison of similar genes across other squamates may reveal the pattern and timing of such a modification. In particular, it would be interesting to obtain full-length sequences of <italic>Amcr-UA</italic> orthologues from the other closely related iguana species used in this study to determine whether they possess the same truncation at the 3′ end of the CDS.</p> <p>However, it is also conceivable that intact exons coding for Tm/Cyt domains do exist downstream of the available genomic sequence and simply are not transcribed, at least not at easily detectable levels in the blood. Therefore, future work should involve the use of large insert vectors to explore adjacent stretches of genomic DNA in order to rule out this possibility. In addition, it is equally important to measure expression levels of Amcr-UA in different tissues, and perform a more exhaustive search for a full-length version of this molecule containing a membrane-spanning domain.</p> <p>Despite the apparent lack of Tm/Cyt domains, which are necessary for classical MHC antigen presentation, the results presented here suggest that the <italic>Amcr</italic>-<italic>UA</italic> sequence represents a functional class I gene and not a pseudogene. The presence of the gene in the cDNA pool shows that it is expressed at some level in the blood. In addition, both the genomic and cDNA sequences do not contain any additional stop codons in the CDS, as might be found in a pseudogene. There are also no stop codons in any of the exon 3 <italic>UA</italic> sequences from other iguanine species.</p> <p>The presence of a polyadenylation signal in the 3′ UTR is essential for the transcription of a functional gene. While neither of the canonical polyadenylation signals (AATAAA or ATTAAA) were found in the <italic>Amcr-UA</italic> cDNA sequence, several alternative sites were identified, including ATAAA, which is known from fish ##REF##11717131##[41]##, and AAGAAA, from humans ##REF##10585854##[31]##. Moreover, Beaudoing et al. ##REF##10899149##[30]## showed that the signal sequence and its position from the poly(A) region vary greatly within the human genome.</p> <p>The higher divergence in intron 3 compared to exon 3 among the iguanine <italic>UA</italic> sequences suggests that this gene is conserved and under some degree of purifying selection, or at least is not evolving neutrally as would be expected for a nonfunctional pseudogene. However, other physiological and molecular data, including expression at the protein level, must be collected in order to confirm the functionality of <italic>Amcr-UA</italic>. But for the purpose of discussion, we will proceed with the assumption that this gene is functional in order to explore its characteristics and relationship with other vertebrate class I genes.</p> <p>Despite the absence of a Tm/Cyt region, the <italic>Amcr-UA</italic> sequence shares many of the conserved peptide-binding residues that are suggestive of classical class I function (##FIG##2##Figure 3##). However, the comparison of <italic>UA</italic> exon 3 sequences between Iguaninae taxa shows that there is little adaptive divergence in the antigen-binding region between species, suggesting that this gene is under purifying selection, rather than balancing-selection, and has a conserved function. Thus, it is possible that the product of this gene is involved in antigen binding, but not in the classical sense.</p> <p>Although Kaufman et al. ##REF##7654997##[24]## interpreted the presence of certain conserved residues at peptide-binding sites as preliminary evidence of classical function, many of these amino acid character states are maintained in nonclassical genes such as the human HLA-E, -F, and -G loci, as well as the mouse H2-M3 gene. None of these loci differ from the classical mammalian type by more than three out of the nine residues ##REF##10380691##[36]##. Thus, it is not at all unprecedented that <italic>Amcr-UA</italic> displays several nonclassical features while still showing signs of peptide-binding ability.</p> <p>Since a lack of polymorphism is an important criterion for describing nonclassical loci, sequence data was collected from exon 3 and intron 3 of <italic>Amcr-UA</italic> for three marine iguana populations. The near absence of exon 3 variability in these samples seems to support the pattern of purifying selection indicating a conserved function. However, the simultaneous lack of diversity in the adjacent intron indicates either strong linkage to the conserved exon or perhaps that an insufficient amount of time has passed for a large number of substitutions to accumulate in the intron. The latter pattern would not be surprising since mitochondrial DNA evidence suggests that marine iguana populations are not genetically diverse and may have recently expanded in the Galápagos archipelago ##UREF##8##[42]##. Therefore, the unclear evolutionary history of this species makes it difficult to attribute the dearth of polymorphism in <italic>Amcr-UA</italic> to its nonclassical function.</p> <p>While rare, several truncated class I molecules are known to lack Tm/Cyt regions in humans and mice. For example, the mouse Q10 gene contains a 13 bp deletion in the exon encoding the Tm domain, causing a frame shift and the introduction of a premature termination signal downstream in the same exon. In addition, the remaining transcribed portion of the Tm domain has numerous polar amino acid residues that would likely prevent its insertion into the cell membrane. Therefore, it seems clear that the Q10 molecule is not involved in classical antigen presentation. It also possesses several other nonclassical features, including a lack of polymorphism and expression that is almost exclusive to mouse liver cells ##REF##6956903##[43]##–##REF##6420070##[45]##. Several studies have shown that Q10 is likely secreted and can bind a wide array of non-self peptides in a similar manner to classical molecules, but its function is still not well understood ##REF##6411351##[44]##, ##REF##10657640##[46]##. While there is no phylogenetic similarity between <italic>Amcr-UA</italic> and the mouse Q10, the characteristics of the latter gene demonstrate that truncated, soluble, secreted class I molecules can exist which lack genetic variation and maintain protein binding capabilities.</p> <p>Numerous other classical and nonclassical class I molecules are known to exist in soluble form for secretion. For example, the human HLA-G locus possesses all of the coding features of a membrane-bound class I receptor, but is sometimes subjected to alternative transcription where the Tm/Cyt domains are deleted. This molecule is expressed specifically in placental tissue and is secreted during pregnancy ##REF##10082427##[47]##. The role of this gene is also not well understood, but is thought to be involved in inducing apoptosis in activated maternal CD8<sup>+</sup> T cells ##REF##10843658##[48]##. Classical human class I genes (HLA-A, -B, and -C) are also known to exist in soluble forms and to play a role in cell death of activated T cells ##REF##12594841##[49]##. Nevertheless, a unique feature of the <italic>Amcr-UA</italic> sequence which distinguishes it from the Q10 and HLA-G loci, as well as other soluble class I molecules, is that there does not appear be any sign of a Tm/Cyt coding sequence in the genomic DNA, regardless of whether it is transcribed or not; but again, additional support for this conclusion must come from further collection of sequence data in the 3′ region of the gene.</p> <p>Some truncated class Ib molecules, such as those in the mouse Qa-2 family, are not necessarily secreted, but are rather linked to the cell membrane through a glycosylphosphatidylinositol (GPI) anchor that is added to the carboxyl terminus of the protein during posttranslational modification ##REF##2953793##[50]##, ##REF##8876057##[51]##. Therefore, the apparent lack of a membrane-spanning domain in <italic>Amcr-UA</italic> doesn't rule out its expression on the cell surface or its involvement in the primary T cell response.</p> <p>Overall, the <italic>Amcr-UA</italic> CDS showed the highest similarity with iguanine <italic>UB</italic> and Ameiva lizard sequences. In addition, phylogenetic reconstruction supported the monophyly of all available squamate class I sequences. The basal position of <italic>Amcr-UA</italic> relative to all other squamate class I sequences suggests that this gene diverged very early in the evolution of this reptilian order. The most recent common ancestor of iguanines and the two groups represented by the Ameiva lizard (Family: Teiidae) and northern water snake (Suborder: Serpentes) is estimated to have existed between 179–206 million years ago ##REF##16286089##[14]##, providing a minimum time for the split of <italic>Amcr-UA</italic> from class Ia genes in squamates. Outside of the squamate grouping, the tree topology is not well supported, suggesting that class I sequences are highly divergent among vertebrate groups.</p> <p>In summary, the <italic>Amcr-UA</italic> sequence possesses several characteristics of a functional, non-classical class I gene with a conserved protein structure. Additional work must be conducted to understand whether it is expressed at the level of the protein and what its position is in the genome relative to published classical loci ##REF##18488213##[16]##. While <italic>Amcr-UA</italic> is most closely related to the other published squamate sequences, it does not cluster with class Ia sequences from the same species, suggesting that it has long been on a separate evolutionary trajectory. Further characterization of <italic>UA</italic>-like sequences from other squamates will reveal whether this gene is part of a lineage that has maintained a non-classical function over the course of squamate evolution.</p>
<p>Conceived and designed the experiments: SG. Performed the experiments: SG. Analyzed the data: SG LdP. Contributed reagents/materials/analysis tools: AC. Wrote the paper: SG.</p> <p>Squamates are a diverse order of vertebrates, representing more than 7,000 species. Yet, descriptions of full-length major histocompatibility complex (MHC) genes in this group are nearly absent from the literature, while the number of MHC studies continues to rise in other vertebrate taxa. The lack of basic information about MHC organization in squamates inhibits investigation into the relationship between MHC polymorphism and disease, and leaves a large taxonomic gap in our understanding of amniote MHC evolution. Here, we use both cDNA and genomic sequence data to characterize a class I MHC gene (<italic>Amcr-UA</italic>) from the Galápagos marine iguana, a member of the squamate subfamily Iguaninae. <italic>Amcr-UA</italic> appears to be functional since it is expressed in the blood and contains many of the conserved peptide-binding residues that are found in classical class I genes of other vertebrates. In addition, comparison of <italic>Amcr-UA</italic> to homologous sequences from other iguanine species shows that the antigen-binding portion of this gene is under purifying selection, rather than balancing selection, and therefore may have a conserved function. A striking feature of <italic>Amcr-UA</italic> is that both the cDNA and genomic sequences lack the transmembrane and cytoplasmic domains that are necessary to anchor the class I receptor molecule into the cell membrane, suggesting that the product of this gene is secreted and consequently not involved in classical class I antigen-presentation. The truncated and conserved character of <italic>Amcr-UA</italic> lead us to define it as a nonclassical gene that is related to the few available squamate class I sequences. However, phylogenetic analysis placed <italic>Amcr-UA</italic> in a basal position relative to other published classical MHC genes from squamates, suggesting that this gene diverged near the beginning of squamate diversification.</p>
<p>Maria Moreno and Stephen Dellaporta (Yale University) provided helpful advice for carrying out RACE protocol. Larry Buckley (Rochester Institute of Technology) contributed genomic DNA from four of the iguanine species used in this study (<italic>C. defensor</italic>, <italic>C. clarki</italic>, C. <italic>rileyi</italic>, and <italic>C. carinata</italic>). The <italic>C. cornuta</italic> specimen was provided by the Peabody Museum of Natural History with the help of Greg Watkins-Colwell. The Charles Darwin Research Station (Santa Cruz, Galápagos, Ecuador) and the Galápagos National Park provided logistical support for the collection and export of marine iguana blood and RNA. Cruz Marquez helped with marine iguana blood sampling. Special thanks to Ylenia Chiari and Hilary Miller for helpful comments on the manuscript.</p>
<fig id="pone-0002859-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002859.g001</object-id><label>Figure 1</label><caption><title>Map of the <italic>Amcr-UA</italic> class I gene based on genomic and cDNA data.</title><p>All exons are scaled relative to each other. Introns are also scaled relative to each other, but are more condensed than exons. Asterisk in the α-3 amino acid sequence represents the termination signal. Numbers above arrows show the position of primers used in this study and correspond to the numbering scheme in ##TAB##0##Table 1## – except for primer 11, which has been published elsewhere ##UREF##4##[18]##.</p></caption></fig> <fig id="pone-0002859-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002859.g002</object-id><label>Figure 2</label><caption><title>3′ sequence map of the <italic>Amcr-UA</italic> fragment generated from cDNA.</title><p>Possible polyadenylation sites are labeled with a solid underline. The asterisk marks the termination signal at the end of exon 4 (α-3 domain). The eleven amino acid residues that mark the unusual extension of the α-3 domain are labeled above their respective coding nucleotide sequences. The sequence of primer MHC17H3-5R-R1 is labeled with a broken underline; this was the reverse primer used to generate the genomic DNA sequence and shows how far the genomic sequence stretches downstream relative to the cDNA.</p></caption></fig> <fig id="pone-0002859-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002859.g003</object-id><label>Figure 3</label><caption><title>Amino acid alignment of <italic>Amcr-UA</italic> with other vertebrate class I sequences.</title><p>Coding domains are separated according to Koller and Orr ##REF##2982951##[52]## and exon/intron information from the <italic>Amcr-UA</italic> genomic sequence. Numerical labels in the alignment refer to amino acid positions in the <italic>Amcr-UA</italic> sequence, not in the alignment itself. Shaded columns indicate amino acid positions that are conserved or have expected functions. These amino acid positions also contain the following additional labels: “asterisks” = conserved peptide-binding residues of antigen N- and C- termini; “diamonds” = salt bridge-forming residues; “circles” = disulfide bridge-forming cysteines; “squares” = N-glycosylation site; “CD8” = expected CD8 binding site. The boxed sequence from <italic>Amcr-</italic>UA represents the 11 amino acid extension that is encoded in exon 4, which corresponds to the α-3 domain. Information and accession numbers of other vertebrate class I sequences can be found in the <xref ref-type="sec" rid="s2">methods</xref> section.</p></caption></fig> <fig id="pone-0002859-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002859.g004</object-id><label>Figure 4</label><caption><title>Class I phylogenetic reconstruction using α-3 (exon 4) sequence data from the major vertebrate groups.</title><p>The top number at each node is the Bayesian posterior distribution value supporting the topology displayed. The middle and bottom numbers are bootstrap values from maximum likelihood and neighbor joining searches respectively. Branch lengths are represented as expected number of substitutions per site and were inferred along with topology in the Bayesian analysis.</p></caption></fig> <fig id="pone-0002859-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002859.g005</object-id><label>Figure 5</label><caption><title>(a) Nucleotide and (b) amino-acid alignments of iguanine exon 3 (α-2 domain) sequences.</title><p>The lower-case letters and shaded boxes represent specific codon positions that show variation in amino acid residues between species. Species abbreviations are described in the Materials and <xref ref-type="sec" rid="s2">methods</xref> section.</p></caption></fig>
<table-wrap id="pone-0002859-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002859.t001</object-id><label>Table 1</label><caption><title>Gene-specific primers used in this study.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Primer Name</td><td align="left" rowspan="1" colspan="1">Sequence (5′ to 3′)</td><td align="left" rowspan="1" colspan="1">Target Region</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">(1) MHCI-R3</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">ASGTAYYTCBBCAGCCACT</named-content>
</td><td align="left" rowspan="1" colspan="1">Initial exon 3 genomic amplification</td></tr><tr><td align="left" rowspan="1" colspan="1">(2) IgNC-F1</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">TGYGAGCTGAGGAAAGATGGGAGCATAG</named-content>
</td><td align="left" rowspan="1" colspan="1">3′ RACE</td></tr><tr><td align="left" rowspan="1" colspan="1">(3) IgNC-F2</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">TTACCARTGTGCTTATGATGGGAGGGAC</named-content>
</td><td align="left" rowspan="1" colspan="1">3′ RACE (Nested)</td></tr><tr><td align="left" rowspan="1" colspan="1">(4) MHC17H3-5R-R1</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">TTCATGCAGTCCAAAGGCAGCAG</named-content>
</td><td align="left" rowspan="1" colspan="1">5′ RACE; Genomic fragment</td></tr><tr><td align="left" rowspan="1" colspan="1">(5) MHC17H3-5R-R2</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">TCTTCTGCCTTGCTTCTGTCAAATATGGAG</named-content>
</td><td align="left" rowspan="1" colspan="1">5′ RACE (Nested)</td></tr><tr><td align="left" rowspan="1" colspan="1">(6) MHC1-Gr2-F1</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">ACCGAGAGGGTTGAGCTGGAGAG</named-content>
</td><td align="left" rowspan="1" colspan="1">Genomic fragment</td></tr><tr><td align="left" rowspan="1" colspan="1">(7) MHCI-Int-R1</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">TGAGGCTGASAGAGTGTAATCCTCCC</named-content>
</td><td align="left" rowspan="1" colspan="1">Iguaninae exon3/intron3 fragment</td></tr><tr><td align="left" rowspan="1" colspan="1">(8) MHC1-Int2-F2</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">AGRTGTRTGTAATATATCATCCAGG</named-content>
</td><td align="left" rowspan="1" colspan="1">Intron3/exon4 fragment (All Iguaninae)</td></tr><tr><td align="left" rowspan="1" colspan="1">(9) MHC1-R4</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">CAGRCYGGYRTGMTCCACRCGSCAC</named-content>
</td><td align="left" rowspan="1" colspan="1">Intron3/exon4 fragment (<italic>Amcr; Cosu</italic>)</td></tr><tr><td align="left" rowspan="1" colspan="1">(10) MHC1-R6</td><td align="left" rowspan="1" colspan="1">
<named-content content-type="gene">TCTCCTTGGGGTAGAAGCCRTC</named-content>
</td><td align="left" rowspan="1" colspan="1">Intron3/exon4 fragment (<italic>Cyco</italic>; <italic>Ctde</italic>; <italic>Cyca</italic>; <italic>Ctcl</italic>; <italic>Cyri</italic>)</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0002859-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002859.t002</object-id><label>Table 2</label><caption><title>Percent amino acid identities between <italic>Amcr-UA</italic> and class I sequences from other vertebrates.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Species</td><td colspan="4" align="left" rowspan="1">Percent Identity with <italic>Amcr-UA</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">α1</td><td align="left" rowspan="1" colspan="1">α2</td><td align="left" rowspan="1" colspan="1">α3</td><td align="left" rowspan="1" colspan="1">α1, α2, α3</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">Marine iguana (<italic>Amcr-UB</italic>)</td><td align="left" rowspan="1" colspan="1">43.5–46.7</td><td align="left" rowspan="1" colspan="1">74.7–75.8</td><td align="left" rowspan="1" colspan="1">52.2</td><td align="left" rowspan="1" colspan="1">56.2–57.8</td></tr><tr><td align="left" rowspan="1" colspan="1">Land iguana (<italic>Cosu-UB</italic>)</td><td align="left" rowspan="1" colspan="1">45.7–50.0</td><td align="left" rowspan="1" colspan="1">68.1–75.8</td><td align="left" rowspan="1" colspan="1">56.5</td><td align="left" rowspan="1" colspan="1">57.8–60.0</td></tr><tr><td align="left" rowspan="1" colspan="1">Green iguana (<italic>Igig-UB</italic>)</td><td align="left" rowspan="1" colspan="1">46.7–48.9</td><td align="left" rowspan="1" colspan="1">73.6–76.7</td><td align="left" rowspan="1" colspan="1">53.3</td><td align="left" rowspan="1" colspan="1">57.8–59.5</td></tr><tr><td align="left" rowspan="1" colspan="1">Ameiva Lizard (LC1)</td><td align="left" rowspan="1" colspan="1">46.7</td><td align="left" rowspan="1" colspan="1">64.8</td><td align="left" rowspan="1" colspan="1">48.4</td><td align="left" rowspan="1" colspan="1">53.3</td></tr><tr><td align="left" rowspan="1" colspan="1">Tuatara (<italic>Sppu-U</italic>)</td><td align="left" rowspan="1" colspan="1">45.5–50.0</td><td align="left" rowspan="1" colspan="1">50.5</td><td align="left" rowspan="1" colspan="1">36.3</td><td align="left" rowspan="1" colspan="1">44.1–45.6</td></tr><tr><td align="left" rowspan="1" colspan="1">Chicken (<italic>Gaga</italic>-B-F10)</td><td align="left" rowspan="1" colspan="1">40.9</td><td align="left" rowspan="1" colspan="1">52.2</td><td align="left" rowspan="1" colspan="1">36.3</td><td align="left" rowspan="1" colspan="1">43.1</td></tr><tr><td align="left" rowspan="1" colspan="1">Mallard (<italic>Anpl</italic>-Du2)</td><td align="left" rowspan="1" colspan="1">45.5</td><td align="left" rowspan="1" colspan="1">56.0</td><td align="left" rowspan="1" colspan="1">40.7</td><td align="left" rowspan="1" colspan="1">47.4</td></tr><tr><td align="left" rowspan="1" colspan="1">Warbler (<italic>Acar</italic>-cN3)</td><td align="left" rowspan="1" colspan="1">43.2</td><td align="left" rowspan="1" colspan="1">46.2</td><td align="left" rowspan="1" colspan="1">44.0</td><td align="left" rowspan="1" colspan="1">44.4</td></tr><tr><td align="left" rowspan="1" colspan="1">Axolotl (<italic>Amme</italic>-3)</td><td align="left" rowspan="1" colspan="1">43.2</td><td align="left" rowspan="1" colspan="1">54.9</td><td align="left" rowspan="1" colspan="1">34.8</td><td align="left" rowspan="1" colspan="1">44.4</td></tr><tr><td align="left" rowspan="1" colspan="1">Xenopus (<italic>Xela</italic>-UAA-1f)</td><td align="left" rowspan="1" colspan="1">44.3</td><td align="left" rowspan="1" colspan="1">44.0</td><td align="left" rowspan="1" colspan="1">37.4</td><td align="left" rowspan="1" colspan="1">41.9</td></tr><tr><td align="left" rowspan="1" colspan="1">Human (<italic>Hosa</italic>-HLA-A)</td><td align="left" rowspan="1" colspan="1">43.7</td><td align="left" rowspan="1" colspan="1">53.8</td><td align="left" rowspan="1" colspan="1">38.0</td><td align="left" rowspan="1" colspan="1">45.2</td></tr><tr><td align="left" rowspan="1" colspan="1">Mouse (<italic>Mumu</italic>-Q10)</td><td align="left" rowspan="1" colspan="1">44.8</td><td align="left" rowspan="1" colspan="1">57.8</td><td align="left" rowspan="1" colspan="1">34.8</td><td align="left" rowspan="1" colspan="1">45.7</td></tr><tr><td align="left" rowspan="1" colspan="1">Mouse (<italic>Mumu</italic>-H2K)</td><td align="left" rowspan="1" colspan="1">43.7</td><td align="left" rowspan="1" colspan="1">52.2</td><td align="left" rowspan="1" colspan="1">31.5</td><td align="left" rowspan="1" colspan="1">42.4</td></tr><tr><td align="left" rowspan="1" colspan="1">Wallaby (<italic>Maru</italic>-UB*01)</td><td align="left" rowspan="1" colspan="1">40.0</td><td align="left" rowspan="1" colspan="1">54.9</td><td align="left" rowspan="1" colspan="1">38.0</td><td align="left" rowspan="1" colspan="1">44.3</td></tr><tr><td align="left" rowspan="1" colspan="1">Platypus (<italic>Oran</italic>-2-1)</td><td align="left" rowspan="1" colspan="1">40.9</td><td align="left" rowspan="1" colspan="1">56.0</td><td align="left" rowspan="1" colspan="1">32.6</td><td align="left" rowspan="1" colspan="1">43.2</td></tr><tr><td align="left" rowspan="1" colspan="1">Possum (<italic>Trvu</italic>-UB)</td><td align="left" rowspan="1" colspan="1">46.7</td><td align="left" rowspan="1" colspan="1">56.0</td><td align="left" rowspan="1" colspan="1">34.8</td><td align="left" rowspan="1" colspan="1">45.8</td></tr><tr><td align="left" rowspan="1" colspan="1">Nurse Shark (<italic>Gici</italic>-UAA)</td><td align="left" rowspan="1" colspan="1">41.4</td><td align="left" rowspan="1" colspan="1">44.0</td><td align="left" rowspan="1" colspan="1">32.2</td><td align="left" rowspan="1" colspan="1">39.2</td></tr><tr><td align="left" rowspan="1" colspan="1">Rainbow Trout (<italic>Onmy</italic>-UBA)</td><td align="left" rowspan="1" colspan="1">33.3</td><td align="left" rowspan="1" colspan="1">40.7</td><td align="left" rowspan="1" colspan="1">24.7</td><td align="left" rowspan="1" colspan="1">33.0</td></tr><tr><td align="left" rowspan="1" colspan="1">Zebrafish (<italic>Dare</italic>-UBA)</td><td align="left" rowspan="1" colspan="1">27.6</td><td align="left" rowspan="1" colspan="1">51.6</td><td align="left" rowspan="1" colspan="1">24.7</td><td align="left" rowspan="1" colspan="1">34.8</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0002859-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002859.t003</object-id><label>Table 3</label><caption><title>Uncorrected pairwise distances (<italic>p-</italic>distance) between <italic>Amcr-UA</italic> and similar sequences from other iguanine species.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">7</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<bold>1</bold>
<italic>Ctde</italic>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.04435</td><td align="left" rowspan="1" colspan="1">0.03226</td><td align="left" rowspan="1" colspan="1">0.02823</td><td align="left" rowspan="1" colspan="1">0.03226</td><td align="left" rowspan="1" colspan="1">0.04435</td><td align="left" rowspan="1" colspan="1">0.04839</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>2</bold>
<italic>Ctcl</italic>
</td><td align="left" rowspan="1" colspan="1">0.04546</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.03629</td><td align="left" rowspan="1" colspan="1">0.03226</td><td align="left" rowspan="1" colspan="1">0.03629</td><td align="left" rowspan="1" colspan="1">
<italic>0.06048</italic>
</td><td align="left" rowspan="1" colspan="1">0.05645</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>3</bold>
<italic>Cyco</italic>
</td><td align="left" rowspan="1" colspan="1">0.04789</td><td align="left" rowspan="1" colspan="1">0.05722</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.00403</td><td align="left" rowspan="1" colspan="1">0.00806</td><td align="left" rowspan="1" colspan="1">0.03226</td><td align="left" rowspan="1" colspan="1">0.02823</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>4</bold>
<italic>Cyca</italic>
</td><td align="left" rowspan="1" colspan="1">0.04693</td><td align="left" rowspan="1" colspan="1">0.05822</td><td align="left" rowspan="1" colspan="1">0.01975</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.00403</td><td align="left" rowspan="1" colspan="1">0.02823</td><td align="left" rowspan="1" colspan="1">0.02419</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>5</bold>
<italic>Cyri</italic>
</td><td align="left" rowspan="1" colspan="1">0.04791</td><td align="left" rowspan="1" colspan="1">0.05430</td><td align="left" rowspan="1" colspan="1">0.02073</td><td align="left" rowspan="1" colspan="1">0.00988</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.03226</td><td align="left" rowspan="1" colspan="1">0.02823</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>6</bold>
<italic>Cosu</italic>
</td><td align="left" rowspan="1" colspan="1">0.04671</td><td align="left" rowspan="1" colspan="1">
<italic>0.05694</italic>
</td><td align="left" rowspan="1" colspan="1">0.03846</td><td align="left" rowspan="1" colspan="1">0.03552</td><td align="left" rowspan="1" colspan="1">0.03850</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.01210</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>7</bold>
<italic>Amcr-UA</italic>
</td><td align="left" rowspan="1" colspan="1">0.05366</td><td align="left" rowspan="1" colspan="1">0.06577</td><td align="left" rowspan="1" colspan="1">0.04924</td><td align="left" rowspan="1" colspan="1">0.04633</td><td align="left" rowspan="1" colspan="1">0.05120</td><td align="left" rowspan="1" colspan="1">0.02650</td><td align="left" rowspan="1" colspan="1"/></tr></tbody></table></alternatives></table-wrap>
<table-wrap-foot><fn id="nt101"><p>Numbers before primer refer to location in the gene as indicated by arrows in ##FIG##0##Figure 1##.</p></fn><fn id="nt102"><p>Footnote: See <xref ref-type="sec" rid="s2">Materials and methods</xref> for Iguaninae taxa abbreviations.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt103"><p>Values were calculated for each extracellular domain separately and for all three domains combined.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt104"><p>Distances were calculated separately for exon 3 (above the diagonal) and intron 3 (below the diagonal). See <xref ref-type="sec" rid="s2">methods</xref> for Iguaninae taxa abbreviations. Pairwise values between <italic>Ctcl</italic> and <italic>Cosu</italic> are italicized because they deviate from the otherwise typical pattern of higher divergence in intron 3 versus exon 3.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>Funding for this project was provided by the Yale Institute for Biospheric Studies (YIBS). Work was also supported by pilot study funds from the YIBS Center for Field Ecology. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p></fn></fn-group>
<graphic xlink:href="pone.0002859.g001"/> <graphic id="pone-0002859-t001-1" xlink:href="pone.0002859.t001"/> <graphic xlink:href="pone.0002859.g002"/> <graphic xlink:href="pone.0002859.g003"/> <graphic id="pone-0002859-t002-2" xlink:href="pone.0002859.t002"/> <graphic xlink:href="pone.0002859.g004"/> <graphic id="pone-0002859-t003-3" xlink:href="pone.0002859.t003"/> <graphic xlink:href="pone.0002859.g005"/>
[{"label": ["1"], "element-citation": ["\n"], "surname": ["Klein"], "given-names": ["J"], "year": ["1986"], "article-title": ["Natural history of the major histocompatibility complex."], "publisher-loc": ["New York"], "publisher-name": ["Wiley"], "fpage": ["xv, 775"]}, {"label": ["11"], "element-citation": ["\n"], "surname": ["Hess", "Edwards"], "given-names": ["CM", "SV"], "year": ["2002"], "article-title": ["The evolution of the major histocompatibility complex in birds."], "source": ["Bioscience"], "volume": ["52"], "fpage": ["423"], "lpage": ["431"]}, {"label": ["12"], "element-citation": ["\n"], "surname": ["Belov", "Deakin", "Papenfuss", "Baker", "Melman"], "given-names": ["K", "JE", "AT", "ML", "SD"], "year": ["2006"], "article-title": ["Reconstructing an ancestral mammalian immune supercomplex from a marsupial major histocompatibility complex."], "source": ["PLoS Biol"], "volume": ["4"], "fpage": ["317"], "lpage": ["328"]}, {"label": ["17"], "element-citation": ["\n"], "surname": ["Flajnik"], "given-names": ["M"], "year": ["2004"], "article-title": ["Comparative genomics of the MHC."], "source": ["Tissue Antigens"], "volume": ["64"], "fpage": ["328"], "lpage": ["328"]}, {"label": ["18"], "element-citation": ["\n"], "surname": ["Radtkey", "Becker", "Miller", "Riblet", "Case"], "given-names": ["RR", "B", "RD", "R", "TJ"], "year": ["1996"], "article-title": ["Variation and evolution of class I Mhc in sexual and parthenogenetic geckos."], "source": ["Proc R Soc London Ser B"], "volume": ["263"], "fpage": ["1023"], "lpage": ["1032"]}, {"label": ["19"], "element-citation": ["\n"], "surname": ["Sambrook", "Russell"], "given-names": ["J", "DW"], "year": ["2001"], "article-title": ["Molecular cloning : a laboratory manual."], "publisher-loc": ["Cold Spring Harbor, N.Y."], "publisher-name": ["Cold Spring Harbor Laboratory Press"]}, {"label": ["25"], "element-citation": ["\n"], "surname": ["Nylander"], "given-names": ["JAA"], "year": ["2004"], "article-title": ["MrModeltest v2. Evolutionary Biology Centre, Uppsala University: Program distributed by the author"]}, {"label": ["39"], "element-citation": ["\n"], "surname": ["Rassmann"], "given-names": ["K"], "year": ["1997"], "article-title": ["Evolutionary age of the Gal\u00e1pagos iguanas predates the age of the present Gal\u00e1pagos Islands."], "source": ["Mol Phylogenet and Evol"], "volume": ["7"], "fpage": ["158"], "lpage": ["172"]}, {"label": ["42"], "element-citation": ["\n"], "surname": ["Rassmann", "Tautz", "Trillmich", "Gliddon"], "given-names": ["K", "D", "F", "C"], "year": ["1997"], "article-title": ["The microevolution of the Gal\u00e1pagos marine iguana "], "italic": ["Amblyrhynchus cristatus"], "source": ["Mol Ecol"], "volume": ["6"], "fpage": ["437"], "lpage": ["452"]}]
PLoS One. 2008 Aug 6; 3(8):e2859