<title>Background</title> <p>The evolution of the vertebrate globin superfamily has been extensively studied for many decades by comparing the structure and function of members of the gene families. These are principally haemoglobin, myoglobin, cytoglobin and neuroglobin and, more recently, globin X (in fish and amphibians [##REF##15356282##1##]) and globin Y (specific to amphibians [##REF##16484786##2##]).</p> <p>Haemoglobin genes (alpha- and beta-globin) are of particular interest because of their critical role in oxygen transportation from the respiratory surfaces to the inner organs, and because of the dire effects of mutations in human globin genes that cause haemoglobinopathies [##UREF##0##3##]. The genes contained in the alpha (α)- and beta (β)-globin clusters are expressed at different stages of development and in different tissues. Together, gene products from both clusters form the functional tetrameric haemoglobin molecules needed to fulfil oxygen requirements.</p> <p>The evolutionary history of α- and β-globin genes can be traced back to the common ancestors of fish, amphibians and amniotes (reptiles, birds and mammals), by comparing gene structure and composition of α- and β-globin clusters across vertebrates. In the amphibians <italic>Xenopus laevis </italic>and <italic>X. tropicalis</italic>, α- and β-globin genes are tightly juxtaposed as 5'-α-β-3' [##REF##16484786##2##,##REF##510792##4##, ####REF##7407927##5##, ##REF##7001356##6####7001356##6##]. In the Antarctic notothenioid fish (<italic>Notothenia coriiceps, N. angustata, Trematomus hansoni, T. pennellii</italic>), there is also a single 5'-α-β-3' locus [##REF##14516071##7##], although in pufferfish (<italic>Fugu rubripes</italic>) there are two globin clusters (one with α-globin genes and the other with both α- and β-globin genes), which are located on different chromosomes [##REF##12517812##8##].</p> <p>In amniotes, α- and β-globin clusters are located on different chromosomes. It was proposed that the ancestral α- and β-globin genes were located together in the common ancestor of amniotes, as they are in fish and amphibians, but became separated, either by chromosome fission or translocation between α- and β-genes, or by chromosome/genome or <italic>in trans </italic>duplication and gene loss [##REF##7407927##5##].</p> <p>Further duplications then occurred in amniote lineages. The ancestral α-globin gene is thought to have duplicated twice before the divergence of the bird-mammalian lineages, to produce progenitors of embryonic globin genes π/ζ, and adult α<sup>D </sup>and α<sup>A</sup>, all of which are present in birds (for example, the chicken <italic>Gallus gallus</italic>) [##REF##6273837##9##, ####REF##711749##10##, ##REF##6248855##11####6248855##11##] and mammals [##REF##16024261##12##,##REF##15855277##13##]. The order and timing of these duplications is still debated, as is their origin: for instance, α<sup>D </sup>may have evolved by duplication either of adult α<sup>A </sup>(see [##REF##16024261##12##]), or of an embryonic α-like gene [##REF##17586601##14##]. After the avian and mammalian lineages diverged, there were further tandem duplications of the π/ζ and α<sup>A </sup>lineages to produce more complex marsupial and eutherian ('placental') mammalian α-globin clusters, 5'-ζ-ψζ'-α<sup>D</sup>-ψα<sup>3</sup>-α<sup>2</sup>-α<sup>1</sup>-θ-3' (see [##REF##16024261##12##,##REF##3952001##15##, ####REF##6446404##16##, ##REF##282616##17##, ##REF##7407925##18####7407925##18##]). The timing of these duplication events is also uncertain, because we do not know whether these seven α-like globin genes all existed at the stem of the mammalian radiation.</p> <p>As for many other gene families [##REF##15892878##19##], comparisons of globin genes between distantly related mammals have provided unique insight into the evolution and function of the mammalian globins. Marsupials diverged from eutherian mammals about 148 million years ago (MYA), and mammalian Subclass Theria that contains these groups diverged from monotremes (Subclass Prototheria) about 166 MYA [##REF##17392779##20##], so comparisons between these major mammal groups provide depth for evolutionary comparisons. Monotremes retain many anatomical and developmental features shared with birds and reptiles. Their small genome, too, and disjunct chromosome size classes are reminiscent of reptile genomes, and the 10 sex chromosomes in a karyotype of 52 chromosomes is unique among mammals [##REF##18185982##21##, ####REF##15502814##22##, ##REF##15534209##23####15534209##23##]. Their importance for comparative studies is now increasingly recognised after the sequencing of the genome of a monotreme, <italic>Ornithorhynchus anatinus </italic>(platypus), to a depth of six to eight times by the Washington University Genome Centre, St Louis [##REF##18464734##24##].</p> <p>Indeed, studies of marsupial globins have clarified the timing of some of the duplications. The finding of single ε- (embryonic) and β-globin (adult) genes together in the marsupial β-globin cluster indicated that a two-gene cluster (ε-β) was present in the common therian ancestor [##REF##8265626##25##, ####REF##8752008##26##, ##REF##3375246##27##, ##REF##10381323##28####10381323##28##]. Genes in the cluster were further duplicated to produce the ancestral eutherian β-globin cluster of 5'-ε-γ-η-δ-β-3' (see [##REF##6527390##29##, ####REF##8383794##30##, ##REF##6599973##31##, ##UREF##1##32####1##32##]), which then underwent further tandem duplication events. In contrast, the bird (<italic>G. gallus</italic>) β-like globin genes (ε-β<sup>H</sup>-β<sup>A</sup>-ρ) show very little homology to the mammalian β-like globin genes [##REF##3454296##33##,##REF##8307571##34##].</p> <p>The discovery of a β-like globin gene (ω <italic>-</italic>globin) adjacent (3') to the α-globin cluster in marsupials led to a re-interpretation of globin evolution in birds and mammals [##REF##11158601##35##,##REF##15461421##36##]. Comparative sequence and phylogenetic analysis suggested that the ω-globin gene was more closely related to bird β-like globin genes than to other mammalian β-like globin genes. The specific function of the ω-globin gene is not yet known, but it is expressed just before birth and in the early stages of pouch young development [##REF##9342240##37##]. In addition, the ω-globin product binds to α-like globin chains to form functional haemoglobin, so it is likely to be involved in oxygen transportation [##REF##11158601##35##, ####REF##15461421##36##, ##REF##9342240##37####9342240##37##].</p> <p>This finding of a remnant β-like globin gene (ω <italic>-</italic>globin) beside the α-globin cluster in marsupials [##REF##11158601##35##,##REF##15461421##36##] provided some support for the alternative hypothesis [##REF##7407927##5##] that the α- and β-globin clusters in birds and mammals arose by <italic>in trans </italic>duplication of a chromosomal region, rather than simply by separation of the ancestral α-β globin cluster by chromosome fission or translocation. Wheeler et al. [##REF##11158601##35##,##REF##15461421##36##] proposed that before the divergence of birds and mammals (&gt;315 MYA), the chromosome region bearing the ancestral α-β clusters duplicated to form two clusters (α1-β1 and α2-β2) on different chromosomes, and their contents diverged independently in mammals and birds by silencing of some genes within each cluster (Figure ##FIG##0##1##). To account for the apparent orthology of the marsupial ω-globin gene and bird β-like globin genes, Wheeler et al. [##REF##11158601##35##,##REF##15461421##36##] suggested that the α1 and β2 were silenced in the eutherian lineage, but β2 was retained in marsupials as the ω-globin. In contrast, α2 and β1 were silenced in the bird lineage (Figure ##FIG##0##1##). On this hypothesis, then, both the α clusters and the β clusters of birds and mammals are paralogous (that is, evolved independently from ancient duplicates in an amniote ancestor) rather than orthologous (that is, diverged from the same ancestral cluster in an amniote ancestor).</p> <p>This paralogy hypothesis (which rests on the rather weak orthology between the chicken β and marsupial ω), as well as the dates and types of other duplications, could be further tested by studying globin genes of monotreme mammals, and using comparative data to infer the ancestral globin gene arrangement of a mammal ancestor 166 MYA. The availability of platypus genomic sequences now provides an efficient way to discover all of the globin genes and regulatory signals, and to understand their function and evolution. Studies of globin genes in monotremes are also interesting because the specialized features and lifestyle of these unique mammals may have given rise to special adaptations of globin genes to fulfil unusual oxygen requirements. These features include the need for oxygen by diffusion through the egg membrane to the embryo after birth and the physiological response to hypoxic conditions during hibernation, burrowing and diving [##UREF##2##38##, ####REF##6067227##39##, ##REF##5119901##40####5119901##40##].</p> <p>Little is known about monotreme α- and β-globin families. More than 30 years ago, studies of adult blood revealed a single adult α and β globin protein in the platypus [##REF##4464824##41##,##REF##1191127##42##] and echidna (<italic>Tachyglossus aculeatus </italic>[##REF##4798231##43##,##REF##4663350##44##]). Lee et al. [##REF##10381323##28##] later isolated an adult β-globin gene in the echidna that encoded a polypeptide identical to the previously isolated echidna β-globin [##REF##4663350##44##]. To date, there is no evidence of any monotreme embryonic ζ- or ε-globin genes.</p> <p>We used platypus genomic sequences from bacterial artificial chromosomes (BACs) to characterise the α- and β-globin gene families of the platypus and investigate their molecular evolution. In particular, we searched for embryonic and ω-globin genes and any novel globin genes that might fulfil the requirements for oxygen transport under hypoxic conditions. We investigated the genome context in order to infer the structure and origin of the ancestral α- and β-globin clusters at the stem of the mammalian radiation. Our results strongly support the hypothesis that the mammalian α- and β-globin clusters are orthologous to the avian α- and β-globin clusters, respectively, and that the β cluster evolved by transposition of a copy of the beta-like ω-globin gene in an amniote ancestor.</p>
<title>Methods</title> <title>Isolation and purification of probes to screen for platypus β-globins</title> <p>At the start of this project there were no trace sequences available for any globin genes in the platypus trace archive. We therefore designed probes to screen the platypus male Oa_Bb BAC library (Clemson University Genomic Institute, USA). The platypus β-globin-specific primers OaBGF (5'-tggacccagaggttctttgac-3') and OaBGR (5'-tgcaattcactcagcttggag-3') were designed from the reference tammar β-globin sequence [GenBank: <ext-link ext-link-type="gen" xlink:href="AY450928">AY450928</ext-link>] using Primer3 [##REF##10547847##64##]. Amplification by PCR was performed in a final volume of 25 μl, with 40 ng genomic DNA, 1× Buffer (Roche, Australia), 0.2 mM dNTPs, 0.05 U Taq (Roche, Australia) and 1 μM each of forward and reverse primers. PCR cycling conditions were: 94°C for 2 minutes, then 35 cycles of 94°C for 30 seconds, 50 to 60°C for 30 seconds, 72°C for 1 minute, followed by 72°C for 10 minutes. The PCR products were sub-cloned according to the TOPO TA cloning<sup>® </sup>Kit Protocol (Invitrogen, Australia) and the resulting plasmids were purified according to the centrifugation protocol of Wizard<sup>® </sup>Plus SV Minipreps DNA Purification System (Promega, Australia). The purified plasmids were confirmed to contain PCR products of a partial platypus β-globin gene (167 bp) by sequencing at the Australian Genome Research Facility (AGRF, Brisbane, Australia) using M13 forward (5'-gtaaaacgacggccag-3') and M13 reverse (5'-caggaaacagctatgac-3') primers. Once confirmed, they were used as probes to screen the platypus BAC library.</p> <title>Screening the platypus BAC library for β-like globin genes</title> <p>The platypus BAC library filters were pre-hybridised at 65°C with Church Buffer (1 mM EDTA, 0.5 M phosphate buffer, 7% (w/v) SDS) including 1% BSA for 4 hours. The platypus β-globin probes (25 ng) were labelled with <sup>32</sup>P-dATP using the Megaprime DNA labelling System (GE Healthcare, Australia) following the manufacturer's instructions. The probes were allowed to hybridise to the filters at stringent conditions (65°C with the above buffer) for 24 hours and then washed twice for 15 minutes each in 2 × SSC/0.1%SDS and 1 × SSC/0.1%SDS. Autoradiography was carried out for 14 days at -80°C with an intensifying cassette.</p> <title>Identification of platypus BAC clones containing α-like globin genes</title> <p>Unlike β, BACs were not screened for α-like globin genes. Instead they were identified directly from the Encyclopaedia Of DNA Elements Project [##REF##15499007##65##], in which the α-globin cluster is one of the targeted regions [##REF##16024261##12##,##UREF##5##66##]. Two platypus BAC clones (Oa_Bb-2L7 and Oa_Bb-131M24), which were sequenced but not yet annotated, were identified by computational analysis (below) to contain parts of the α-globin cluster and a ω-globin gene.</p> <title>Isolation and purification of DNA from BAC clones</title> <p>DNA from the identified BAC clones (including those that were screened) was extracted using Wizard<sup>® </sup>Plus SV Minipreps DNA Purification System (Promega, Australia). The purified BAC clones were then subjected to Dot or Southern Blot to confirm the presence of α- or β-globin genes respectively.</p> <title>Confirmation of BACs containing globin genes</title> <p>Dot blot methods were used to verify the presence of the α-like globin genes. In a plate containing Luria broth agar with chloramphenicol, a Hybond N<sup>+ </sup>(GE Healthcare, Australia) filter was placed and multiple 1 μl of liquid culture BAC clones were spotted onto the filter. The plate was incubated at 37°C overnight and then the filter was soaked in Denaturation Solution (0.5 M NaOH and 1.5 M NaCl) for 5 minutes, followed by soaking twice in Neutralisation Solution (0.5 M Tris-Cl pH 7.4 and 1.5 M NaOH) for 5 minutes each. The filter was then rinsed in 2 × SSC, soaked in 0.4 M NaOH for 20 minutes and washed with 6 × SSC to remove all cellular debris. The filters were then screened with the platypus α-globin probes using the standard library screening procedure (above).</p> <p>Southern blotting was used to verify the presence of the β-like globin genes. In a 40 μl reaction, 20 to 40 ng BAC DNA was digested with 10 U of restriction enzyme,<italic>HIND III </italic>(Roche, Australia). The reaction was incubated at 37°C for at least 4 hours and separated by electrophoresis on a 0.8% agarose gel overnight at 40 V. The DNA fragments were transferred onto a Hybond N<sup>+ </sup>(GE Healthcare, Australia) nylon filter overnight by capillary action following the manufacturer's instructions, and cross-linked in 0.4 M NaOH for 20 minutes. These filters were then screened with the platypus β-globin probes using the standard library screening procedure (above).</p> <title>Fluorescence <italic>in situ</italic> hybridisation (FISH)</title> <p>Male platypus metaphase spreads were prepared and <italic>in situ </italic>mapping was performed using two-colour FISH as described previously by McMillan et al. [##REF##18185982##21##]. The verified BACs containing the α-like globin genes (ζ and ζ': Oa_Bb-2L7) and β-like globin genes (ε and β: Oa_Bb-484F22) were labelled with different fluorochromes and then hybridised to the chromosomes. The signals were detected by fluorescent microscopy, where at least twenty metaphase images were captured and analysed.</p> <title>Sequence data of the platypus BAC clones containing the α- and β-globin clusters</title> <p>Information about the platypus BAC clones containing the α-like globin genes along with the ω-globin gene was obtained directly from the ENCODE Project [##UREF##5##66##]. Their sequence information was obtained from GenBank; accession numbers: <ext-link ext-link-type="gen" xlink:href="AC195438">AC195438</ext-link> (Oa_Bb-2L7) and <ext-link ext-link-type="gen" xlink:href="AC203513">AC203513</ext-link> (Oa_Bb-131M24).</p> <p>The BAC clone containing the β-like globin genes that were found from the library screening procedure were sequenced at the Washington University Genome Sequencing Centre (St Louis, USA). The sequence information for this BAC clone was obtained from GenBank: <ext-link ext-link-type="gen" xlink:href="AC192436">AC192436</ext-link> (Oa_Bb-484F22).</p> <title>Computational characterisation of the α- and β-globin clusters in the platypus</title> <p>Using sequence information of <ext-link ext-link-type="gen" xlink:href="AC195438">AC195438</ext-link>, <ext-link ext-link-type="gen" xlink:href="AC203513">AC203513</ext-link> and <ext-link ext-link-type="gen" xlink:href="AC192436">AC192436</ext-link>, genes were predicted by GENSCAN [##REF##9149143##46##] and GenomeScan [##REF##11337476##47##] using default settings. All predicted gene sequences were then subjected to a BLAST search of the translated nucleotide acid (BlastX) and protein (BlastP) databases to confirm their identities.</p> <title>Promoter analyses</title> <p>Transcription factor binding motifs were predicted in the 200 bp promoter region located 5' to the predicted platypus α- and β-like genes and <italic>GBY </italic>by rVista 2.0 [##REF##15215384##67##] using user-defined consensus sequences for 'CACCC', 'CAAT', 'TATA', GATA1 ('WGATAR' [##REF##8262042##51##]) and EKLF ('NGNGTGGGN' [##REF##8262042##51##]). The same criteria were used to predict the same motifs in marsupials (<italic>Didelphis virginiana </italic>[ζ and ψζ': <ext-link ext-link-type="gen" xlink:href="AC139599">AC139599</ext-link>] and <italic>Sminthopsis macroura </italic>[α<sup>D</sup>, ψα<sup>3</sup>, α<sup>2</sup>, α<sup>1</sup>, ω: <ext-link ext-link-type="gen" xlink:href="AC146781">AC146781</ext-link>; and ε, β: <ext-link ext-link-type="gen" xlink:href="AC148754">AC148754</ext-link>]) and in humancs [ζ, ψζ', α<sup>D</sup>, ψα<sup>3</sup>, α<sup>2</sup>, α<sup>1</sup>: <ext-link ext-link-type="gen" xlink:href="NG_000006">NG_000006</ext-link>; and ε, β: <ext-link ext-link-type="gen" xlink:href="NG_000007">NG_000007</ext-link>] for consistency in comparison.</p> <title>Confirmation of α1 and α3 by BLAST search and Southern blot</title> <p>To confirm that the presence of two almost identical genes (α<sup>1 </sup>and α<sup>3</sup>) was real rather than an assembly error, the boundaries (~300 bp) of the homologous regions were investigated by a BLAST search against the platypus WGS database. The raw sequences of best hits were extracted from NCBI GenBank, cleaned and aligned in Sequencher v4.8 (Gene Codes Corporation, Michigan) using default settings.</p> <p>Southern blotting was also used to verify the presence of α<sup>1 </sup>and α<sup>3 </sup>genes. In a 30 μl reaction, 100 μg BAC DNA (Oa_Bb: 131M24, 130N2, 150K14 and 223I12) was digested with 10 U of restriction enzyme,<italic>EcoRV </italic>(Roche, Australia). The reaction was incubated at 37°C for at least 4 hours and separated by electrophoresis on a 0.8% agarose gel overnight at 40 V. The DNA fragments were transferred onto a Hybond N<sup>+ </sup>(GE Healthcare, Australia) nylon filter overnight by capillary action following the manufacturer's instructions, and cross-linked in 0.4 M NaOH for 20 minutes. These filters were then screened with the platypus α<sup>1</sup>/α<sup>3 </sup>(test) and α<sup>2 </sup>(control) probes using the standard library screening procedure (above).</p> <title>RT-PCR analyses</title> <p>To remove DNA contamination, RNAs derived from adult male platypus liver, kidney, spleen, testis, brain and lungs were DNase treated using a DNA-<italic>free</italic>™ kit according to the manufacturer's instructions (Applied Biosystems, Australia). Treated RNAs were then reverse transcribed using Superscript III (Invitrogen, Australia) following the manufacturer's instructions. Primers were designed against predicted α- and β- like and <italic>GBY </italic>globin gene sequences using Primer3 [##REF##10547847##64##]. In each case, the region amplified spanned an intron so that the origin of the template (gDNA or cDNA) was immediately obvious. Primer sequences and the expected sizes of amplified cDNA and gDNA bands are shown in Table ##TAB##1##2##. PCR reactions and cycling conditions were the same as for screening for the platypus β-globin genes, above. The positive bands were directly sequenced by AGRF (Brisbane, Australia) to confirm their identities. The blood contamination in the tested samples had minimal effect on the observed expression pattern, as some tissues (for example, lung and liver) showed no amplification despite containing large quantities of blood.</p> <title>Phylogenetic analyses</title> <p>Phylogenetic analyses were employed to verify the identities of the platypus globin genes and study the evolutionary relationships of the different members of the α- and β-globin gene families. This study was restricted to the coding domains of the α- and β-globin members and the accession numbers of the sequences used are given in the legends of Figures ##FIG##1##2## and ##FIG##2##3##. Phylogenetic analyses were conducted using MP in PAUP* v.4.0b10 [##REF##12504223##68##], and a BI approach using MrBayes v.3.1.2 [##REF##11524383##69##]. Concordance of trees from each of the different methods, bootstrap proportions and posterior probability estimates were used to examine the robustness of nodes.</p> <p>MP analyses were conducted for the entire coding sequence matrix and after excluding third codon positions using a heuristic search option and default options (TBR branch swapping), with the exception of using random stepwise addition repeated 100 times. Character state optimisation for MP trees used the DELTRAN option. MP bootstrap analyses [##UREF##6##70##] were carried out using 1000 bootstrap pseudoreplicates, employing a heuristic search option with random stepwise addition.</p> <p>The program MODELTEST [##REF##9918953##71##] and the Akaike information criterion (AIC) were used to assess the most appropriate model for BI analyses. The MODELTEST analyses were facilitated using the program MrMTgui v1.0 [##UREF##7##72##]. The MODELTEST analysis was carried out on separate codon positions for α- and β-globin data sets. For α-globin sequences, a general time reversible (GTR) model [##REF##2338834##73##], with a proportion of invariant sites (I) and unequal rates among sites [##UREF##8##74##], modelled with a gamma distribution (G) was found to be the most appropriate model to use for first and second codon positions, and a GTR+G model was appropriate for third codon positions under the AIC. For β-globin sequences a GTR+I+G model was considered appropriate for first positions, and a GTR+G model was found to be appropriate for second and third codon positions. The MrBayes analysis was carried out applying these different models to each codon position using an unlinked analysis, with default uninformative priors. Four chains were run simultaneously for 2 million generations in two independent runs, sampling trees every 100 generations. The program TRACER (version 1.3; [##UREF##9##75##]) was used to assess tree and parameter convergence. For both the α-globin and β-globin analyses all effective sample sizes for all parameters were &gt;1297, indicating a sufficient sample of the parameter space had been taken. A burn-in of 2000 trees (equivalent to 200,000 generations) was chosen for each independent run of MrBayes, with a &gt;50% posterior probability consensus tree constructed from the remaining 36,002 trees (18,001 trees each run).</p> <p>A BI analysis using MrBayes (version 3.1.2) was also carried out using protein sequence data from β-globin genes. A mixed protein model was used, allowing the optimum model of protein evolution to be assessed from a selection of nine fixed-rate models. The optimum model was found to be the Dayhoff model with a posterior probability of 1.0. The analyses were conducted using two million generations in two independent runs, sampling trees every 100 generations. A burn-in of 2,000 trees was used for each run with a 50% consensus tree constructed from the remaining 36,002 trees.</p>
<title>Results</title> <title>Identification of BAC clones containing the α- and β-globin clusters</title> <p>The draft sequence assembly of platypus [##REF##18464734##24##] is readily available on the University of California Santa Cruz (UCSC) Genome Browser [##UREF##3##45##]. However, currently the assembly is incomplete for the α- and β-globin clusters, as individual globin genes appear on different contigs. There are also sequences of the platypus BAC clones available in NCBI GenBank that are not yet annotated and assembled, nor is part of the platypus genome assembly. Two of these are Oa_Bb-2L7 [GenBank:<ext-link ext-link-type="gen" xlink:href="AC195438">AC195438</ext-link>] and Oa_Bb-131M24 [<ext-link ext-link-type="gen" xlink:href="AC203513">AC203513</ext-link>], which were identified from the Encyclopaedia of DNA Elements Project to contain parts of the α-globin cluster (see Methods). The BAC clone Oa_Bb-484F22 [GenBank: <ext-link ext-link-type="gen" xlink:href="AC192436">AC192436</ext-link>] containing the β-globin cluster was obtained by screening a male platypus BAC library (Clemson University Genomic Institute, USA) and was subsequently fully sequenced and assembled by the Washington University Genome Sequencing Centre (St Louis, USA). These sequences were therefore used in this study to characterise the whole α- and β-globin clusters in the platypus.</p> <p>Genes in these sequenced BAC clones were predicted by programs GENSCAN [##REF##9149143##46##] and GenomeScan [##REF##11337476##47##]. Many genes were predicted, which were then used for BLAST searches of nucleotide (BlastN) and amino acid (BlastP) databases to help identify them (data not shown). Phylogenetic analyses were also conducted for the platypus α- and β-like globin genes to further verify the identity of each gene (see below and also Figures ##FIG##1##2##, ##FIG##2##3## and ##FIG##3##4## below). With only one exception (platypus ε-globin, see below), the identities of all of the genes inferred by BLAST analyses were supported by phylogenetic analyses with high posterior probabilities and bootstrap support values.</p> <title>Predictions and characterisation of genes in the platypus α-globin cluster</title> <p>One BAC (Oa_Bb-2L7) contained two embryonic α-like globin genes, and a second BAC (Oa_Bb-131M24) contained six α-like globin genes and a β-like globin gene (see Additional file ##SUPPL##0##1##). These two BACs were found to overlap by 10,066 base pairs (bp), resulting in a contig of 330,126 bp that contained the entire platypus α-globin cluster and flanking genes.</p> <p>The 330,126 bp α-globin contig was found to contain six α-like globin genes, a β-like globin gene, and a gene that bore little similarity to α- and β-like globin genes but some similarity to cytoglobins (Figure ##FIG##4##5A##). These six α-like globin genes have a three-exon/two-intron structure and conserved donor/acceptor splice sites (GT/AG) typical of all vertebrate α-like globin genes. They are separated from each other by 2 to 6 kilobase pairs (kb). Full details of the exon/intron lengths, location of the putative poly-A addition site (AATAAA) and the lengths of the coding domains with the predicted encoded polypeptide for each predicted gene are given in Table ##TAB##0##1##. Figure ##FIG##4##5B## shows the predictions for some of the well-characterised protein-binding sites in the 5' promoter region (about 200 bp 5' to the cap site of each gene). These include CACCC [##REF##8428663##48##], CAAT [##REF##2376599##49##], TATA [##REF##9642099##50##], GATA 1 [##REF##8262042##51##], EKLF (Erythroid Krüppel-like Factor; [##REF##9744863##52##]) and have been experimentally shown to control the stage- and tissue-specific expression of α- and β-like globin genes in other mammals [##REF##9642099##50##,##REF##17042767##53##, ####REF##8161819##54##, ##REF##9461381##55####9461381##55##].</p> <p>Two genes at the 5' end of the α-globin cluster were both identified as ζ-like (referred to here as ζ and ζ') and predicted to encode polypeptides of 142 amino acids (aa), which are typical of known functional mammalian α-like globin genes. The amino acid sequence alignment of ζ and ζ' shows 95% identity. In the promoter region of both genes, CACCC and CAAT consensus boxes are conserved at similar positions, and in comparable order to that of human ζ and ζ' (Figure ##FIG##4##5B##).</p> <p>Adjoining the two ζ-like globin genes, four other α-like globin genes were identified. One was an orthologue of bird and reptilian α<sup>D</sup>, and the other three were orthologues of adult α genes (here called α<sup>3</sup>, α<sup>2 </sup>and α<sup>1</sup>). The long and uninterrupted open reading frame (ORF) of α<sup>D </sup>strongly suggests that it encodes a functional polypeptide of 141 aa, typical of known functional α<sup>D </sup>globin genes. The platypus α<sup>D </sup>globin gene contains introns of 1450 bp (intron 1) and 1610 bp (intron 2) that are very large compared with those of other α-like globins, which are usually less than 1000 bp.</p> <p>Analyses of the platypus adult α-like globin genes reveal three adult (α<sup>3</sup>, α<sup>2 </sup>and α<sup>1</sup>) globin genes in the α-globin cluster. The sequence of α<sup>3 </sup>(the most 5' gene, adjacent to α<sup>D</sup>) was found to be almost identical to α<sup>1 </sup>(the most 3' gene) in their exon and intron regions, as well as in flanking regions of about 130 bp on both sides. The coding region was 100% identical, and just two sites in intron 1 were found to be different between the two genes. In order to confirm that identification of these two identical genes was not due to an error in the assembly of the original sequence data, the boundaries of the region containing the homology between α<sup>1 </sup>and α<sup>3 </sup>was further analysed by a BLAST search of the platypus whole-genome shotgun (WGS) database (data not shown). Two contigs were identified with homology to α<sup>1 </sup>and α<sup>3</sup>; these had identical sequences on one side of the boundary but different sequences on the other, confirming the presence of two separate genes. Further confirmation was obtained by performing a Southern blot on the α-globin-containing BACs, digested with an enzyme (<italic>EcoRV</italic>) that does not cut within the α<sup>1</sup>, α<sup>2</sup>and α<sup>3 </sup>(data not shown). Probing with α<sup>1</sup>/α<sup>3 </sup>revealed two bright bands, corresponding to α<sup>1 </sup>and α<sup>3</sup>, and one fainter band between them, corresponding to α<sup>2</sup>. Probing with α<sup>2 </sup>produced the same three bands, but in this case the middle one was brighter, corresponding to α<sup>2</sup>, and the outer bands were fainter, corresponding to α<sup>1 </sup>and α<sup>3</sup>. These analyses confirmed the existence of separate genes α<sup>1 </sup>and α<sup>3 </sup>in the platypus α-globin cluster. The α<sup>2 </sup>gene, located between α<sup>1 </sup>and α<sup>3</sup>, was distinct from both genes in the coding sequence (with 83% homology), in intron lengths (intron 1: 405 bp in α<sup>1</sup>/α<sup>3 </sup>and 720 bp in α<sup>2</sup>; intron 2: 151 bp in α<sup>1</sup>/α<sup>3 </sup>and 155 bp in α<sup>2</sup>) and in the promoter region (Figure ##FIG##4##5B##).</p> <p>The amino acid sequence encoded by α<sup>1 </sup>and α<sup>3 </sup>was identical to the platypus adult α-chain previously identified by Whittaker and Thompson [##REF##4464824##41##], implying that at least one of these genes is expressed in the adult platypus. The coding domain of α<sup>1 </sup>and α<sup>3 </sup>is shorter (426 bp) than that of α<sup>2 </sup>(429 bp), because it lacks the first three nucleotides of exon 1. The ORF of α<sup>2 </sup>gives a strong indication that it is translated into a functional polypeptide of 142 aa, typical of known functional mammalian α-like globin genes.</p> <p>On the 3' side of the six α-like globin genes, a β-like globin gene was predicted, which was identified as the orthologue of the marsupial ω-globin gene. This platypus ω-globin gene has a typical three-exon/two-intron structure, conserved donor/acceptor splice sites, and encodes a polypeptide of 146 aa, typical of all vertebrate β-like globin genes (Table ##TAB##0##1##). The promoter region located 5' of the ω-globin initiation codon contains conserved sites for CAAT-EKLF-CACCC in an order identical to that of marsupial ω-globin gene.</p> <p>Unexpectedly, GenomeScan predicted a gene based on the protein similarities with the α- and β-polypeptide chains, approximately 1.5 kb 3' of the ω-globin gene. Like other α- and β-globins, this gene also has a three-exon/two-intron structure and conserved donor/acceptor splice sites (Table ##TAB##0##1##). The lengths of its exons 1, 2, and 3 are 98, 223 and 144 bp, respectively, compared with 92, 223 and 129 bp in other β-like globin genes. However, it has much larger introns of 3364 bp (intron 1) and 3053 bp (intron 2). The long and uninterrupted ORF of this gene can be translated into a polypeptide of 154 aa, which is atypical of any known α- or β-like globin genes. A BLAST search of the amino acid sequence of this gene obtained the best hit with Globin Y (<italic>gby</italic>) of the amphibian <italic>X. laevis </italic>(identity score of 39%), and weaker identity scores with Cytoglobins (<italic>cygb</italic>) of other species, such as the fish <italic>Danio rerio </italic>(27%), <italic>X. tropicalis </italic>(26%), chicken (28%) and human (25%) at the protein level. We designated this gene '<italic>GBY</italic>' based on similarities with <italic>X. laevis gby</italic>, and its similar position adjoining the globin cluster [##REF##16484786##2##]. The predicted polypeptide of platypus <italic>GBY </italic>(154 aa) was shorter than <italic>X. laevis gby </italic>(156 aa), and quite different from <italic>X. laevis cygb </italic>(179 aa), <italic>D. rerio cygb1 </italic>(174 aa) and <italic>cygb2 </italic>(179 aa), and human <italic>CYGB </italic>(190 aa). Using the Expressed Sequence Tag (EST) database, a BLAST search of the platypus <italic>GBY </italic>also obtained an identity score of 38% with <italic>X. tropicalis gby </italic>that was expressed in both tadpoles and adults, but produced no significant matches with any other mammalian genes. The present work was the first opportunity to analyse the promoter region of any <italic>GBY </italic>gene (Figure ##FIG##4##5B##).</p> <title>Predictions and characterisation of genes in the platypus β-globin cluster</title> <p>In the platypus, only two β-like globin genes were predicted within the 129,521 bp BAC clone (Oa_Bb-484F22) by GENSCAN and GenomeScan (see Additional file ##SUPPL##0##1##). When the predicted amino acid sequences were subjected to BLAST search, the 5' gene had best hits with mammalian embryonic ε-globin genes. Although the phylogenetic analyses using Bayesian inference (BI; see below) indicated that this gene was more closely related to the platypus and echidna adult β-globin genes than to therian ε-globin genes, the position of this gene on the 5' end of the β-globin cluster and expression data (see below) supports its orthology with mammalian embryonic ε-globin genes, and is henceforth referred to as ε. The 3' gene encoded a protein identical to the previously identified platypus adult β-chain [##REF##1191127##42##], and is henceforth referred to as β.</p> <p>Both genes encode polypeptides of 146 aa, typical of known functional mammalian β-like globin genes. The promoter region of the platypus β has conserved sites of CACCC and CAAT in all three extant of mammals. However, the promoter region of the platypus ε appears to be quite different from other mammalian ε-globin genes and even from the platypus β (Figure ##FIG##4##5B##). The promoter of platypus ε contains only one predicted motif (CAAT), whereas the promoters of other mammalian ε, β and the platypus β contain many predicted motifs.</p> <title>Expression studies of the platypus α- and β-like globin genes</title> <p>Transcription studies were performed to gain insight into the expression and function of all of the predicted platypus globin genes. Adult liver, kidney, spleen, testis, lung and brain were obtained for this project: no embryonic samples were available (or are ever likely to be available) for this vulnerable and iconic species. Observation of the expression of any of the predicted genes would constitute a good indication that the gene is transcriptionally active and functional.</p> <p>Reverse-transcriptase polymerase chain reaction (RT-PCR) of all predicted platypus genes showed that they are all expressed in at least some of these adult platypus tissues (Figure ##FIG##5##6##). Platypus genes α<sup>1</sup>/α<sup>3</sup>, α<sup>2 </sup>and β, whose orthologues are usually expressed in the bone marrow of an adult human, were expressed in almost all platypus tissues tested, suggesting a broader expression of these genes in the monotreme lineage. Surprisingly, the genes ζ, ζ' and ε, whose therian orthologues are expressed only at embryonic stages of development, were expressed in adult spleen and testis, but not in the other tissues of adult platypus. This suggests that persistent expression of these genes in some adult tissues was selected for in the platypus, perhaps in response to its aquatic lifestyle and the hypoxic conditions of a confined burrow. Also, the expression pattern of platypus ε is similar to embryonic α-like ζ and ζ' but different from that of adult globin genes (α<sup>1</sup>/α<sup>3</sup>, α<sup>2 </sup>and β). The ω and α<sup>D </sup>globin genes, whose functions are unknown, were also expressed mainly in the spleen. <italic>GBY </italic>was expressed in all adult platypus tissues, most strongly in testis.</p> <title>Phylogenetic analyses</title> <p>Phylogenetic analyses of the α-like globin genes using BI and maximum parsimony (MP) produced several noteworthy results. The platypus adult α globin genes (α<sup>1</sup>/α<sup>3 </sup>and α<sup>2</sup>) grouped closely together to the exclusion of eutherian and marsupial α- and θ-globin genes for all analyses, although posterior probability (69%) and bootstrap support (66%) for this arrangement were relatively weak (Figure ##FIG##1##2##). This finding suggests that the duplication leading to the marsupial and eutherian θ-globin lineage occurred after the divergence of the monotreme and therian lineages. This is consistent with the absence of a θ-globin gene from the region between platypus α<sup>1</sup>- and ω-globin, its expected location based on its position in marsupial α-globin clusters [##REF##16024261##12##,##REF##15983873##56##].</p> <p>Both platypus ζ-globin genes grouped closely together and formed a sister group relationship with chicken π, supported by a high posterior probability of 97% (Figure ##FIG##1##2##). A sister group relationship was also found in MP trees for analyses of the entire platypus coding region (bootstrap support &lt;50%), and when third positions in the codon were excluded, was supported by 73% bootstrap pseudoreplicates (data not shown). This differs from the expectation that platypus ζ-globin genes would group with other mammalian ζ-globin genes to the exclusion of chicken π, suggesting that other factors (for example, purifying selection) operated to maintain a similar sequence in birds and monotremes.</p> <p>There is still considerable uncertainty in the phylogenetic position of the α<sup>D</sup>-globin clade. It has recently been proposed that the α<sup>D </sup>globin lineage resulted from duplication of the embryonic α-globin lineage, with phylogenetic analyses supporting a sister lineage relationship of these lineages to the exclusion of the adult α-globin lineage [##REF##17586601##14##]. However, this arrangement was not supported in BI analyses of the data set used here, and the position of the α<sup>D </sup>lineage was different in the different analyses. Analyses using BI (Figure ##FIG##1##2##) supported the sister lineage relationship of the α<sup>D </sup>and adult α-globin lineages (as proposed by Cooper et al. [##REF##16024261##12##]), with 87% posterior probability support. In contrast, all MP analyses supported the sister lineage status of α<sup>D </sup>and embryonic α-globin genes, indicating an uncertainty in the phylogenetic position of the α<sup>D</sup>-globin clade.</p> <p>Phylogenetic analyses of the β-globin genes provided results similar to recently reported phylogenetic analyses [##REF##11158601##35##,##REF##15461421##36##], with one notable exception. The BI analyses of coding sequence data (Figure ##FIG##2##3##) provided strong support (99% posterior probability) for the sister relationship of bird and mammalian β-like globin genes, contradicting previously published phylogenies of mammalian β-globin genes showing a sister relationship of marsupial ω-globin and bird β-like globin genes [##REF##11158601##35##,##REF##15461421##36##]. MP analyses (Figure ##FIG##3##4##), excluding third position in the codon, gave a similar tree arrangement, albeit with very low bootstrap support (&lt;50%). In marked contrast to the BI analyses of DNA sequence data, BI protein analyses (data not shown) supported the sister relationship of bird β-like globin and mammal ω-globin lineages with a high posterior probability (99%).</p> <p>Lastly, phylogenetic analyses using BI indicated that the platypus ε gene was more closely related to the platypus and echidna adult β-globin genes than to therian ε-globin genes, suggesting it may not be orthologous to marsupial and eutherian ε-globin (Figure ##FIG##2##3##). BI analyses of β-globin protein data and MP analyses of the coding sequence data, with third codon positions excluded, grouped the gene as an ancestral lineage to eutherian and monotreme adult β-globin genes (see Figure ##FIG##3##4##). This ancestral position suggests that the lineage evolved following duplication of an ancestral β-globin gene prior to the divergence of monotremes and therians.</p> <title>Location of the α- and β-globin clusters in the platypus</title> <p>The location of the verified BAC clones containing the α- (Oa_Bb-2L7) and β-globin (Oa_Bb-484F22) clusters in the platypus was determined by fluorescence <italic>in situ </italic>hybridisation (FISH) (Figure ##FIG##6##7##). The β-globin cluster localised to one of the largest autosomes, giving unambiguous signals on the long arm of chromosome 2 (2q5.1). The α-globin cluster localised to the smallest autosome, 21, whose two arms are not distinguishable by size or DAPI banding pattern [##REF##18185982##21##]. This is the first gene that has been localised on the platypus chromosome 21.</p> <title>Loci flanking the α- and β-globin clusters in the platypus and other vertebrates</title> <p>To explore the genome context of the α- and β-globin clusters in the platypus and other vertebrates, the platypus BAC sequences and the genomes of other sequenced species were searched for loci residing beside the α- and β-globin clusters.</p> <p>As well as globin genes, GENSCAN predicted within the platypus α-globin 330,126 bp contig many genes that flank the platypus α-globin cluster (Figure ##FIG##4##5A##), which were identified by BLAST analyses. These include <italic>IL9RP3-POLR3K</italic>-<italic>C16orf33</italic>-<italic>C16orf8-MPG-C16orf35 </italic>upstream (5') of the α-globin cluster, and, <italic>LUC7L-ITFG3-RGS11</italic>-<italic>ARHGDIG</italic>-<italic>PDIA2</italic>-<italic>AXIN1 </italic>downstream (3') of the α-globin cluster (Figure ##FIG##4##5A##).</p> <p>To compare the α-globin flanking loci of the platypus and other vertebrates, the genes closest to the α-globin cluster, <italic>MPG</italic>, <italic>C16orf35 </italic>and <italic>LUC7L </italic>were searched for in the human, opossum (<italic>Monodelphis domestica</italic>), chicken, frog (<italic>X. tropicalis</italic>) and zebrafish (<italic>D. rerio</italic>) genomes that were accessible from Ensembl [##UREF##4##57##]. Figure ##FIG##7##8A## shows that the locations of <italic>MPG</italic>, <italic>C16orf35 </italic>and <italic>LUC7L </italic>are conserved adjacent to the α-globin cluster of birds and mammals, and in the same position adjacent to the α-β cluster of amphibians, and all but <italic>LUC7L </italic>were also present in fish. These results are consistent with the previous analyses of Flint et al. [##REF##11157800##58##] and Hughes et al. [##REF##15998734##59##]. Thus the flanking loci analyses reveal that the genome context of the platypus α-globin cluster is the same as the α-globin clusters in therian mammals and birds, and this is the same as for the α-β cluster of fish and frogs.</p> <p>GENSCAN also predicted numerous genes other than globin genes in the platypus β-globin BAC (484F22). These were identified by a BLAST search as members of the olfactory receptor gene (<italic>ORG</italic>) family that are responsible for odour detection. Three conserved <italic>ORG </italic>members were identified at the 5' end of the platypus β-globin cluster and one conserved <italic>ORG </italic>member at the 3' end (Figure ##FIG##4##5A##).</p> <p>To compare β-globin flanking loci, <italic>ORG </italic>genes, as well as other genes that are closest to the β-globin cluster in other species, <italic>RRM1</italic>, <italic>CCKBR </italic>and <italic>ILK </italic>were searched for in the human, opossum, chicken and zebrafish genomes that were accessible from Ensembl [##UREF##4##57##]. Data from frog (<italic>X. tropicalis</italic>) was not useful since all of these loci lie on different contigs or scaffolds due to assembly problems. The locations of multiple <italic>ORG </italic>genes, <italic>RRM1</italic>, <italic>CCKBR </italic>and <italic>ILK </italic>were found to be conserved adjacent to β-globin cluster of birds and mammals [##REF##16554811##60##,##REF##10220430##61##], but not for the α-β cluster of fish and frogs, nor beside the second α-β cluster of zebrafish and pufferfish (Figure ##FIG##7##8B##). Thus the genome context of the platypus β-globin cluster is the same as in therian mammals and birds, but this is different from the α-β cluster of fish and frogs.</p>
<title>Discussion</title> <p>The phylogenetic position of monotremes makes comparisons with platypus of special value for exploring the organization, function and evolution of mammalian genes and genomes. The availability of platypus genome sequence data now makes many such studies possible, and have been used here to characterise the platypus α- and β-globin gene clusters and explore their evolutionary history.</p> <title>The platypus α-globin gene cluster</title> <p>The platypus α-globin cluster contains at least eight genes within more than 40 kb, including six α-like globin genes (including the identical α<sup>1 </sup>and α<sup>3</sup>), one β-like globin gene (ω-globin) and a gene belonging to another member of the globin super-family (<italic>GBY</italic>) arranged in the order 5'-ζ-ζ'-α<sup>D</sup>-α<sup>3</sup>-α<sup>2</sup>-α<sup>1</sup>-ω-<italic>GBY</italic>-3' (Figure ##FIG##4##5A##). The cluster maps to chromosome 21, the smallest autosome in platypus. All eight genes are likely to be functional since their expression was detected in tissues of an adult platypus.</p> <p>The platypus α-globin cluster is almost identical to the arrangement of α-like globin genes in the ancestral therian cluster reported by Cooper et al. [##REF##16024261##12##]. The one exception is the absence of a θ-globin gene from the platypus cluster. Phylogenetic analyses support the basal position of the monotreme adult α-globin lineage relative to marsupial and eutherian α- and θ-globin lineages, implying that the duplication of an adult α-globin to produce θ-globin occurred in the therian lineage after its divergence from the monotreme lineage (Figure ##FIG##8##9B##). However, although the numbers and arrangements of genes is so similar in platypus and therians, the presence of three adult α-globin genes and two embryonic ζ-globin genes in their common ancestor was not supported by phylogenetic analyses, which showed independent groupings of the three adult and embryonic genes within each separate mammalian lineage (Figure ##FIG##1##2## and see Cooper et al. [##REF##16024261##12##]). This result can be interpreted literally as resulting from independent duplications in each mammalian lineage to produce three adult and two embryonic genes in each. However, this seems unlikely to explain the convergence in gene number of the α-globin cluster in these distantly related mammalian lineages. We suggest that a more parsimonious explanation is that the common ancestor of monotremes and therians contained three adult α-globin genes and two ζ-globin genes, which were homogenised by ongoing gene conversion events, leading to the gene tree that does not match the duplication history of the individual genes. The close similarity of the platypus α<sup>3 </sup>and α<sup>1 </sup>loci suggests a very recent gene conversion event that homogenised their sequences. Therefore, we propose that the platypus α-globin cluster of eight genes (ζ-ζ'-α<sup>D</sup>-α<sup>3</sup>-α<sup>2</sup>-α<sup>1</sup>-ω-<italic>GBY</italic>) represents the ancestral mammalian α-globin cluster arrangement (Figure ##FIG##8##9B##), in which all genes were transcriptionally active.</p> <p>Importantly, the platypus α-globin cluster contains a copy of the β-like ω-globin gene, also found in the marsupial α-globin cluster, but absent in humans, supporting the hypothesis that ω-globin was present in the common ancestor of all mammals. Phylogenetic analyses also confirm the ancient ancestry of the ω-globin gene, as concluded by Wheeler et al. [##REF##11158601##35##,##REF##15461421##36##]. Among adult platypus tissues this gene was expressed only in the spleen. In marsupials, expression of the ω-globin gene was detected just prior to birth and during early pouch young development [##REF##9342240##37##], although the site of expression was not studied, and there was no evidence of adult expression in blood cells.</p> <title>Discovery of a mammalian GBY globin gene adjoining the α-globin cluster</title> <p>We discovered a globin gene <italic>GBY </italic>in the platypus that is adjacent (3') to ω in the α-globin cluster. It has a typical three-exon/two-intron structure like other α/β-globin genes, contains an ORF encoding a polypeptide chain of 154 aa, and is expressed in almost all adult tissues, most strongly in testis. The amino acid sequence is unrelated to any of the other globin genes in the cluster, so it is unlikely to be derived by duplication of α- or ω-globin within the monotreme lineage. Rather, it shows sequence similarity to <italic>gby </italic>of <italic>X. tropicalis </italic>and <italic>X. laevis</italic>, a gene thought to be related to cytoglobins [##REF##16484786##2##].</p> <p>Little is known of the function of amphibian <italic>gby</italic>, or its relationship with other globins. Fuchs et al. [##REF##16484786##2##] reported that amphibian <italic>gby </italic>encodes a <italic>bona fide </italic>globin of 156 aa, having all of the sequence features of a functional respiratory protein. <italic>gby </italic>was expressed in all adult tissues tested in <italic>X. laevis</italic>, most strongly in ovary, kidney and eye, and was present in 20 expressed sequence tag clones from different stages of <italic>X. laevis </italic>and <italic>X. tropicalis </italic>embryonic and adult development [##REF##16484786##2##], suggesting that it is expressed in embryonic as well as adult stages. Phylogenetic analysis of all vertebrate globins [##REF##16484786##2##] showed that the <italic>gby </italic>lineage diverged at the base of two separate clades, one comprising all vertebrate cytoglobins, myoglobins, agnathan globins and bird globin E, and the other comprising the haemoglobin α- and β-chains.</p> <p>The position of platypus <italic>GBY </italic>adjacent to the α-globin cluster and flanked by <italic>LUC7L </italic>mirrors its position in <italic>X. tropicalis </italic>between the main α-β cluster and <italic>LUC7L </italic>[##REF##16484786##2##]. Another common feature of both was strong expression in gonads (ovary in <italic>X. laevis </italic>[##REF##16484786##2##] and testis in platypus), so <italic>GBY </italic>has sex-related expression in both lineages. Thus <italic>GBY </italic>is not specific to amphibians, as was thought, but was a component of the cluster in an ancient tetrapod, and has been lost, or has diverged beyond recognition, in birds and therian mammals.</p> <title>The platypus β-globin gene cluster</title> <p>Characterisation of the platypus β-globin cluster revealed two β-like globin genes over about 13.2 kb that are arranged in the same order as marsupials, 5'-ε-β-3' (Figure ##FIG##4##5A##). This cluster is located on platypus chromosome 2q5.1. Both genes appear to be transcriptionally active and are likely to be functional.</p> <p>At the time of revising this paper, an independent paper on monotreme β-like globin genes was published by Opazo et al. [##REF##18216242##62##] in which they reported the presence of ω, ε<sup>P </sup>and β<sup>P </sup>in the platypus. Largely on the basis of phylogenetic analyses of flanking and coding sequence data, they proposed that platypus ε<sup>P </sup>and β<sup>P </sup>were not 1:1 orthologues of therian ε and β, respectively, and arose by independent duplication of an ancestral β-globin gene in the monotreme lineage, with a separate duplication event, just prior to the divergence of therians, producing the progenitors of ε and β of therians. This hypothesis was strongly supported by our BI phylogenetic (Figure ##FIG##2##3##) analyses, but not by MP analyses of coding sequence data, with third codon sites excluded (Figure ##FIG##3##4##), or BI analyses of protein sequence data (not shown). These contradictory analyses highlight the difficulty in resolving deep relationships among globin genes, particularly when the time periods between duplication and speciation events are relatively small, the phylogenetic signal at third codon positions is potentially saturated, and non-synonymous sites may be subjected to purifying or positive selection. Despite a very high posterior probability (100%) for the grouping of platypus ε with monotreme β, this value is a Bayesian probability and depends on the model adequately representing the evolution of the gene. Furthermore, although it was reported [##REF##18216242##62##] that the 5' flanking sequences of platypus ε and β were similar, we found no evidence for similarity of the promoter signals of these two genes (Figure ##FIG##4##5B##).</p> <p>We consider that a more parsimonious explanation is that the platypus ε is orthologous to the marsupial and eutherian embryonic β-like globin lineages (ε and γ), and arose by duplication of an ancestral β-globin gene prior to the mammalian radiation (166 MYA; Figure ##FIG##8##9B##). The sequence of platypus ε may have been homogenised by some gene conversion events, leading it to group with other monotreme adult β-like globin genes. In addition to the MP analyses reported above, this explanation is further supported by the conserved position of ε to the 5' side of the adult β-globin gene in the platypus cluster, which is similar to that found in other therian β-globin gene clusters [##REF##8752008##26##]; see also [##REF##6527390##29##]). Amino acid sequence analyses (BlastP) also provided additional support for the orthology of platypus ε to other mammalian ε-globin genes. Although we were unable to examine the expression of the genes in embryonic tissues, it was found that the expression profile of the platypus ε was similar to the embryonic α-like globins ζ and ζ' of the platypus, but not to the adult β-globin gene, supporting its potential role as an embryonic β-like globin gene.</p> <title>The ω-globin gene and the evolution of the β-globin cluster</title> <p>The discovery of the marsupial ω-globin gene in the α-globin cluster [##REF##11158601##35##,##REF##15461421##36##] was critical in re-interpreting the relationships of the α- and β-globin clusters in amniotes (reptiles, birds and mammals) to favour the hypothesis that these clusters in birds and mammals are paralogous, having diverged independently from different ancestral copies of the vertebrate α-β-globin locus [##REF##11171947##63##].</p> <p>Our observation of an ω-globin gene in the α-globin cluster in the platypus, as well as in the marsupials, confirms that the ancestral mammal α-globin cluster contained a β-like globin gene that was lost in eutherians, as proposed by Wheeler et al. [##REF##11158601##35##,##REF##15461421##36##]. However, the position of monotreme and marsupial ω in the phylogeny (Figure ##FIG##2##3##) is more consistent with the original hypothesis [##REF##7407927##5##] that mammal and bird β-globin are orthologous, having descended from the same β-globin progenitor in an amniote ancestor, and this is strongly supported by flanking sequence data (see below). Our data support the proposition that the ω <italic>-</italic>globin gene represents an ancient β-like globin gene lineage that is ancestral to a group containing both mammalian and bird β-globins with a high posterior probability (99%). This arrangement, however, was not supported by analyses of amino acid sequence data, indicating that there is uncertainty in the phylogenetic position of ω-globin relative to bird β-globins, or that convergent evolution of bird β-globin genes and ω-globin resulted in their similarity at the protein level. To further resolve the key question of whether bird and mammal β-globin gene clusters are orthologous we carried out comparative analyses of flanking loci of the α- and β-globin clusters.</p> <title>Genome context of vertebrate α- and β-globin clusters</title> <p>We found that the platypus α-globin cluster is flanked by <italic>MPG</italic>, <italic>C16orf35</italic>, <italic>GBY </italic>and <italic>LUC7L</italic>, and that the same genes (except <italic>GBY</italic>) flank the α-globin cluster in mammals and birds [##REF##11157800##58##,##REF##15998734##59##]. The same genes flank the α-β cluster of frog, and even zebrafish and the α-cluster of pufferfish [##REF##12517812##8##] (except <italic>GBY </italic>and <italic>LUC7L</italic>), implying that a very ancient region containing these genes (5'-<italic>MPG</italic>-<italic>C16orf35</italic>-α-β-<italic>GBY</italic>-<italic>LUC7L</italic>-3'), or perhaps an even larger region, was present in their common ancestor and has been conserved since the evolution of jawed vertebrates more than 450 MYA.</p> <p>In contrast, the amniote β-globin clusters reside in a very different genome, sharing none of the flanking loci with the mammal and bird α-globin clusters, or the α-β cluster of frogs and fish. In platypus, as well as in therian mammals [##REF##16554811##60##,##REF##10220430##61##], the β-globin clusters are flanked by numerous <italic>ORG </italic>genes on both sides. In birds, also, the β-globin cluster is embedded in <italic>ORG </italic>genes [##REF##16554811##60##]. Even the outside loci <italic>RRM1</italic>, <italic>CCKBR </italic>and <italic>ILK </italic>lie in the same orientation with respect to the bird and mammalian β-globin clusters [##REF##16554811##60##], suggesting that the 5'-<italic>RRM1-ORG-β </italic>(cluster)<italic>-ORG-CCKBR-ILK</italic>-3' arrangement has been conserved since before the divergence of birds and mammals, more than 315 MYA. Therefore, the bird β-globin cluster is orthologous to the β-globin clusters of mammals.</p>
<title>Conclusion</title> <title>New model for the evolution of α- and β-globin clusters in amniotes</title> <p>This analysis of flanking loci, in addition to the phylogenetic analyses reported above, refutes the prevailing hypothesis that mammal and bird α- and β-globin clusters evolved from different (paralogous) copies of an ancestral α-β-globin region containing <italic>MPG-C16orf35-α </italic>(cluster)-β (cluster)<italic>-GBY-LUC7L</italic>. Rather, the context of β-globin clusters within olfactory receptor genes in birds as well as mammals suggests that a copy of a β-globin locus was moved into a region replete with <italic>ORG </italic>genes before the divergence of birds and mammals 315 MYA. The precise mechanism for this translocation is unknown, but is likely to be either by transposition of a tandem duplicate of an ancestral β-globin gene, or retrotransposition of an intron-containing primary transcript. Phylogenetic analyses suggest that this ancestral β-globin gene within the α-globin cluster is represented by the platypus and marsupial ω-globin gene. The transposed β-globin gene then independently duplicated several times within the avian and mammalian lineages to form the different clusters of differentially expressed β-globin genes. Full details of this new model are given in Figure ##FIG##8##9A## and ##FIG##8##9B##.</p> <p>This hypothesis could be further tested by investigating the gene organization of the α- and β-globin clusters in reptiles such as lizards and snakes, which form a sister group to birds. Our hypothesis predicts that reptiles should possess a <italic>MPG-C16orf35</italic>-α (cluster)-β (cluster)<italic>-GBY-LUC7L </italic>cluster, and an unlinked <italic>RRM1-ORG-β </italic>(cluster)<italic>-ORG-CCBKR-ILK </italic>cluster like birds and mammals. The full genome sequence of the first reptilian species,<italic>Anolis carolinensis</italic>, will provide an opportunity to test this hypothesis.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Vertebrate alpha (α)- and beta (β)-globin gene families exemplify the way in which genomes evolve to produce functional complexity. From tandem duplication of a single globin locus, the α- and β-globin clusters expanded, and then were separated onto different chromosomes. The previous finding of a fossil β-globin gene (ω) in the marsupial α-cluster, however, suggested that duplication of the α-β cluster onto two chromosomes, followed by lineage-specific gene loss and duplication, produced paralogous α- and β-globin clusters in birds and mammals. Here we analyse genomic data from an egg-laying monotreme mammal, the platypus (<italic>Ornithorhynchus anatinus</italic>), to explore haemoglobin evolution at the stem of the mammalian radiation.</p> <title>Results</title> <p>The platypus α-globin cluster (chromosome 21) contains embryonic and adult α- globin genes, a β-like ω-globin gene, and the <italic>GBY </italic>globin gene with homology to cytoglobin, arranged as 5'-ζ-ζ'-α<sup>D</sup>-α<sup>3</sup>-α<sup>2</sup>-α<sup>1</sup>-ω-<italic>GBY</italic>-3'. The platypus β-globin cluster (chromosome 2) contains single embryonic and adult globin genes arranged as 5'-ε-β-3'. Surprisingly, all of these globin genes were expressed in some adult tissues. Comparison of flanking sequences revealed that all jawed vertebrate α-globin clusters are flanked by <italic>MPG</italic>-<italic>C16orf35 </italic>and <italic>LUC7L</italic>, whereas all bird and mammal β-globin clusters are embedded in olfactory genes. Thus, the mammalian α- and β-globin clusters are orthologous to the bird α- and β-globin clusters respectively.</p> <title>Conclusion</title> <p>We propose that α- and β-globin clusters evolved from an ancient <italic>MPG-C16orf35-α</italic>-β-<italic>GBY-LUC7L </italic>arrangement 410 million years ago. A copy of the original β (represented by ω in marsupials and monotremes) was inserted into an array of olfactory genes before the amniote radiation (&gt;315 million years ago), then duplicated and diverged to form orthologous clusters of β-globin genes with different expression profiles in different lineages.</p>
<title>Authors' contributions</title> <p>VSP designed and performed most of the experiments and analysed the data. VSP also drafted the main manuscript. SJBC conducted phylogenetic analyses and contributed to the writing of the manuscript. JED helped in designing the experiments and trouble-shooting experiments. BF, TG, WCW and RKW were involved in sequencing the platypus BAC clone (Oa_Bb-484F22). JAMG conceived and supervised the research and contributed to the writing of the manuscript. All authors read and approved the final manuscript.</p> <title>Competing Interests</title> <p>The authors declare that they have no competing interests.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We thank Dr Frank Grützner (The University of Adelaide) and Tim Hore (The Australian National University) for providing platypus tissues and RNA preparations. We also thank Dr Paul Waters for providing computational help with the identification and confirmation of two identical but different genes α<sup>1 </sup>and α<sup>3</sup>.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Current proposed model for the evolution of α- and β-globin clusters from paralogous clusters in different lineages</bold>. The unlinked α- and β-globin clusters in birds and mammals evolved from an ancient <italic>in trans </italic>duplication of the ancestral linked α-β cluster, followed by differential gene silencing (marked with X). This resulted in bird β-like globin genes (β2) orthologous to the marsupial ω-globin gene (β2 beside the α-globin cluster) but paralogous to mammalian β-like globin genes (β1). Adapted from Wheeler et al. [##REF##15461421##36##].</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Evolutionary relationships among vertebrate α-like globin genes using a 50% majority rule consensus phylogram from an analysis using Bayesian Inference</bold>. The tree was constructed using mixed models of evolution for each codon position (see methods) and estimated base frequencies in an unlinked analysis using MrBayes (v. 3.1.2). Numbers adjacent to branches refer to % posterior probabilities. GenBank accession numbers for sequences are: Virginian Opossum (<italic>Didelphis virginiana</italic>) ζ<sup>1</sup>, ζ<sup>2</sup>, α<sup>1</sup>, α<sup>2</sup>, θ [<ext-link ext-link-type="gen" xlink:href="AC139599.2">AC139599.2</ext-link>, <ext-link ext-link-type="gen" xlink:href="AC148752.1">AC148752.1</ext-link>]; Stripe-faced Dunnart (<italic>Sminthopsis macroura</italic>) α<sup>D</sup>, α<sup>2</sup>, θ [<ext-link ext-link-type="gen" xlink:href="AC146781">AC146781</ext-link>]; Brazilian Opossum (<italic>Monodelphis domestica</italic>) α [TI# 453585430]; Tammar wallaby (<italic>Macropus eugenii</italic>) θ [<ext-link ext-link-type="gen" xlink:href="AY459590">AY459590</ext-link>], α [<ext-link ext-link-type="gen" xlink:href="AY459589">AY459589</ext-link>]; ζ [<ext-link ext-link-type="gen" xlink:href="AY789121">AY789121</ext-link>], ζ' [<ext-link ext-link-type="gen" xlink:href="AY789122">AY789122</ext-link>]; Horse (<italic>Equus caballus</italic>) θ (ψ α) [<ext-link ext-link-type="gen" xlink:href="Y00284">Y00284</ext-link>], α<sup>1 </sup>[<ext-link ext-link-type="gen" xlink:href="M17902">M17902</ext-link>], ζ [<ext-link ext-link-type="gen" xlink:href="X07051">X07051</ext-link>]; pig (<italic>Sus scrofa</italic>) α<sup>D </sup>[<ext-link ext-link-type="gen" xlink:href="AC145444">AC145444</ext-link>]; cat (<italic>Felis catus</italic>) α<sup>D </sup>[<ext-link ext-link-type="gen" xlink:href="AC130194">AC130194</ext-link>]; cow (<italic>Bos taurus</italic>) α<sup>D </sup>[<ext-link ext-link-type="gen" xlink:href="AC150547">AC150547</ext-link>]; Goat (<italic>Capra hircus</italic>) α [<ext-link ext-link-type="gen" xlink:href="J00043">J00043</ext-link>]; Human (<italic>Homo sapiens</italic>) α1 [<ext-link ext-link-type="gen" xlink:href="V00491">V00491</ext-link>], θ [<ext-link ext-link-type="gen" xlink:href="X06482">X06482</ext-link>], ζ [<ext-link ext-link-type="gen" xlink:href="NM_005332">NM_005332</ext-link>]; <italic>mu</italic>/α<sup>D </sup>chain [<ext-link ext-link-type="gen" xlink:href="AY698022">AY698022</ext-link>]; Mouse (<italic>Mus muscularis</italic>) α<sup>1 </sup>[<ext-link ext-link-type="gen" xlink:href="NM_008218">NM_008218</ext-link>], ζ [<ext-link ext-link-type="gen" xlink:href="X62302">X62302</ext-link>]; Rabbit (<italic>Oryctolagus cuniculus</italic>)α [<ext-link ext-link-type="gen" xlink:href="X04751">X04751</ext-link>]; Eastern Quoll (<italic>Dasyurus viverrinus</italic>) α [<ext-link ext-link-type="gen" xlink:href="M14567">M14567</ext-link>]; Chicken (<italic>Gallus gallus) </italic>α<sup>A</sup>, π, α<sup>D </sup>[<ext-link ext-link-type="gen" xlink:href="AF098919">AF098919</ext-link>]; Duck (<italic>Cairina moschata</italic>) α<sup>D </sup>[<ext-link ext-link-type="gen" xlink:href="X01831">X01831</ext-link>]; Pigeon (<italic>Columba livia</italic>) α<sup>D </sup>[<ext-link ext-link-type="gen" xlink:href="AB001981">AB001981</ext-link>]; Turtle (<italic>Geochelone nigra</italic>) α<sup>D </sup>[SEG# <ext-link ext-link-type="gen" xlink:href="AB1165195">AB1165195</ext-link>]; Zebrafish (<italic>Danio rerio</italic>) α<sup>1 </sup>[<ext-link ext-link-type="gen" xlink:href="NM_131257">NM_131257</ext-link>]; Salamander (<italic>Hynobius retardatus</italic>) larval α [<ext-link ext-link-type="gen" xlink:href="AB034756">AB034756</ext-link>]; Salamander (<italic>Pleurodeles waltlii</italic>) α [<italic>M</italic>13365]; Frog (<italic>Xenopus laevis</italic>) α <italic>I </italic>[<ext-link ext-link-type="gen" xlink:href="X14259">X14259</ext-link>], larval (tadpole) α T5 [<ext-link ext-link-type="gen" xlink:href="X02798">X02798</ext-link>]; Yellowtail (<italic>Seriola quinqueradiata</italic>) α<sup>A </sup>[<ext-link ext-link-type="gen" xlink:href="AB034639">AB034639</ext-link>]; Salmon (<italic>Salmo salar</italic>) α [<ext-link ext-link-type="gen" xlink:href="X97289">X97289</ext-link>]; Southern Puffer (<italic>Sphoeroides nephelus</italic>) α<sup>2 </sup>[<ext-link ext-link-type="gen" xlink:href="AY016023">AY016023</ext-link>]; Platypus (<italic>Ornithorhynchus anatinus</italic>) ζ, ζ', α<sup>D</sup>, α<sup>3</sup>, α<sup>2</sup>, α<sup>1 </sup>[<ext-link ext-link-type="gen" xlink:href="AC203513">AC203513</ext-link>].</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Evolutionary relationships among vertebrate β-like globin genes using a 50% majority rule consensus phylogram from an analysis using Bayesian inference</bold>. The tree was constructed using mixed models of evolution for each codon position (see methods) and estimated base frequencies in an unlinked analysis using MrBayes (v. 3.1.2). Numbers adjacent to branches refer to % posterior probabilities. GenBank accession numbers for sequences are: Fat-tailed Dunnart (<italic>Sminthopsis crassicaudata</italic>) β [<ext-link ext-link-type="gen" xlink:href="Z69592">Z69592</ext-link>], ε [<ext-link ext-link-type="gen" xlink:href="Z48632">Z48632</ext-link>], ω [<ext-link ext-link-type="gen" xlink:href="AY014770">AY014770</ext-link>]; Stripe-faced Dunnart (<italic>S. macroura</italic>) β, ε [<ext-link ext-link-type="gen" xlink:href="AC148754">AC148754</ext-link>]; Virginian Opossum (<italic>Didelphis virginiana</italic>) β [<ext-link ext-link-type="gen" xlink:href="J03643">J03643</ext-link>], ε [<ext-link ext-link-type="gen" xlink:href="J03642">J03642</ext-link>]; Brazilian Opossum (<italic>Monodelphis domestica</italic>) β [<ext-link ext-link-type="gen" xlink:href="XM_001365299">XM_001365299</ext-link>], ε [<ext-link ext-link-type="gen" xlink:href="XM_001364448">XM_001364448</ext-link>], ω [<ext-link ext-link-type="gen" xlink:href="XM_001364828">XM_001364828</ext-link>]; Tammar Wallaby (<italic>Macropus eugenii</italic>) β [<ext-link ext-link-type="gen" xlink:href="AY450928">AY450928</ext-link>], ε [<ext-link ext-link-type="gen" xlink:href="AY450927">AY450927</ext-link>], ω [<ext-link ext-link-type="gen" xlink:href="AY014769">AY014769</ext-link>]; African clawed frog (<italic>Xenopus laevis</italic>) larval β <italic>I </italic>[<ext-link ext-link-type="gen" xlink:href="NM_001086273">NM_001086273</ext-link>], larval βII [<ext-link ext-link-type="gen" xlink:href="NM_001088028">NM_001088028</ext-link>]; Western clawed frog (<italic>X. tropicalis</italic>) β [<ext-link ext-link-type="gen" xlink:href="NM_203528">NM_203528</ext-link>], larval ε1 [<ext-link ext-link-type="gen" xlink:href="NM_001016495">NM_001016495</ext-link>]; Chicken (<italic>Gallus gallus</italic>) β (β<sup>A</sup>) [<ext-link ext-link-type="gen" xlink:href="NM_205489">NM_205489</ext-link>], ε [<ext-link ext-link-type="gen" xlink:href="NM_001004390">NM_001004390</ext-link>], γ (β<sup>A</sup>) [<ext-link ext-link-type="gen" xlink:href="NM_001031489">NM_001031489</ext-link>]; Duck (<italic>Cairina moschata</italic>) β [<ext-link ext-link-type="gen" xlink:href="J00926">J00926</ext-link>], ε [<ext-link ext-link-type="gen" xlink:href="X15740">X15740</ext-link>]; Human (<italic>Homo sapiens</italic>) β [<ext-link ext-link-type="gen" xlink:href="NM_000518">NM_000518</ext-link>], γ [<ext-link ext-link-type="gen" xlink:href="BC130459">BC130459</ext-link>], ε [<ext-link ext-link-type="gen" xlink:href="NM_005330">NM_005330</ext-link>]; Mouse (<italic>Mus musculus</italic>) β (β1) [<ext-link ext-link-type="gen" xlink:href="NM_008220">NM_008220</ext-link>], γ (β h0) [<ext-link ext-link-type="gen" xlink:href="NW_001030869">NW_001030869</ext-link>], ε (ε<sup>y</sup>) [<ext-link ext-link-type="gen" xlink:href="M26897">M26897</ext-link>]; Goat (<italic>Capra hirus</italic>)<italic>β </italic>(β<sup>A</sup>) [<ext-link ext-link-type="gen" xlink:href="DQ350619">DQ350619</ext-link>], ε (ε<sup>I</sup>) [<ext-link ext-link-type="gen" xlink:href="X01912">X01912</ext-link>], γ [<ext-link ext-link-type="gen" xlink:href="M15388">M15388</ext-link>]; Rabbit (<italic>Oryctolagus cuniculus</italic>) β, γ, ε [<ext-link ext-link-type="gen" xlink:href="M18818">M18818</ext-link>]; Echidna (<italic>Tachyglossus aculeatus) </italic>β [<ext-link ext-link-type="gen" xlink:href="L23800">L23800</ext-link>]; Pufferfish (<italic>Fugu rubripes</italic>) β [<ext-link ext-link-type="gen" xlink:href="AY170464">AY170464</ext-link>]; Zebrafish (<italic>Danio rerio</italic>) ε1 [<ext-link ext-link-type="gen" xlink:href="NM_001103130">NM_001103130</ext-link>]; Platypus β, ε [<ext-link ext-link-type="gen" xlink:href="AC192436">AC192436</ext-link>], ω [<ext-link ext-link-type="gen" xlink:href="AC203513">AC203513</ext-link>].</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Evolutionary relationships among vertebrate β-like globin genes analysed by maximum parsimony (MP) trees of length 926 (one of eight trees)</bold>. Third position in codons were excluded in the MP analyses, which were conducted using a heuristic search in PAUP* v.4.0b10 [##REF##15499007##65##]. The tree is rooted using pufferfish β-globin. Numbers adjacent to branches represent % bootstrap values (&gt;50%) from MP heuristic analyses of 1000 pseudoreplicates. Accession numbers for sequences are given in the caption of Figure 3.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>Gene structure of the platypus α- and β-globin clusters and flanking loci, and comparisons of their promoter regions with other mammals</bold>. (A) The platypus α-globin cluster contains six α-like globin genes (red), a β-like (ω) globin gene (blue) and a distantly related globin gene, <italic>GBY </italic>(green), which are flanked by <italic>IL9RP3-POLR3K</italic>-<italic>C16orf33</italic>-<italic>C16orf8-MPG-C16orf35 </italic>on the 5' end and <italic>LUC7L-ITFG3-RGS11</italic>-<italic>ARHGDIG</italic>-<italic>PDIA2</italic>-<italic>AXIN1 </italic>on the 3' end (black). The platypus β-globin cluster contains only two genes, ε and β (blue), which are flanked on both sides by <italic>ORG </italic>genes (black). (B) Relative positions of the putative transcription factor binding sites in the 200 bp promoter region located upstream of the predicted platypus, marsupial (<italic>Didelphis virginiana </italic>ζ and ψζ', and <italic>Sminthopsis macroura </italic>α<sup>D</sup>, ψα<sup>3</sup>, α<sup>2</sup>, α<sup>1</sup>, ω, ε and β) and human α- and β-like globin genes. For the platypus <italic>GBY </italic>no data was available from other species, including <italic>Xenopus tropicalis</italic>, for comparisons.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p><bold>Expression of all predicted α- and β-like globin genes including <italic>GBY </italic>in an adult platypus</bold>. For each of the platypus predicted genes, expression was investigated by reverse transcriptase polymerase chain reaction in adult liver, kidney, spleen, testis, brain and lung. Primers for each gene were designed between two exons so that it would result in a product distinguishable from genomic contamination of cDNA. The negative control (last lane) contained no cDNA. All genes were expressed in one or more tissues, indicating that they are transcriptionally active and might be functional.</p></caption></fig> <fig position="float" id="F7"><label>Figure 7</label><caption><p><bold>Chromosomal location of the platypus α- and β-globin clusters</bold>. Two-colour fluorescence in situ hybridisation showing the location of the α-globin cluster on chromosome 21 (green) and the β-globin cluster on chromosome 2q5.1 (red). The chromosomes are counterstained with DAPI (blue).</p></caption></fig> <fig position="float" id="F8"><label>Figure 8</label><caption><p><bold>Loci flanking vertebrate α-globin (A) and β-globin (B) clusters</bold>. The relative locations of flanking loci (A) <italic>MPG</italic>, <italic>C16orf35</italic>, <italic>LUC7L </italic>and <italic>GBY </italic>and (B) <italic>RRM1</italic>, <italic>CCKBR</italic>, <italic>ILK </italic>and <italic>ORG </italic>genes were searched for beside the α-β globin cluster in zebrafish (<italic>Danio rerio</italic>) and frog (<italic>Xenopus tropicalis</italic>), and beside the separate α-globin and β-globin clusters in chicken (<italic>Gallus gallus</italic>), opossum (<italic>Monodelphis domestica</italic>) and human (<italic>Homo sapiens) </italic>from Ensembl [##UREF##4##57##]. The pufferfish (<italic>Fugu rubripes</italic>) flanking loci shown here were adapted from Gillemans et al. [##REF##12517812##8##]. For the platypus, the α-globin flanking loci were characterised in this study, and <italic>ORG </italic>genes surrounding the platypus β-globin cluster were discovered: however, the BAC clone (484F22) was too small to cover the region containing the loci <italic>RRM1</italic>, <italic>CCKBR </italic>and <italic>ILK</italic>. In <italic>X. tropicalis LUC7L </italic>was found on another scaffold (466 from Ensembl) but sequence analyses by Fuchs et al. [##REF##16484786##2##] suggested that <italic>LUC7L </italic>resides 3' to the frog α-β-<italic>GBY </italic>cluster. The flanking loci as well as the α- and β-globin clusters are differentiated by colour.</p></caption></fig> <fig position="float" id="F9"><label>Figure 9</label><caption><p><bold>Proposed model for the evolution of the α- and β-globin clusters in vertebrate lineages</bold>. (A) A region containing <italic>MPG-C16orf35</italic>-α-β-<italic>GBY</italic>-<italic>LUC7L </italic>represented the ancient α-β globin cluster of jawed vertebrates (&gt;450 MYA), which is seen in the amphibian lineage. This region further duplicated and underwent some gene silencing in teleost fish. In an amniote ancestor of reptiles, birds and mammals (&gt;315 MYA), a copy of an ancestral β-globin gene from this region was inserted into a different chromosome within a region replete with multiple copies of <italic>ORG </italic>genes. The original amniote β-globin gene survives as the ω-globin gene (β1) in the α-globin cluster of marsupials and monotremes, whereas the transposed β-globin gene (β2) duplicated several times to form different clusters in the different lineages. (B) Tandem duplications of the ancestral amniote α-globin gene produced a three-gene (π-α<sup>D</sup>-α<sup>A</sup>) cluster in the avian lineage. In the mammalian lineage, further duplications gave rise to a six-gene (ζ-ζ'-α<sup>D</sup>-α<sup>3</sup>-α<sup>2</sup>-α<sup>1</sup>) cluster with ongoing gene conversion events homogenising the embryonic and adult genes. In monotremes, the ancestral ω (β1) and <italic>GBY </italic>are retained. After the divergence of monotreme and therian mammals, there was an additional duplication of α<sup>2 </sup>to form θ, giving rise to the seven-gene cluster (ζ-ζ'-α<sup>D</sup>-α<sup>3</sup>-α<sup>2</sup>-α<sup>1</sup>-θ) in marsupials and eutherians. Marsupials also retain the ancestral ω but may have lost <italic>GBY </italic>gene; eutherians retain no identifiable remnant of either gene. Furthermore, the ancestral transposed β2-globin gene duplicated independently in birds and mammals. Before the mammalian radiation, we propose that the ancestral β2 gene duplicated to form a two-gene β-globin cluster (ε-β) as seen in monotremes and marsupials, except that ongoing gene conversion events homogenised platypus ε to group with monotreme β genes. After the divergence of marsupial and eutherian mammals, there were further tandem duplications of these two genes to produce complex β-globin cluster (ε-γ-η-δ-β) in eutherians.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Gene-structure of the predicted platypus α- and β-like globin genes and <italic>GBY</italic></p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Length/</bold><break/><bold>Genes</bold><break/></td><td align="left"><bold>Exon 1 </bold><break/><bold>(bp)</bold><break/></td><td align="left"><bold>Intron 1 </bold><break/><bold>(bp)</bold><break/></td><td align="left"><bold>Exon 2 </bold><break/><bold>(bp)</bold><break/></td><td align="left"><bold>Intron 2 </bold><break/><bold>(bp)</bold><break/></td><td align="left"><bold>Exon 3 </bold><break/><bold>(bp)</bold><break/></td><td align="left"><bold>Position of </bold><break/><bold>Poly-A</bold><break/></td><td align="left"><bold>CDS </bold><break/><bold>(bp)</bold><break/></td><td align="left"><bold>Poly- </bold><break/><bold>peptide </bold><break/><bold>(aa)</bold></td></tr></thead><tbody><tr><td align="left"><bold>ζ</bold></td><td align="left">95</td><td align="left">337</td><td align="left">205</td><td align="left">114</td><td align="left">129</td><td align="left">+119</td><td align="left">429</td><td align="left">142</td></tr><tr><td align="left"><bold>ζ<sup>1</sup></bold></td><td align="left">95</td><td align="left">336</td><td align="left">205</td><td align="left">102</td><td align="left">129</td><td align="left">+133</td><td align="left">429</td><td align="left">142</td></tr><tr><td align="left"><bold>α<sup>D</sup></bold></td><td align="left">92</td><td align="left">1450</td><td align="left">205</td><td align="left">1610</td><td align="left">129</td><td align="left">+77</td><td align="left">426</td><td align="left">141</td></tr><tr><td align="left"><bold>α<sup>1</sup></bold>/<bold>α</bold><sup><bold>3</bold></sup></td><td align="left">92</td><td align="left">405</td><td align="left">205</td><td align="left">151</td><td align="left">129</td><td align="left">+94</td><td align="left">426</td><td align="left">141</td></tr><tr><td align="left"><bold>α<sup>2</sup></bold></td><td align="left">95</td><td align="left">720</td><td align="left">205</td><td align="left">155</td><td align="left">129</td><td align="left">+115</td><td align="left">429</td><td align="left">142</td></tr><tr><td align="left"><bold>ω</bold></td><td align="left">92</td><td align="left">256</td><td align="left">223</td><td align="left">111</td><td align="left">129</td><td align="left">+69</td><td align="left">444</td><td align="left">147</td></tr><tr><td align="left"><bold><italic>GBY</italic></bold></td><td align="left">98</td><td align="left">3364</td><td align="left">223</td><td align="left">3053</td><td align="left">144</td><td align="left">+141</td><td align="left">465</td><td align="left">154</td></tr><tr><td align="left"><bold>ε</bold></td><td align="left">92</td><td align="left">143</td><td align="left">223</td><td align="left">474</td><td align="left">129</td><td align="left">+96</td><td align="left">444</td><td align="left">147</td></tr><tr><td align="left"><bold>β</bold></td><td align="left">92</td><td align="left">153</td><td align="left">223</td><td align="left">438</td><td align="left">129</td><td align="left">+71</td><td align="left">444</td><td align="left">147</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>PCR primers used for amplification of the α- and β-like globin genes including <italic>GBY </italic>from the platypus gDNA and cDNA</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Gene</bold><break/></td><td align="left"><bold>Forward Primer</bold><break/></td><td align="left"><bold>Reverse Primer</bold><break/></td><td align="left"><bold>gDNA</bold><break/><bold>(bp)</bold></td><td align="left"><bold>cDNA </bold><break/><bold>(bp)</bold></td></tr></thead><tbody><tr><td align="left"><bold>ζ</bold></td><td align="left">GGCCGACAAGACCGCAGTCATCTCCC</td><td align="left">CCCGATGGCGCTGATGACT</td><td align="left">527</td><td align="left">190</td></tr><tr><td align="left"><bold>ζ<sup>1</sup></bold></td><td align="left">TGACCAAAGGCGACAAGACCT</td><td align="left">CCCCGATGGCACCGATGACC</td><td align="left">534</td><td align="left">198</td></tr><tr><td align="left"><bold>α<sup>D</sup></bold></td><td align="left">GAGGCTGTGAAGAACCTGGA</td><td align="left">GGTGTACTCCCCTTGCAGAT</td><td align="left">1793</td><td align="left">153</td></tr><tr><td align="left"><bold>α<sup>2</sup></bold></td><td align="left">TGGCCCACCTCGATGACCTGG</td><td align="left">GGGAAGGTGTCTGGCCACC</td><td align="left">289</td><td align="left">134</td></tr><tr><td align="left"><bold>α<sup>1</sup></bold>/<bold>α</bold><sup><bold>3</bold></sup></td><td align="left">GCAAGGCCGCCGGTCACGGC</td><td align="left">CGCTGTCCATGTCATCGAAGTGCC</td><td align="left">597</td><td align="left">192</td></tr><tr><td align="left"><bold>ω</bold></td><td align="left">ATTGTGTCCATCTGGGGAAA</td><td align="left">GCTTGGCAAAGTTGCTCTTC</td><td align="left">488</td><td align="left">232</td></tr><tr><td align="left"><bold><italic>GBY</italic></bold></td><td align="left">CTGGAAACAGGTGTGCAAGA</td><td align="left">CTATCTCCGGGGTGTAGCAG</td><td align="left">3202</td><td align="left">149</td></tr><tr><td align="left"><bold>ε</bold></td><td align="left">ATCTGAGCGCTGAGGAGAAG</td><td align="left">GACAGGTTGCCGAAGGAGTCA</td><td align="left">285</td><td align="left">142</td></tr><tr><td align="left"><bold>β</bold></td><td align="left">CTGTGGGGGAAAGTGAACAT</td><td align="left">GGTCAGCACCTTAGCACCAT</td><td align="left">321</td><td align="left">168</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p><bold>Annotation of the platypus α- and β-like and <italic>GBY </italic>globin genes</bold>. This table shows the predicted positions of six α-like, ω and <italic>GBY </italic>globin genes in the platypus BAC clone Oa_Bb-131M24 [GenBank: <ext-link ext-link-type="gen" xlink:href="AC203513">AC203513</ext-link>], two α-like globin genes in BAC clone Oa_Bb-2L7 [GenBank: <ext-link ext-link-type="gen" xlink:href="AC195438">AC195438</ext-link>], and two β-like globin genes in BAC clone Oa_Bb-484F22 [GenBank: <ext-link ext-link-type="gen" xlink:href="AC192436">AC192436</ext-link> reverse direction].</p></caption></supplementary-material>
<table-wrap-foot><p>For each predicted gene, the length of the exons and introns, position of poly-A addition site (AATAAA) from the stop codon, and the length of their putative coding domain (CDS) and encoded polypeptide are shown. All genes contained consensus splice sites (GT/AG) in both introns.</p></table-wrap-foot>
<graphic xlink:href="1741-7007-6-34-1"/> <graphic xlink:href="1741-7007-6-34-2"/> <graphic xlink:href="1741-7007-6-34-3"/> <graphic xlink:href="1741-7007-6-34-4"/> <graphic xlink:href="1741-7007-6-34-5"/> <graphic xlink:href="1741-7007-6-34-6"/> <graphic xlink:href="1741-7007-6-34-7"/> <graphic xlink:href="1741-7007-6-34-8"/> <graphic xlink:href="1741-7007-6-34-9"/>
<media xlink:href="1741-7007-6-34-S1.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Dickerson", "Geis"], "given-names": ["RE", "I"], "source": ["Hemoglobin: Structure, Function, Evolution and Pathology"], "year": ["1983"], "publisher-name": ["Menlo Park, CA: Benjamin Cummings"], "fpage": ["65"], "lpage": ["116"]}, {"surname": ["Weatherall"], "given-names": ["DJ"], "source": ["The Structure, Organization, and Regulation of Human Genes"], "year": ["1991"], "edition": ["3"], "publisher-name": ["Oxford: Oxford University Press"]}, {"surname": ["Andersen", "Mesch", "Lovell", "Nicol"], "given-names": ["NA", "U", "DJ", "SC"], "article-title": ["The effects of sex, season and hybernation on haematology and blood viscosity of free-ranging echidnas ("], "italic": ["Tachyglossus aculeatus"], "source": ["Can J Zool"], "year": ["2000"], "volume": ["78"], "fpage": ["174"], "lpage": ["181"], "pub-id": ["10.1139/cjz-78-2-174"]}, {"article-title": ["University of California Santa Cruz (UCSC) Genome Browser"]}, {"article-title": ["Ensembl Genome Browser"]}, {"article-title": ["NIH Intramural Sequencing Center (NISC)"]}, {"surname": ["Felsenstein"], "given-names": ["J"], "article-title": ["Confidence limits on phylogenies: an approach using the bootstrap"], "source": ["Evolution"], "year": ["1985"], "volume": ["39"], "fpage": ["783"], "lpage": ["791"], "pub-id": ["10.2307/2408678"]}, {"article-title": ["MrMTgui v1.0"]}, {"surname": ["Yang"], "given-names": ["Z"], "article-title": ["Among-site rate variation and its impact on phylogenetic analyses"], "source": ["Trends Ecol Evol"], "year": ["1996"], "volume": ["11"], "fpage": ["367"], "lpage": ["372"], "pub-id": ["10.1016/0169-5347(96)10041-0"]}, {"surname": ["Rambaut", "Drummond"], "given-names": ["A", "AJ"], "article-title": ["Traver v1.3"], "year": ["2004"]}]
BMC Biol. 2008 Jul 25; 6:34