Disruption and pseudoautosomal localization of the major histocompatibility complex in monotremes
© Dohm et al.; licensee BioMed Central Ltd. 2007
Received: 28 May 2007
Accepted: 29 August 2007
Published: 29 August 2007
The monotremes, represented by the duck-billed platypus and the echidnas, are the most divergent species within mammals, featuring a flamboyant mix of reptilian, mammalian and specialized characteristics. To understand the evolution of the mammalian major histocompatibility complex (MHC), the analysis of the monotreme genome is vital.
We characterized several MHC containing bacterial artificial chromosome clones from platypus (Ornithorhynchus anatinus) and the short-beaked echidna (Tachyglossus aculeatus) and mapped them onto chromosomes. We discovered that the MHC of monotremes is not contiguous and locates within pseudoautosomal regions of two pairs of their sex chromosomes. The analysis revealed an MHC core region with class I and class II genes on platypus and echidna X3/Y3. Echidna X4/Y4 and platypus Y4/X5 showed synteny to the human distal class III region and beyond. We discovered an intron-containing class I pseudogene on platypus Y4/X5 at a genomic location equivalent to the human HLA-B,C region, suggesting ancestral synteny of the monotreme MHC. Analysis of male meioses from platypus and echidna showed that MHC chromosomes occupy different positions in the meiotic chains of either species.
Molecular and cytogenetic analyses reveal new insights into the evolution of the mammalian MHC and the multiple sex chromosome system of monotremes. In addition, our data establish the first homology link between chicken microchromosomes and the smallest chromosomes in the monotreme karyotype. Our results further suggest that segments of the monotreme MHC that now reside on separate chromosomes must once have been syntenic and that the complex sex chromosome system of monotremes is dynamic and still evolving.
The major histocompatibility complex (MHC) is of central importance for adaptive and innate immunity in vertebrates . Sequencing MHCs from several species of eutherian mammals and human has led to the identification of approximately 220 genes located within an interval of 3.5-4 Mbp [2, 3]. The MHC region contains genes encoding class I and class II receptors that are involved in peptide display, genes that are responsible for peptide generation and transport, as well as genes encoding complement factors or cytokines. Many other genes with functions that are not related to immunity and defense are also located within the MHC.
In contrast to the conserved MHC architecture, much plasticity is observed within class I gene containing segments and, to a lesser extent, in the class II region. Orthologous class I genes are found in closely related species, for example, human and chimpanzees. Class I gene orthology becomes less apparent when more distantly related mammals are compared. This is attributed to the ongoing process of new class I genes arising by gene duplication and, at the same time, their disappearance from the genome by mutational inactivation or by deletion. For example, the sequences of the classical class I genes from rat form a separate clade when compared to class I genes from mouse, indicating species-specific evolution since separation of the two lineages 20 million years (Myr) ago .
In addition to comparative studies on eutherian mammals, other species have been studied and compared to the eutherians: recently, the sequence of a metatherian mammal (marsupial) MHC from the opossum (Monodelphis domestica) was assembled, based on whole genome shotgun sequence data . Its annotation revealed a core MHC region, resembling the eutherian class II region, that, unlike in eutherian mammals, also encodes the class I genes. The opossum MHC was found to be syntenic to a human MHC region of several megabase pairs. Also in Xenopus tropicalis extensive synteny with the eutherian MHC has been found. Ohta et al.  analyzed the scaffolds that were generated from assembling Xenopus whole genome shotgun data for the presence of MHC-encoded genes, including framework genes. They found eight 200-900 kbp scaffolds that contained orthologs to many of the genes in the eutherian MHC. The interpretation is that the syntenies of the MHC region are ancient and predate the amphibian-mammalian split 350 Myr ago even though it has not been shown that these scaffolds reside on a single chromosome. In contrast, the MHC of the chicken, representing a lineage that is closer to mammals than Xenopus (distance aves-mammalia of 310 Myr ), is highly derived and does not align well with the human MHC . The chicken MHC is remarkably different and contains only 19 genes within 92 kbp on chicken chromosome 16.
Monotremes are prototherian mammals, possessing fur and mammary glands, but uniquely combining reptilian features, for example, egg-laying with mammalian characteristics. They are the earliest offshoot of the mammalian clade that separated from theria (consisting of metatheria and eutheria, or marsupials and placental mammals, respectively) over 160 Myr ago . The only extant prototheria are the duck-billed platypus (Ornithorhynchus anatinus), the short-beaked echidna (Tachyglossus aculeatus) and three species of long-beaked echidnas (Zaglossus sp.) . The geographic range of current monotremes is restricted to Australia and nearby islands, including Papua New Guinea. In the following, when mentioning echidna, we refer to the short-beaked echidna, also known as the spiny anteater.
Because of their unique phylogenetic position, monotremes are ideal species to unravel the characteristics of the ancestral mammalian MHC. So far, little is known about MHC genes in monotremes and there are only two reports where individual cDNA clones were identified as MHC class I or class II genes [16, 17].
The general importance of monotremes for understanding the evolution and function of the mammalian genome is now generally accepted and is reflected by the presently ongoing efforts to sequence the platypus genome, with funding from the National Institute of Health (NIH) .
The platypus karyotype consists of 21 pairs of autosomes and 10 sex chromosomes (X1Y1-X5Y5 in males and X1X1-X5X5 in females). In addition, the platypus karyotype contains a number of very small chromosomes that have been proposed to be bird-like microchromosomes; however, the homology between these chromosomes and chicken microchromosomes has been controversial and no homologies have been identified to date [19, 20]. The sex chromosomes assemble as a chain at prophase I during male meiosis where they adopt an XY alternating pattern, which ensures the segregation into X1-X5 and Y1-Y5 bearing sperm [21, 22]. Although XY shared regions are required to ensure meiotic pairing, chromosome painting could not identify pseudoautosomal regions on some of the sex chromosomes so far. The echidna karyotype contains 27 autosomes and 10 sex chromosomes in females and 9 sex chromosomes in males, which where determined by counting the number of elements of the meiotic chain [23, 24]. Very little is known about the gene content of both sex chromosome systems.
Here we describe the identification and characterization of genomic clones encompassing parts of the monotreme MHC. The unexpected localization of MHC genes on sex chromosomes of platypus and echidna provides valuable information about chromosomal homologies between aves, monotremes and eutheria, as well as homologies between platypus and echidna sex chromosomes.
Identification of BAC clones from the platypus and echidna MHC
We identified bacterial artificial chromosome (BAC) clones from the platypus MHC by hybridization with oligonucleotide probes designed from the platypus class I gene Oran2-1 , the class II gene DZB , and from the MHC framework genes Bat1 and Pou5f1 encoding a DEAD-Box helicase and a stem cell transcription factor, respectively. Probes for the framework genes were designed by aligning the chicken Bat1 and Pou5f1 mRNA sequences with the rat MHC . A 36-mer oligonucleotide probe for Bat1 corresponded to rat BAT1 amino acid residues 265-276, which are 100% conserved between rat, chicken and zebrafish. A 37-mer oligonucleotide probe was generated matching a segment of Pou5f1 whose translation had 100% amino acid identity in rat, chicken and zebrafish (rat POU5F1 residues 174-185).
Using probes based on platypus MHC sequences, we searched for BACs from the corresponding regions in the echidna genome. A platypus Tap2 probe identified two positive clones in echidna (BACs 48g5 and 287o10; Additional data file 1a). Separate Southern blots prepared from these clones were hybridized with the Oran2-1 probe  and with a pool of two different class II gene oligonucleotide probes (Additional data file 1b,c). Echidna BAC 48g5 was found to contain class I and class II gene sequences, while BAC 287o10 contained class I sequences but no class II genes. Shared restriction fragments indicated that these two clones were from the same locus and overlapped.
Further echidna BAC colony screens were done using oligonucleotide probes designed on the platypus genes Bat1 and Pou5f1. For each gene, three probes were synthesized and pooled before labeling (one pool per gene). Echidna BAC clones 107j4, 129j16, 152g23 and 268a21 could be confirmed on Southern blots to contain both the Bat1 and Pou5f1 genes (Additional data file 1d,e).
BAC sequence annotation
We generated reference gene models for the genes identified within platypus BACs 1-3, combining the gene prediction programs Genscan  and Genewise  as tools and utilized rat MHC protein sequences  as input for Genewise predictions. The translation products from predictions with either program were compared to each other and to the rat protein as a reference by pairwise blast (Additional data file 7). In the absence of monotreme expressed sequence tag (EST) sequences for the majority of genes, the annotations covered open reading frames (ORFs) and did not include untranslated regions (UTRs). Genscan and Genewise complement each other: The accuracy of Genewise is very good in coding regions that are well conserved, while theoretical gene models predicted by Genscan are exclusively based on signals and content. In total, we annotated 25 genes and two pseudogenes in the three BAC sequences. Within 320,837 sequenced bases in BAC1 and BAC2, one gene could be annotated with full ORF plus UTRs (DZB gene), 20 genes were annotated with their full ORFs only, two genes had missing 5' ends because of low conservation and two genes were only partially contained in the sequenced interval. One annotated feature was a non-processed pseudogene fragment that lacked an intact ORF, though partial reconstruction of the protein sequence was possible. Only one single gene feature was discovered and annotated as a pseudogene within BAC3 (see below).
Characterization of a processed class I pseudogene containing platypus BAC clone
Platypus BAC3 was assembled to a finished sequence of 146,150 bp. Fluorescence in situ hybridization (FISH) mapping located BAC3 on the long arm of platypus chromosome 3 (data not shown). During annotation, no genes were identified with the exception of an intronless class I gene that contained a non-interrupted ORF of 359 amino acids. This gene showed 99% sequence identity to OranPS1-2, a class I sequence previously described as a pseudogene . The class I homology was considered as the reason why the clone had been positive in the BAC colony screen, even though the clone was not MHC-derived.
Analysis of the platypus core MHC
The three framework genes were discovered as the platypus orthologs of the immunoproteasome subunit Psmb8, as well as Tap1 (only partially contained in the BAC sequence) and Tap2, which encode for subunits of the antigen peptide transporter. These genes are typically found in the class II region of mammalian MHCs, even though they are involved in the generation and transport of peptides for display by class I molecules. However, in contrast to the MHCs of eutheria and similar to the opossum MHC, platypus possesses an MHC with class II and class I genes being co-localized (Figure 1a,b).
Analysis of the platypus MHC Bat4-Cdsninterval
The sequence of BAC2 was assembled to a finished length of 194,463 bp. It contained 20 genes in a region that extended from Bat4 to Cdsn (Figures 1c and 2b). Thus, this sequence covers part of the class III region up to Bat1/Mccd1, spans a segment that in the human MHC encodes the HLA-B and HLA-C genes, and contains part of the adjacent framework region distal to Pou5f1. The gene order and transcriptional orientation of genes is the same as in other mammalian MHCs. However, we could not identify Ncr3 in the platypus sequence. Ncr3 encodes the natural cytotoxicity triggering receptor 3 and has not been annotated in the opossum MHC but was discovered in both the human HLA and rat RT1. Human and rat NCR3 proteins evolve rapidly and share only 64% amino acid identity (Additional data file 7), which is much lower than the average identity of 88% that has been observed for protein orthologs in these species . Thus, similarity to a putative platypus Ncr3 gene may have been too low for detection. Between Aif1 and Lst1 where Ncr3 is located in human and rat, Genscan predicts a platypus gene with correct position and orientation. This prediction could not be verified by homology searches in the NCBI databases, but might be a very diverged version of the Ncr3 gene. Alternatively, Ncr3 might have been deleted from the platypus genome or might be eutheria specific.
Between platypus proteins and their orthologs in rat or human we found an average of 61-62% amino acid identity (Additional data file 7) and observed much variation in the degree of sequence conservation (range: 27-99%). For instance, BAT1, which is involved in mRNA splicing and transport, was 99% conserved, showing that little innovation in such processes has occurred over the past 160 Myr. In contrast, the framework gene Lst1, encoding a protein of unknown function, had an identity of less than 30%.
The region between Bat1 and Pou5f1 harbors a class I gene block in the MHC of placental mammals. In rat, this region spans 340 kbp and contains 16 class I genes. In platypus, the corresponding region is only 19 kbp in size and contains a class I gene fragment (Oran-ps1) that consists of two exon fragments encoding α2 and α3 domains that are separated by a 283 bp intron. This suggests that the interval once encoded at least one active class I gene in platypus. The sequence identity of the Oran-ps1 translation product to known and new monotreme class I sequences is only 30-40%. We could not find the active counterpart of Oran-ps1 by searching the ENSEMBL platypus whole-genome assembly Oana-5.0 , suggesting that Oran-ps1 was inactivated a long time ago. The 283 bp intron within Oran-ps1 is not conserved in other class I loci and cannot, therefore, be used for phylogenetic analysis to date the age of the gene.
Chromosomal mapping of platypus MHC BAC clones
FISH mapping of echidna MHC BACs
The platypus MHC is unusual when compared to MHCs from eutherian mammals, because it contains a core region that encodes both class I and class II genes, though such a common class I/II region has also been described in the opossum . The coincidence of compact class I/II regions in both a monotreme and a marsupial supports the hypothesis that this is the ancestral organization of the mammalian MHC. However, a platypus class I gene fragment, Oran-ps1, exists at a position that contains a class I gene block in eutheria. Also, the opossum MHC contains a class I-like MIC gene at this position (Figure 1). Both findings contradict the strict concept of a core MHC and indicate that early mammals already encoded class I genes outside a common class I/II region at positions equivalent to class I gene blocks in eutheria, embedded between framework genes.
The intron-containing class I pseudogene Oran-ps1 and the surrounding region that corresponds to the distal class III region and beyond map to platypus chromosome Y4/X5. In contrast, the core MHC locates to X3/Y3. In the eutherian MHC, class I genes are distributed among several class I gene clusters, separated by hundreds of kilobase-pairs of intervening regions. Despite such distances, concerted evolution takes place between different class I blocks within a species. For instance, the rat RT1-A and RT1-CE regions are more than 1.2 Mbp apart and have a history of transposition events that led to the presence of similar class I genes and shared non-class I gene sequences (gene fragments) in both intervals. Similar but unrelated transposition events were observed within the mouse MHC region, showing that reorganization across large genomic distances is fairly common . These distances seem to be large when considering linear DNA but the regions may be actually juxtaposed in chromatin, facilitating transposition. Oran-ps1 is flanked by Bat1 and Pou5f1 and these framework genes also flank a class I gene cluster in the eutherian MHC. Thus, Oran-ps1 may have been generated during concerted evolution of class I genes in the platypus MHC at a time when the distal class III region and the core MHC were still on the same chromosome. However, concerted evolution between MHC segments on different platypus chromosomes could be still ongoing, by means of trans-interaction of separate chromosomes .
We localized the monotreme MHC to PARs of sex chromosomes. This is the first time MHC genes have been mapped onto sex chromosomes in any mammal or vertebrate. In all mammals investigated so far, the MHC region is located on a single autosome, and in most cases without interruptions. Exceptions are the swine leukocyte antigen (SLA) complex, where the centromere of chromosome 7 separates the class II region from the other parts of the MHC  and the recently described class I genes of the tammar wallaby, which were found to be dispersed in the wallaby genome .
The multiple sex chromosomes in platypus and echidna pair at the first meiotic division and form a sex chromosome chain during prophase I. PARs are necessary to ensure pairing and segregation of most mammalian sex chromosomes. The ten sex chromosomes of platypus, therefore, would require nine PARs to connect the sex chromosomes at meiosis. The extent of XY homology (that is, the size of PARs; see below) is clearly different among the sex chromosomes. While some PARs are large enough and could easily be identified by chromosome painting, other PARs seem to be much smaller and did not show XY homology using whole chromosome paints [21, 22]. We found that MHC genes reside on two of these small PARs in platypus. BAC2 was mapped to the PARs of Y4 and the short arm of X5. This fits with the orientation of X5 in the chain where the short arm lies towards Y4 . Together with the fact that chromosome painting could not detect these PARs, it shows that the Y4X5 shared region is much smaller than the extensive homology found between, for example, X1 and Y1.
PAR sizes vary widely. For instance, the human PAR regions PAR1 and PAR2 differ eight-fold with respect to their sizes, and in the mouse, PAR1 is almost four-fold smaller than in human . The small size of PARs and an obligatory crossover to ensure proper segregation at meiosis results in a high recombination frequency in mammalian PARs . It is therefore interesting to note that the MHC, which is one of the recombination hotspots in the human genome, resides within the PAR, another recombination hotspot.
Another special feature of the monotreme karyotype is the presence of several tiny chromosomes. Initially they where regarded similarly to reptilian microchromosomes [38, 39], which was later disputed on the basis of even distribution of chromosome size in the monotreme karyotype . Microchromosomes have long been regarded as heterochromatic elements but have been shown recently to represent the gene rich regions of the chicken genome . The identification of the ten sex chromosomes in platypus revealed that some distinctively small chromosomes are not autosomes, but part of the sex chromosome system (Y3, X4, Y4). The mapping of MHC genes (class I, class II, Tap1, Tap2) on X3Y3 provides the first evidence for homology to the MHC bearing chicken microchromosome 16. Likewise, five of the MHC framework genes that map on Y4X5 (Bat2, Aif1, Tnf/Lta, Atp6v1g2) locate to chicken microchromosome 17, establishing another novel synteny link between the genomes of monotremes and birds.
The comparative mapping of the MHC in echidna and platypus provides evidence for homology between the sex chromosome systems in platypus and echidna. Our results support the idea that an originally syntenic MHC has been disrupted as a result of the evolution of the platypus and echidna meiotic sex chromosome chain. Following an initial X-autosome translocation in a common platypus/echidna ancestor, subsequent translocations resulted in the differences observed in echidna and platypus sex chromosome systems.
Materials and methods
BAC colony arrays, probe design and labeling, and hybridization
Macroarrays and clones were obtained from the CUGI BAC/EST resource centre (Clemson, SC, USA) for the platypus BAC library OA_B and from the AGI BAC/EST resource centre (Tuscon, AZ, USA) for the echidna BAC library TA_Ba. Platypus colony arrays were hybridized with the following oligonucleotide probes: Bat1-Gga/Rno, 5'-TACTACGTGAAACTGAAGGACAACGAGAAGAACCGG-3'; Pou5f1-Gga/Rno, 5'-CAGACAACCATCTGCCGCTTCGAGGCCCTGCAGCTCA-3'; probe Oan-class I, 5'-AGGACGTCAGCCTGCGGTGCCGGGCCCTCGGCTTC-3', designed from the α3-domain of platypus class I gene Oran2-1 (AY112715); and platypus class II gene probe Oan-class II, 5'-CAGGGAATGGAGCGGCCCATCCTCCCCATGTTGAA-3', designed from the 3' UTR of the DZB gene (AY288074).
A Tap2 probe was amplified from platypus BAC 462c1 using primers: Tap2-Orna-For, 5'-CACCGTGACCCCTCCATCTTCCTCA-3'; and Tap2-Orna-Rev, 5'-GTGGGGGAGAGCTCTGGGGGCTGGACG-3'. Tap2-positive echidna BAC clones were characterized with the platypus Oan-class I probe (see above) and the two class II gene probes 5'-CGGGAGACATCCAAGTCCAGTGGTTGCGGAATGGA-3' and 5'-TGATGCGTCATGGAGACTGGACCTTCCAGGTCCTG-3'. Pools of three oligonucleotides based on platypus Bat1 and Pou5f1 were used to identify echidna BACs. Oligonucleotide sequences were: Oran-Bat1-1, 5'-GCCTCAACCTGAAGCACATTAAACACTTCATCCTG-3'; Oran-Bat1-2, 5'-CGACGTGCAGGATCGCTTCGAGGTCAACATCAGCG-3'; Oran-Bat1-3, 5'-CAGGTGGTGATCTTCGTGAAGTCAGTGCAGCGCTG-3'; Oran-Pou5f1-1, 5'-ATGGCCGGACACCTGGGTCCCGACTTCGCCTTCTC-3'; Oran-Pou5f1-2, 5'-GACCACCATCTGCCGCTTCGAGGCCCAGCAGCTGA-3'; and Oran-Pou5f1-3, 5'-CCGCGTCTGGTTCTGTAACCGCCGGCAGAAAGGCA-3'.
We labeled 25 pmol of an oligonucleotide using 25 μCi of γ32PdATP (Amersham, Little Chalfont, UK) and 10 u T4-Polynucleotide kinase (New England Biolabs, Ipswich, MA, USA) in a volume of 15 μl at 37°C for 60 minutes. Unincorporated γ32PdATP was removed using G-50 MicroColumns (Amersham). Filters were hybridized in Church buffer (0.25 M Na2HPO4, 1 mM EDTA, 5% SDS) at 58°C over night. Thereafter, filters were washed twice. The first washing step was carried out using 2 × SSC/0.1% SDS wash solution (wash I). The second wash was done using 0.5 × SSC/0.1% SDS (wash II). Each incubation lasted 20 minutes. Filters were exposed on Kodak XAR X-ray film (Sigma-Aldrich, Taufkirchen, Germany) for up to 14 days at -80°C. Procedures for PCR fragment labeling were similar, except that labeling was done with the random primed labeling kit (Roche Diagnostics, Mannheim, Germany), following the manufacturer's instructions, under the addition of 25 μCi of α32PdATP (Amersham). Hybridizations and washing steps were carried out at 65°C. The second washing step was done using 0.1 × SSC/0.1% SDS.
Shotgun library construction and sequencing
BAC clones were grown in 1,000 ml cultures at 37°C with agitation (200 rpm) in a shaking incubator over night. DNA was prepared using the alkaline lysis method and purification on a CsCl gradient in a Beckman centrifuge (VTi65i2 rotor) at 45,000 rpm for 20 h, followed by DNA extraction . Finally, DNA was dissolved in 100 μl 1 × TE. Using a Branson Sonifier 250 (Branson Ultrasonic, Danbury, CT, USA), 10 μg BAC DNA was sheared in a volume of 200 μl to an average size of 1.5 kbp, and after precipitation blunt-ended with Klenow fragment (USB, Cleveland, OH, USA) and T4 DNA polymerase (New England Biolabs). End-repaired fragments were size selected by electrophoresis on a gel of 0.8% SeaPlaque GTG agarose (Biozym, Hess. Oldendorf, Germany) in 1 × TAE buffer (40 mM Tris base, 1 mM EDTA, 20 mM acetic acid). A gel slice that contained fragments of 1.2-2 kbp was cut out, molten at 55°C and incubated at 37°C for 3 h in the presence of 10 u agarase (Sigma-Aldrich) per 100 μl of solution. Following chloroform-phenol extraction, the DNA was precipitated and dissolved in 10 μl 1 × TE buffer. Ligation with 30 ng of dephosphorylated, SmaI cut pUC19 and 200-300 ng size-selected DNA was carried out at 16°C overnight, using 400 u T4 DNA ligase (New England Biolabs). A 1 μl ligation was used for transforming 40 μl of ElectroTen Blue Escherichia coli cells (Stratagene, La Jolla, CA, USA), following the manufacturer's recommendations. Cells were grown on 22 cm × 22 cm Nunc dishes on LB agar supplemented with 1 ml 2% X-Gal, 250 μl 1 M Isopropyl-β-D-thiogalactopyranosid (IPTG) and 400 μl 50 mg/ml ampicillin per 400 ml medium. Colonies were arrayed into 384-well plates (Genetix, Dornach, Germany) using in-house custom-built robotic devices. Plates contained 2YT medium supplemented with ampicillin (50 μg/ml) and 1 × HMFM freezing mix (0.4 mM MgSO4 × 7H2O, 1.6 mM Tri-sodium citrate × 2H2O, 6.8 mM (NH4)2SO4, 3.6% glycerol, 13.2 mM KH2PO2 and 26 mM K2HPO4). Plates were incubated at 37°C for at least 12 h and subsequently shock-frozen on dry ice. Sequencing was carried out on PCR amplified pUC19 inserts, using the BigDye terminator chemistry on ABI3700 automated sequencers. Initial sequence coverage of BACs was up to ten-fold. Thereafter, gaps were closed directionally. Analyzed regions were manually edited in GAP4 . Annotation of assembled BACs was carried out using Genscan and Genewise [25, 26]. Exon coordinates in PIPmaker format and predicted translations are available (Additional data files 8 and 9).
We searched the NCBI protein database for MHC class I genes and class II α chain genes in all species groups ranging from amphibians to mammals, based on the NCBI TaxBrowser classification and recorded all the species found. Opossum class I sequences were taken from . We selected only sequences that, by their length, indicated that they were full-coding or nearly so, thereby excluding gene fragments from further processing. The panel of sequences was further reduced by manual inspection of alignments; that is, if several similar sequences representing a gene were found in the data set, only one was kept. We tried to maximize the number of taxonomic entities while keeping the dataset manageable at the same time. In order to achieve this, some taxonomic groups were left represented by a single species only, for example, rat for Glires and human for primates. Also, the number of paralogs per species was reduced in most cases, except for marsupials and monotremes. Multiple alignments were carried out with clustalw  using default settings. We calculated NJ trees using Poisson correction as substitution model with homogeneous pattern among lineages, pairwise deletion, and gamma distributed rates among sites with α = 2. We compared the NJ branching to MP and ML trees (Figures 3 and 4; Additional data files 3-6). ML trees were calculated using TreePuzzle  and visualized with TreeView . NJ and MP tree calculation and visualization were carried out with MEGA 3.1 .
Preparation of chromosomes and meiotic cells
Mitotic metaphase chromosomes were prepared from established platypus and echidna fibroblast cell lines. Primary cultures were set up from toe web from animals captured at the Upper Barnard River, New South Wales, Australia, during breeding season (AEC permit no. S-49-2006 to FG). The captured animals were killed with an intraperitoneal injection of pentobarbitone sodium (Nembutal, Boehringer Ingelheim, North Ryde, NSW, Austrialia) at a dose of 0.1 mg/g body weight. Meiotic cells and sperm were obtained by crushing the testis. The material was either directly fixed in methanol/acetic acid (3:1) or incubated in 0.075 M KCl at 37°C as hypotonic treatment and then fixed.
Fluorescence in situhybridization of BAC clones
DNA (1 μg) from BAC clones was directly labeled with spectrum orange or spectrum green (Vysis, Abbot Molecular, Des Plaines, IL, USA) using random primers and Klenow polymerase. Hybridization by FISH to platypus or echidna metaphase chromosomes and echidna meiotic chromosomes was performed under standard conditions. Briefly, the slides were treated with 100 μg/ml RNase A/2 × SSC 37°C for 30 minutes and with 0.01% pepsin in 10 mM HCl at 37°C for 10 minutes. After re-fixing for 10 minutes in 1 × PBS, 50 mM MgCl2, 1% formaldehyde, the preparations were dehydrated in an ethanol series. Slides were denatured for 2.5 minutes at 75°C in 70% formamide, 2 × SSC, pH 7.0 and again dehydrated. For hybridization of one half slide, 10 μl of probe DNA was co-precipitated with 10-20 μg of boiled genomic platypus or echidna DNA as competitor, and 50 μg salmon sperm DNA as carrier, and dissolved in 50% formamide, 10% dextran sulfate, 2 × SSC. The hybridization mixture was denatured for 10 minutes at 80°C. Pre-annealing of repetitive DNA sequences was carried out for 30 minutes at 37°C. The slides were hybridized overnight in a moist chamber at 37°C and, thereafter, washed three times for 5 minutes in 50% formamide, 2 × SSC at 42°C and once for 5 minutes in 0.1 × SSC. Chromosomes and cell nuclei were counterstained with 1 μg/ml DAPI in 2 × SSC for 1 minute and mounted in 90% glycerol, 0.1 M Tris-HCl, pH 8.0 and 2.3% DABCO. Images were taken with a Zeiss AxioImager Z.1 epifluorescence microscope equipped with a CCD camera and Zeiss Axiovision software.
The sequence data from this study have been submitted to GenBank under accession numbers EU030442-EU030444.
Additional data files
The following additional data are available with the online version of this paper. Additional data file 1 is a figure showing the characterization of echidna BAC clones on Southern blots. Additional data file 2 is a figure showing a phylogenetic tree of MHC class II genes shown in Figure 1. Additional data file 3 is a figure showing an MHC class I gene maximum likelihood phylogenetic tree. Additional data file 4 is a figure showing a MHC class I gene maximum parsimony phylogenetic tree. Additional data file 5 is a figure showing a MHC class II gene maximum likelihood phylogenetic tree. Additional data file 6 is a figure showing a MHC class II gene maximum parsimony phylogenetic tree. Additional data file 7 is a table listing the evolutionary conservation of MHC non-class I, non-class II proteins. Additional data file 8 is a table listing gene models in PIPMaker format and ORF translation for platypus BAC 462c1. Additional data file 9 is a table listing gene models in PIPMaker format and ORF translation for platypus BAC 466a15. Additional data file 10 contains Additional data file 10 is a table listing accession numbers for MHC class II genes used for phylogenetic analysis in Figure 3. Additional data file 11 is a table listing accession numbers for MHC class I genes used for phylogenetic analysis in Figure 4.
bacterial artificial chromosome
expressed sequence tag
fluorescence in situ hybridization
human leukocyte antigen
major histocompatibility complex
open reading frame
swine leukocyte antigen
We are grateful to Stefanie Palczewski and Marion Klein for their excellent technical assistance. Lisa Königsmayr and Sebastian Dohm characterized BAC clones during their summer internships in Berlin. We thank Kathy Belov and Jennifer A Marshall Graves for inspiring discussions at the beginning of the project. We are grateful to Russell Jones and Heath Ecroyd for help with the sample collection. Facilities were provided by Macquarie Generation and Glenrock Station, NSW. Approval to collect animals was granted by the New South Wales National Parks and Wildlife Services, New South Wales Fisheries, Environment and the Animal Experimentation and Ethics committee of the University of Adelaide. Cornelia Lange contributed text for the preparation of shotgun libraries. This work was supported by the Max-Planck-Gesellschaft and the Australian Research Council (Discovery Project and Fellowship schemes). Frank Grützner is an ARC Research Fellow and Enkhjargal Tsend-Ayush is an ARC Postdoctoral Research Fellow.
- Horton R, Wilming L, Rand V, Lovering RC, Bruford EA, Khodiyar VK, Lush MJ, Povey S, Talbot CC, Wright MW, et al: Gene map of the extended human MHC. Nat Rev Genet. 2004, 5: 889-899. 10.1038/nrg1489.PubMedView ArticleGoogle Scholar
- Kulski JK, Shiina T, Anzai T, Kohara S, Inoko H: Comparative genomic analysis of the MHC: the evolution of class I duplication blocks, diversity and complexity from shark to man. Immunol Rev. 2002, 190: 95-122. 10.1034/j.1600-065X.2002.19008.x.PubMedView ArticleGoogle Scholar
- Kelley J, Walter L, Trowsdale J: Comparative genomics of major histocompatibility complexes. Immunogenetics. 2005, 56: 683-695. 10.1007/s00251-004-0717-7.PubMedView ArticleGoogle Scholar
- Hurt P, Walter L, Sudbrak R, Klages S, Müller I, Shiina T, Inoko H, Lehrach H, Günther E, Reinhardt R, et al: The genomic sequence and comparative analysis of the rat major histocompatibility complex. Genome Res. 2004, 14: 631-639. 10.1101/gr.1987704.PubMedPubMed CentralView ArticleGoogle Scholar
- MHC sequencing consortium: Complete sequence and gene map of a human major histocompatibility complex. Nature. 1999, 401: 921-923. 10.1038/44853.View ArticleGoogle Scholar
- Kumanovics A, Takada T, Lindahl KF: Genomic organization of the mammalian MHC. Annu Rev Immunol. 2003, 21: 629-657. 10.1146/annurev.immunol.21.090501.080116.PubMedView ArticleGoogle Scholar
- Renard C, Hart E, Sehra H, Beasley H, Coggill P, Howe K, Harrow J, Gilbert J, Sims S, Rogers J, et al: The genomic sequence and analysis of the swine major histocompatibility complex. Genomics. 2006, 88: 96-110. 10.1016/j.ygeno.2006.01.004.PubMedView ArticleGoogle Scholar
- Liu H, Liu K, Wang J, Ma RZ: A BAC clone-based physical map of ovine major histocompatibility complex. Genomics. 2006, 88: 88-95. 10.1016/j.ygeno.2006.02.006.PubMedView ArticleGoogle Scholar
- Amadou C: Evolution of the Mhc class I region: the framework hypothesis. Immunogenetics. 1999, 49: 362-367. 10.1007/s002510050507.PubMedView ArticleGoogle Scholar
- Belov K, Deakin JE, Papenfuss AT, Baker ML, Melman SD, Siddle HV, Gouin N, Goode DL, Sargeant TJ, Robinson MD, et al: Reconstructing an ancestral mammalian immune supercomplex from a marsupial major histocompatibility complex. PLoS Biol. 2006, 4: e46-10.1371/journal.pbio.0040046.PubMedPubMed CentralView ArticleGoogle Scholar
- Ohta Y, Goetz W, Hossain MZ, Nonaka M, Flajnik MF: Ancestral organization of the MHC revealed in the amphibian Xenopus. J Immunol. 2006, 176: 3674-3685.PubMedView ArticleGoogle Scholar
- Kumar S, Hedges SB: A molecular timescale for vertebrate evolution. Nature. 1998, 392: 917-920. 10.1038/31927.PubMedView ArticleGoogle Scholar
- Kaufman J, Milne S, Gobel TW, Walker BA, Jacob JP, Auffray C, Zoorob R, Beck S: The chicken B locus is a minimal essential major histocompatibility complex. Nature. 1999, 401: 923-925. 10.1038/44856.PubMedView ArticleGoogle Scholar
- Bininda-Emonds OR, Cardillo M, Jones KE, MacPhee RD, Beck RM, Grenyer R, Price SA, Vos RA, Gittleman JL, Purvis A: The delayed rise of present-day mammals. Nature. 2007, 446: 507-512. 10.1038/nature05634.PubMedView ArticleGoogle Scholar
- Flannery TF, Groves CP: A revision of the genus Zaglossus (Monotremata, Tachyglossidae) with description of new species and subspecies. Mammalia. 1998, 62: 367-396.View ArticleGoogle Scholar
- Miska KB, Harrison GA, Hellman L, Miller RD: The major histocompatibility complex in monotremes: an analysis of the evolution of Mhc class I genes across all three mammalian subclasses. Immunogenetics. 2002, 54: 381-393. 10.1007/s00251-002-0484-2.PubMedView ArticleGoogle Scholar
- Belov K, Lam MK, Hellman L, Colgan DJ: Evolution of the major histocompatibility complex: Isolation of class II beta cDNAs from two monotremes, the platypus and the short-beaked echidna. Immunogenetics. 2003, 55: 402-411. 10.1007/s00251-003-0598-1.PubMedView ArticleGoogle Scholar
- Platypus Genome Sequencing Proposal. [http://www.genome.gov/12512287]
- Van Brink J: L'expression morphologique de la diagametie chez les sauropsids et les monotremes. Chromosoma. 1959, 10: 1-72. 10.1007/BF00396564.PubMedView ArticleGoogle Scholar
- Grützner F, Graves JA: A platypus' eye view of the mammalian genome. Curr Opin Genet Dev. 2004, 14: 642-649. 10.1016/j.gde.2004.09.006.PubMedView ArticleGoogle Scholar
- Grützner F, Rens W, Tsend-Ayush E, El-Mogharbel N, O'Brien PC, Jones RC, Ferguson-Smith MA, Marshall Graves JA: In the platypus a meiotic chain of ten sex chromosomes shares genes with the bird Z and mammal X chromosomes. Nature. 2004, 432: 913-917. 10.1038/nature03021.PubMedView ArticleGoogle Scholar
- Rens W, Grützner F, O'Brien PC, Fairclough H, Ferguson-Smith MA: Resolution and evolution of the duck-billed platypus karyotype with an X1Y1X2Y2X3Y3X4Y4X5Y5 male sex chromosome constitution. Proc Natl Acad Sci USA. 2004, 101: 16257-16261. 10.1073/pnas.0405702101.PubMedPubMed CentralView ArticleGoogle Scholar
- Bick YAE: The meiotic chain of chromosomes of Monotremata. Platypus and Echidnas. Edited by: Augee ML. 1992, Sydney: Royal Zoological Society of New South Wales, 64-68.Google Scholar
- Watson JM, Meyne J, Graves JAM: Studies of the chromosomes of the echidna meiotic translocation chain. Platypus and Echidnas. Edited by: Augee ML. 1992, Sydney: Royal Zoological Society of New South Wales, 53-63.Google Scholar
- Burge C, Karlin S: Prediction of complete gene structures in human genomic DNA. J Mol Biol. 1997, 268: 78-94. 10.1006/jmbi.1997.0951.PubMedView ArticleGoogle Scholar
- Birney E, Durbin R: Using GeneWise in the Drosophila annotation experiment. Genome Res. 2000, 10: 547-548. 10.1101/gr.10.4.547.PubMedPubMed CentralView ArticleGoogle Scholar
- Westerman M, Edwards D: DNA hybridization and the phylogeny of monotremes. Platypus and Echidnas. Edited by: Augee ML. 1992, Sydney: Royal Zoological Society of New South Wales, 28-34.Google Scholar
- Kirsch JA, Mayer GC: The platypus is not a rodent: DNA hybridization, amniote phylogeny and the palimpsest theory. Philos Trans R Soc Lond B Biol Sci. 1998, 353: 1221-1237. 10.1098/rstb.1998.0278.PubMedPubMed CentralView ArticleGoogle Scholar
- Belov K, Hellman L: Platypus IgM and the divergence of the two extant monotreme lineages. Aust Mammalogy. 2003, 25: 87-94.Google Scholar
- Springer MS, Murphy WJ, Eizirik E, O'Brien SJ: Placental mammal diversification and the Cretaceous-Tertiary boundary. Proc Natl Acad Sci USA. 2003, 100: 1056-1061. 10.1073/pnas.0334222100.PubMedPubMed CentralView ArticleGoogle Scholar
- Gibbs RA, Weinstock GM, Metzker ML, Muzny DM, Sodergren EJ, Scherer S, Scott G, Steffen D, Worley KC, Burch PE, et al: Genome sequence of the Brown Norway rat yields insights into mammalian evolution. Nature. 2004, 428: 493-521. 10.1038/nature02426.PubMedView ArticleGoogle Scholar
- Explore the Platypus Genome. [http://www.ensembl.org/Ornithorhynchus_anatinus/index.html]
- Krueger C, Osborne CS: Raising the curtains on interchromosomal interactions. Trends Genet. 2006, 22: 637-639. 10.1016/j.tig.2006.09.008.PubMedView ArticleGoogle Scholar
- Smith TP, Rohrer GA, Alexander LJ, Troyer DL, Kirby-Dobbels KR, Janzen MA, Cornwell DL, Louis CF, Schook LB, Beattie CW: Directed integration of the physical and genetic linkage maps of swine chromosome 7 reveals that the SLA spans the centromere. Genome Res. 1995, 5: 259-271. 10.1101/gr.5.3.259.PubMedView ArticleGoogle Scholar
- Deakin JE, Siddle HV, Cross JG, Belov K, Graves JA: Class I genes have split from the MHC in the tammar wallaby. Cytogenet Genome Res. 2007, 116: 205-211. 10.1159/000098188.PubMedView ArticleGoogle Scholar
- Blaschke RJ, Rappold G: The pseudoautosomal regions, SHOX and disease. Curr Opin Genet Dev. 2006, 16: 233-239. 10.1016/j.gde.2006.04.004.PubMedView ArticleGoogle Scholar
- Jeffreys AJ, May CA: Intense and highly localized gene conversion activity in human meiotic crossover hot spots. Nat Genet. 2004, 36: 151-156. 10.1038/ng1287.PubMedView ArticleGoogle Scholar
- Matthey R: Les Chromosomes des Vertébrés. 1949, Lausanne: F RougeGoogle Scholar
- White MJD: Animal Cytology and Evolution. 1973, Cambridge: Cambridge University Press, 3Google Scholar
- McQueen HA, Fantes J, Cross SH, Clark VH, Archibald AL, Bird AP: CpG islands of chicken are concentrated on microchromosomes. Nat Genet. 1996, 12: 321-324. 10.1038/ng0396-321.PubMedView ArticleGoogle Scholar
- Sambrook J, Russell DW: Molecular Cloning: A Laboratory Manual. 2001, Cold Spring Harbor: Cold Spring Harbor Laboratory Press, 3Google Scholar
- GAP4 Sequence Assembly and Contig Editing Program. [http://www.sanger.ac.uk/Software/production/staden/]
- Monodelphis domestica Resources. [http://bioinf.wehi.edu.au/opossum]
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMedPubMed CentralView ArticleGoogle Scholar
- Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18: 502-504. 10.1093/bioinformatics/18.3.502.PubMedView ArticleGoogle Scholar
- Page RDM: TREEVIEW: An application to display phylogenetic trees on personal computers. Comput Appl Biosci. 1996, 12: 357-358.PubMedGoogle Scholar
- Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform. 2004, 5: 150-163. 10.1093/bib/5.2.150.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.