<title>Introduction</title> <p>CCCTC-binding factor (CTCF) is a ubiquitously expressed protein that binds to more than 20,000 sites within the human genome ##REF##2284094##[1]##–##REF##17512414##[3]##. The distribution of these binding sites, along with experimental data from several well-characterised loci (reviewed ##REF##16953792##[4]##) indicates that CTCF acts as an insulator protein genome-wide, defining boundaries for gene clusters or segregating alternative promoters. This can affect gene expression, for instance at the well-studied chicken <italic>ß–globin</italic> locus, where CTCF binding to the FII insulator leads to transcriptional silencing by blocking the effects of a nearby enhancer ##REF##10458613##[5]##.</p> <p>CTCF is also required for inter-chromosomal interactions such as pairing of the X chromosomes during initiation of X chromosome inactivation ##REF##17952071##[6]## and even co-localisation of non-homologous chromosomes ##REF##16614224##[7]##. It is now considered that CTCF contributes more broadly to the establishment of nuclear compartments where transcription is enhanced or repressed ##REF##17913488##[8]##,##REF##17950379##[9]##, rather than functioning only to insulate neighbouring regions of the genome from each other. Given these diverse and significant roles, it is not surprising that CTCF is essential for life (reviewed ##REF##17950379##[9]##). Furthermore, point mutation and loss of heterozygosity of <italic>CTCF</italic> is associated with human cancer, identifying <italic>CTCF</italic> as an important candidate tumour-suppressor gene ##REF##11782357##[10]##.</p> <p>The CTCF protein, and the nucleotide sequence that encodes it, can conceptually be divided into three separate domains (##FIG##0##Figure 1A##). The central (ZF) domain contains ten Cys<sub>2</sub>His<sub>2</sub> zinc-fingers (ZFs), and one Cys<sub>2</sub>HisCys ZF, combinations of which are used to bind various DNA sequences ##REF##8649389##[11]##. Flanking the ZF domain are the N- and C-terminal domains, which interact with other DNA-binding proteins, histones and histone modifying proteins, and the large subunit of polymerase II (reviewed ##REF##17913488##[8]##). In all three of its domains, CTCF shows extraordinary conservation throughout vertebrates ##REF##8649389##[11]##–##REF##16647825##[14]##, and even non-vertebrates ##REF##15678159##[15]##, reflecting the considerable functional constraint CTCF must face due to its multiple essential roles and many interacting partners.</p> <p>In humans and mice, a paralogue of <italic>CTCF</italic> has been identified known as CTCF-like (<italic>CTCFL</italic>), or as it was originally named (and how we will refer to it hereafter), Brother Of Regulator of Imprinted Sites (<italic>BORIS</italic>) ##REF##12011441##[16]##. Human and mouse BORIS posses a suite of ZFs with binding capability, sequence and underlying gene structure that is extremely similar to CTCF (##FIG##0##Figure 1A##). However, the N- and C-terminal domains of human and mouse BORIS show almost no similarity to CTCF, implying that although they can bind the same DNA, they are likely to act differently at these sites.</p> <p>One example of how CTCF and BORIS may function differently comes from their effects on the regulation of genomic imprinting, which is responsible for parent-of-origin specific, mono-allelic gene expression in about 100 mammalian genes ##REF##15990197##[17]##. The most extensively studied imprinted gene, insulin-like growth factor 2 (<italic>IGF2</italic>), is expressed exclusively from the paternally-derived chromosome in eutherian (‘placental’) mammals ##REF##1997210##[18]##–##REF##15044262##[21]## and marsupial mammals ##REF##10603082##[22]##,##REF##15652708##[23]##. Located downstream of <italic>IGF2</italic> is the untranslated RNA <italic>H19</italic>, which is expressed solely from the maternally derived chromosome ##REF##1709450##[24]##,##REF##1380925##[25]##. Biallelic expression of <italic>IGF2</italic> was discovered in the egg-laying monotreme mammals ##REF##11479919##[26]##, birds ##REF##10603082##[22]##,##REF##11455432##[27]## and fish ##REF##15654625##[28]##, implying that imprinting of this region evolved at the same time as viviparity, 180-210MYA ##REF##12878472##[29]##.</p> <p>In mice, imprinted expression of <italic>Igf2/H19</italic> depends on the imprint control region (ICR), an insulator element located between these two genes. The ICR is methylated during spermatogenesis, specifically marking the paternally-derived chromosome ##REF##9234689##[30]##. CTCF binds to the ICR, but only on the unmethylated, maternally-derived chromosome. When bound to the maternally-derived ICR, CTCF performs many functions including protecting the ICR from methylation ##REF##12461525##[31]##–##REF##16928784##[33]##, blocking <italic>Igf2</italic> access to a downstream enhancer (resulting in <italic>Igf2</italic> silencing <italic>in cis</italic>
##REF##10839547##[34]##,##REF##10839546##[35]##), and simultaneously activating <italic>H19</italic> expression ##REF##16928784##[33]##. CTCF is thought to orchestrate these events through the formation of maternal-specific chromosomal loops ##REF##16815976##[36]## and the establishment of local chromatin modifications ##REF##18039862##[37]##. Thus, CTCF acts somatically to ‘<italic>interpret</italic>’ the differential methylation mark of the ICR acquired during gametogenesis, resulting in imprinted expression of <italic>Igf2/H19</italic>.</p> <p>In contrast, BORIS appears to be essential for the <italic>establishment</italic> of differential methylation at the <italic>IGF2/H19</italic> ICR ##REF##17048991##[38]##. In mouse testes, BORIS is bound to the <italic>Igf2/H19</italic> ICR during the time when the ICR becomes methylated. Methylation is accomplished by members of the <italic>de novo</italic> methyltransferase 3 family, of which DNMT3L is essential to this process ##REF##15215868##[39]##,##REF##15318244##[40]## and DNMT3A/3B are partially redundant ##REF##17616512##[41]##. Transgenes containing the mouse ICR were methylated in <italic>Xenopus</italic> oocytes only when co-injected with BORIS, DNMT3L, one of DNMT3A/3B and a histone modifier called protein arginine methyltransferase 7 (PRMT7) ##REF##17048991##[38]##. Thus, BORIS and CTCF both bind to the ICR through their common ZF domain, yet appear to act differently at this site. BORIS <italic>establishes</italic> differential methylation of the ICR and later CTCF <italic>interprets</italic> this mark, resulting in imprinted expression of <italic>Igf2/H19</italic>.</p> <p>Significantly, in humans and mice <italic>CTCF</italic> and <italic>BORIS</italic> show mutually exclusive expression; <italic>BORIS</italic> is transcribed only in certain parts of the developing and adult testes, whereas <italic>CTCF</italic> is expressed in all other regions tested ##REF##12011441##[16]##,##REF##17048991##[38]##. The only reported instances of <italic>BORIS</italic> expression outside of the testes is in various types of cancers ##REF##16140943##[42]##–##REF##17957795##[47]##. This mutually exclusive expression pattern could be explained in part by the recent discovery that CTCF actually binds to the promoter of <italic>BORIS</italic> and negatively regulates its expression ##REF##17962299##[44]##.</p> <p>
<italic>BORIS</italic> is associated with a large group of potentially oncogenic “cancer-testis” (CT) genes, which also show testis-specific, or gonad-specific, expression in healthy individuals, but are highly expressed in cancers ##REF##16034368##[48]##. CTCF binds to the promoter of many CT-genes in healthy somatic tissue where these genes are silenced ##REF##16140943##[42]##,##REF##16140944##[43]##,##REF##17260018##[49]##,##REF##17406683##[50]##. However, this repression is disrupted by conditional expression of BORIS, which replaces CTCF binding at the promoter and subsequently causes local demethylation and gene activation ##REF##16140943##[42]##,##REF##16140944##[43]##,##REF##17260018##[49]##. Similarly, CTCF-binding has a repressive effect on the promoter of the anti-apoptotic gene <italic>BAG1</italic>, whereas BORIS performs oppositely, altering histone methylation and upregulating <italic>BAG1</italic> expression ##REF##18413740##[51]##. The discovery that <italic>CTCF</italic> and <italic>BORIS</italic> have opposite effects on transcription of <italic>BAG1</italic>, some CT-genes, and on the epigenetic status of the <italic>IGF2/H19</italic> ICR, has lead to the (albeit controversial ##REF##17957795##[47]##) hypothesis that CTCF and BORIS are antagonistic regulators of the common loci to which they bind, and that inappropriate interactions between them is cancer promoting ##REF##12011441##[16]##,##REF##12191639##[52]##.</p> <p>Comparisons between the genomes of mammals and other vertebrates are powerful tools in understanding how human genes and their products are regulated, what their function is and how and why they evolved ##REF##10521336##[53]##,##UREF##0##[54]##. Indeed, much of CTCF function has been characterised in chicken, including its capacity as an insulator protein ##REF##10458613##[5]## and recent studies have revealed the extreme conservation of <italic>CTCF</italic> sequence and function in amphibians ##REF##11900981##[13]##, fish ##REF##16647825##[14]## and even invertertebrates such as <italic>Drosophila</italic>
##REF##15678159##[15]##.</p> <p>Despite this, <italic>CTCF</italic> has not been characterised in non-eutherian mammals or reptiles and whether <italic>BORIS</italic> exists outside humans and mice is not even known. From reported failures to find <italic>BORIS</italic> sequence in chicken and fish ##REF##16647825##[14]##,##REF##12011441##[16]## it has been proposed that <italic>BORIS</italic> arose recently from duplication of <italic>CTCF</italic> in an early mammal ##REF##12011441##[16]##. However, here we report that <italic>BORIS</italic> orthologues are present in all major mammalian groups and at least two reptilian species, proving that <italic>BORIS</italic> evolution occurred much earlier than has been recognised. We examined the expression pattern of <italic>CTCF</italic> and <italic>BORIS</italic> in the three major mammalian clades and a reptile, discovering that although <italic>CTCF</italic> is ubiquitously expressed in all species, <italic>BORIS</italic> became progressively specialised to testis throughout amniote evolution. We consider these new data with respect to current theories regarding <italic>CTCF</italic> and <italic>BORIS</italic> as antagonistic epigenetic regulators and their roles in governing genomic imprinting at the <italic>IGF2/H19</italic> locus.</p>
<title>Materials and Methods</title> <title>Tissue</title> <p>Adult cattle tissue was sourced from commercial abattoirs processing farmed animals from New South Wales, Australia. Tissue from adult wallaby and platypus were sourced from a captive breeding colony of wallabies and a platypus tissue collection, both held at the Research School of Biological Sciences, Australian National University, Canberra, Australia. Juvenile central bearded dragon tissues samples were sourced from a captive breeding colony held at the University of Canberra, Australia. All tissue (excluding testes samples) was from females, except for platypus tissue which was male. The captivity and sacrifice of all animals was approved by the Australian National University (wallaby and platypus) and University of Canberra (bearded dragon) Animal Experimentation Ethics Committees (AEECP R.CG.08.03, R.CG.02.00 and CEAE 04/04 respectively). Sourcing of cattle tissue was exempt from AEEC approval, as these animals were not sacrificed primarily for research purposes (Simon Bain, ANU AEEC).</p> <title>Nucleic Acid Extraction, Amplification, and Sequencing</title> <p>Genomic DNA extraction was performed on liver tissue samples following the standard protocol for mammalian tissue ##UREF##2##[64]##. Total RNA was extracted using the GenElute Mammalian Total RNA Miniprep Kit (Sigma-Aldrich) according to the manufacturer's instructions. Eluted RNA was treated by DNAse digestion, using the DNA-free Dnase kit (Ambion) as recommended by the manufacturer. All samples were checked for quality and purity on a 1.2% denaturing formaldehyde agarose gel ##UREF##2##[64]##. RNA was tested for genomic DNA contamination by PCR prior to first strand synthesis of cDNA. Approximately 800 ng of purified RNA was used to create cDNA using the SuperScript III Reverse Transcriptase system (Invitrogen) according to manufacturer's instructions. All first strand synthesis reactions were undertaken using random hexamer primers except for Rapid Amplification of cDNA Ends (RACE) experiments, where the GeneRacer Oligo dT primer (Invitrogen) was used. Conventional PCR amplifications were performed in a 50 µL reaction, including either 1 µL of undiluted cDNA or 200 ng of genomic DNA as a template, 0.2 µM of each primer (##SUPPL##2##Table S1##) and the following reagents from Invitrogen; 1X PCR Buffer, 0.8 mM dNTP mixture (0.2 mM each), 1.5 mM MgCl<sub>2</sub> and 0.2 µL of Platinum Taq DNA Polymerase. Cycling conditions used were as follows; 94°C, 2 min; 34×(94°C, 30 sec; 61°C, 30 sec; 72°C, 1 min); 72°C, 10 min. When amplifications over 1000 bp were performed, extension times were increased by 1 min/kb. Nested PCR amplifications for 3′ RACE were also undertaken using this protocol, except with reduced cycle numbers, modified primer concentration and increased annealing temperatures as stipulated in the GeneRacer kit (Invitrogen) protocol.</p> <p>For initial gene expression studies 5 µL of <italic>CTCF</italic> and <italic>BORIS</italic> amplified products were combined together with 6 µL of loading buffer (30% glycerol, with light Bromophenol blue staining) and subjected to electrophoresis for 40min at 7.6 V/cm on a 1% agarose gel with TAE buffer and SYBR Safe DNA gel stain (Invitrogen). Gels photographs were illuminated with blue light and exposed using the Gel Logic 100 Imaging System (Kodak). Other than cropping, no alterations to these images were performed.</p> <p>Full length <italic>CTCF</italic> and <italic>BORIS</italic> cDNAs were amplified from liver and testes samples respectively and cloned using the TOPO TA Cloning Kit. Recombinant plasmid DNA was purified using the Wizard Plus SV Miniprep System and then combined with relevant primers (##SUPPL##2##Table S1##) for sequencing at the Australian Genome Research Facility.</p> <title>Real-Time PCR</title> <p>Real-time PCR was performed in 20 µL reactions using the QuantiTect SYBR Green PCR Kit (Qiagen) according to manufacturer's instructions. Amplifications were performed and detected with a Rotorgene 3000 cycler (Corbett Research) using the following cycling conditions; 95°C, 15 min; 45×(94°C, 30 sec; 58°C, 30 sec; 72°C, 20 sec); 72°C, 10 min. All experimental amplifications were performed in triplicate and averaged over two or three concordant results which varied by C<sub>t</sub> values of less than 0.7. Levels of <italic>CTCF</italic> and <italic>BORIS</italic> relative to <italic>GAPDH</italic> in each tissue and species were calculated using the comparative quantitation software supplied by Rotorgene. All products were checked for specificity by melt-curve analysis and electrophoresis.</p> <p>Primers used in this analysis were designed for each species from similar intron-spanning regions of <italic>CTCF</italic>, <italic>BORIS</italic> and <italic>GAPDH</italic> (##SUPPL##2##Table S1##). These primers were selected for high amplification efficiency (&gt;1.65) and low primer-dimer. A 10-fold serial dilution of testis cDNA was undertaken to determine the amplification range and performance of <italic>BORIS</italic> primers at low template concentrations, because <italic>BORIS</italic> (unlike <italic>CTCF</italic> and <italic>GAPDH</italic>) is known to have low or undetectable expression in many tissues ##REF##12011441##[16]##,##REF##17957795##[47]##. We found that <italic>BORIS</italic> transcripts could be detected reliably down to the 10<sup>−3</sup> dilution. Primers for the positive control gene <italic>GAPDH</italic> (##SUPPL##2##Table S1##) were designed from sequence deposited on NCBI for cattle (NM_001034034.1), platypus (EH003224) and wallaby (EF654515 and trace archive data). For bearded dragon, <italic>GAPDH</italic> primers were designed from sequence we determined ourselves by PCR amplification and sequencing (EU784660).</p> <title>Bioinformatic Analysis</title> <p>Homology searches were performed using BLASTn and tBLASTn ##REF##2231712##[65]## against the non-redundant, expressed sequence tag and trace archive databases at the NCBI website (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov">http://www.ncbi.nlm.nih.gov</ext-link>) or release 46 of the Ensembl website (<ext-link ext-link-type="uri" xlink:href="http://www.ensembl.org">http://www.ensembl.org</ext-link>). For species in which gene prediction was not available, or was unrealistic, we performed our own gene predictions using Genomescan ##REF##11337476##[66]## and local alignment. A multiple alignment of the resulting set of predicted and experimentally determined cDNA sequence (##SUPPL##3##Table S2##) was produced using ClustalW2 with default parameters (<ext-link ext-link-type="uri" xlink:href="http://www.ebi.ac.uk/Tools/clustalw2">http://www.ebi.ac.uk/Tools/clustalw2</ext-link>). Phylogenetic analysis was performed on aligned cDNA sequences by the neighbour-joining method with uncorrected distance measure, using the phylogenetic program PAUP* version 4.0 b 10 ##UREF##3##[67]##. 1,000 replications were performed for bootstrap analysis. Protein coding predictions of these cDNA sequences were also made and aligned using ClustalW2. This alignment was then used to calculate pairwise identity between selected orthologues using MacVector v9.5.2.</p> <p>The conserved block of genes orthologous to the region surrounding human <italic>BORIS</italic> was identified in amniote species by BLAST with the criteria of unique reciprocal best-hits back to the query sequence in the human genome. Genomic sequence from these orthologous blocks was extracted from Ensembl and aligned using the LAGAN algorithm ##REF##12654723##[68]## available on the mVISTA website with default parameters (<ext-link ext-link-type="uri" xlink:href="http://genome.lbl.gov/vista/mvista/submit.shtml">http://genome.lbl.gov/vista/mvista/submit.shtml</ext-link>).</p>
<title>Results</title> <p>We isolated, sequenced and characterised <italic>CTCF</italic> and <italic>BORIS</italic> homologues in eutherians, marsupials, monotremes and reptiles, and studied their expression profiles in one species from each of these vertebrate groups.</p> <title>Cloning and Characterisation of <italic>CTCF</italic> and <italic>BORIS</italic> Orthologues in Vertebrates</title> <p>Homologues of <italic>CTCF</italic> and <italic>BORIS</italic> were amplified from a range of amniotes by reverse-transcriptase PCR (RT-PCR) and rapid amplification of cDNA ends, using primers designed from sequenced genomic data or evolutionarily conserved regions (##SUPPL##2##Table S1##). Full-length or near full-length protein coding cDNA sequences were retrieved in this way from our model eutherian, marsupial, monotreme and reptilian species; domestic cattle (<italic>Bos taurus</italic>), tammar wallaby (<italic>Macropus eugenii</italic>), duck-billed platypus (<italic>Ornithorhynchus anatinus</italic>) and central bearded dragon (<italic>Pogona vitticeps</italic>) respectively (accession numbers EU527852-EU527858). Similarity searches, using these sequences and other annotated <italic>CTCF</italic> and <italic>BORIS</italic> sequences as queries, were conducted in a variety of databases hosted at NCBI (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov">http://www.ncbi.nlm.nih.gov</ext-link>) and Ensembl (<ext-link ext-link-type="uri" xlink:href="http://www.ensembl.org">http://www.ensembl.org</ext-link>). This approach identified a further 37 homologues of these genes in vertebrates (##SUPPL##3##Table S2##).</p> <p>From the largest region of common overlap between these homologues, a neighbour joining tree was constructed, revealing two distinct clusters of sequence (##FIG##1##Figure 2##). One of these clusters contained previously annotated copies of <italic>CTCF</italic> from human (NM_006565), mouse (NM_007794), rat (NM_031824), cattle (NM_001075748), chicken (NM_205332) and zebrafish (NM_001001844). The other cluster contained annotated <italic>BORIS</italic> sequence from human (NM_080618) and mouse (NM_001081387). The branch separating these two clusters was supported by a 100% bootstrap value. This unambiguously defined which sequences were <italic>CTCF</italic> orthologues and which were <italic>BORIS</italic> orthologues.</p> <p>In line with previous studies, we detected <italic>CTCF</italic> orthologues in all major vertebrate groups ##REF##9756895##[12]##–##REF##16647825##[14]##. Included in this cluster were closely related duplicate <italic>CTCF</italic> sequences (designated <italic>CTCF</italic>1 and <italic>CTCF</italic>2) from stickleback and medaka. The <italic>BORIS</italic> cluster included, as well as orthologues from many eutherian species, clear orthologues in two marsupials (Gray short-tailed opossum, <italic>Monodelphis domestica</italic>, and wallaby), a monotreme (platypus) and two reptiles (bearded dragon and green anole, <italic>Anolis carolinensis</italic>). No orthologues of <italic>BORIS</italic> could be detected using nucleotide BLAST, or translated BLAST searches in genomes of any avian (chicken, <italic>Gallus gallus</italic>; and zebra finch, <italic>Taeniopygia guttata</italic>), amphibian (Western-clawed frog, <italic>Xenopus tropicalis</italic>), teleost fish (puffer fish, <italic>Takifugu rubripes</italic> and <italic>Tetraodon nigroviridis</italic>; zebrafish, <italic>Danio rerio</italic>; stickleback, <italic>Gasterosteus aculeatus</italic>; and medaka <italic>Oryzias latipes</italic>) or primitive vertebrate (sea lamprey, <italic>Petromyzon marinus</italic>).</p> <p>Although BLAST searches failed to identify any sequence orthologous to <italic>BORIS</italic> in bird, amphibian and fish genomes, it remained possible that <italic>BORIS</italic> is present in these genomes but is too diverged to detect using standard alignment methods. This seemed particularly likely for chicken, as birds are a sister taxon to the reptiles, in which we discovered <italic>BORIS</italic> orthologues. Applying a strategy previously used in the search for divergent genes ##REF##17333539##[55]##, we sought orthologues of markers on either side of human <italic>BORIS.</italic> We located such sequences in multiple species, and searched the dividing spaces for <italic>BORIS</italic>-like sequence.</p> <p>We found that genes flanking <italic>BORIS</italic> in humans were part of a single large block of genes (<italic>TMEPAI-BMP7</italic>) clustered together in the same orientation in all tetrapods (data not shown). Genes from this block were either not clustered together, or were not present in sea lamprey and teleost fish genomes.</p> <p>We aligned the regions containing genes immediately adjacent to <italic>BORIS</italic> (<italic>PCK1</italic> and <italic>RBM38</italic>) between human, mouse, dog, opossum, platypus, chicken, green anole and frog (##FIG##2##Figure 3##). As before, we could detect no <italic>BORIS</italic> orthologues in frog, but found some similarity between the first zinc finger of <italic>BORIS</italic> and a 108-bp region of the chicken <italic>PCK1-RBM38</italic> intergenic sequence. When this sequence was used a query for reciprocal BLAST against the entire human genome, the best alignments were to the first zinc finger of <italic>CTCF</italic> and <italic>BORIS</italic>, indicating that these sequences were homologous. We could uncover no evidence for this sequence being part of an active gene other than finding that it overlaps with an Ensembl <italic>ab-initio</italic> gene prediction (GENSCAN00000030237). We therefore conclude that the sequence is a degraded relic of <italic>BORIS</italic>.</p> <title>Sequence Analysis and Gene Structure of <italic>CTCF</italic> and <italic>BORIS</italic> in Vertebrates</title> <p>We aligned predicted CTCF and BORIS proteins and found that, like previously annotated versions of these proteins, all possessed eleven ZFs, ten of which belong to the Cys<sub>2</sub>His<sub>2</sub> class, and one which belongs to the Cys<sub>2</sub>HisCys class (##SUPPL##2##Figure S1##). As previously reported for human, mouse, chicken, zebrafish and frog ##REF##9756895##[12]##–##REF##16647825##[14]##,##REF##12011441##[16]##, we found that all vertebrate CTCF orthologues are extremely highly conserved throughout their entire length. From pairwise alignments over the entire length of its sequence, we found 92% average identity between human CTCF and other selected vertebrate CTCF sequences (##TAB##0##Table 1##). In comparison, similarity of BORIS orthologues, to each other and to CTCF, was largely restricted to the region encoding the ZFs. When human BORIS was compared to other BORIS sequences the average identity was 80.4% within the ZF domain, but less than 35% similar in the other two regions. Moreover, comparisons of human BORIS with CTCF sequences produced an average identity of 74.1% within the ZF domain, but less than 15% conservation within the other regions.</p> <p>
<italic>CTCF</italic> genomic sequence from human, mouse, zebrafish and frog were all reported to have ten protein-encoding exons when they were first characterised ##REF##11900981##[13]##,##REF##16647825##[14]##,##REF##12011441##[16]##. In contrast, chicken <italic>CTCF</italic> was reported by Klenova <italic>et al</italic>. ##REF##9756895##[12]## to only have seven protein coding exons, four of which contained all eleven zinc fingers. We analysed the gene structure of <italic>CTCF</italic> in all species from which there was full genomic sequence and found that all sequences, including chicken <italic>CTCF</italic>, contained ten exons in total, with seven ZF exons (##FIG##0##Figure 1B##). <italic>BORIS</italic> orthologues were also found to have a very similar structure, especially within the ZF domain where intron-exon boundaries were identical.</p> <title>Gene Expression Analysis</title> <p>One of the most remarkable characteristics of <italic>CTCF</italic> and <italic>BORIS</italic> is that in humans and mice they show apparently mutually exclusive expression. <italic>BORIS</italic> is transcribed only in specific parts of the testis, while <italic>CTCF</italic> is expressed in all tissues except those expressing <italic>BORIS</italic>
##REF##12011441##[16]##,##REF##17048991##[38]##. This expression pattern underpins the hypothesis that <italic>BORIS</italic> is the key regulator establishing the male germline imprint of <italic>IGF2</italic>, that it acts antagonistically to <italic>CTCF</italic> and defines its inclusion within the cancer-testis group of genes.</p> <p>To determine if this expression pattern is conserved more widely in vertebrates, we examined the transcription of <italic>CTCF</italic> and <italic>BORIS</italic> in cattle, wallaby, platypus and bearded dragon. Initially, we performed 35 cycles of RT-PCR on a series of tissues using <italic>CTCF/BORIS</italic> primers anchored within at least one of the ZFs and a surrounding non-zinc finger region (##SUPPL##2##Table S1##). <italic>CTCF</italic> transcripts were detected in this way for all tissues and animals tested (##FIG##3##Figure 4A##). <italic>BORIS</italic> transcripts were detected only in the gonads of cattle and wallaby; strongly in testes, and weakly in ovarian samples. In contrast, <italic>BORIS</italic> was amplified from a much wider set of somatic and reproductive tissues in platypus (brain, heart, liver, kidney and testis); and bearded dragon (brain, lung, liver, kidney, spleen, testis and ovary). To minimise the possibility we were observing tissue specific splice variants, we repeated our RT-PCR experiments using primers from different regions of <italic>BORIS</italic> (##SUPPL##2##Table S1##), and found similar results (data not shown).</p> <p>Despite these discoveries, due to the nature of conventional ‘end-point observed’ RT-PCR our initial experiments were semi-quantitative at best. Thus, we were unsure if the expression we were observing was at a level which was biologically relevant. A recent publication using the quantitative real-time PCR technique found that although <italic>BORIS</italic> expression is considered to be restricted to the testis and some tumours, <italic>BORIS</italic> transcripts could be detected in other tissues up to 0.3% of the level of <italic>BORIS</italic> in the testis ##REF##17957795##[47]##. The authors concluded from this that expression of <italic>BORIS</italic> less than 0.3% of the level in testis was not biologically relevant.</p> <p>We performed real-time PCR amplifications of <italic>CTCF</italic> and <italic>BORIS</italic> on all our available tissues in triplicate and comparatively quantified their respective levels using the Corbett Research Rotorgene system, with SYBR Green as the fluorescent DNA-binding dye. Differences in template concentration within a species were taken into consideration by normalising our results to the housekeeping gene glyceraldehyde-3-phosphate dehydrogenase (<italic>GAPDH</italic>). In agreement with our initial RT-PCR experiments, amplification of <italic>CTCF</italic> occurred in all tissues and species reproducibly, but with up to 50-fold variation between tissues (##SUPPL##1##Figure S2##) not unlike that seen previously in developing zebrafish and frog ##REF##11900981##[13]##,##REF##16647825##[14]##.</p> <p>Like previous experiments, we found that <italic>BORIS</italic> amplifications by real-time PCR were predominantly from the testis, with consistently high expression between 10% and 100% of the level of <italic>GAPDH</italic> (##SUPPL##1##Figure S2##). As expected, <italic>BORIS</italic> amplification was also detected on multiple occasions outside of the testis, particularly within platypus and bearded dragon, and at levels well within the expected limitations of our assay (see <xref ref-type="sec" rid="s4">methods</xref>). When we applied the 0.3% cut-off defined by Kholmanskikh <italic>et al</italic>., ##REF##17957795##[47]## to our results, we found that as in humans, levels of <italic>BORIS</italic> in somatic tissues were below this threshold in cattle and wallaby (&lt;0.2% of testis expression), while ovarian <italic>BORIS</italic> levels were just on (cattle) or above (wallaby) this threshold (##FIG##3##Figure 4B##).</p> <p>Expression of <italic>BORIS</italic> outside of the testis in platypus and bearded dragon was much higher. <italic>BORIS</italic> transcripts in the liver and kidney of platypus was within 6–10% of that found in platypus testis (##FIG##3##Figure 4B##), and was at a level comparable to <italic>CTCF</italic> expression found in these tissues (##SUPPL##1##Figure S2##). Likewise, levels of <italic>BORIS</italic> transcripts in bearded dragon brain, kidney and ovary were 2–5% of the level of <italic>BORIS</italic> in testis. <italic>BORIS</italic> transcripts were detected in other five other somatic tissues of platypus (brain, heart and spleen) and bearded dragon (lung and spleen) at levels just on or above the 0.3% threshold.</p>
<title>Discussion</title> <p>We examined the sequence and expression of key epigenetic regulators <italic>CTCF</italic> and <italic>BORIS</italic> in vertebrates through cloning, sequencing, bioinformatic analysis and quantitative gene expression experiments. Our results are at variance with the hypothesis that <italic>BORIS</italic> arose recently by duplication of <italic>CTCF</italic> in mammals, and was quickly specialised for a role in germ cell imprinting that was complementary, or even antagonistic, to the role of <italic>CTCF</italic>.</p> <title>
<italic>BORIS</italic> First Arose in Early Amniotes</title> <p>Previous studies established that <italic>CTCF</italic> is a highly conserved and ubiquitous gene in humans and mice, as well as other vertebrates including birds, fish and amphibians ##REF##9756895##[12]##–##REF##16647825##[14]##. Our studies on cattle, wallaby, platypus and dragon lizard confirm the expectation that <italic>CTCF</italic> is highly conserved in all vertebrate groups, and is expressed to varying degrees in all tissues of eutherian, marsupial and monotreme mammals, as well as reptiles (##FIG##3##Figures 4A## and ##SUPPL##1##S2##).</p> <p>In contrast to the well-studied <italic>CTCF</italic> gene, much less is known about the evolutionary history and function of <italic>BORIS</italic>. <italic>BORIS</italic> sequence was previously determined only in humans and mice ##REF##12011441##[16]##, and no orthologue was detected in chicken. This gave rise to the speculation that <italic>BORIS</italic> duplicated from <italic>CTCF</italic> only recently in the mammal lineage. In addition, chicken <italic>CTCF</italic> was reported to have a gene structure significantly different from that of mammal <italic>CTCF</italic> and <italic>BORIS</italic>. Chicken was therefore considered to represent the ancestral gene structure, and an alteration of <italic>CTCF</italic> gene structure was proposed to have occurred in the mammalian ancestor, followed by a duplication to give rise to <italic>BORIS</italic>.</p> <p>We found that chicken <italic>CTCF</italic> was not, after all, different in structure from mammal <italic>CTCF</italic> as was previously reported ##REF##9756895##[12]##,##REF##12011441##[16]## (##FIG##0##Figure 1B##). The chicken genome project had not been undertaken when chicken <italic>CTCF</italic> was initially sequenced, so sequence coverage from this early study may not have been sufficient to build a reliable assembly of the region. Alternatively, the <italic>CTCF</italic> clone that was sequenced may have been a cDNA and genomic DNA chimaera.</p> <p>Unexpectedly, we found orthologues of <italic>BORIS</italic> in at least two reptilian species (bearded dragon and green anole). This means that the duplication of <italic>CTCF</italic> which gave rise to <italic>BORIS</italic> must have occurred prior to the divergence of sauropsids (birds and reptiles) and mammals 210–310 million years ago (##FIG##4##Figure 5##). In agreement with Loukinov <italic>et al</italic>. ##REF##12011441##[16]## we could find no full orthologue of <italic>BORIS</italic> within the chicken genome. However, by analysing the intergenic region between markers flanking the expected site of <italic>BORIS</italic> in chicken, we did discover a small 108-bp segment of DNA homologous to the first zinc finger of <italic>BORIS</italic> (##FIG##2##Figure 3##). Although this region of DNA may be part of another functional gene, we consider that it is unlikely to be functionally related to other <italic>BORIS</italic> orthologues, given that no other regions showed conservation, even the usually well-conserved zinc fingers. We conclude that either <italic>BORIS</italic> succumbed to pseudogenisation in birds some time after they diverged from reptiles, or underwent a rapid functional change leaving behind only small traces of its evolutionary past.</p> <p>In an extension of previous studies ##REF##12011441##[16]##, we found extremely high conservation between vertebrate CTCF orthologues, but observed that BORIS homologues were similar to each other, and to CTCF, only within the ZF domains (##TAB##0##Table 1##). These observations support the prediction of Loukinov <italic>et al</italic>. that any major differences between CTCF and BORIS function are probably attributable to the N- and C-terminal domains, given these are the most divergent. In fact, the N- and C-terminal domains of CTCF and BORIS contained only small pockets of sequence that were obviously alignable (##SUPPL##0##Figure S1##). Interestingly, two of these conserved regions overlapped the start and end of these proteins, implying that the duplication that gave rise to <italic>BORIS</italic> must have involved the entire <italic>CTCF</italic> sequence.</p> <p>The ZF domain of CTCF orthologues we examined showed an almost perfect (99.5% average) identity with human CTCF. In comparison, the average conservation between the ZF domain of human BORIS and other BORIS orthologues was much lower (80.4%). This suggests that BORIS experienced a decrease in functional constraint relative to CTCF, initially because it was a duplicate, and presently because it only binds a subset of the sites bound by CTCF. Alternatively, BORIS may bind some sequences not recognised by CTCF ##REF##17048991##[38]##. In support of this, we found that although all of the amino acids thought to perform protein-DNA interactions ##REF##11395410##[56]## were 100% conserved for vertebrate CTCF, many were not conserved in some, or all BORIS orthologues (##SUPPL##0##Figure S1##). It would be interesting to investigate this further by mapping BORIS binding sites in the genome relative to the published CTCF binding sites ##REF##17382889##[2]##,##REF##17512414##[3]##.</p> <title>Expression of Ancestral <italic>BORIS</italic> in Somatic Tissue</title> <p>Our RT-PCR experiments showed two main patterns of <italic>BORIS</italic> expression (##FIG##3##Figure 4##). In the marsupial and the eutherian (wallaby and cattle respectively) we found predominantly testis-specific expression with some ovarian expression, whereas in the reptile (bearded dragon) and the monotreme (platypus) we detected expression of <italic>BORIS</italic> in multiple somatic tissues as well as the gonads. When these experiments were repeated using quantitative real-time PCR, we discovered that after 45 cycles of PCR, some <italic>BORIS</italic> transcripts could be detected outside of the germline in cattle and wallaby. However, we found that these levels of <italic>BORIS</italic> were extremely low, approaching the limits of detection and falling under a previously defined threshold for meaningful expression of <italic>BORIS</italic>
##REF##17957795##[47]##. Thus, we expect that <italic>BORIS</italic> function in cattle and wallaby is absent from somatic tissues, just as is predicted in humans and mice. More experiments will be required to determine if the ovarian expression of <italic>BORIS</italic> in cattle and particularly wallaby is functionally significant.</p> <p>The highest levels of <italic>BORIS</italic> expression outside of the testis were found in platypus and bearded dragon. The most striking examples of these came from the liver and kidney of platypus and the brain, kidney and ovary of bearded dragon, which were at levels 2-10% of <italic>BORIS</italic> expression in the testis (##FIG##3##Figure 4B##). In two cases (platypus liver and kidney) this level of expression was close to the level of <italic>CTCF</italic> expression within the same tissues (##SUPPL##1##Figure S2##). These results strongly suggest that <italic>BORIS</italic> in platypus and bearded dragon functions outside of the testes, including in the ovary and multiple somatic tissues. This finding is of significance because it indicates that <italic>BORIS</italic> had wide expression in an ancestral amniote, similar to that of <italic>CTCF,</italic> the gene from which it arose by duplication.</p> <p>The question then arises, why was the <italic>CTCF</italic> duplicate (or “proto-<italic>BORIS</italic>”) initially retained and why did it succumb to evolutionary change? The high degree of <italic>CTCF</italic> conservation throughout vertebrates implies that it is a gene under extreme functional constraint. Accordingly, perhaps it is not surprising that a <italic>CTCF</italic> duplicate in a new genomic environment would be retained and undergo sub-functionalisation, alleviating some mutational load upon CTCF. Subfunctionalisation may also explain why duplicate copies of <italic>CTCF</italic> have been retained in the genomes of medaka and stickleback (##FIG##1##Figure 2##), following whole-genome duplication of early teleost fish ##REF##12618368##[57]##.</p> <p>Although we observed <italic>CTCF</italic> and <italic>BORIS</italic> expression alongside each other in some tissues of monotremes and reptiles, these genes are apparently not co-expressed in humans and mice and may even be antagonistic. CTCF and BORIS bind competitively to common sites and display opposing effects on the epigenetic status of the <italic>Igf2/H19</italic> ICR and transcription of <italic>BAG1</italic> and the CT-genes ##REF##12011441##[16]##,##REF##16140943##[42]##,##REF##16140944##[43]##,##REF##17260018##[49]##,##REF##12191639##[52]##. Thus, at some stage during the evolution of therian mammals, <italic>CTCF</italic> and <italic>BORIS</italic> evolved mutually exclusive expression and potential antagonism. As our studies were performed on whole tissues, we could not resolve whether <italic>CTCF</italic> and <italic>BORIS</italic> show mutually exclusive expression amongst the many discrete cell-types in testis and ovary in wallaby and cattle, so we cannot pinpoint when mutually exclusive expression arose in therian mammals after their divergence from the monotremes.</p> <title>
<italic>BORIS</italic> Specialisation Correlates with the Evolution of Imprinting</title> <p>To date, the only non-pathological function proposed for <italic>BORIS</italic> is the establishment of paternal-specific methylation at the <italic>Igf2/H19</italic> ICR in mice ##REF##12011441##[16]##,##REF##17048991##[38]##. If found to be true for mice, it seems likely that this function is conserved in humans, since they also possess a paternally methylated CTCF-dependent insulator (the ICR) ##REF##10839547##[34]##,##REF##10839546##[35]## and testis-specific <italic>BORIS</italic> expression which is exclusive of <italic>CTCF</italic>
##REF##12011441##[16]##. Moreover, differential methylation of CTCF/BORIS binding sites upstream of a maternally-expressed <italic>H19</italic> orthologue has been discovered in sheep and wallaby ##REF##14654216##[58]##-##UREF##1##[60]##, suggesting that the mouse model of <italic>Igf2/H19</italic> imprinted regulation and <italic>BORIS</italic> function may be conserved throughout all therians. Yet, <italic>BORIS</italic> is not expected to have this function in reptiles and monotremes, or the amniotic ancestor from which <italic>BORIS</italic> first arose, as <italic>IGF2</italic> imprinting evolved after the divergence of monotremes from therian mammals. Our finding that <italic>BORIS</italic> expression is gonad-specific in wallaby and cattle, both of which possess imprinting of <italic>IGF2</italic>
##REF##15044262##[21]##,##REF##15652708##[23]##, implies that restriction of <italic>BORIS</italic> expression to the germline correlates with the evolution of genomic imprinting at <italic>IGF2/H19</italic> and other loci (reviewed ##REF##17683825##[61]##).</p> <p>In support of this, the evolution of another essential regulator of the <italic>Igf2/H19</italic> ICR is also strongly correlated with the evolution of imprinting. Orthologues of the <italic>de novo</italic> methyltransferase family member <italic>DNMT3L</italic> are present in eutherians and marsupials (which posses imprinting), but apparently not in chicken, fish ##REF##16575165##[62]## or platypus (T.H., unpublished data) which are thought to lack imprinting.</p> <title>Model of CTCF and BORIS Evolution in Amniotes</title> <p>We propose that a duplication of <italic>CTCF</italic> occurred in a common ancestor of all amniotes, probably some time after their divergence from amphibians 350-310MYA (##FIG##4##Figure 5##). We predict that originally this ‘proto-<italic>BORIS</italic>’ functioned alongside <italic>CTCF</italic>, perhaps subfunctionalising to take on tissue-specific roles from the highly conserved and functionally constrained CTCF protein. When genomic imprinting arose in early therian mammals 210-180MYA, <italic>BORIS</italic> was recruited to perform imprint establishment in germ cells, and <italic>CTCF</italic> imprint interpretation at <italic>IGF2/H19</italic> and potentially other imprinted genes. We speculate that this specialisation marked the start of antagonism between <italic>BORIS</italic> and <italic>CTCF</italic>, through the development of opposing epigenetic effects at the common loci to which they bound. The result of this was restriction of <italic>BORIS</italic> expression to the gonads of early therian mammals, and later restriction to the testes in the ancestor of humans and mice.</p> <p>The divergent nature and proposed clash of function between <italic>CTCF</italic> and <italic>BORIS</italic> has often been described as ‘sibling-rivalry’ ##REF##12011441##[16]##,##REF##12191639##[52]##. Our results show that this rivalry did not always exist, and ironically may have evolved in response to the evolution of genomic imprinting, which is in turn thought to have evolved from other conflicts in the family ##REF##2035190##[63]##.</p>
<p>Conceived and designed the experiments: TAH. Performed the experiments: TAH. Analyzed the data: TAH JED JAMG. Wrote the paper: TAH JAMG.</p> <p>CTCF is an essential, ubiquitously expressed DNA-binding protein responsible for insulator function, nuclear architecture, and transcriptional control within vertebrates. The gene <italic>CTCF</italic> was proposed to have duplicated in early mammals, giving rise to a paralogue called “brother of regulator of imprinted sites” (<italic>BORIS</italic> or <italic>CTCFL</italic>) with DNA binding capabilities similar to <italic>CTCF</italic>, but testis-specific expression in humans and mice. <italic>CTCF</italic> and <italic>BORIS</italic> have opposite regulatory effects on human cancer-testis genes, the anti-apoptotic <italic>BAG1</italic> gene, the insulin-like growth factor 2/H19 imprint control region (<italic>IGF2/H19</italic> ICR), and show mutually exclusive expression in humans and mice, suggesting that they are antagonistic epigenetic regulators. We discovered orthologues of <italic>BORIS</italic> in at least two reptilian species and found traces of its sequence in the chicken genome, implying that the duplication giving rise to <italic>BORIS</italic> occurred much earlier than previously thought. We analysed the expression of <italic>CTCF</italic> and <italic>BORIS</italic> in a range of amniotes by conventional and quantitative PCR. <italic>BORIS</italic>, as well as <italic>CTCF</italic>, was found widely expressed in monotremes (platypus) and reptiles (bearded dragon), suggesting redundancy or cooperation between these genes in a common amniote ancestor. However, we discovered that <italic>BORIS</italic> expression was gonad-specific in marsupials (tammar wallaby) and eutherians (cattle), implying that a functional change occurred in <italic>BORIS</italic> during the early evolution of therian mammals. Since therians show imprinting of <italic>IGF2</italic> but other vertebrate taxa do not, we speculate that <italic>CTCF</italic> and <italic>BORIS</italic> evolved specialised functions along with the evolution of imprinting at this and other loci, coinciding with the restriction of <italic>BORIS</italic> expression to the germline and potential antagonism with <italic>CTCF</italic>.</p> <title>Author Summary</title> <p>Epigenetic mechanisms heritably change gene expression without altering DNA sequence. Currently, little is known about the evolution of epigenetic traits, and the genes that control them. <italic>CTCF</italic> is an essential epigenetic regulator that is expressed widely in the tissues of vertebrates and modifies the transcription of genes by altering their location within the nucleus. <italic>CTCF</italic> duplicated at some time during vertebrate evolution, giving rise to a similar gene called <italic>BORIS</italic> with expression that is limited to parts of the testes of humans and mice, but whose function is largely unknown. <italic>BORIS</italic> may contribute to the regulation of genomic imprinting, a form of epigenetic control specific to live-bearing mammals. We discovered <italic>BORIS</italic> in all mammal groups and reptiles, implying that its genesis from <italic>CTCF</italic> occurred much earlier than previously thought, preceding genomic imprinting by over 100 million years. <italic>CTCF</italic> and <italic>BORIS</italic> have not previously been found expressed together except in tumours, leading to the hypothesis that <italic>CTCF</italic> and <italic>BORIS</italic> have conflicting functions that cause cancer when allowed to overlap. We found that <italic>CTCF</italic> and <italic>BORIS</italic> are expressed alongside each other in multiple somatic and reproductive tissues of a reptile (bearded dragon) and an egg-laying mammal (platypus), but that <italic>BORIS</italic> is restricted to the gonads of live-bearing mammals (cattle and wallaby). This indicates that <italic>BORIS</italic> specialised during mammalian evolution, in concert with the evolution of genomic imprinting.</p>
<title>Supporting Information</title>
<p>The authors would like to acknowledge the following people for assisting in the supply of tissue for this study: Mr. Arthur Cullen and Dr. Scott Hammond (cattle), Dr. Lauren Marotte (wallaby), Dr. Frank Grützner (platypus) and Dr. Tariq Ezaz and Dr. Alex Quinn (bearded dragon). The authors would also like to thank Dr. Paul Waters and Miss Donna Bond for assistance with phylogenetic analysis and real-time PCR respectively.</p>
<fig id="pgen-1000169-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000169.g001</object-id><label>Figure 1</label><caption><title>Gene structure of <italic>CTCF</italic> and <italic>BORIS</italic>.</title><p>(A) <italic>CTCF</italic> and <italic>BORIS</italic> share a similar ZF domain, but different N- and C-terminal domains. (B) All vertebrate <italic>CTCF</italic> orthologues posses ten exons. Intron-exon boundaries are identical between all <italic>CTCF</italic> and <italic>BORIS</italic> orthologues within the ZF domain (grey). Note, genomic coverage of platypus <italic>BORIS</italic> is incomplete at the 5′ end (*).</p></caption></fig> <fig id="pgen-1000169-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000169.g002</object-id><label>Figure 2</label><caption><title>Neighbour-joining tree showing relationships between members of the <italic>CTCF</italic> and <italic>BORIS</italic> gene family.</title><p>Sequence we determined experimentally (blue) and discovered by <italic>in silico</italic> similarity searches (black) form two distinct clusters with previously annotated <italic>CTCF</italic> and <italic>BORIS</italic> orthologues (bold). These clusters are separated from each other by a branch with 100% bootstrap value (thick line). Accession numbers and scientific names for these sequences and species are shown in ##SUPPL##3##Table S2##.</p></caption></fig> <fig id="pgen-1000169-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000169.g003</object-id><label>Figure 3</label><caption><title>Human genomic sequence encompassing <italic>PCK1-BORIS-RBM38</italic> (top) compared to the orthologous regions in other amniotes.</title><p>High similarity over a 100bp window is seen for most exonic sequence (blue) and some untranslated regions (UTR, light blue) or non coding regions (NCR, pink). Despite no similarity to any other region of human <italic>BORIS</italic> to the chicken <italic>PCK1-RBM38</italic> region, the peak labelled with a star is homologous to the first ZF of <italic>BORIS</italic>.</p></caption></fig> <fig id="pgen-1000169-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000169.g004</object-id><label>Figure 4</label><caption><title>Expression analysis of <italic>CTCF</italic> and <italic>BORIS</italic>.</title><p>(A) Conventional RT-PCR of <italic>CTCF</italic> and <italic>BORIS</italic> after 35 cycles. Note, cattle brain, platypus ovary and bearded dragon muscle could not be tested due to tissue unavailability or poor RNA quality. (B) <italic>BORIS</italic> transcript levels relative to the positive control gene <italic>GAPDH</italic> as quantified by real-time PCR. To assist comparisons between species, we set <italic>BORIS</italic> expression in the testis to 1 and adjusted all other values within the same species proportionally. As found in humans, expression of <italic>BORIS</italic> in somatic tissues of cattle and wallaby did not exceed 0.3% of that found in testis, although ovarian expression is just on or above this threshold. In contrast, levels of <italic>BORIS</italic> expression in platypus (liver and kidney) and bearded dragon (brain, kidney and ovary) are well above this threshold and in some cases match <italic>CTCF</italic> expression (##SUPPL##1##Figure S2##).</p></caption></fig> <fig id="pgen-1000169-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000169.g005</object-id><label>Figure 5</label><caption><title>Proposed model of <italic>CTCF</italic> and <italic>BORIS</italic> evolution in amniotes.</title><p>The expression of <italic>CTCF</italic> and <italic>BORIS</italic> is indicated (black = expressed, white = not expressed) within various taxa. The ancestral expression pattern of <italic>BORIS</italic> is wide, including multiple somatic tissues (reptiles and monotremes), but becomes progressively restricted in therian mammals with gonad specific expression in marsupials and cattle (##FIG##3##Figure 4##) and testis-specific expression in humans and mice ##REF##12011441##[16]##,##REF##17048991##[38]##,##REF##17957795##[47]##. Significant events in the evolution of <italic>CTCF</italic> and <italic>BORIS</italic> are marked with respect to the phylogenetic tree.</p></caption></fig>
<table-wrap id="pgen-1000169-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000169.t001</object-id><label>Table 1</label><caption><title>Average pairwise similarity (%) between regions of human CTCF/BORIS and other vertebrate orthologues.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td colspan="2" align="left" rowspan="1"/><td colspan="4" align="left" rowspan="1">Human CTCF</td><td colspan="4" align="left" rowspan="1">Human BORIS</td></tr><tr><td colspan="2" align="left" rowspan="1"/><td align="left" rowspan="1" colspan="1">N-term</td><td align="left" rowspan="1" colspan="1">ZF</td><td align="left" rowspan="1" colspan="1">C-term</td><td align="left" rowspan="1" colspan="1">Total</td><td align="left" rowspan="1" colspan="1">N-term</td><td align="left" rowspan="1" colspan="1">ZF</td><td align="left" rowspan="1" colspan="1">C-term</td><td align="left" rowspan="1" colspan="1">Total</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<bold>Average pairwise identity with vertebrate</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>CTCF orthologues</bold>
</td><td align="left" rowspan="1" colspan="1">90.1</td><td align="left" rowspan="1" colspan="1">99.5</td><td align="left" rowspan="1" colspan="1">80.7</td><td align="left" rowspan="1" colspan="1">92.0</td><td align="left" rowspan="1" colspan="1">14.2</td><td align="left" rowspan="1" colspan="1">74.1</td><td align="left" rowspan="1" colspan="1">10.2</td><td align="left" rowspan="1" colspan="1">38.0</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>BORIS orthologues</bold>
</td><td align="left" rowspan="1" colspan="1">9.7</td><td align="left" rowspan="1" colspan="1">72.8</td><td align="left" rowspan="1" colspan="1">8.6</td><td align="left" rowspan="1" colspan="1">35.8</td><td align="left" rowspan="1" colspan="1">32.3</td><td align="left" rowspan="1" colspan="1">80.4</td><td align="left" rowspan="1" colspan="1">23.7</td><td align="left" rowspan="1" colspan="1">53.1</td></tr></tbody></table></alternatives></table-wrap>
<supplementary-material content-type="local-data" id="pgen.1000169.s001"><label>Figure S1</label><caption><p>Alignment of predicted CTCF and BORIS sequence in a range of vertebrates. Amino acids from the 11 zinc fingers (blue), likely to interact with DNA are shown indicated with arrows, and regions of similarity at the start and end of CTCF and BORIS are highlighted in red. Note, this figure is not intended for printing. Rather, it should be viewed in a program such as Adobe Reader, which has zoom function.</p><p>(12.39 MB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000169.s002"><label>Figure S2</label><caption><p>Quantification of (A) <italic>CTCF</italic> and (B) <italic>BORIS</italic> transcripts in various tissues and species relative to <italic>GAPDH</italic> as determined by real-time PCR.</p><p>(0.47 MB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000169.s003"><label>Table S1</label><caption><p>Primers used in this study.</p><p>(0.19 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000169.s004"><label>Table S2</label><caption><p>Vertebrate homologues of <italic>CTCF</italic> and <italic>BORIS</italic>.</p><p>(0.11 MB DOC)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><label/><p>Note, as not all sequences we discovered (##SUPPL##3##Table S2##) were of ideal length for pairwise comparisons, identities were calculated using human, mouse, dog, cattle, elephant, opossum, wallaby, platypus, chicken, bearded dragon, green anole, frog (<italic>X. tropicalis</italic> and <italic>X. laevis</italic>) and zebrafish sequences.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p>This work was supported by funding from the Australian Research Council and a block grant from the Research School of Biological Sciences.</p></fn></fn-group>
<graphic xlink:href="pgen.1000169.g001"/> <graphic xlink:href="pgen.1000169.g002"/> <graphic xlink:href="pgen.1000169.g003"/> <graphic id="pgen-1000169-t001-1" xlink:href="pgen.1000169.t001"/> <graphic xlink:href="pgen.1000169.g004"/> <graphic xlink:href="pgen.1000169.g005"/>
<media xlink:href="pgen.1000169.s001.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000169.s002.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000169.s003.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000169.s004.doc"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["54"], "element-citation": ["\n"], "surname": ["Volff"], "given-names": ["J-N"], "year": ["2006"], "article-title": ["Vertebrate genomes."], "publisher-loc": ["Basel; New York"], "publisher-name": ["Karger"], "fpage": ["vii, 216"]}, {"label": ["60"], "element-citation": ["\n"], "surname": ["Smits", "Mungall", "Griffiths-Jones", "Smith", "Beury"], "given-names": ["G", "AJ", "S", "P", "D"], "year": ["2008"], "article-title": ["Conservation of the H19 noncoding RNA and H19-IGF2 imprinting mechanism in therians."], "source": ["Nat Genet. In press"], "comment": ["doi:10.1038/ng.168"]}, {"label": ["64"], "element-citation": ["\n"], "surname": ["Sambrook", "Maniatis", "Fritsch"], "given-names": ["J", "T", "EF"], "year": ["1989"], "article-title": ["Molecular cloning : a laboratory manual."], "publisher-loc": ["Cold Spring Harbor, N.Y"], "publisher-name": ["Cold Spring Harbor Laboratory"]}, {"label": ["67"], "element-citation": ["\n"], "surname": ["Swofford"], "given-names": ["DL"], "year": ["2003"], "article-title": ["PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). 4 ed."], "publisher-loc": ["Sunderland, Massachusetts"], "publisher-name": ["Sinauer Associates"]}]
PLoS Genet. 2008 Aug 29; 4(8):e1000169